One-dimensional, time dependent simulation of the planetary boundary layer over a 48-hour period

International Nuclear Information System (INIS)

Haschke, D.; Gassmann, F.; Rudin, F.

1978-05-01

Results of a one-dimensional, time dependent simulation of the planetary boundary layer are given. First, a description of the mathematical model used is given and its approximations are discussed. Then a description of the initial and boundary conditions used for the simulation is given. Results are discussed with respect to their agreement with observed data and their precision. It can be demonstrated that a simulation of the planetary boundary layer is possible with satisfactory precision. The incompleteness of observed data gives, however, problems with their use and thus introduces uncertainties into the simulation. As a consequence, the report tries to point to the inherent limitations of such a simulation. (Auth.)

A note on the physical interpretation of frequency dependent boundary conditions in a digital waveguide mesh

DEFF Research Database (Denmark)

Escolano-Carrasco, José; Jacobsen, Finn

2007-01-01

Digital waveguide mesh (DWM) is a popular method for time domain modelling of sound fields. DWM consists of a recursive digital filter structure where a D'Alembert solution of the wave equation is propagated. One of the attractive characteristics of this method is related to the simplicity...... model of the boundary does not agree with the behaviour of a locally reacting surface, and this can give rise to contradictions in the physical interpretation of the reflected sound field. This paper analyses the behaviour of frequency dependent boundary conditions in DWM in order to obtain a physical...

Analytical solution for beam with time-dependent boundary conditions versus response spectrum

International Nuclear Information System (INIS)

Gou, P.F.; Panahi, K.K.

2001-01-01

This paper studies the responses of a uniform simple beam for which the supports are subjected to time-dependent conditions. Analytical solution in terms of series was presented for two cases: (1) Two supports of a simple beam are subjected to a harmonic motion, and (2) One of the two supports is stationary while the other is subjected to a harmonic motion. The results of the analytical solution were investigated and compared with the results of conventional response spectrum method using the beam finite element model. One of the applications of the results presented in this paper can be used to assess the adequacy and accuracy of the engineering approaches such as response spectra methods. It has been found that, when the excitation frequency equals the fundamental frequency of the beam, the results from response spectrum method are in good agreement with the exact calculation. The effects of initial conditions on the responses are also examined. It seems that the non-zero initial velocity has pronounced effects on the displacement time histories but it has no effect on the maximum accelerations. (author)

Time-dependent embedding

OpenAIRE

Inglesfield, J. E.

2007-01-01

A method of solving the time-dependent Schr\\"odinger equation is presented, in which a finite region of space is treated explicitly, with the boundary conditions for matching the wave-functions on to the rest of the system replaced by an embedding term added on to the Hamiltonian. This time-dependent embedding term is derived from the Fourier transform of the energy-dependent embedding potential, which embeds the time-independent Schr\\"odinger equation. Results are presented for a one-dimensi...

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.

Weak solutions to the full Navier-Stokes-Fourier system with slip boundary conditions in time dependent domains

Czech Academy of Sciences Publication Activity Database

Kreml, Ondřej; Mácha, Václav; Nečasová, Šárka; Wróblewska-Kamińska, A.

2018-01-01

Roč. 109, January (2018), s. 67-92 ISSN 0021-7824 R&D Projects: GA ČR GA13-00522S; GA MŠk(CZ) 7AMB16PL060 Institutional support: RVO:67985840 Keywords : compressible Navier–Stokes–Fourier equations * time-varying domain * slip boundary conditions Subject RIV: BA - General Mathematics OBOR OECD: Pure mathematics Impact factor: 1.802, year: 2016 http://www.sciencedirect.com/science/article/pii/S0021782417301381?via%3Dihub

Weak solutions to the full Navier-Stokes-Fourier system with slip boundary conditions in time dependent domains

Czech Academy of Sciences Publication Activity Database

Kreml, Ondřej; Mácha, Václav; Nečasová, Šárka; Wróblewska-Kamińska, A.

2018-01-01

Roč. 109, January (2018), s. 67-92 ISSN 0021-7824 R&D Projects: GA ČR GA13-00522S; GA MŠk(CZ) 7AMB16PL060 Institutional support: RVO:67985840 Keywords : compressible Navier–Stokes–Fourier equations * time-varying domain * slip boundary conditions Subject RIV: BA - General Mathematics OBOR OECD: Pure mathematics Impact factor: 1.802, year: 2016 http://www. science direct.com/ science /article/pii/S0021782417301381?via%3Dihub

Stability properties of a heat equation with state-dependent parameters and asymmetric boundary conditions

DEFF Research Database (Denmark)

Backi, Christoph Josef; Bendtsen, Jan Dimon; Leth, John

2015-01-01

In this work the stability properties of a partial differential equation (PDE) with state-dependent parameters and asymmetric boundary conditions are investigated. The PDE describes the temperature distribution inside foodstuff, but can also hold for other applications and phenomena. We show...

A Time-dependent Heliospheric Model Driven by Empirical Boundary Conditions

Science.gov (United States)

Kim, T. K.; Arge, C. N.; Pogorelov, N. V.

2017-12-01

Consisting of charged particles originating from the Sun, the solar wind carries the Sun's energy and magnetic field outward through interplanetary space. The solar wind is the predominant source of space weather events, and modeling the solar wind propagation to Earth is a critical component of space weather research. Solar wind models are typically separated into coronal and heliospheric parts to account for the different physical processes and scales characterizing each region. Coronal models are often coupled with heliospheric models to propagate the solar wind out to Earth's orbit and beyond. The Wang-Sheeley-Arge (WSA) model is a semi-empirical coronal model consisting of a potential field source surface model and a current sheet model that takes synoptic magnetograms as input to estimate the magnetic field and solar wind speed at any distance above the coronal region. The current version of the WSA model takes the Air Force Data Assimilative Photospheric Flux Transport (ADAPT) model as input to provide improved time-varying solutions for the ambient solar wind structure. When heliospheric MHD models are coupled with the WSA model, density and temperature at the inner boundary are treated as free parameters that are tuned to optimal values. For example, the WSA-ENLIL model prescribes density and temperature assuming momentum flux and thermal pressure balance across the inner boundary of the ENLIL heliospheric MHD model. We consider an alternative approach of prescribing density and temperature using empirical correlations derived from Ulysses and OMNI data. We use our own modeling software (Multi-scale Fluid-kinetic Simulation Suite) to drive a heliospheric MHD model with ADAPT-WSA input. The modeling results using the two different approaches of density and temperature prescription suggest that the use of empirical correlations may be a more straightforward, consistent method.

Incorporation of exact boundary conditions into a discontinuous galerkin finite element method for accurately solving 2d time-dependent maxwell equations

KAUST Repository

Sirenko, Kostyantyn

2013-01-01

A scheme that discretizes exact absorbing boundary conditions (EACs) to incorporate them into a time-domain discontinuous Galerkin finite element method (TD-DG-FEM) is described. The proposed TD-DG-FEM with EACs is used for accurately characterizing transient electromagnetic wave interactions on two-dimensional waveguides. Numerical results demonstrate the proposed method\\'s superiority over the TD-DG-FEM that employs approximate boundary conditions and perfectly matched layers. Additionally, it is shown that the proposed method can produce the solution with ten-eleven digit accuracy when high-order spatial basis functions are used to discretize the Maxwell equations as well as the EACs. © 1963-2012 IEEE.

Radiating star with a time-dependent Karmarkar condition

Energy Technology Data Exchange (ETDEWEB)

Naidu, Nolene Ferrari; Maharaj, Sunil D. [University of KwaZulu-Natal, Astrophysics and Cosmology Research Unit, School of Mathematics, Statistics and Computer Science, Durban (South Africa); Govender, Megandhren [Durban University of Technology, Department of Mathematics, Faculty of Applied Sciences, Durban (South Africa)

2018-01-15

In this paper we employ the Karmarkar condition (Proc Indian Acad Sci A 27:56, 1948) to model a spherically symmetric radiating star undergoing dissipative gravitational collapse in the form of a radial heat flux. A particular solution of the boundary condition renders the Karmarkar condition independent of time which allows us to fully specify the spatial behaviour of the gravitational potentials. The interior solution is smoothly matched to Vaidya's outgoing solution across a time-like hypersurface which yields the temporal behaviour of the model. Physical analysis of the matter and thermodynamical variables show that the model is well-behaved. (orig.)

Incorporation of exact boundary conditions into a discontinuous galerkin finite element method for accurately solving 2d time-dependent maxwell equations

KAUST Repository

Sirenko, Kostyantyn; Liu, Meilin; Bagci, Hakan

2013-01-01

A scheme that discretizes exact absorbing boundary conditions (EACs) to incorporate them into a time-domain discontinuous Galerkin finite element method (TD-DG-FEM) is described. The proposed TD-DG-FEM with EACs is used for accurately characterizing

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

Asymptotics of the Eigenvalues of a Self-Adjoint Fourth Order Boundary Value Problem with Four Eigenvalue Parameter Dependent Boundary Conditions

Directory of Open Access Journals (Sweden)

Manfred Möller

2013-01-01

Full Text Available Considered is a regular fourth order ordinary differential equation which depends quadratically on the eigenvalue parameter λ and which has separable boundary conditions depending linearly on λ. It is shown that the eigenvalues lie in the closed upper half plane or on the imaginary axis and are symmetric with respect to the imaginary axis. The first four terms in the asymptotic expansion of the eigenvalues are provided.

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

Directory of Open Access Journals (Sweden)

Cornelis van der Mee

2005-01-01

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

Stokes system with solution-dependent threshold slip boundary conditions: Analysis, approximation and implementation

Czech Academy of Sciences Publication Activity Database

Haslinger, Jaroslav; Kučera, R.; Šátek, V.

2017-01-01

Roč. 22, October 2017 (2017), s. 1-14 ISSN 1081-2865 R&D Projects: GA MŠk LQ1602; GA ČR(CZ) GA17-01747S Institutional support: RVO:68145535 Keywords : Stokes system * threshold slip boundary conditions * solution dependent slip function Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 2.953, year: 2016 http://journals.sagepub.com/doi/abs/10.1177/1081286517716222

Stokes system with solution-dependent threshold slip boundary conditions: Analysis, approximation and implementation

Czech Academy of Sciences Publication Activity Database

Haslinger, Jaroslav; Kučera, R.; Šátek, V.

2017-01-01

Roč. 22, October 2017 (2017), s. 1-14 ISSN 1081-2865 R&D Projects: GA MŠk LQ1602; GA ČR(CZ) GA17-01747S Institutional support: RVO:68145535 Keywords : Stokes system * threshold slip boundary conditions * solution dependent slip function Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 2.953, year: 2016 http:// journals .sagepub.com/doi/abs/10.1177/1081286517716222

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.

Artificial Boundary Conditions for the Numerical Simulation of Unsteady Acoustic Waves

National Research Council Canada - National Science Library

Tsynkov, S. V

2003-01-01

We construct non-local artificial boundary conditions (ABCs) for the numerical simulation of genuinely time-dependent acoustic waves that propagate from a compact source in an unbounded unobstructed space...

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.

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.

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

Boundary-integral equation formulation for time-dependent inelastic deformation in metals

Energy Technology Data Exchange (ETDEWEB)

Kumar, V; Mukherjee, S

1977-01-01

The mathematical structure of various constitutive relations proposed in recent years for representing time-dependent inelastic deformation behavior of metals at elevated temperatues has certain features which permit a simple formulation of the three-dimensional inelasticity problem in terms of real time rates. A direct formulation of the boundary-integral equation method in terms of rates is discussed for the analysis of time-dependent inelastic deformation of arbitrarily shaped three-dimensional metallic bodies subjected to arbitrary mechanical and thermal loading histories and obeying constitutive relations of the kind mentioned above. The formulation is based on the assumption of infinitesimal deformations. Several illustrative examples involving creep of thick-walled spheres, long thick-walled cylinders, and rotating discs are discussed. The implementation of the method appears to be far easier than analogous BIE formulations that have been suggested for elastoplastic problems.

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.

Hybrid Fourier pseudospectral/discontinuous Galerkin time-domain method for arbitrary boundary conditions

NARCIS (Netherlands)

Pagan Munoz, R.; Hornikx, M.C.J.

The wave-based Fourier Pseudospectral time-domain (Fourier-PSTD) method was shown to be an effective way of modeling outdoor acoustic propagation problems as described by the linearized Euler equations (LEE), but is limited to real-valued frequency independent boundary conditions and predominantly

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.

Stabilization of time domain acoustic boundary element method for the interior problem with impedance boundary conditions.

Science.gov (United States)

Jang, Hae-Won; Ih, Jeong-Guon

2012-04-01

The time domain boundary element method (BEM) is associated with numerical instability that typically stems from the time marching scheme. In this work, a formulation of time domain BEM is derived to deal with all types of boundary conditions adopting a multi-input, multi-output, infinite impulse response structure. The fitted frequency domain impedance data are converted into a time domain expression as a form of an infinite impulse response filter, which can also invoke a modeling error. In the calculation, the response at each time step is projected onto the wave vector space of natural radiation modes, which can be obtained from the eigensolutions of the single iterative matrix. To stabilize the computation, unstable oscillatory modes are nullified, and the same decay rate is used for two nonoscillatory modes. As a test example, a transient sound field within a partially lined, parallelepiped box is used, within which a point source is excited by an octave band impulse. In comparison with the results of the inverse Fourier transform of a frequency domain BEM, the average of relative difference norm in the stabilized time response is found to be 4.4%.

Implementation of Unsplit Perfectly Matched Layer Absorbing Boundary Condition in 3 Dimensional Finite Difference Time Domain Method

Directory of Open Access Journals (Sweden)

B. U. Musa

2017-04-01

Full Text Available The C++ programming language was used to implement three-dimensional (3-D finite-difference time-domain (FDTD technique to simulate radiation of high frequency electromagnetic waves in free space. To achieve any meaningful results the computational domain of interest should have to be truncated in some way and this is achieved by applying absorbing boundary conditions. A uniaxial perfectly matched layer (UPML absorbing boundary condition is used in this work. The discretised equations of the UPML in FDTD time stepping scheme were derived and has been successfully implemented using the computer program. Simulation results showed that the UPML behaves as an absorber. This was confirmed by comparing the results with another boundary condition, the Mur ABC.

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.

Accurate characterization of 3D diffraction gratings using time domain discontinuous Galerkin method with exact absorbing boundary conditions

KAUST Repository

Sirenko, Kostyantyn

2013-07-01

Exact absorbing and periodic boundary conditions allow to truncate grating problems\\' infinite physical domains without introducing any errors. This work presents exact absorbing boundary conditions for 3D diffraction gratings and describes their discretization within a high-order time-domain discontinuous Galerkin finite element method (TD-DG-FEM). The error introduced by the boundary condition discretization matches that of the TD-DG-FEM; this results in an optimal solver in terms of accuracy and computation time. Numerical results demonstrate the superiority of this solver over TD-DG-FEM with perfectly matched layers (PML)-based domain truncation. © 2013 IEEE.

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

Discontinuous Sturm-Liouville Problems with Eigenvalue Dependent Boundary Condition

Energy Technology Data Exchange (ETDEWEB)

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

2014-12-15

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

Aspect of ECMWF downscaled Regional Climate Modeling in simulating Indian summer monsoon rainfall and dependencies on lateral boundary conditions

Science.gov (United States)

Ghosh, Soumik; Bhatla, R.; Mall, R. K.; Srivastava, Prashant K.; Sahai, A. K.

2018-03-01

Climate model faces considerable difficulties in simulating the rainfall characteristics of southwest summer monsoon. In this study, the dynamical downscaling of European Centre for Medium-Range Weather Forecast's (ECMWF's) ERA-Interim (EIN15) has been utilized for the simulation of Indian summer monsoon (ISM) through the Regional Climate Model version 4.3 (RegCM-4.3) over the South Asia Co-Ordinated Regional Climate Downscaling EXperiment (CORDEX) domain. The complexities of model simulation over a particular terrain are generally influenced by factors such as complex topography, coastal boundary, and lack of unbiased initial and lateral boundary conditions. In order to overcome some of these limitations, the RegCM-4.3 is employed for simulating the rainfall characteristics over the complex topographical conditions. For reliable rainfall simulation, implementations of numerous lower boundary conditions are forced in the RegCM-4.3 with specific horizontal grid resolution of 50 km over South Asia CORDEX domain. The analysis is considered for 30 years of climatological simulation of rainfall, outgoing longwave radiation (OLR), mean sea level pressure (MSLP), and wind with different vertical levels over the specified region. The dependency of model simulation with the forcing of EIN15 initial and lateral boundary conditions is used to understand the impact of simulated rainfall characteristics during different phases of summer monsoon. The results obtained from this study are used to evaluate the activity of initial conditions of zonal wind circulation speed, which causes an increase in the uncertainty of regional model output over the region under investigation. Further, the results showed that the EIN15 zonal wind circulation lacks sufficient speed over the specified region in a particular time, which was carried forward by the RegCM output and leads to a disrupted regional simulation in the climate model.

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.

Effect of real-time boundary wind conditions on the air flow and pollutant dispersion in an urban street canyon—Large eddy simulations

Science.gov (United States)

Zhang, Yun-Wei; Gu, Zhao-Lin; Cheng, Yan; Lee, Shun-Cheng

2011-07-01

Air flow and pollutant dispersion characteristics in an urban street canyon are studied under the real-time boundary conditions. A new scheme for realizing real-time boundary conditions in simulations is proposed, to keep the upper boundary wind conditions consistent with the measured time series of wind data. The air flow structure and its evolution under real-time boundary wind conditions are simulated by using this new scheme. The induced effect of time series of ambient wind conditions on the flow structures inside and above the street canyon is investigated. The flow shows an obvious intermittent feature in the street canyon and the flapping of the shear layer forms near the roof layer under real-time wind conditions, resulting in the expansion or compression of the air mass in the canyon. The simulations of pollutant dispersion show that the pollutants inside and above the street canyon are transported by different dispersion mechanisms, relying on the time series of air flow structures. Large scale air movements in the processes of the air mass expansion or compression in the canyon exhibit obvious effects on pollutant dispersion. The simulations of pollutant dispersion also show that the transport of pollutants from the canyon to the upper air flow is dominated by the shear layer turbulence near the roof level and the expansion or compression of the air mass in street canyon under real-time boundary wind conditions. Especially, the expansion of the air mass, which features the large scale air movement of the air mass, makes more contribution to the pollutant dispersion in this study. Comparisons of simulated results under different boundary wind conditions indicate that real-time boundary wind conditions produces better condition for pollutant dispersion than the artificially-designed steady boundary wind conditions.

Modelling of pressurized water reactor fuel, rod time dependent radial heat flow with boundary element method; Modeliranje spremenljivega radijalnega toplotnega toka tlacnovodne gorivne palice z metodo robnih elementov

Energy Technology Data Exchange (ETDEWEB)

Sarler, B [Institut Jozef Stefan, Ljubljana (Yugoslavia)

1987-07-01

The basic principles of the boundary element method numerical treatment of the radial flow heat diffusion equation are presented. The algorithm copes the time dependent Dirichlet and Neumann boundary conditions, temperature dependent material properties and regions from different materials in thermal contact. It is verified on the several analytically obtained test cases. The developed method is used for the modelling of unsteady radial heat flow in pressurized water reactor fuel rod. (author)

Numerical method for time-dependent localized corrosion analysis with moving boundaries by combining the finite volume method and voxel method

International Nuclear Information System (INIS)

Onishi, Yuki; Takiyasu, Jumpei; Amaya, Kenji; Yakuwa, Hiroshi; Hayabusa, Keisuke

2012-01-01

Highlights: ► A novel numerical method to analyze time dependent localized corrosion is developed. ► It takes electromigration, mass diffusion, chemical reactions, and moving boundaries. ► Our method perfectly satisfies the conservation of mass and electroneutrality. ► The behavior of typical crevice corrosion is successfully simulated. ► Both verification and validation of our method are carried out. - Abstract: A novel numerical method for time-dependent localized corrosion analysis is presented. Electromigration, mass diffusion, chemical reactions, and moving boundaries are considered in the numerical simulation of localized corrosion of engineering alloys in an underwater environment. Our method combines the finite volume method (FVM) and the voxel method. The FVM is adopted in the corrosion rate calculation so that the conservation of mass is satisfied. A newly developed decoupled algorithm with a projection method is introduced in the FVM to decouple the multiphysics problem into the electrostatic, mass transport, and chemical reaction analyses with electroneutrality maintained. The polarization curves for the corroding metal are used as boundary conditions for the metal surfaces to calculate the corrosion rates. The voxel method is adopted in updating the moving boundaries of cavities without remeshing and mesh-to-mesh solution mapping. Some modifications of the standard voxel method, which represents the boundaries as zigzag-shaped surfaces, are introduced to generate smooth surfaces. Our method successfully reproduces the numerical and experimental results of a capillary electrophoresis problem. Furthermore, the numerical results are qualitatively consistent with the experimental results for several examples of crevice corrosion.

Reflected forward-backward SDEs and obstacle problems with boundary conditions

Directory of Open Access Journals (Sweden)

Jin Ma

2001-01-01

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

Time-dependent free boundary equilibrium and resistive diffusion in a tokamak plasma

International Nuclear Information System (INIS)

Selig, G.

2012-12-01

In a Tokamak, in order to create the necessary conditions for nuclear fusion to occur, a plasma is maintained by applying magnetic fields. Under the hypothesis of an axial symmetry of the tokamak, the study of the magnetic configuration at equilibrium is done in two dimensions, and is deduced from the poloidal flux function. This function is solution of a non linear partial differential equation system, known as equilibrium problem. This thesis presents the time dependent free boundary equilibrium problem, where the circuit equations in the tokamak coils and passive conductors are solved together with the Grad-Shafranov equation to produce a dynamic simulation of the plasma. In this framework, the Finite Element equilibrium code CEDRES has been improved in order to solve the aforementioned dynamic problem. Consistency tests and comparisons with the DINA-CH code on an ITER vertical instability case have validated the results. Then, the resistive diffusion of the plasma current density has been simulated using a coupling between CEDRES and the averaged one-dimensional diffusion equation, and it has been successfully compared with the integrated modeling code CRONOS. (author)

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.

Element stacking method for topology optimization with material-dependent boundary and loading conditions

DEFF Research Database (Denmark)

Yoon, Gil Ho; Park, Y.K.; Kim, Y.Y.

2007-01-01

A new topology optimization scheme, called the element stacking method, is developed to better handle design optimization involving material-dependent boundary conditions and selection of elements of different types. If these problems are solved by existing standard approaches, complicated finite...... element models or topology optimization reformulation may be necessary. The key idea of the proposed method is to stack multiple elements on the same discretization pixel and select a single or no element. In this method, stacked elements on the same pixel have the same coordinates but may have...... independent degrees of freedom. Some test problems are considered to check the effectiveness of the proposed stacking method....

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

Bayesian Inference for Linear Parabolic PDEs with Noisy Boundary Conditions

KAUST Repository

Ruggeri, Fabrizio; Sawlan, Zaid A; Scavino, Marco; Tempone, Raul

2016-01-01

In this work we develop a hierarchical Bayesian setting to infer unknown parameters in initial-boundary value problems (IBVPs) for one-dimensional linear parabolic partial differential equations. Noisy boundary data and known initial condition are assumed. We derive the likelihood function associated with the forward problem, given some measurements of the solution field subject to Gaussian noise. Such function is then analytically marginalized using the linearity of the equation. Gaussian priors have been assumed for the time-dependent Dirichlet boundary values. Our approach is applied to synthetic data for the one-dimensional heat equation model, where the thermal diffusivity is the unknown parameter. We show how to infer the thermal diffusivity parameter when its prior distribution is lognormal or modeled by means of a space-dependent stationary lognormal random field. We use the Laplace method to provide approximated Gaussian posterior distributions for the thermal diffusivity. Expected information gains and predictive posterior densities for observable quantities are numerically estimated for different experimental setups.

Bayesian Inference for Linear Parabolic PDEs with Noisy Boundary Conditions

KAUST Repository

Ruggeri, Fabrizio

2015-01-07

In this work we develop a hierarchical Bayesian setting to infer unknown parameters in initial-boundary value problems (IBVPs) for one-dimensional linear parabolic partial differential equations. Noisy boundary data and known initial condition are assumed. We derive the likelihood function associated with the forward problem, given some measurements of the solution field subject to Gaussian noise. Such function is then analytically marginalized using the linearity of the equation. Gaussian priors have been assumed for the time-dependent Dirichlet boundary values. Our approach is applied to synthetic data for the one-dimensional heat equation model, where the thermal diffusivity is the unknown parameter. We show how to infer the thermal diffusivity parameter when its prior distribution is lognormal or modeled by means of a space-dependent stationary lognormal random field. We use the Laplace method to provide approximated Gaussian posterior distributions for the thermal diffusivity. Expected information gains and predictive posterior densities for observable quantities are numerically estimated for different experimental setups.

Bayesian Inference for Linear Parabolic PDEs with Noisy Boundary Conditions

KAUST Repository

Ruggeri, Fabrizio

2016-01-06

In this work we develop a hierarchical Bayesian setting to infer unknown parameters in initial-boundary value problems (IBVPs) for one-dimensional linear parabolic partial differential equations. Noisy boundary data and known initial condition are assumed. We derive the likelihood function associated with the forward problem, given some measurements of the solution field subject to Gaussian noise. Such function is then analytically marginalized using the linearity of the equation. Gaussian priors have been assumed for the time-dependent Dirichlet boundary values. Our approach is applied to synthetic data for the one-dimensional heat equation model, where the thermal diffusivity is the unknown parameter. We show how to infer the thermal diffusivity parameter when its prior distribution is lognormal or modeled by means of a space-dependent stationary lognormal random field. We use the Laplace method to provide approximated Gaussian posterior distributions for the thermal diffusivity. Expected information gains and predictive posterior densities for observable quantities are numerically estimated for different experimental setups.

Subgeometric Ergodicity Analysis of Continuous-Time Markov Chains under Random-Time State-Dependent Lyapunov Drift Conditions

Directory of Open Access Journals (Sweden)

Mokaedi V. Lekgari

2014-01-01

Full Text Available We investigate random-time state-dependent Foster-Lyapunov analysis on subgeometric rate ergodicity of continuous-time Markov chains (CTMCs. We are mainly concerned with making use of the available results on deterministic state-dependent drift conditions for CTMCs and on random-time state-dependent drift conditions for discrete-time Markov chains and transferring them to CTMCs.

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

Temperature dependence of grain boundary free energy and elastic constants

International Nuclear Information System (INIS)

Foiles, Stephen M.

2010-01-01

This work explores the suggestion that the temperature dependence of the grain boundary free energy can be estimated from the temperature dependence of the elastic constants. The temperature-dependent elastic constants and free energy of a symmetric Σ79 tilt boundary are computed for an embedded atom method model of Ni. The grain boundary free energy scales with the product of the shear modulus times the lattice constant for temperatures up to about 0.75 the melting temperature.

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.

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

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.

Finite element time domain modeling of controlled-Source electromagnetic data with a hybrid boundary condition

DEFF Research Database (Denmark)

Cai, Hongzhu; Hu, Xiangyun; Xiong, Bin

2017-01-01

method which is unconditionally stable. We solve the diffusion equation for the electric field with a total field formulation. The finite element system of equation is solved using the direct method. The solutions of electric field, at different time, can be obtained using the effective time stepping...... method with trivial computation cost once the matrix is factorized. We try to keep the same time step size for a fixed number of steps using an adaptive time step doubling (ATSD) method. The finite element modeling domain is also truncated using a semi-adaptive method. We proposed a new boundary...... condition based on approximating the total field on the modeling boundary using the primary field corresponding to a layered background model. We validate our algorithm using several synthetic model studies....

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.

Storm-time ring current: model-dependent results

Directory of Open Access Journals (Sweden)

N. Yu. Ganushkina

2012-01-01

Full Text Available The main point of the paper is to investigate how much the modeled ring current depends on the representations of magnetic and electric fields and boundary conditions used in simulations. Two storm events, one moderate (SymH minimum of −120 nT on 6–7 November 1997 and one intense (SymH minimum of −230 nT on 21–22 October 1999, are modeled. A rather simple ring current model is employed, namely, the Inner Magnetosphere Particle Transport and Acceleration model (IMPTAM, in order to make the results most evident. Four different magnetic field and two electric field representations and four boundary conditions are used. We find that different combinations of the magnetic and electric field configurations and boundary conditions result in very different modeled ring current, and, therefore, the physical conclusions based on simulation results can differ significantly. A time-dependent boundary outside of 6.6 RE gives a possibility to take into account the particles in the transition region (between dipole and stretched field lines forming partial ring current and near-Earth tail current in that region. Calculating the model SymH* by Biot-Savart's law instead of the widely used Dessler-Parker-Sckopke (DPS relation gives larger and more realistic values, since the currents are calculated in the regions with nondipolar magnetic field. Therefore, the boundary location and the method of SymH* calculation are of key importance for ring current data-model comparisons to be correctly interpreted.

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

OpenAIRE

Li, Zhiyuan; Huang, Xinchi; Yamamoto, Masahiro

2018-01-01

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

Modeling Space-Time Dependent Helium Bubble Evolution in Tungsten Armor under IFE Conditions

International Nuclear Information System (INIS)

Qiyang Hu; Shahram Sharafat; Nasr Ghoniem

2006-01-01

The High Average Power Laser (HAPL) program is a coordinated effort to develop Laser Inertial Fusion Energy. The implosion of the D-T target produces a spectrum of neutrons, X-rays, and charged particles, which arrive at the first wall (FW) at different times within about 2.5 μs at a frequency of 5 to 10 Hz. Helium is one of several high-energy charged particle constituents impinging on the candidate tungsten armored low activation ferritic steel First Wall. The spread of the implanted debris and burn helium energies results in a unique space-time dependent implantation profile that spans about 10 μm in tungsten. Co-implantation of X-rays and other ions results in spatially dependent damage profiles and rapid space-time dependent temperature spikes and gradients. The rate of helium transport and helium bubble formation will vary significantly throughout the implanted region. Furthermore, helium will also be transported via the migration of helium bubbles and non-equilibrium helium-vacancy clusters. The HEROS code was developed at UCLA to model the spatial and time-dependent helium bubble nucleation, growth, coalescence, and migration under transient damage rates and transient temperature gradients. The HEROS code is based on kinetic rate theory, which includes clustering of helium and vacancies, helium mobility, helium-vacancy cluster stability, cavity nucleation and growth and other microstructural features such as interstitial loop evolution, grain boundaries, and precipitates. The HEROS code is based on space-time discretization of reaction-diffusion type equations to account for migration of mobile species between neighboring bins as single atoms, clusters, or bubbles. HAPL chamber FW implantation conditions are used to model helium bubble evolution in the implanted tungsten. Helium recycling rate predictions are compared with experimental results of helium ion implantation experiments. (author)

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

Finite-difference time-domain modeling of curved material interfaces by using boundary condition equations method

International Nuclear Information System (INIS)

Lu Jia; Zhou Huaichun

2016-01-01

To deal with the staircase approximation problem in the standard finite-difference time-domain (FDTD) simulation, the two-dimensional boundary condition equations (BCE) method is proposed in this paper. In the BCE method, the standard FDTD algorithm can be used as usual, and the curved surface is treated by adding the boundary condition equations. Thus, while maintaining the simplicity and computational efficiency of the standard FDTD algorithm, the BCE method can solve the staircase approximation problem. The BCE method is validated by analyzing near field and far field scattering properties of the PEC and dielectric cylinders. The results show that the BCE method can maintain a second-order accuracy by eliminating the staircase approximation errors. Moreover, the results of the BCE method show good accuracy for cylinder scattering cases with different permittivities. (paper)

Eulerian and Lagrangian views of a turbulent boundary layer flow using time-resolved tomographic PIV

NARCIS (Netherlands)

Schröder, A.; Geisler, R.; Staack, K.; Elsinga, G.E.; Scarano, F.; Wieneke, B.; Henning, A.; Poelma, C.; Westerweel, J.

2010-01-01

Coherent structures and their time evolution in the logarithmic region of a turbulent boundary layer investigated by means of 3D space–time correlations and time-dependent conditional averaging techniques are the focuses of the present paper. Experiments have been performed in the water tunnel at TU

Brain response to prosodic boundary cues depends on boundary position

Directory of Open Access Journals (Sweden)

Julia eHolzgrefe

2013-07-01

Full Text Available Prosodic information is crucial for spoken language comprehension and especially for syntactic parsing, because prosodic cues guide the hearer’s syntactic analysis. The time course and mechanisms of this interplay of prosody and syntax are not yet well understood. In particular, there is an ongoing debate whether local prosodic cues are taken into account automatically or whether they are processed in relation to the global prosodic context in which they appear. The present study explores whether the perception of a prosodic boundary is affected by its position within an utterance. In an event-related potential (ERP study we tested if the brain response evoked by the prosodic boundary differs when the boundary occurs early in a list of three names connected by conjunctions (i.e., after the first name as compared to later in the utterance (i.e., after the second name. A closure positive shift (CPS — marking the processing of a prosodic phrase boundary — was elicited only for stimuli with a late boundary, but not for stimuli with an early boundary. This result is further evidence for an immediate integration of prosodic information into the parsing of an utterance. In addition, it shows that the processing of prosodic boundary cues depends on the previously processed information from the preceding prosodic context.

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.

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.

Pomarning-eddington approximation for time-dependent radiation transfer in finite slab media

International Nuclear Information System (INIS)

El-Wakil, S.A.; Degheidy, A.R.; Sallah, M.

2005-01-01

The time-dependent monoenergetic radiation transfer equation with linear anisotropic scattering is proposed. Pomraning-Eddington approximation is used to calculate the radiation intensity in finite plane-parallel media. Numerical results are done for the isotropic media. Shielding calculations are shown for reflectivity and transmissivity at different times. The medium is assumed to have specular-reflecting boundaries. Two different weight functions are introduced to force the boundary conditions to be fulfilled

Response moderation models for conditional dependence between response time and response accuracy.

Science.gov (United States)

Bolsinova, Maria; Tijmstra, Jesper; Molenaar, Dylan

2017-05-01

It is becoming more feasible and common to register response times in the application of psychometric tests. Researchers thus have the opportunity to jointly model response accuracy and response time, which provides users with more relevant information. The most common choice is to use the hierarchical model (van der Linden, 2007, Psychometrika, 72, 287), which assumes conditional independence between response time and accuracy, given a person's speed and ability. However, this assumption may be violated in practice if, for example, persons vary their speed or differ in their response strategies, leading to conditional dependence between response time and accuracy and confounding measurement. We propose six nested hierarchical models for response time and accuracy that allow for conditional dependence, and discuss their relationship to existing models. Unlike existing approaches, the proposed hierarchical models allow for various forms of conditional dependence in the model and allow the effect of continuous residual response time on response accuracy to be item-specific, person-specific, or both. Estimation procedures for the models are proposed, as well as two information criteria that can be used for model selection. Parameter recovery and usefulness of the information criteria are investigated using simulation, indicating that the procedure works well and is likely to select the appropriate model. Two empirical applications are discussed to illustrate the different types of conditional dependence that may occur in practice and how these can be captured using the proposed hierarchical models. © 2016 The British Psychological Society.

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.

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

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.

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

Modelling of Edge Insulation Depending on Boundary Conditions for the Ground Level

Science.gov (United States)

Stolarska, Agata; Strzałkowski, Jarosław

2017-10-01

The article presents results of CFD software aided simulations of a thermal bridge, existing at the wall-slab on ground connection. Calculations were made for different variants of the edge insulation location. Schemes without any edge insulation, with some vertical insulation, horizontal, diagonal, and diagonal combined with insulation used as formwork under the slab on ground were analysed. Each variant was differentiated with boundary conditions for the ground. Vertical borders of the model in the ground, as well as the lower border were described in the first solution as adiabatic, while in the second case, a variable temperature value, depending on the ground depth, was set. For comparison, additional calculations were conducted for non-stationary conditions, in which the initial temperature of the ground was set to the average annual temperature of air. The calculations were based on the location of Szczecin, for which the outside air temperature was set to -16.0°C. Results obtained from the simulation were then used to determine the thermal bridge parameters, in particular, thermal coupling coefficient and linear thermal transmittance. The effect of the set of boundary conditions is clearly seen. In general, for all the five variants, lower values of heat fluxes and linear thermal transmittances were obtained, when variable temperature in the ground was assumed. From the point of view of energy balance, it is more favourable to use the values of ψg obtained when the ground temperature is taken into account. The data breakdown shows that application of the actual temperature distribution in the ground to a model has a strong effect on distribution of the 0.0°C isotherm. The adiabatic model indicates that the ground under the slab freezes, while the model, which takes into account the temperature of the ground, shows that the ground under the floor has positive temperatures and the 0.0°C isotherm reaches only the edge of the outer wall. Moreover, the

Solving Two -Dimensional Diffusion Equations with Nonlocal Boundary Conditions by a Special Class of Padé Approximants

Directory of Open Access Journals (Sweden)

Mohammad Siddique

2010-08-01

Full Text Available Parabolic partial differential equations with nonlocal boundary conditions arise in modeling of a wide range of important application areas such as chemical diffusion, thermoelasticity, heat conduction process, control theory and medicine science. In this paper, we present the implementation of positivity- preserving Padé numerical schemes to the two-dimensional diffusion equation with nonlocal time dependent boundary condition. We successfully implemented these numerical schemes for both Homogeneous and Inhomogeneous cases. The numerical results show that these Padé approximation based numerical schemes are quite accurate and easily implemented.

Cellular Neural Network Method for Critical Slab with Albedo Boundary Condition

International Nuclear Information System (INIS)

Pirouzmanda, A.; Hadada, K.; Suh, K. Y.

2010-01-01

The neutron transport problems have been studied theoretically and numerically for years. A number of researchers have studied the criticality problems of one-speed neutrons in homogeneous slabs and spheres using various methods. The Chebyshev polynomial approximation method (T N method) has lately been developed and improved for the neutron transport equation in slab geometry. The one-speed time-dependent neutron transport equation using the Cellular Neural Network (CNN) for the vacuum boundary condition has previously been solved. In this paper, we demonstrate the capacity of CNN in calculating the critical slab thickness for different boundary conditions and its variation with moments N. The architecture of the CNN has already been dealt with thoroughly. Essentially, the CNN is used to model a first-order system of the partial differential equations (PDEs). The original equations in the T N approximation are also a set of PDEs. The CNN approach lends itself to analog VLSI implementation. In this study, the CNN model is implemented using the HSpice software package

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.

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.

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

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.

Numerical Solution of Time-Dependent Problems with a Fractional-Power Elliptic Operator

Science.gov (United States)

Vabishchevich, P. N.

2018-03-01

A time-dependent problem in a bounded domain for a fractional diffusion equation is considered. The first-order evolution equation involves a fractional-power second-order elliptic operator with Robin boundary conditions. A finite-element spatial approximation with an additive approximation of the operator of the problem is used. The time approximation is based on a vector scheme. The transition to a new time level is ensured by solving a sequence of standard elliptic boundary value problems. Numerical results obtained for a two-dimensional model problem are presented.

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.

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.

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)

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

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.

Time reversal method with stabilizing boundary conditions for Photoacoustic tomography

International Nuclear Information System (INIS)

Chervova, Olga; Oksanen, Lauri

2016-01-01

We study an inverse initial source problem that models photoacoustic tomography measurements with array detectors, and introduce a method that can be viewed as a modification of the so called back and forth nudging method. We show that the method converges at an exponential rate under a natural visibility condition, with data given only on a part of the boundary of the domain of wave propagation. In this paper we consider the case of noiseless measurements. (paper)

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)

Generic short-time propagation of sharp-boundaries wave packets

Science.gov (United States)

Granot, E.; Marchewka, A.

2005-11-01

A general solution to the "shutter" problem is presented. The propagation of an arbitrary initially bounded wave function is investigated, and the general solution for any such function is formulated. It is shown that the exact solution can be written as an expression that depends only on the values of the function (and its derivatives) at the boundaries. In particular, it is shown that at short times (t << 2mx2/hbar, where x is the distance to the boundaries) the wave function propagation depends only on the wave function's values (or its derivatives) at the boundaries of the region. Finally, we generalize these findings to a non-singular wave function (i.e., for wave packets with finite-width boundaries) and suggest an experimental verification.

Time evolution of some quantum-mechanical systems. Wavefunction cloning in evolving rotating systems. Finite range boundary conditions for time dependent Schroedinger Equation; Evolution temporelle de quelques systemes quantiques. Le clonage de la fonction d`onde dans l`evolution au cours du temps de systemes tournants. Formulation de conditions aux limites a distance finie pour l`equation de Schroedinger dependante du temps

Energy Technology Data Exchange (ETDEWEB)

Arvieu, R.; Carbonell, J.; Gignoux, C.; Mangin-Brinet, M. [Inst. des Sciences Nucleaires, Grenoble-1 Univ., 38 (France); Rozmej, P. [Uniwersytet Marii Curie-Sklodowskiej, Lublin (Poland)

1997-12-31

The time evolution of coherent rotational wave packets associated to a diatomic molecule or to a deformed nucleus has been studied. Assuming a rigid body dynamics the J(J+1) law leads to a mechanism of cloning: the way function is divided into wave packets identical to the initial one at specific time. Applications are studied for a nuclear wave packed formed by Coulomb excitation. Exact boundary conditions at finite distance for the solution of the time-dependent Schroedinger equation are derived. A numerical scheme based on Crank-Nicholson method is proposed to illustrate its applicability in several examples. (authors) 3 refs.

Simulation of time-dependent free-surface Navier-Stokes flows

International Nuclear Information System (INIS)

Muldowney, G.P.

1989-01-01

Two numerical methods for simulation of time-dependent free-surface Navier-Stokes flows are developed. Both techniques are based on semi-implicit time advancement of the momentum equations, integral formulation of the spatial problem at each timestep, and spectral-element discretization to solve the resulting integral equation. Central to each algorithm is a boundary-specific solution step which permits the spatial treatment in two dimensions to be performed in O(N 3 ) operations per timestep despite the presence of deforming geometry. The first approach is a domain-integral formulation involving integrals over the entire flow domain of kernel functions which arise in time-differencing the Navier-Stokes equations. The second is a particular-solution formulation which replaces domain integration with an iterative scheme to generate particular velocity and pressure fields on individual elements, followed by a patching step to produce a particular solution continuous over the full domain. Two of the most difficult aspects of viscous free-surface flow simulations, namely time-dependent geometry and nontrivial boundary conditions, are well accommodated by these integral equation techniques. In addition the methods offer spectral accuracy in space and admit arbitrarily high-order discretization in time. For large-scale computations and/or long-term time advancement the domain-integral algorithm must be executed on a supercomputer to deliver results in reasonable processing time. A detailed simulation of gas liquid flow with full resolution of the free phase boundary requires approximately five CPU hours at 80 megaflops

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

CFD Modeling of Non-Neutral Atmospheric Boundary Layer Conditions

DEFF Research Database (Denmark)

Koblitz, Tilman

model results. A method is developed how to simulate the time-dependant non-neutral ABL flow over complex terrain: a precursor simulation is used to specify unsteady inlet boundary conditions on complex terrain domains. The advantage of the developed RANS model framework is its general applicability...... characteristics of neutral and non-neutral ABL flow. The developed ABL model significantly improves the predicted flow fields over both flat and complex terrain, when compared against neutral models and measurements....... cost than e.g. using large-eddy simulations. The developed ABL model is successfully validated using a range of different test cases with increasing complexity. Data from several large scale field campaigns, wind tunnel experiments, and previous numerical simulations is presented and compared against...

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.

Ensemble-marginalized Kalman filter for linear time-dependent PDEs with noisy boundary conditions: Application to heat transfer in building walls

KAUST Repository

Iglesias, Marco

2017-11-26

In this work, we present the ensemble-marginalized Kalman filter (EnMKF), a sequential algorithm analogous to our previously proposed approach [1,2], for estimating the state and parameters of linear parabolic partial differential equations in initial-boundary value problems when the boundary data are noisy. We apply EnMKF to infer the thermal properties of building walls and to estimate the corresponding heat flux from real and synthetic data. Compared with a modified Ensemble Kalman Filter (EnKF) that is not marginalized, EnMKF reduces the bias error, avoids the collapse of the ensemble without needing to add inflation, and converges to the mean field posterior using $50\\\\%$ or less of the ensemble size required by EnKF. According to our results, the marginalization technique in EnMKF is key to performance improvement with smaller ensembles at any fixed time.

Effect of the time spent by the photon in the absorbed state on the time-dependent transfer of radiation

International Nuclear Information System (INIS)

Rao, D.M.; Rangarajan, K.E.; Peraiah, A.

1990-01-01

The time-dependent transfer equation is derived for a two-level atomic model which takes both bound-bound and bound-free transitions into account. A numerical scheme is proposed for solving the monochromatic time-dependent transfer equation when the time spent by the photon in the absorbed state is significant. The method can be easily extended to solve the problem of time-dependent line formation of the bound-free continuum. It is used here to study three types of boundary conditions of the incident radiation incident on a scattering atmosphere. The quantitative results show that the relaxation of the radiation field depends on the optical depth of the medium and on the ray's angle of emergence. 21 refs

Time-dependent behavior of rough discontinuities under shearing conditions

Science.gov (United States)

Wang, Zhen; Shen, Mingrong; Ding, Wenqi; Jang, Boan; Zhang, Qingzhao

2018-02-01

The mechanical properties of rocks are generally controlled by their discontinuities. In this study, the time-dependent behavior of rough artificial joints under shearing conditions was investigated. Based on Barton’s standard profile lines, samples with artificial joint surfaces were prepared and used to conduct the shear and creep tests. The test results showed that the shear strength of discontinuity was linearly related to roughness, and subsequently an empirical equation was established. The long-term strength of discontinuity can be identified using the inflection point of the isocreep-rate curve, and it was linearly related to roughness. Furthermore, the ratio of long-term and instantaneous strength decreased with the increase of roughness. The shear-stiffness coefficient increased with the increase of shear rate, and the influence of shear rate on the shear stiffness coefficient decreased with the decrease of roughness. Further study of the mechanism revealed that these results could be attributed to the different time-dependent behavior of intact and joint rocks.

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.

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”

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.

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

Directory of Open Access Journals (Sweden)

Zhang Xuemei

2009-01-01

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

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 Elliptic Boundary Value Problems at Resonance with Nonlinear Wentzell Boundary Conditions

Directory of Open Access Journals (Sweden)

Ciprian G. Gal

2017-01-01

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

Relativistic time-dependent Fermion-mass renormalization using statistical regularization

Science.gov (United States)

Kutnink, Timothy; McMurray, Christian; Santrach, Amelia; Hockett, Sarah; Barcus, Scott; Petridis, Athanasios

2017-09-01

The time-dependent electromagnetically self-coupled Dirac equation is solved numerically by means of the staggered-leap-frog algorithm with reflecting boundary conditions. The stability region of the method versus the interaction strength and the spatial-grid size over time-step ratio is established. The expectation values of several dynamic operators are then evaluated as functions of time. These include the fermion and electromagnetic energies and the fermion dynamic mass. There is a characteristic, non-exponential, oscillatory dependence leading to asymptotic constants of these expectation values. In the case of the fermion mass this amounts to renormalization. The dependence of the expectation values on the spatial-grid size is evaluated in detail. Furthermore, the contribution of positive and negative energy states to the asymptotic values and the gauge fields is analyzed. Statistical regularization, employing a canonical ensemble whose temperature is the inverse of the grid size, is used to remove the grid-size and momentum-dependence and produce a finite result in the continuum limit.

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.

Reactive scattering theory for molecular transitions in time-dependent fields

International Nuclear Information System (INIS)

Peskin, U.; Miller, W.H.

1995-01-01

A new approach is introduced for computing probabilities of molecular transitions in time-dependent fields. The method is based on the stationary (t,t') representation of the Schroedinger equation and is shown to be equivalent to infinite order time-dependent perturbation theory. Bound-to-bound (i.e., photoexcitation) and bound-to-continuum (i.e., photoreaction) transitions are regarded as reactive collisions with the ''time coordinate'' as the reaction coordinate in an extended Hilbert space. A numerical method based on imposing absorbing boundary conditions for the time coordinate in a discrete variable representation framework is introduced. A single operation of the Green's operator provides all the state-specific transition probabilities as well as partial state-resolved (inclusive) reaction probabilities. Illustrative numerical applications are given for model systems

Greenland Regional and Ice Sheet-wide Geometry Sensitivity to Boundary and Initial conditions

Science.gov (United States)

Logan, L. C.; Narayanan, S. H. K.; Greve, R.; Heimbach, P.

2017-12-01

Ice sheet and glacier model outputs require inputs from uncertainly known initial and boundary conditions, and other parameters. Conservation and constitutive equations formalize the relationship between model inputs and outputs, and the sensitivity of model-derived quantities of interest (e.g., ice sheet volume above floatation) to model variables can be obtained via the adjoint model of an ice sheet. We show how one particular ice sheet model, SICOPOLIS (SImulation COde for POLythermal Ice Sheets), depends on these inputs through comprehensive adjoint-based sensitivity analyses. SICOPOLIS discretizes the shallow-ice and shallow-shelf approximations for ice flow, and is well-suited for paleo-studies of Greenland and Antarctica, among other computational domains. The adjoint model of SICOPOLIS was developed via algorithmic differentiation, facilitated by the source transformation tool OpenAD (developed at Argonne National Lab). While model sensitivity to various inputs can be computed by costly methods involving input perturbation simulations, the time-dependent adjoint model of SICOPOLIS delivers model sensitivities to initial and boundary conditions throughout time at lower cost. Here, we explore both the sensitivities of the Greenland Ice Sheet's entire and regional volumes to: initial ice thickness, precipitation, basal sliding, and geothermal flux over the Holocene epoch. Sensitivity studies such as described here are now accessible to the modeling community, based on the latest version of SICOPOLIS that has been adapted for OpenAD to generate correct and efficient adjoint code.

Bayesian Inference for Linear Parabolic PDEs with Noisy Boundary Conditions

KAUST Repository

Ruggeri, Fabrizio; Sawlan, Zaid A; Scavino, Marco; Tempone, Raul

2015-01-01

have been assumed for the time-dependent Dirichlet boundary values. Our approach is applied to synthetic data for the one-dimensional heat equation model, where the thermal diffusivity is the unknown parameter. We show how to infer the thermal

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.

Effects of Time-Dependent Inflow Perturbations on Turbulent Flow in a Street Canyon

Science.gov (United States)

Duan, G.; Ngan, K.

2017-12-01

Urban flow and turbulence are driven by atmospheric flows with larger horizontal scales. Since building-resolving computational fluid dynamics models typically employ steady Dirichlet boundary conditions or forcing, the accuracy of numerical simulations may be limited by the neglect of perturbations. We investigate the sensitivity of flow within a unit-aspect-ratio street canyon to time-dependent perturbations near the inflow boundary. Using large-eddy simulation, time-periodic perturbations to the streamwise velocity component are incorporated via the nudging technique. Spatial averages of pointwise differences between unperturbed and perturbed velocity fields (i.e., the error kinetic energy) show a clear dependence on the perturbation period, though spatial structures are largely insensitive to the time-dependent forcing. The response of the error kinetic energy is maximized for perturbation periods comparable to the time scale of the mean canyon circulation. Frequency spectra indicate that this behaviour arises from a resonance between the inflow forcing and the mean motion around closed streamlines. The robustness of the results is confirmed using perturbations derived from measurements of roof-level wind speed.

Metric space construction for the boundary of space-time

International Nuclear Information System (INIS)

Meyer, D.A.

1986-01-01

A distance function between points in space-time is defined and used to consider the manifold as a topological metric space. The properties of the distance function are investigated: conditions under which the metric and manifold topologies agree, the relationship with the causal structure of the space-time and with the maximum lifetime function of Wald and Yip, and in terms of the space of causal curves. The space-time is then completed as a topological metric space; the resultant boundary is compared with the causal boundary and is also calculated for some pertinent examples

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

Similarity solutions of time-dependent relativistic radiation-hydrodynamical plane-parallel flows

Science.gov (United States)

Fukue, Jun

2018-04-01

Similarity solutions are examined for the frequency-integrated relativistic radiation-hydrodynamical flows, which are described by the comoving quantities. The flows are vertical plane-parallel time-dependent ones with a gray opacity coefficient. For adequate boundary conditions, the flows are accelerated in a somewhat homologous manner, but terminate at some singular locus, which originates from the pathological behavior in relativistic radiation moment equations truncated in finite orders.

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.

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)

Impact of chemical lateral boundary conditions in a regional air quality forecast model on surface ozone predictions during stratospheric intrusions

Science.gov (United States)

Pendlebury, Diane; Gravel, Sylvie; Moran, Michael D.; Lupu, Alexandru

2018-02-01

A regional air quality forecast model, GEM-MACH, is used to examine the conditions under which a limited-area air quality model can accurately forecast near-surface ozone concentrations during stratospheric intrusions. Periods in 2010 and 2014 with known stratospheric intrusions over North America were modelled using four different ozone lateral boundary conditions obtained from a seasonal climatology, a dynamically-interpolated monthly climatology, global air quality forecasts, and global air quality reanalyses. It is shown that the mean bias and correlation in surface ozone over the course of a season can be improved by using time-varying ozone lateral boundary conditions, particularly through the correct assignment of stratospheric vs. tropospheric ozone along the western lateral boundary (for North America). Part of the improvement in surface ozone forecasts results from improvements in the characterization of near-surface ozone along the lateral boundaries that then directly impact surface locations near the boundaries. However, there is an additional benefit from the correct characterization of the location of the tropopause along the western lateral boundary such that the model can correctly simulate stratospheric intrusions and their associated exchange of ozone from stratosphere to troposphere. Over a three-month period in spring 2010, the mean bias was seen to improve by as much as 5 ppbv and the correlation by 0.1 depending on location, and on the form of the chemical lateral boundary condition.

Modified Pressure-Correction Projection Methods: Open Boundary and Variable Time Stepping

KAUST Repository

Bonito, Andrea

2014-10-31

© Springer International Publishing Switzerland 2015. In this paper, we design and study two modifications of the first order standard pressure increment projection scheme for the Stokes system. The first scheme improves the existing schemes in the case of open boundary condition by modifying the pressure increment boundary condition, thereby minimizing the pressure boundary layer and recovering the optimal first order decay. The second scheme allows for variable time stepping. It turns out that the straightforward modification to variable time stepping leads to unstable schemes. The proposed scheme is not only stable but also exhibits the optimal first order decay. Numerical computations illustrating the theoretical estimates are provided for both new schemes.

Modified Pressure-Correction Projection Methods: Open Boundary and Variable Time Stepping

KAUST Repository

Bonito, Andrea; Guermond, Jean-Luc; Lee, Sanghyun

2014-01-01

© Springer International Publishing Switzerland 2015. In this paper, we design and study two modifications of the first order standard pressure increment projection scheme for the Stokes system. The first scheme improves the existing schemes in the case of open boundary condition by modifying the pressure increment boundary condition, thereby minimizing the pressure boundary layer and recovering the optimal first order decay. The second scheme allows for variable time stepping. It turns out that the straightforward modification to variable time stepping leads to unstable schemes. The proposed scheme is not only stable but also exhibits the optimal first order decay. Numerical computations illustrating the theoretical estimates are provided for both new schemes.

Time integrated spectroscopy of turbid media based on the microscopic beer-lambert law: application to small-size phantoms having different boundary conditions.

Science.gov (United States)

Zhang, H; Urakami, T; Tsuchiya, Y; Lu, Z; Hiruma, T

1999-01-01

Continued work on time-integrated spectroscopy (TIS) is presented to quantify absorber concentrations in turbid media. We investigated the applicability of the TIS method to small-size media that have different boundary conditions by measuring two 20×20×50 mm3 cuboid liquid tissue-like phantoms at various absorption levels (absorption coefficients of the phantom from 2.5×10-3 to 4.4×10-2 mm-1 at 782 nm and from 3.1×10-3 to 2.7×10-2 mm-1 at 831 nm). The scattering and absorbing solution was filled into ordinary and black-anodized aluminum containers to provide different boundary conditions. By means of a single equation, the absorber concentrations have been recovered within errors of a few percent in both cases. This demonstrates that the TIS method can quantify absorbers in small-size media having different boundary conditions. © 1999 Society of Photo-Optical Instrumentation Engineers.

Numerical simulation of time-dependent deformations under hygral and thermal transient conditions

International Nuclear Information System (INIS)

Roelfstra, P.E.

1987-01-01

Some basic concepts of numerical simulation of the formation of the microstructure of HCP are outlined. The aim is to replace arbitrary terms like aging by more realistic terms like bond density in the xerogel and bonds between hydrating particles of HCP. Actual state parameters such as temperature, humidity and degree of hydration can be determined under transient hygral and thermal conditions by solving numerically a series of appropriate coupled differential equations with given boundary conditions. Shrinkage of a composite structure without crack formation, based on calculated moisture distributions, has been determined with numerical concrete codes. The influence of crack formation, tensile strain-hardening and softening on the total deformation of a quasi-homogeneous drying material has been studied by means of model based on FEM. The difference between shrinkage without crack formation and shrinkage with crack formation can be quantified. Drying shrinkage and creep of concrete cannot be separated. The total deformation depends on the superimposed stress fields. Transient hygral deformation can be realistically predicted if the concept of point properties is applied rigorously. Transient thermal deformation has to be dealt with in the same way. (orig./HP)

The Dependence of Chimera States on Initial Conditions

International Nuclear Information System (INIS)

Feng Yue-E; Li Hai-Hong

2015-01-01

A chimera state consisting of both coherent and incoherent groups is a fascinating spatial pattern in non-locally coupled identical oscillators. It is thought that random initial conditions hardly evolve to chimera states. In this work, we study the dependence of chimera states on initial conditions. We show that random initial conditions may lead to chimera states and the chance of realizing chimera states becomes increasing when the model parameters are moving away from the boundary of their stable regime. (paper)

Bound states on the lattice with partially twisted boundary conditions

International Nuclear Information System (INIS)

Agadjanov, D.; Guo, F.-K.; Ríos, G.; Rusetsky, A.

2015-01-01

We propose a method to study the nature of exotic hadrons by determining the wave function renormalization constant Z from lattice simulations. It is shown that, instead of studying the volume-dependence of the spectrum, one may investigate the dependence of the spectrum on the twisting angle, imposing twisted boundary conditions on the fermion fields on the lattice. In certain cases, e.g., the case of the DK bound state which is addressed in detail, it is demonstrated that the partial twisting is equivalent to the full twisting up to exponentially small corrections.

Postural orientation : Age-related changes in variability and time-to-boundary

NARCIS (Netherlands)

Van Wegen, E. E.H.; Van Emmerik, R. E.A.; Riccio, G. E.

2002-01-01

The relation between age-specific postural instability and the detection of stability boundaries was examined. Balance control was investigated under different visual conditions (eyes open/closed) and postural orientations·(forward/backward lean) while standing on a force platform. Dependent

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

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

Directory of Open Access Journals (Sweden)

Meiqiang Feng

2009-01-01

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

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.

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.

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

Science.gov (United States)

Bardhan, Jaydeep P; Knepley, Matthew G

2014-10-07

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

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

International Nuclear Information System (INIS)

Bardhan, Jaydeep P.; Knepley, Matthew G.

2014-01-01

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

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

Energy Technology Data Exchange (ETDEWEB)

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

2014-10-07

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

Time-dependent radiation transfer with rayleigh scattering in finite plane-parallel media using pomraning-eddington approximation

International Nuclear Information System (INIS)

El-Wakil, S.A.; Sallah, M.; Degheidy, A.R.

2005-01-01

The time-dependent radiation transfer equation in plane geometry with Rayleigh scattering is studied. The traveling wave transformation is used to obtain the corresponding stationary-like equation. Pomraning-Eddington approximation is then used to calculate the radiation intensity in finite plane-parallel media. Numerical results and shielding calculations are shown for reflectivity and transmissivity at different times. The medium is assumed to have specular-reflecting boundaries. For the sake of comparison, two different weight functions are introduced and to force the boundary conditions to be fulfilled

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)

Global existence of strong solutions to the three- dimensional incompressible Navier-Stokes equations with special boundary conditions

Science.gov (United States)

Riley, Douglas A.

We study the three-dimensional incompressible Navier- Stokes equations in a domain of the form W'×(0,e) . First, we assume W' is a C3 bounded domain and impose no-slip boundary conditions on 6W'×(0,e ) , and periodic conditions on W'×0,e . Physically, this models fluid flow through a pipe with cross-section W' where the inlet and outlet conditions are assumed periodic. Secondly, we assume W'=(0,l4) ×(0,l5) and impose periodic boundary conditions. This problem is of interest mathematically, and has been more widely considered than the pipe flow problem. For both sets of boundary conditions, we show that a strong solution exists for all time with conditions on the initial data and forcing. We start by recalling that if the forcing function and initial condition do not depend on x3, then a global strong solution exists which also does not depend on x3. Here (x1,x2,x3) ∈W≡W'×( 0,e) . With this observation as motivation, and using an additive decomposition introduced by Raugel and Sell, we split the initial data and forcing into a portion independent of x3 and a remainder. In our first result, we impose a smallness condition on the remainder and assume the forcing function is square- integrable in time as a function into L2(W) . With these assumptions, we prove a global existence theorem that does not require a smallness condition on e or on the portion of the initial condition and forcing independent of x3. However, these quantities do affect the allowable size of the remainder. For our second result, we assume the forcing is only bounded in time as a function into L2(W) . In this case, we need a smallness condition on the initial data, the forcing, and e to obtain global existence. The interesting observation is that the allowable sizes for the initial data and forcing grow as e-->0 . Thus, we obtain a `thin-domain' result as originally obtained by Raugel and Sell. In fact, our results allow the portion of the initial data and forcing independent of x3 to

First-passage time asymptotics over moving boundaries for random walk bridges

OpenAIRE

Sloothaak, F.; Zwart, B.; Wachtel, V.

2017-01-01

We study the asymptotic tail probability of the first-passage time over a moving boundary for a random walk conditioned to return to zero, where the increments of the random walk have finite variance. Typically, the asymptotic tail behavior may be described through a regularly varying function with exponent -1/2, where the impact of the boundary is captured by the slowly varying function. Yet, the moving boundary may have a stronger effect when the tail is considered at a time close to the re...

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.

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

Direct measurement of time dependent diffusion for Ag and Au under ambient conditions

Energy Technology Data Exchange (ETDEWEB)

Yoo, Pil Sun; Jo, Han Yeol; Kim, Tae Kyeong [Hankuk University of Foreign Studies, Yongin (Korea, Republic of)

2014-12-15

Time-dependent diffusion for Ag and Au metal atoms was measured using the scanning tunneling microscope break-junction technique in ambient conditions. We observed that Ag contacts do not form long single-atomic chains compared to Au contacts during the elongation of each metal electrode, and Ag atoms diffuse more quickly than Au atoms after metal contact rupture. This is consistent with previous results of molecular dynamic simulations. Further, we found a correlation between diffusion length and the evolution time on an atomic scale to reveal the time-dependent diffusion for Ag and Au metal atoms.

Analysis of the dual phase lag bio-heat transfer equation with constant and time-dependent heat flux conditions on skin surface

Directory of Open Access Journals (Sweden)

Ziaei Poor Hamed

2016-01-01

Full Text Available This article focuses on temperature response of skin tissue due to time-dependent surface heat fluxes. Analytical solution is constructed for DPL bio-heat transfer equation with constant, periodic and pulse train heat flux conditions on skin surface. Separation of variables and Duhamel’s theorem for a skin tissue as a finite domain are employed. The transient temperature responses for constant and time-dependent boundary conditions are obtained and discussed. The results show that there is major discrepancy between the predicted temperature of parabolic (Pennes bio-heat transfer, hyperbolic (thermal wave and DPL bio-heat transfer models when high heat flux accidents on the skin surface with a short duration or propagation speed of thermal wave is finite. The results illustrate that the DPL model reduces to the hyperbolic model when τT approaches zero and the classic Fourier model when both thermal relaxations approach zero. However for τq = τT the DPL model anticipates different temperature distribution with that predicted by the Pennes model. Such discrepancy is due to the blood perfusion term in energy equation. It is in contrast to results from the literature for pure conduction material, where the DPL model approaches the Fourier heat conduction model when τq = τT . The burn injury is also investigated.

Asymmetric simple exclusion process with position-dependent hopping rates: Phase diagram from boundary-layer analysis.

Science.gov (United States)

Mukherji, Sutapa

2018-03-01

In this paper, we study a one-dimensional totally asymmetric simple exclusion process with position-dependent hopping rates. Under open boundary conditions, this system exhibits boundary-induced phase transitions in the steady state. Similarly to totally asymmetric simple exclusion processes with uniform hopping, the phase diagram consists of low-density, high-density, and maximal-current phases. In various phases, the shape of the average particle density profile across the lattice including its boundary-layer parts changes significantly. Using the tools of boundary-layer analysis, we obtain explicit solutions for the density profile in different phases. A detailed analysis of these solutions under different boundary conditions helps us obtain the equations for various phase boundaries. Next, we show how the shape of the entire density profile including the location of the boundary layers can be predicted from the fixed points of the differential equation describing the boundary layers. We discuss this in detail through several examples of density profiles in various phases. The maximal-current phase appears to be an especially interesting phase where the boundary layer flows to a bifurcation point on the fixed-point diagram.

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

OpenCL-Based FPGA Accelerator for 3D FDTD with Periodic and Absorbing Boundary Conditions

Directory of Open Access Journals (Sweden)

Hasitha Muthumala Waidyasooriya

2017-01-01

Full Text Available Finite difference time domain (FDTD method is a very poplar way of numerically solving partial differential equations. FDTD has a low operational intensity so that the performances in CPUs and GPUs are often restricted by the memory bandwidth. Recently, deeply pipelined FPGA accelerators have shown a lot of success by exploiting streaming data flows in FDTD computation. In spite of this success, many FPGA accelerators are not suitable for real-world applications that contain complex boundary conditions. Boundary conditions break the regularity of the data flow, so that the performances are significantly reduced. This paper proposes an FPGA accelerator that computes commonly used absorbing and periodic boundary conditions in many 3D FDTD applications. Accelerator is designed using a “C-like” programming language called OpenCL (open computing language. As a result, the proposed accelerator can be customized easily by changing the software code. According to the experimental results, we achieved over 3.3 times and 1.5 times higher processing speed compared to the CPUs and GPUs, respectively. Moreover, the proposed accelerator is more than 14 times faster compared to the recently proposed FPGA accelerators that are capable of handling complex boundary conditions.

Role of Grain Boundaries under Long-Time Radiation

Science.gov (United States)

Zhu, Yichao; Luo, Jing; Guo, Xu; Xiang, Yang; Chapman, Stephen Jonathan

2018-06-01

Materials containing a high proportion of grain boundaries offer significant potential for the development of radiation-resistant structural materials. However, a proper understanding of the connection between the radiation-induced microstructural behavior of a grain boundary and its impact at long natural time scales is still missing. In this Letter, point defect absorption at interfaces is summarized by a jump Robin-type condition at a coarse-grained level, wherein the role of interface microstructure is effectively taken into account. Then a concise formula linking the sink strength of a polycrystalline aggregate with its grain size is introduced and is well compared with experimental observation. Based on the derived model, a coarse-grained formulation incorporating the coupled evolution of grain boundaries and point defects is proposed, so as to underpin the study of long-time morphological evolution of grains induced by irradiation. Our simulation results suggest that the presence of point defect sources within a grain further accelerates its shrinking process, and radiation tends to trigger the extension of twin boundary sections.

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

Effect of velocity boundary conditions on the heat transfer and flow topology in two-dimensional Rayleigh-Bénard convection.

Science.gov (United States)

van der Poel, Erwin P; Ostilla-Mónico, Rodolfo; Verzicco, Roberto; Lohse, Detlef

2014-07-01

The effect of various velocity boundary condition is studied in two-dimensional Rayleigh-Bénard convection. Combinations of no-slip, stress-free, and periodic boundary conditions are used on both the sidewalls and the horizontal plates. For the studied Rayleigh numbers Ra between 10(8) and 10(11) the heat transport is lower for Γ=0.33 than for Γ=1 in case of no-slip sidewalls. This is, surprisingly, the opposite for stress-free sidewalls, where the heat transport increases for a lower aspect ratio. In wider cells the aspect-ratio dependence is observed to disappear for Ra ≥ 10(10). Two distinct flow types with very different dynamics can be seen, mostly dependent on the plate velocity boundary condition, namely roll-like flow and zonal flow, which have a substantial effect on the dynamics and heat transport in the system. The predominantly horizontal zonal flow suppresses heat flux and is observed for stress-free and asymmetric plates. Low aspect-ratio periodic sidewall simulations with a no-slip boundary condition on the plates also exhibit zonal flow. In all the other cases, the flow is roll like. In two-dimensional Rayleigh-Bénard convection, the velocity boundary conditions thus have large implications on both roll-like and zonal flow that have to be taken into consideration before the boundary conditions are imposed.

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

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.

Numerical tests of evolution systems, gauge conditions, and boundary conditions for 1D colliding gravitational plane waves

International Nuclear Information System (INIS)

Bardeen, J.M.; Buchman, L.T.

2002-01-01

We investigate how the accuracy and stability of numerical relativity simulations of 1D colliding plane waves depends on choices of equation formulations, gauge conditions, boundary conditions, and numerical methods, all in the context of a first-order 3+1 approach to the Einstein equations, with basic variables some combination of first derivatives of the spatial metric and components of the extrinsic curvature tensor. Hyperbolic schemes, specifically variations on schemes proposed by Bona and Masso and Anderson and York, are compared with variations of the Arnowitt-Deser-Misner formulation. Modifications of the three basic schemes include raising one index in the metric derivative and extrinsic curvature variables and adding a multiple of the energy constraint to the extrinsic curvature evolution equations. Redundant variables in the Bona-Masso formulation may be reset frequently or allowed to evolve freely. Gauge conditions which simplify the dynamical structure of the system are imposed during each time step, but the lapse and shift are reset periodically to control the evolution of the spacetime slicing and the longitudinal part of the metric. We show that physically correct boundary conditions, satisfying the energy and momentum constraint equations, generically require the presence of some ingoing eigenmodes of the characteristic matrix. Numerical methods are developed for the hyperbolic systems based on decomposing flux differences into linear combinations of eigenvectors of the characteristic matrix. These methods are shown to be second-order accurate, and in practice second-order convergent, for smooth solutions, even when the eigenvectors and eigenvalues of the characteristic matrix are spatially varying

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

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

Science.gov (United States)

Bardhan, Jaydeep P.; Knepley, Matthew G.

2014-01-01

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

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)

Causal boundary for stably causal space-times

International Nuclear Information System (INIS)

Racz, I.

1987-12-01

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

A simple method to calculate first-passage time densities with arbitrary initial conditions

Science.gov (United States)

Nyberg, Markus; Ambjörnsson, Tobias; Lizana, Ludvig

2016-06-01

Numerous applications all the way from biology and physics to economics depend on the density of first crossings over a boundary. Motivated by the lack of general purpose analytical tools for computing first-passage time densities (FPTDs) for complex problems, we propose a new simple method based on the independent interval approximation (IIA). We generalise previous formulations of the IIA to include arbitrary initial conditions as well as to deal with discrete time and non-smooth continuous time processes. We derive a closed form expression for the FPTD in z and Laplace-transform space to a boundary in one dimension. Two classes of problems are analysed in detail: discrete time symmetric random walks (Markovian) and continuous time Gaussian stationary processes (Markovian and non-Markovian). Our results are in good agreement with Langevin dynamics simulations.

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.

Experimental Investigation of Separated and Transitional Boundary Layers Under Low-Pressure Turbine Airfoil Conditions

Science.gov (United States)

Hultgren, Lennart S.; Volino, Ralph J.

2002-01-01

Modern low-pressure turbine airfoils are subject to increasingly stronger pressure gradients as designers impose higher loading in an effort to improve efficiency and to reduce part count. The adverse pressure gradients on the suction side of these airfoils can lead to boundary-layer separation, particularly under cruise conditions. Separation bubbles, notably those which fail to reattach, can result in a significant degradation of engine efficiency. Accurate prediction of separation and reattachment is hence crucial to improved turbine design. This requires an improved understanding of the transition flow physics. Transition may begin before or after separation, depending on the Reynolds number and other flow conditions, has a strong influence on subsequent reattachment, and may even eliminate separation. Further complicating the problem are the high free-stream turbulence levels in a real engine environment, the strong pressure gradients along the airfoils, the curvature of the airfoils, and the unsteadiness associated with wake passing from upstream stages. Because of the complicated flow situation, transition in these devices can take many paths that can coexist, vary in importance, and possibly also interact, at different locations and instances in time. The present work was carried out in an attempt to systematically sort out some of these issues. Detailed velocity measurements were made along a flat plate subject to the same nominal dimensionless pressure gradient as the suction side of a modern low-pressure turbine airfoil ('Pak-B'). The Reynolds number based on wetted plate length and nominal exit velocity, Re, was varied from 50;000 to 300; 000, covering cruise to takeoff conditions. Low, 0.2%, and high, 7%, inlet free-stream turbulence intensities were set using passive grids. These turbulence levels correspond to about 0.2% and 2.5% turbulence intensity in the test section when normalized with the exit velocity. The Reynolds number and free

Methodology for determining time-dependent mechanical properties of tuff subjected to near-field repository conditions

International Nuclear Information System (INIS)

Blacic, J.D.; Andersen, R.

1983-01-01

We have established a methodology to determine the time dependence of strength and transport properties of tuff under conditions appropriate to a nuclear waste repository. Exploratory tests to determine the approximate magnitudes of thermomechanical property changes are nearly complete. In this report we describe the capabilities of an apparatus designed to precisely measure the time-dependent deformation and permeability of tuff at simulated repository conditions. Preliminary tests with this new apparatus indicate that microclastic creep failure of tuff occurs over a narrow strain range with little precursory Tertiary creep behavior. In one test, deformation under conditions of slowly decreasing effective pressure resulted in failure, whereas some strain indicators showed a decreasing rate of strain

Methodology for determining time-dependent mechanical properties of tuff subjected to near-field repository conditions

Energy Technology Data Exchange (ETDEWEB)

Blacic, J.D.; Andersen, R.

1983-01-01

We have established a methodology to determine the time dependence of strength and transport properties of tuff under conditions appropriate to a nuclear waste repository. Exploratory tests to determine the approximate magnitudes of thermomechanical property changes are nearly complete. In this report we describe the capabilities of an apparatus designed to precisely measure the time-dependent deformation and permeability of tuff at simulated repository conditions. Preliminary tests with this new apparatus indicate that microclastic creep failure of tuff occurs over a narrow strain range with little precursory Tertiary creep behavior. In one test, deformation under conditions of slowly decreasing effective pressure resulted in failure, whereas some strain indicators showed a decreasing rate of strain.

Comparative study of the free-surface boundary condition in two-dimensional finite-difference elastic wave field simulation

International Nuclear Information System (INIS)

Lan, Haiqiang; Zhang, Zhongjie

2011-01-01

The finite-difference (FD) method is a powerful tool in seismic wave field modelling for understanding seismic wave propagation in the Earth's interior and interpreting the real seismic data. The accuracy of FD modelling partly depends on the implementation of the free-surface (i.e. traction-free) condition. In the past 40 years, at least six kinds of free-surface boundary condition approximate schemes (such as one-sided, centred finite-difference, composed, new composed, implicit and boundary-modified approximations) have been developed in FD second-order elastodynamic simulation. Herein we simulate seismic wave fields in homogeneous and lateral heterogeneous models using these free-surface boundary condition approximate schemes and evaluate their stability and applicability by comparing with corresponding analytical solutions, and then quantitatively evaluate the accuracies of different approximate schemes from the misfit of the amplitude and phase between the numerical and analytical results. Our results confirm that the composed scheme becomes unstable for the V s /V p ratio less than 0.57, and suggest that (1) the one-sided scheme is only accurate to first order and therefore introduces serious errors for the shorter wavelengths, other schemes are all of second-order precision; (2) the new composed, implicit and boundary-modified schemes are stable even when the V s /V p ratio is less than 0.2; (3) the implicit and boundary-modified schemes are able to deal with laterally varying (heterogeneous) free surface; (4) in the corresponding stability range, the one-sided scheme shows remarkable errors in both phase and amplitude compared to analytical solution (which means larger errors in travel-time and reflection strength), the other five approximate schemes show better performance in travel-time (phase) than strength (amplitude)

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

Science.gov (United States)

Azarnavid, Babak; Parand, Kourosh; Abbasbandy, Saeid

2018-06-01

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

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.

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.

Improving Delay-Range-Dependent Stability Condition for Systems with Interval Time-Varying Delay

Directory of Open Access Journals (Sweden)

Wei Qian

2013-01-01

Full Text Available This paper discusses the delay-range-dependent stability for systems with interval time-varying delay. Through defining the new Lyapunov-Krasovskii functional and estimating the derivative of the LKF by introducing new vectors, using free matrices and reciprocally convex approach, the new delay-range-dependent stability conditions are obtained. Two well-known examples are given to illustrate the less conservatism of the proposed theoretical results.

Monopole Giant Resonances and TDHF boundary conditions

International Nuclear Information System (INIS)

Stevenson, P.D.; Almehed, D.; Reinhard, P.-G.; Maruhn, J.A.

2007-01-01

Using time-dependent Hartree-Fock, we induce isoscalar and isovector monopole vibrations and follow the subsequent vibrations of both the same and opposite isospin nature in the N Z nucleus 132 Sn. By suitable scaling of the proton and neutron parts of the excitation operators, the coupling between the modes is studied, and the approximate normal modes found. Chaotic dynamics are then analysed in the isoscalar giant monopole resonance by using reflecting boundaries in a large space to build up a large number of 0 + states whose spacings are then analysed. A Wigner-like distribution is found

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.

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.

Exact solutions to the supply chain equations for arbitrary, time-dependent demands

DEFF Research Database (Denmark)

Warburton, Roger D.H.; Hodgson, J.P.E.; Nielsen, Erland Hejn

2014-01-01

, so users can determine the inventory behavior to any desired precision. To illustrate, we solve the equations for a non-linear, quadratic time-dependence in the demand. For practical use, only a few terms in the series are required, a proposition illustrated by the For All Practical Purposes (FAPP......We study the impact on inventory of an unexpected, non-linear, time-dependent demand and present the exact solutions over time to the supply chain equations without requiring any approximations. We begin by imposing a boundary condition of stability at infinity, from which we derive expressions...... for the estimated demand and the target work in progress when the demand is time-dependent. The resulting inventory equation is solved in terms of the Lambert modes with all of the demand non-linearities confined to the pre-shape function. The series solution is exact, and all terms are reasonably easy to calculate...

Analytical solution for multi-species contaminant transport in finite media with time-varying boundary conditions

Science.gov (United States)

Most analytical solutions available for the equations governing the advective-dispersive transport of multiple solutes undergoing sequential first-order decay reactions have been developed for infinite or semi-infinite spatial domains and steady-state boundary conditions. In this work we present an ...

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.

78 FR 15883 - Standard Time Zone Boundaries

Science.gov (United States)

2013-03-13

...] RIN 2105-AE20 Standard Time Zone Boundaries AGENCY: Office of the Secretary (OST), Department of... time zone boundaries regulations to reflect changes that Congress made to the Uniform Time Act. The... regulations on standard time zone boundaries, 49 CFR Part 71, need to be updated in order to ensure their...

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.

A Two-Dimensional Transverse Magnetic Propagation Model of a Sine Wave Using Mur Boundary Conditions

National Research Council Canada - National Science Library

Korjack, T

1997-01-01

.... The nonreflecting boundary conditions due to Mur were used at the boundary surfaces. Electric field intensity distributions resulted over a progressive time expansion to illustrate the propagation effect over the entire 2-D mesh...

Time-dependent reliability analysis and condition assessment of structures

International Nuclear Information System (INIS)

Ellingwood, B.R.

1997-01-01

Structures generally play a passive role in assurance of safety in nuclear plant operation, but are important if the plant is to withstand the effect of extreme environmental or abnormal events. Relative to mechanical and electrical components, structural systems and components would be difficult and costly to replace. While the performance of steel or reinforced concrete structures in service generally has been very good, their strengths may deteriorate during an extended service life as a result of changes brought on by an aggressive environment, excessive loading, or accidental loading. Quantitative tools for condition assessment of aging structures can be developed using time-dependent structural reliability analysis methods. Such methods provide a framework for addressing the uncertainties attendant to aging in the decision process

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.

Artificial boundary conditions for the numerical simulation of unsteady acoustic waves

International Nuclear Information System (INIS)

Tsynkov, S.V.

2003-01-01

We construct non-local artificial boundary conditions (ABCs) for the numerical simulation of genuinely time-dependent acoustic waves that propagate from a compact source in an unbounded unobstructed space. The key property used for obtaining the ABCs is the presence of lacunae, i.e., sharp aft fronts of the waves, in wave-type solutions in odd-dimension spaces. This property can be considered a manifestation of the Huygens' principle. The ABCs are obtained directly for the discrete formulation of the problem. They truncate the original unbounded domain and guarantee the complete transparency of the new outer boundary for all the outgoing waves. A central feature of the proposed ABCs is that the extent of their temporal non-locality is fixed and limited, and it does not come at the expense of simplifying the original model. It is rather a natural consequence of the existence of lacunae, which is a fundamental property of the corresponding solutions. The proposed ABCs can be built for any consistent and stable finite-difference scheme. Their accuracy can always be made as high as that of the interior approximation, and it will not deteriorate even when integrating over long time intervals. Besides, the ABCs are most flexible from the standpoint of geometry and can handle irregular boundaries on regular grids with no fitting/adaptation needed and no accuracy loss induced. Finally, they allow for a wide range of model settings. In particular, not only one can analyze the simplest advective acoustics case with the uniform background flow, but also the case when the waves' source (or scatterer) is engaged in an accelerated motion

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.

Non-equilibrium effects of core-cooling and time-dependent internal heating on mantle flush events

Directory of Open Access Journals (Sweden)

D. A. Yuen

1995-01-01

Full Text Available We have examined the non-equilibrium effects of core-cooling and time-dependent internal-heating on the thermal evolution of the Earth's mantle and on mantle flush events caused by the two major phase transitions. Both two- and three-dimensional models have been employed. The mantle viscosity responds to the secular cooling through changes in the averaged temperature field. A viscosity which decreases algebraically with the average temperature has been considered. The time-dependent internal-heating is prescribed to decrease exponentially with a single decay time. We have studied the thermal histories with initial Rayleigh numbers between 2 x 107 and 108 . Flush events, driven by the non-equilibrium forcings, are much more dramatic than those produced by the equilibrium boundary conditions and constant internal heating. Multiple flush events are found under non-equilibrium conditions in which there is very little internal heating or very fast decay rates of internal-heating. Otherwise, the flush events take place in a relatively continuous fashion. Prior to massive flush events small-scale percolative structures appear in the 3D temperature fields. Time-dependent signatures, such as the surface heat flux, also exhibits high frequency oscillatory patterns prior to massive flush events. These two observations suggest that the flush event may be a self-organized critical phenomenon. The Nusselt number as a function of the time-varying Ra does not follow the Nusselt vs. Rayleigh number power-law relationship based on equilibrium (constant temperature boundary conditions. Instead Nu(t may vary non-monotonically with time because of the mantle flush events. Convective processes in the mantle operate quite differently under non-equilibrium conditions from its behaviour under the usual equilibrium situations.

Time asymmetries in quantum cosmology and the search for boundary conditions to the Wheeler-DeWitt equation

Energy Technology Data Exchange (ETDEWEB)

Castagnino, Mario [Instituto de Astronomia y Fisica del Espacio, Casilla de Correo 67, Sucursal 28, 1428 Buenos Aires (Argentina); Catren, Gabriel [Instituto de Astronomia y Fisica del Espacio, Casilla de Correo 67, Sucursal 28, 1428 Buenos Aires (Argentina); Ferraro, Rafael [Instituto de Astronomia y Fisica del Espacio, Casilla de Correo 67, Sucursal 28, 1428 Buenos Aires (Argentina)

2002-09-21

This paper addresses the quantization of minisuperspace cosmological models by studying a possible solution to the problem of time and time asymmetries in quantum cosmology. Since general relativity does not have a privileged time variable of Newtonian type, it is necessary, in order to have a dynamical evolution, to select a physical clock. This choice yields, in the proposed approach, to the breaking of the so-called clock-reversal invariance of the theory which is clearly distinguished from the well-known motion-reversal invariance of both classical and quantum mechanics. In light of this new perspective, the problem of imposing proper boundary conditions on the space of solutions of the Wheeler-DeWitt equation is reformulated. The symmetry-breaking formalism of previous papers is analyzed and a clarification of it is proposed in order to satisfy the requirements of the new interpretation.

Time asymmetries in quantum cosmology and the search for boundary conditions to the Wheeler-DeWitt equation

International Nuclear Information System (INIS)

Castagnino, Mario; Catren, Gabriel; Ferraro, Rafael

2002-01-01

This paper addresses the quantization of minisuperspace cosmological models by studying a possible solution to the problem of time and time asymmetries in quantum cosmology. Since general relativity does not have a privileged time variable of Newtonian type, it is necessary, in order to have a dynamical evolution, to select a physical clock. This choice yields, in the proposed approach, to the breaking of the so-called clock-reversal invariance of the theory which is clearly distinguished from the well-known motion-reversal invariance of both classical and quantum mechanics. In light of this new perspective, the problem of imposing proper boundary conditions on the space of solutions of the Wheeler-DeWitt equation is reformulated. The symmetry-breaking formalism of previous papers is analyzed and a clarification of it is proposed in order to satisfy the requirements of the new interpretation

Time-dependent theory of double ionization of helium under XUV radiation

International Nuclear Information System (INIS)

Nikolopoulos, L A A; Lambropoulos, P

2007-01-01

We present non-perturbative time-dependent calculations of single and double ionization of helium, under XUV radiation of photon energy ranging from 40 to 45 eV, through the direct propagation of the time-dependent Schroedinger equation. The time-dependent wavefunction of the atom under the field is expanded in terms of correlated multichannel states normalized with incoming-wave boundary conditions. In addition to presenting a new non-perturbative approach to the three-body problem, in a fully correlated scheme, capable of providing in the same calculation photoelectron energy and angularly resolved spectra, as well as cross sections through the lowest non-vanishing order transition amplitude, we also present a detailed comparison of the values of certain key quantities that have been obtained through a variety of other methods. The degree of agreement we find, while lending credence to the approach and its versatility, also highlights the remaining open questions in this novel context of double ionization

Time-Dependent Close-Coupling Methods for Electron-Atom/Molecule Scattering

International Nuclear Information System (INIS)

Colgan, James

2014-01-01

The time-dependent close-coupling (TDCC) method centers on an accurate representation of the interaction between two outgoing electrons moving in the presence of a Coulomb field. It has been extensively applied to many problems of electrons, photons, and ions scattering from light atomic targets. Theoretical Description: The TDCC method centers on a solution of the time-dependent Schrödinger equation for two interacting electrons. The advantages of a time-dependent approach are two-fold; one treats the electron-electron interaction essentially in an exact manner (within numerical accuracy) and a time-dependent approach avoids the difficult boundary condition encountered when two free electrons move in a Coulomb field (the classic three-body Coulomb problem). The TDCC method has been applied to many fundamental atomic collision processes, including photon-, electron- and ion-impact ionization of light atoms. For application to electron-impact ionization of atomic systems, one decomposes the two-electron wavefunction in a partial wave expansion and represents the subsequent two-electron radial wavefunctions on a numerical lattice. The number of partial waves required to converge the ionization process depends on the energy of the incoming electron wavepacket and on the ionization threshold of the target atom or ion.

Time-Dependent Expression of Arc and Zif268 after Acquisition of Fear Conditioning

Directory of Open Access Journals (Sweden)

Mary E. Lonergan

2010-01-01

Full Text Available Memory consolidation requires transcription and translation of new protein. Arc, an effector immediate early gene, and zif268, a regulatory transcription factor, have been implicated in synaptic plasticity underlying learning and memory. This study explored the temporal expression profiles of these proteins in the rat hippocampus following fear conditioning. We observed a time-dependent increase of Arc protein in the dorsal hippocampus 30-to-90-minute post training, returning to basal levels at 4 h. Zif268 protein levels, however, gradually increased at 30-minute post training before peaking in expression at 60 minute. The timing of hippocampal Arc and zif268 expression coincides with the critical period for protein synthesis-dependent memory consolidation following fear conditioning. However, the expression of Arc protein appears to be driven by context exploration, whereas, zif268 expression may be more specifically related to associative learning. These findings suggest that altered Arc and zif268 expression are related to neural plasticity during the formation of fear memory.

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

KAUST Repository

Ruggeri, Fabrizio

2016-05-12

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

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.

Immersed Boundary-Lattice Boltzmann Method Using Two Relaxation Times

Directory of Open Access Journals (Sweden)

Kosuke Hayashi

2012-06-01

Full Text Available An immersed boundary-lattice Boltzmann method (IB-LBM using a two-relaxation time model (TRT is proposed. The collision operator in the lattice Boltzmann equation is modeled using two relaxation times. One of them is used to set the fluid viscosity and the other is for numerical stability and accuracy. A direct-forcing method is utilized for treatment of immersed boundary. A multi-direct forcing method is also implemented to precisely satisfy the boundary conditions at the immersed boundary. Circular Couette flows between a stationary cylinder and a rotating cylinder are simulated for validation of the proposed method. The method is also validated through simulations of circular and spherical falling particles. Effects of the functional forms of the direct-forcing term and the smoothed-delta function, which interpolates the fluid velocity to the immersed boundary and distributes the forcing term to fixed Eulerian grid points, are also examined. As a result, the following conclusions are obtained: (1 the proposed method does not cause non-physical velocity distribution in circular Couette flows even at high relaxation times, whereas the single-relaxation time (SRT model causes a large non-physical velocity distortion at a high relaxation time, (2 the multi-direct forcing reduces the errors in the velocity profile of a circular Couette flow at a high relaxation time, (3 the two-point delta function is better than the four-point delta function at low relaxation times, but worse at high relaxation times, (4 the functional form of the direct-forcing term does not affect predictions, and (5 circular and spherical particles falling in liquids are well predicted by using the proposed method both for two-dimensional and three-dimensional cases.

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.

Nonsteady heat conduction code with radiation boundary conditions

International Nuclear Information System (INIS)

Fillo, J.A.; Benenati, R.; Powell, J.

1975-01-01

A heat-transfer model for studying the temperature build-up in graphite blankets for fusion reactors is presented. In essence, the computer code developed is for two-dimensional, nonsteady heat conduction in heterogeneous, anisotropic solids with nonuniform internal heating. Thermal radiation as well as bremsstrahlung radiation boundary conditions are included. Numerical calculations are performed for two design options by varying the wall loading, bremsstrahlung, surface layer thickness and thermal conductivity, blanket dimensions, time step and grid size. (auth)

A variable-order time-dependent neutron transport method for nuclear reactor kinetics using analytically-integrated space-time characteristics

International Nuclear Information System (INIS)

Hoffman, A. J.; Lee, J. C.

2013-01-01

A new time-dependent neutron transport method based on the method of characteristics (MOC) has been developed. Whereas most spatial kinetics methods treat time dependence through temporal discretization, this new method treats time dependence by defining the characteristics to span space and time. In this implementation regions are defined in space-time where the thickness of the region in time fulfills an analogous role to the time step in discretized methods. The time dependence of the local source is approximated using a truncated Taylor series expansion with high order derivatives approximated using backward differences, permitting the solution of the resulting space-time characteristic equation. To avoid a drastic increase in computational expense and memory requirements due to solving many discrete characteristics in the space-time planes, the temporal variation of the boundary source is similarly approximated. This allows the characteristics in the space-time plane to be represented analytically rather than discretely, resulting in an algorithm comparable in implementation and expense to one that arises from conventional time integration techniques. Furthermore, by defining the boundary flux time derivative in terms of the preceding local source time derivative and boundary flux time derivative, the need to store angularly-dependent data is avoided without approximating the angular dependence of the angular flux time derivative. The accuracy of this method is assessed through implementation in the neutron transport code DeCART. The method is employed with variable-order local source representation to model a TWIGL transient. The results demonstrate that this method is accurate and more efficient than the discretized method. (authors)

The splitting in potential Crank-Nicolson scheme with discrete transparent boundary conditions for the Schroedinger equation on a semi-infinite strip

International Nuclear Information System (INIS)

Ducomet, Bernard; Zlotnik, Alexander; Zlotnik, Ilya

2014-01-01

We consider an initial-boundary value problem for a generalized 2D time-dependent Schroedinger equation (with variable coefficients) on a semi-infinite strip. For the Crank-Nicolson-type finite-difference scheme with approximate or discrete transparent boundary conditions (TBCs), the Strang-type splitting with respect to the potential is applied. For the resulting method, the unconditional uniform in time L2-stability is proved. Due to the splitting, an effective direct algorithm using FFT is developed now to implement the method with the discrete TBC for general potential. Numerical results on the tunnel effect for rectangular barriers are included together with the detailed practical error analysis confirming nice properties of the method. (authors)

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

Homogenized boundary conditions and resonance effects in Faraday cages

Science.gov (United States)

Hewitt, I. J.

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 effect’). Taking the limit as the number of wires in the cage tends to infinity, we use the asymptotic method of multiple scales to derive continuum models for the shielding, involving homogenized boundary conditions on an effective cage boundary. We show how the resulting models depend on key cage parameters such as the size and shape of the wires, and, in the electromagnetic case, on the frequency and polarization of the incident field. In the electromagnetic case, there are resonance effects, whereby at frequencies close to the natural frequencies of the equivalent solid shell, the presence of the cage actually amplifies the incident field, rather than shielding it. By appropriately modifying the continuum model, we calculate the modified resonant frequencies, and their associated peak amplitudes. We discuss applications to radiation containment in microwave ovens and acoustic scattering by perforated shells. PMID:27279775

Homogenized boundary conditions and resonance effects in Faraday cages

Science.gov (United States)

Hewett, D. P.; Hewitt, I. J.

2016-05-01

We present a mathematical study of two-dimensional electrostatic and electromagnetic shielding by a cage of conducting wires (the so-called `Faraday cage effect'). Taking the limit as the number of wires in the cage tends to infinity, we use the asymptotic method of multiple scales to derive continuum models for the shielding, involving homogenized boundary conditions on an effective cage boundary. We show how the resulting models depend on key cage parameters such as the size and shape of the wires, and, in the electromagnetic case, on the frequency and polarization of the incident field. In the electromagnetic case, there are resonance effects, whereby at frequencies close to the natural frequencies of the equivalent solid shell, the presence of the cage actually amplifies the incident field, rather than shielding it. By appropriately modifying the continuum model, we calculate the modified resonant frequencies, and their associated peak amplitudes. We discuss applications to radiation containment in microwave ovens and acoustic scattering by perforated shells.

Providing thermal-hydraulic boundary conditions to the reactor code TINTE through a Flownex-TINTE coupling - HTR2008-58110

International Nuclear Information System (INIS)

Marais, D.; Greyvenstein, G. P.

2008-01-01

TINTE is a well established reactor analysis code which models the transient behaviour of pebble bed reactor cores but it does not include the capabilities to model a power conversion unit (PCU). This raises the issue that TINTE cannot model full system transients. One way to overcome this problem is to supply TINTE with time-dependant thermal-hydraulic boundary conditions which are obtained from PCU simulations. This study investigates a method to provide boundary conditions for the nuclear code TINTE during full system transients. This was accomplished by creating a high level interface between the systems CFD code Flownex and TINTE. An indirect coupling method is explored whereby characteristics of the PCU are matched to characteristics of the nuclear core. This method eliminates the need to iterate between the two codes. A number of transients are simulated using the coupled code and then compared against stand-alone Flownex simulations. The coupling method introduces relatively small errors when reproducing mass flow, temperature and pressure in steady state analysis, but become more pronounced when dealing with fast thermal-hydraulic transients. Decreasing the maximum time step length of TINTE reduces this problem, but increases the computational time. Copyright ASME 2008. (authors)

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

Simulating transient dynamics of the time-dependent time fractional Fokker-Planck systems

Science.gov (United States)

Kang, Yan-Mei

2016-09-01

For a physically realistic type of time-dependent time fractional Fokker-Planck (FP) equation, derived as the continuous limit of the continuous time random walk with time-modulated Boltzmann jumping weight, a semi-analytic iteration scheme based on the truncated (generalized) Fourier series is presented to simulate the resultant transient dynamics when the external time modulation is a piece-wise constant signal. At first, the iteration scheme is demonstrated with a simple time-dependent time fractional FP equation on finite interval with two absorbing boundaries, and then it is generalized to the more general time-dependent Smoluchowski-type time fractional Fokker-Planck equation. The numerical examples verify the efficiency and accuracy of the iteration method, and some novel dynamical phenomena including polarized motion orientations and periodic response death are discussed.

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.

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.

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.

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.

Effect of reactor finiteness on the boundary condition at the surface of a booster section

International Nuclear Information System (INIS)

Wassef, W.A.

1982-01-01

Effect of reactor finiteness on the boundary condition at the surface of an absorbing booster embedded in the reactor core is studied and formulated. The model used in these calculations depends on the Pl-Transport coupling technique. This method takes into consideration the rigorous neutron transport behavior inside the booster medium, while the Pl-approximation in the bulk of the scattering medium surrounding the booster which can be considered infinite in most practical applications. The neutron flux gradient parallel to the surface of the booster is considered. The geometrical configuration of the reactor core cross section is circular or rectangular. Finiteness of the reactor is introduced in the general formulation through its dimensions or buckling. Extensive numerical results are given to demonstrate the dependence of the boundary condition at the surface of the booster section on the reactor finiteness and the different physical parameters

Attainable conditions and exact invariant for the time-dependent harmonic oscillator

Energy Technology Data Exchange (ETDEWEB)

Guasti, Manuel Fernandez [Lab. de Optica Cuantica, Dep. de Fisica, Universidad A. Metropolitana, Unidad Iztapalapa, Mexico DF, Ap. Post. 55-534 (Mexico)

2006-09-22

The time-dependent oscillator equation is solved numerically for various trajectories in amplitude and phase variables. The solutions exhibit a finite time-dependent parameter whenever the squared amplitude times the derivative of the phase is invariant. If the invariant relationship does not hold, the time-dependent parameter has divergent singularities. These observations lead to the proposition that the harmonic oscillator equation with finite time-dependent parameter must have amplitude and phase solutions fulfilling the invariant relationship. Since the time-dependent parameter or the potential must be finite for any real oscillator implementation, the invariant must hold for any such physically realizable system.

Attainable conditions and exact invariant for the time-dependent harmonic oscillator

International Nuclear Information System (INIS)

Guasti, Manuel Fernandez

2006-01-01

The time-dependent oscillator equation is solved numerically for various trajectories in amplitude and phase variables. The solutions exhibit a finite time-dependent parameter whenever the squared amplitude times the derivative of the phase is invariant. If the invariant relationship does not hold, the time-dependent parameter has divergent singularities. These observations lead to the proposition that the harmonic oscillator equation with finite time-dependent parameter must have amplitude and phase solutions fulfilling the invariant relationship. Since the time-dependent parameter or the potential must be finite for any real oscillator implementation, the invariant must hold for any such physically realizable system

A Resistive Boundary Condition Enhanced DGTD Scheme for the Transient Analysis of Graphene

KAUST Repository

Li, Ping; Jiang, Li; Bagci, Hakan

2015-01-01

In this paper, the electromagnetic (EM) features of graphene are characterized by a discontinuous Galerkin timedomain (DGTD) algorithm with a resistive boundary condition (RBC). The atomically thick graphene is equivalently modeled using a RBC by regarding the graphene as an infinitesimally thin conductive sheet. To incorporate RBC into the DGTD analysis, the surface conductivity of the graphene composed of contributions from both intraband and interband terms is firstly approximated by rational basis functions using the fastrelaxation vector-fitting (FRVF) method in the Laplace-domain. Next, through the inverse Laplace transform, the corresponding time-domain matrix equations in integral can be obtained. Finally, these matrix equations are solved by time-domain finite integral technique (FIT). For elements not touching the graphene sheet, however, the well-known Runge-Kutta (RK) method is employed to solve the two first-order time-derivative Maxwell’s equations. The application of the surface boundary condition significantly alleviates the memory consuming and the limitation of time step size required by Courant-Friedrichs-Lewy (CFL) condition. To validate the proposed algorithm, various numerical examples are presented and compared with available references.

A Resistive Boundary Condition Enhanced DGTD Scheme for the Transient Analysis of Graphene

KAUST Repository

Li, Ping

2015-04-24

In this paper, the electromagnetic (EM) features of graphene are characterized by a discontinuous Galerkin timedomain (DGTD) algorithm with a resistive boundary condition (RBC). The atomically thick graphene is equivalently modeled using a RBC by regarding the graphene as an infinitesimally thin conductive sheet. To incorporate RBC into the DGTD analysis, the surface conductivity of the graphene composed of contributions from both intraband and interband terms is firstly approximated by rational basis functions using the fastrelaxation vector-fitting (FRVF) method in the Laplace-domain. Next, through the inverse Laplace transform, the corresponding time-domain matrix equations in integral can be obtained. Finally, these matrix equations are solved by time-domain finite integral technique (FIT). For elements not touching the graphene sheet, however, the well-known Runge-Kutta (RK) method is employed to solve the two first-order time-derivative Maxwell’s equations. The application of the surface boundary condition significantly alleviates the memory consuming and the limitation of time step size required by Courant-Friedrichs-Lewy (CFL) condition. To validate the proposed algorithm, various numerical examples are presented and compared with available references.

Specific features of time-dependent Psub(N) approximations in spherical geometry

International Nuclear Information System (INIS)

Peltzer, P.; Pucker, N.

1979-01-01

Approximations to the time-dependent linear transport equation can result in more serious distortions in the description of the actual physical situation than in the stationary problem. This is demonstrated in detail for the case of a neutron pulse in spherical geometry, treated within a P 1 approximation. One has to pay special attention to the singularity at r = 0 and to the effect of the boundary conditions. Effects similar to those shown here are also to be expected in connection with Psub(N) approximations of higher order. (Auth.)

Determination of a source term for a time fractional diffusion equation with an integral type over-determining condition

Directory of Open Access Journals (Sweden)

Timurkhan S. Aleroev

2013-12-01

Full Text Available We consider a linear heat equation involving a fractional derivative in time, with a nonlocal boundary condition. We determine a source term independent of the space variable, and the temperature distribution for a problem with an over-determining condition of integral type. We prove the existence and uniqueness of the solution, and its continuous dependence on the data.

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.

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

A hybrid time-domain discontinuous galerkin-boundary integral method for electromagnetic scattering analysis

KAUST Repository

Li, Ping; Shi, Yifei; Jiang, Lijun; Bagci, Hakan

2014-01-01

A scheme hybridizing discontinuous Galerkin time-domain (DGTD) and time-domain boundary integral (TDBI) methods for accurately analyzing transient electromagnetic scattering is proposed. Radiation condition is enforced using the numerical flux on the truncation boundary. The fields required by the flux are computed using the TDBI from equivalent currents introduced on a Huygens' surface enclosing the scatterer. The hybrid DGTDBI ensures that the radiation condition is mathematically exact and the resulting computation domain is as small as possible since the truncation boundary conforms to scatterer's shape and is located very close to its surface. Locally truncated domains can also be defined around each disconnected scatterer additionally reducing the size of the overall computation domain. Numerical examples demonstrating the accuracy and versatility of the proposed method are presented. © 2014 IEEE.

A hybrid time-domain discontinuous galerkin-boundary integral method for electromagnetic scattering analysis

KAUST Repository

Li, Ping

2014-05-01

A scheme hybridizing discontinuous Galerkin time-domain (DGTD) and time-domain boundary integral (TDBI) methods for accurately analyzing transient electromagnetic scattering is proposed. Radiation condition is enforced using the numerical flux on the truncation boundary. The fields required by the flux are computed using the TDBI from equivalent currents introduced on a Huygens\\' surface enclosing the scatterer. The hybrid DGTDBI ensures that the radiation condition is mathematically exact and the resulting computation domain is as small as possible since the truncation boundary conforms to scatterer\\'s shape and is located very close to its surface. Locally truncated domains can also be defined around each disconnected scatterer additionally reducing the size of the overall computation domain. Numerical examples demonstrating the accuracy and versatility of the proposed method are presented. © 2014 IEEE.

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

The fermion boundary condition and the THETA-angle in QED2

International Nuclear Information System (INIS)

Hrasko, P.

1983-09-01

The order parameter of the Schwinger model is calculated in the Euclidean functional integral approach. It is shown that the symmetry breaking angle THETA is intimately connected to the boundary condition imposed on the fermions. The transition to the Euclidean description involves both imaginary time and imaginary THETA. (author)

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

Effects of Heterogeneity and Uncertainties in Sources and Initial and Boundary Conditions on Spatiotemporal Variations of Groundwater Levels

Science.gov (United States)

Zhang, Y. K.; Liang, X.

2014-12-01

Effects of aquifer heterogeneity and uncertainties in source/sink, and initial and boundary conditions in a groundwater flow model on the spatiotemporal variations of groundwater level, h(x,t), were investigated. Analytical solutions for the variance and covariance of h(x, t) in an unconfined aquifer described by a linearized Boussinesq equation with a white noise source/sink and a random transmissivity field were derived. It was found that in a typical aquifer the error in h(x,t) in early time is mainly caused by the random initial condition and the error reduces as time goes to reach a constant error in later time. The duration during which the effect of the random initial condition is significant may last a few hundred days in most aquifers. The constant error in groundwater in later time is due to the combined effects of the uncertain source/sink and flux boundary: the closer to the flux boundary, the larger the error. The error caused by the uncertain head boundary is limited in a narrow zone near the boundary but it remains more or less constant over time. The effect of the heterogeneity is to increase the variation of groundwater level and the maximum effect occurs close to the constant head boundary because of the linear mean hydraulic gradient. The correlation of groundwater level decreases with temporal interval and spatial distance. In addition, the heterogeneity enhances the correlation of groundwater level, especially at larger time intervals and small spatial distances.

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.

Magnetospheric conditions near the equatorial footpoints of proton isotropy boundaries

Directory of Open Access Journals (Sweden)

V. A. Sergeev

2015-12-01

Full Text Available Data from a cluster of three THEMIS (Time History of Events and Macroscale Interactions during Substorms spacecraft during February–March 2009 frequently provide an opportunity to construct local data-adaptive magnetospheric models, which are suitable for the accurate mapping along the magnetic field lines at distances of 6–9 Re in the nightside magnetosphere. This allows us to map the isotropy boundaries (IBs of 30 and 80 keV protons observed by low-altitude NOAA POES (Polar Orbiting Environmental Satellites to the equatorial magnetosphere (to find the projected isotropy boundary, PIB and study the magnetospheric conditions, particularly to evaluate the ratio KIB (Rc/rc; the magnetic field curvature radius to the particle gyroradius in the neutral sheet at that point. Special care is taken to control the factors which influence the accuracy of the adaptive models and mapping. Data indicate that better accuracy of an adaptive model is achieved when the PIB distance from the closest spacecraft is as small as 1–2 Re. For this group of most accurate predictions, the spread of KIB values is still large (from 4 to 32, with the median value KIB ~13 being larger than the critical value Kcr ~ 8 expected at the inner boundary of nonadiabatic angular scattering in the current sheet. It appears that two different mechanisms may contribute to form the isotropy boundary. The group with K ~ [4,12] is most likely formed by current sheet scattering, whereas the group having KIB ~ [12,32] could be formed by the resonant scattering of low-energy protons by the electromagnetic ion-cyclotron (EMIC waves. The energy dependence of the upper K limit and close proximity of the latter event to the plasmapause locations support this conclusion. We also discuss other reasons why the K ~ 8 criterion for isotropization may fail to work, as well as a possible relationship between the two scattering mechanisms.

Moment-based boundary conditions for lattice Boltzmann simulations of natural convection in cavities

KAUST Repository

Allen, Rebecca

2016-06-29

We study a multiple relaxation time lattice Boltzmann model for natural convection with moment-based boundary conditions. The unknown primary variables of the algorithm at a boundary are found by imposing conditions directly upon hydrodynamic moments, which are then translated into conditions for the discrete velocity distribution functions. The method is formulated so that it is consistent with the second order implementation of the discrete velocity Boltzmann equations for fluid flow and temperature. Natural convection in square cavities is studied for Rayleigh numbers ranging from 103 to 108. An excellent agreement with benchmark data is observed and the flow fields are shown to converge with second order accuracy. Copyright © 2016 Inderscience Enterprises Ltd.

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

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

Directory of Open Access Journals (Sweden)

Qingkai Kong

2012-02-01

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

On the boundary conditions and validity of the neutral shielding model of a refuelling pellet

International Nuclear Information System (INIS)

Chang, C.T.

1982-02-01

By comparing the ablation time of a hydrogen pellet in a tokamak discharge with the time required for the sublimation process, the vaporization of the pellet is shown to be a dynamic phase transition - i.e. the transport of heat is due to the propagation of an evaporation front. Based on this finding, an alternative boundary condition, consistent with the energy conservation law, is formulated. Computational results utilizing the new boundary condition indicate that the ablatant near the pellet surface is hotter and less dense compared with the results which make use of the previous condition of the vanishing flux. The discrepancy between the two solutions becomes less significant once the ablatant reaches the sonic radius. The scaling law of the pellet ablation rate is unaffected by this change of boundary condition. The present analysis shows that the validity of the neutral shielding model is based mainly on the existence of a thin envelope around the pellet where strong energy absorption occurs and is insensitive to the actual vaporization process occuring at the pellet surface. (Auth.)

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

Time-dependent simulations of disk-embedded planetary atmospheres

Science.gov (United States)

Stökl, A.; Dorfi, E. A.

2014-03-01

At the early stages of evolution of planetary systems, young Earth-like planets still embedded in the protoplanetary disk accumulate disk gas gravitationally into planetary atmospheres. The established way to study such atmospheres are hydrostatic models, even though in many cases the assumption of stationarity is unlikely to be fulfilled. Furthermore, such models rely on the specification of a planetary luminosity, attributed to a continuous, highly uncertain accretion of planetesimals onto the surface of the solid core. We present for the first time time-dependent, dynamic simulations of the accretion of nebula gas into an atmosphere around a proto-planet and the evolution of such embedded atmospheres while integrating the thermal energy budget of the solid core. The spherical symmetric models computed with the TAPIR-Code (short for The adaptive, implicit RHD-Code) range from the surface of the rocky core up to the Hill radius where the surrounding protoplanetary disk provides the boundary conditions. The TAPIR-Code includes the hydrodynamics equations, gray radiative transport and convective energy transport. The results indicate that diskembedded planetary atmospheres evolve along comparatively simple outlines and in particular settle, dependent on the mass of the solid core, at characteristic surface temperatures and planetary luminosities, quite independent on numerical parameters and initial conditions. For sufficiently massive cores, this evolution ultimately also leads to runaway accretion and the formation of a gas planet.

Sparse grid spectral methods for the numerical solution of partial differential equations with periodic boundary conditions

International Nuclear Information System (INIS)

Kupka, F.

1997-11-01

This thesis deals with the extension of sparse grid techniques to spectral methods for the solution of partial differential equations with periodic boundary conditions. A review on boundary and initial-boundary value problems and a discussion on numerical resolution is used to motivate this research. Spectral methods are introduced by projection techniques, and by three model problems: the stationary and the transient Helmholtz equations, and the linear advection equation. The approximation theory on the hyperbolic cross is reviewed and its close relation to sparse grids is demonstrated. This approach extends to non-periodic problems. Various Sobolev spaces with dominant mixed derivative are introduced to provide error estimates for Fourier approximation and interpolation on the hyperbolic cross and on sparse grids by means of Sobolev norms. The theorems are immediately applicable to the stability and convergence analysis of sparse grid spectral methods. This is explicitly demonstrated for the three model problems. A variant of the von Neumann condition is introduced to simplify the stability analysis of the time-dependent model problems. The discrete Fourier transformation on sparse grids is discussed together with its software implementation. Results on numerical experiments are used to illustrate the performance of the new method with respect to the smoothness properties of each example. The potential of the method in mathematical modelling is estimated and generalizations to other sparse grid methods are suggested. The appendix includes a complete Fortran90 program to solve the linear advection equation by the sparse grid Fourier collocation method and a third-order Runge-Kutta routine for integration in time. (author)

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

An Evaluation of a Phase-Lag Boundary Condition for Francis Hydroturbine Simulations Using a Pressure-Based Solver

Science.gov (United States)

Wouden, Alex; Cimbala, John; Lewis, Bryan

2014-11-01

While the periodic boundary condition is useful for handling rotational symmetry in many axisymmetric geometries, its application fails for analysis of rotor-stator interaction (RSI) in multi-stage turbomachinery flow. The inadequacy arises from the underlying geometry where the blade counts per row differ, since the blade counts are crafted to deter the destructive harmonic forces of synchronous blade passing. Therefore, to achieve the computational advantage of modeling a single blade passage per row while preserving the integrity of the RSI, a phase-lag boundary condition is adapted to OpenFOAM® software's incompressible pressure-based solver. The phase-lag construct is accomplished through restating the implicit periodic boundary condition as a constant boundary condition that is updated at each time step with phase-shifted data from the coupled cells adjacent to the boundary. Its effectiveness is demonstrated using a typical Francis hydroturbine modeled as single- and double-passages with phase-lag boundary conditions. The evaluation of the phase-lag condition is based on the correspondence of the overall computational performance and the calculated flow parameters of the phase-lag simulations with those of a baseline full-wheel simulation. Funded in part by DOE Award Number: DE-EE0002667.

Automatic generation of boundary conditions using demons nonrigid image registration for use in 3-D modality-independent elastography.

Science.gov (United States)

Pheiffer, Thomas S; Ou, Jao J; Ong, Rowena E; Miga, Michael I

2011-09-01

Modality-independent elastography (MIE) is a method of elastography that reconstructs the elastic properties of tissue using images acquired under different loading conditions and a biomechanical model. Boundary conditions are a critical input to the algorithm and are often determined by time-consuming point correspondence methods requiring manual user input. This study presents a novel method of automatically generating boundary conditions by nonrigidly registering two image sets with a demons diffusion-based registration algorithm. The use of this method was successfully performed in silico using magnetic resonance and X-ray-computed tomography image data with known boundary conditions. These preliminary results produced boundary conditions with an accuracy of up to 80% compared to the known conditions. Demons-based boundary conditions were utilized within a 3-D MIE reconstruction to determine an elasticity contrast ratio between tumor and normal tissue. Two phantom experiments were then conducted to further test the accuracy of the demons boundary conditions and the MIE reconstruction arising from the use of these conditions. Preliminary results show a reasonable characterization of the material properties on this first attempt and a significant improvement in the automation level and viability of the method.

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.

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.

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)

Vacuum quantum effect for curved boundaries in static Robertson-Walker space-time

International Nuclear Information System (INIS)

Setare, M.R.; Sadeghi, J.

2009-01-01

The energy-momentum tensor for a massless conformally coupled scalar field in the region between two curved boundaries in k=-1 static Robertson-Walker space-time is investigated. We assume that the scalar field satisfies the Dirichlet boundary condition on the boundaries. k=-1 Robertson-Walker space is conformally related to the Rindler space, as a result we can obtain vacuum expectation values of energy-momentum tensor for conformally invariant field in Robertson-Walker space from the corresponding Rindler counterpart by the conformal transformation.

A boundary integral approach to unstable solidification

International Nuclear Information System (INIS)

Strain, J.

1989-01-01

We consider the supercooled Stefan problem with a general anisotropic curvature- and velocity-dependent boundary condition on the moving interface. We present numerical methods, based on an integral equation formulation and including a new algorithm for moving curves with curvature-dependent velocity. These methods compute a periodic interface with O(Δt) accuracy, where Δt is the time step. Previous work has been limited to short time spans and achieved slightly less than O(Δt 1/2 ) accuracy. Accurate numerical results are seen to agree with the predictions of linear stability theory. This agreement has eluded previous authors, because their numerical methods suffered from grid effects and their linear stability theory was incorrect. We study the long-time evolution of an unstable interface. Our computations exhibit the beginnings of a sidebranching instability when the boundary condition includes anisotropy and tip-splitting in the isotropic case. copyright 1989 Academic Press, Inc

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

Ozone and meteorological boundary-layer conditions at Summit, Greenland, during 3-21 June 2000

Energy Technology Data Exchange (ETDEWEB)

Helmig, D.; Boulter, J.; David, D.; Birks, J.W.; Cullen, N.J.; Steffen, K. [University of Colorado, Boulder, CO (United States). Cooperative Institute for Research in Environmental Sciences; Johnson, B.J.; Oltmans, S.J. [National Oceanic and Atmospheric Administration, Boulder, CO (United States). Climate Monitoring and Diagnostics Laboratory

2002-06-01

The temporal and spatial distributions of boundary-layer ozone were studied during June 2000 at Summit, Greenland, using surface-level measurements and vertical profiling from a tethered balloon platform. Three weeks of continuous ozone surface data, 133 meteorological vertical profile data and 82 ozone vertical profile data sets were collected from the surface to a maximum altitude of 1400 m above ground. The lower atmosphere at Summit was characterized by the prevalence of strong stable conditions with strong surface temperature inversions. These inversions reversed to neutral to slightly unstable conditions between {approx} 9.00 and 18.00 h local time with the formation of shallow mixing heights of {approx} 70-250 m above the surface. The surface ozone mixing ratio ranged from 39 to 68 ppbv and occasionally had rapid changes of up to 20 ppb in 12 h. The diurnal mean ozone mixing ratio showed diurnal trends indicating meteorological and photochemical controls of surface ozone. Vertical profiles were within the range of 37-76 ppb and showed strong stratification in the lower troposphere. A high correlation of high ozone/low water vapor air masses indicated the transport of high tropospheric/low stratospheric air into the lower boundary layer. An {approx} 0.1-3 ppb decline of the ozone mixing ratio towards the surface was frequently observed within the neutrally stable mixed layer during midday hours. These data suggest that the boundary-layer ozone mixing ratio and ozone depletion and deposition to the snowpack are influenced by the boundary-layer ozone mixing ratio and ozone depletion and deposition to the snowpack are influenced by photochemical processes and/or transport phenomena that follow diurnal dependencies. With 37 ppb of ozone being the lowest mixing ratio measured in all data no evidence was seen for the occurrence of ozone depletion episodes similar to those that have been reported within the boundary layer at coastal Arctic sites during springtime

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.

CME Simulations with Boundary Conditions Derived from Multiple Viewpoints of STEREO

Science.gov (United States)

Singh, T.; Yalim, M. S.; Pogorelov, N. V.

2017-12-01

Coronal Mass Ejections (CMEs) are major drivers of extreme space weather conditions, which is a matter of huge concern for our modern technologically dependent society. Development of numerical approaches that would reproduce CME propagation through the interplanetary space is an important step towards our capability to predict CME arrival time at Earth and their geo-effectiveness. It is also important that CMEs are propagating through a realistic, data-driven background solar wind (SW). In this study, we use a version of the flux-rope-driven Gibson-Low (GL) model to simulate CMEs. We derive inner boundary conditions for the GL flux rope model using the Graduate Cylindrical Shell (GCS) method. This method uses viewpoints from STEREO A and B, and SOHO/LASCO coronagraphs to determine the size and orientation of a CME flux rope as it starts to erupt from Sun. A flux rope created this way is inserted into an SDO/HMI vector magnetogram driven SW background obtained with the Multi-Scale Fluid-Kinetic Simulation Suite (MS-FLUKSS). Numerical results are compared with STEREO, SDO/AIA and SOHO/LASCO observations in particular in terms of the CME speed, acceleration and magnetic field structure.

Robust routing for hazardous materials transportation with conditional value-at-risk on time-dependent networks.

Science.gov (United States)

2012-11-01

New methods are proposed for mitigating risk in hazardous materials (hazmat) transportation, based on Conditional : Value-at-Risk (CVaR) measure, on time-dependent vehicular networks. While the CVaR risk measure has been : popularly used in financial...

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

DEFF Research Database (Denmark)

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

1995-01-01

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

Analysis on Forced Vibration of Thin-Wall Cylindrical Shell with Nonlinear Boundary Condition

Directory of Open Access Journals (Sweden)

Qiansheng Tang

2016-01-01

Full Text Available Forced vibration of thin-wall cylindrical shell under nonlinear boundary condition was discussed in this paper. The nonlinear boundary was modeled as supported clearance in one end of shell and the restraint was assumed as linearly elastic in the radial direction. Based on Sanders’ shell theory, Lagrange equation was utilized to derive the nonlinear governing equations of cylindrical shell. The displacements in three directions were represented by beam functions and trigonometric functions. In the study of nonlinear dynamic responses of thin-wall cylindrical shell with supported clearance under external loads, the Newmark method is used to obtain time history, frequency spectrum plot, phase portraits, Poincare section, bifurcation diagrams, and three-dimensional spectrum plot with different parameters. The effects of external loads, supported clearance, and support stiffness on nonlinear dynamics behaviors of cylindrical shell with nonlinear boundary condition were discussed.

A quasi-Lagrangian finite element method for the Navier-Stokes equations in a time-dependent domain

Science.gov (United States)

Lozovskiy, Alexander; Olshanskii, Maxim A.; Vassilevski, Yuri V.

2018-05-01

The paper develops a finite element method for the Navier-Stokes equations of incompressible viscous fluid in a time-dependent domain. The method builds on a quasi-Lagrangian formulation of the problem. The paper provides stability and convergence analysis of the fully discrete (finite-difference in time and finite-element in space) method. The analysis does not assume any CFL time-step restriction, it rather needs mild conditions of the form $\\Delta t\\le C$, where $C$ depends only on problem data, and $h^{2m_u+2}\\le c\\,\\Delta t$, $m_u$ is polynomial degree of velocity finite element space. Both conditions result from a numerical treatment of practically important non-homogeneous boundary conditions. The theoretically predicted convergence rate is confirmed by a set of numerical experiments. Further we apply the method to simulate a flow in a simplified model of the left ventricle of a human heart, where the ventricle wall dynamics is reconstructed from a sequence of contrast enhanced Computed Tomography images.

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.

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

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.

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

International Nuclear Information System (INIS)

Mohammed, M.H.H.

2012-01-01

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

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

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.

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

Imposition of Dirichlet Boundary Conditions in Element Free Galerkin Method through an Object-Oriented Implementation

Directory of Open Access Journals (Sweden)

Samira Hosseini

Full Text Available Abstract One of the main drawbacks of Element Free Galerkin (EFG method is its dependence on moving least square shape functions which don’t satisfy the Kronecker Delta property, so in this method it’s not possible to apply Dirichlet boundary conditions directly. The aim of the present paper is to discuss different aspects of three widely used methods of applying Dirichlet boundary conditions in EFG method, called Lagrange multipliers, penalty method, and coupling with finite element method. Numerical simulations are presented to compare the results of these methods form the perspective of accuracy, convergence and computational expense. These methods have been implemented in an object oriented programing environment, called INSANE, and the results are presented and compared with the analytical solutions.

Breaking generalized covariance, classical renormalization, and boundary conditions from superpotentials

International Nuclear Information System (INIS)

Livshits, Gideon I.

2014-01-01

Superpotentials offer a direct means of calculating conserved charges associated with the asymptotic symmetries of space-time. Yet superpotentials have been plagued with inconsistencies, resulting in nonphysical or incongruent values for the mass, angular momentum, and energy loss due to radiation. The approach of Regge and Teitelboim, aimed at a clear Hamiltonian formulation with a boundary, and its extension to the Lagrangian formulation by Julia and Silva have resolved these issues, and have resulted in a consistent, well-defined and unique variational equation for the superpotential, thereby placing it on a firm footing. A hallmark solution of this equation is the KBL superpotential obtained from the first-order Lovelock Lagrangian. Nevertheless, here we show that these formulations are still insufficient for Lovelock Lagrangians of higher orders. We present a paradox, whereby the choice of fields affects the superpotential for equivalent on-shell dynamics. We offer two solutions to this paradox: either the original Lagrangian must be effectively renormalized, or that boundary conditions must be imposed, so that space-time be asymptotically maximally symmetric. Non-metricity is central to this paradox, and we show how quadratic non-metricity in the bulk of space-time contributes to the conserved charges on the boundary, where it vanishes identically. This is a realization of the gravitational Higgs mechanism, proposed by Percacci, where the non-metricity is the analogue of the Goldstone boson

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.

Assessment of a PML Boundary Condition for Simulating an MRI Radio Frequency Coil

Directory of Open Access Journals (Sweden)

Yunsuo Duan

2008-01-01

Full Text Available Computational methods such as the finite difference time domain (FDTD play an important role in simulating radiofrequency (RF coils used in magnetic resonance imaging (MRI. The choice of absorbing boundary conditions affects the final outcome of such studies. We have used FDTD to assess the Berenger's perfectly matched layer (PML as an absorbing boundary condition for computation of the resonance patterns and electromagnetic fields of RF coils. We first experimentally constructed a high-pass birdcage head coil, measured its resonance pattern, and used it to acquire proton (1H phantom MRI images. We then computed the resonance pattern and B1 field of the coil using FDTD with a PML as an absorbing boundary condition. We assessed the accuracy and efficiency of PML by adjusting the parameters of the PML and comparing the calculated results with measured ones. The optimal PML parameters that produce accurate (comparable to the experimental findings FDTD calculations are then provided for the birdcage head coil operating at 127.72 MHz, the Larmor frequency of 1H at 3 Tesla (T.

Time-dependent deformation of concrete under multiaxial stress conditions. Final report

International Nuclear Information System (INIS)

McDonald, J.E.

1975-10-01

An investigation was made on the time-dependent deformation behavior of concrete in the presence of temperature, moisture, and loading conditions similar to those encountered in a prestressed concrete reactor vessel (PCRV). This investigation encompassed one concrete strength (6000 psi at 28 days), three aggregate types (chert, limestone, and graywacke), one cement (Type II), two types of specimens (as-cast and air-dried), two levels of temperature during test (73 and 150 0 F), and four types of loading (uniaxial, hydrostatic, biaxial, and triaxial). This effort was intended primarily as a data report; the experimental procedures and results are presented in detail. A comprehensive evaluation of the effects of various parameters and their interactions on the behavior of concrete is not included. However, a number of general comparisons were made concerning the effect of the various test conditions on concrete behavior. Based on this limited evaluation of the data, general conclusions and recommendations for additional work were formulated

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

Directory of Open Access Journals (Sweden)

Long Yuhua

2017-12-01

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

Transient hydrodynamic finite-size effects in simulations under periodic boundary conditions

Science.gov (United States)

Asta, Adelchi J.; Levesque, Maximilien; Vuilleumier, Rodolphe; Rotenberg, Benjamin

2017-06-01

We use lattice-Boltzmann and analytical calculations to investigate transient hydrodynamic finite-size effects induced by the use of periodic boundary conditions. These effects are inevitable in simulations at the molecular, mesoscopic, or continuum levels of description. We analyze the transient response to a local perturbation in the fluid and obtain the local velocity correlation function via linear response theory. This approach is validated by comparing the finite-size effects on the steady-state velocity with the known results for the diffusion coefficient. We next investigate the full time dependence of the local velocity autocorrelation function. We find at long times a crossover between the expected t-3 /2 hydrodynamic tail and an oscillatory exponential decay, and study the scaling with the system size of the crossover time, exponential rate and amplitude, and oscillation frequency. We interpret these results from the analytic solution of the compressible Navier-Stokes equation for the slowest modes, which are set by the system size. The present work not only provides a comprehensive analysis of hydrodynamic finite-size effects in bulk fluids, which arise regardless of the level of description and simulation algorithm, but also establishes the lattice-Boltzmann method as a suitable tool to investigate such effects in general.

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

Integrable Time-Dependent Quantum Hamiltonians

Science.gov (United States)

Sinitsyn, Nikolai A.; Yuzbashyan, Emil A.; Chernyak, Vladimir Y.; Patra, Aniket; Sun, Chen

2018-05-01

We formulate a set of conditions under which the nonstationary Schrödinger equation with a time-dependent Hamiltonian is exactly solvable analytically. The main requirement is the existence of a non-Abelian gauge field with zero curvature in the space of system parameters. Known solvable multistate Landau-Zener models satisfy these conditions. Our method provides a strategy to incorporate time dependence into various quantum integrable models while maintaining their integrability. We also validate some prior conjectures, including the solution of the driven generalized Tavis-Cummings model.

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.

Validation of a CFD code for Unsteady Flows with cyclic boundary Conditions

International Nuclear Information System (INIS)

Kim, Jong-Tae; Kim, Sang-Baik; Lee, Won-Jae

2006-01-01

Currently Lilac code is under development to analyze thermo-hydraulics of a high-temperature gas-cooled reactor (GCR). Interesting thermo-hydraulic phenomena in a nuclear reactor are usually unsteady and turbulent. The analysis of the unsteady flows by using a three dimension CFD code is time-consuming if the flow domain is very large. Hopefully, flow domains commonly encountered in the nuclear thermo-hydraulics is periodic. So it is better to use the geometrical characteristics in order to reduce the computational resources. To get the benefits from reducing the computation domains especially for the calculations of unsteady flows, the cyclic boundary conditions are implemented in the parallelized CFD code LILAC. In this study, the parallelized cyclic boundary conditions are validated by solving unsteady laminar and turbulent flows past a circular cylinder

Performance advantages of CPML over UPML absorbing boundary conditions in FDTD algorithm

Science.gov (United States)

Gvozdic, Branko D.; Djurdjevic, Dusan Z.

2017-01-01

Implementation of absorbing boundary condition (ABC) has a very important role in simulation performance and accuracy in finite difference time domain (FDTD) method. The perfectly matched layer (PML) is the most efficient type of ABC. The aim of this paper is to give detailed insight in and discussion of boundary conditions and hence to simplify the choice of PML used for termination of computational domain in FDTD method. In particular, we demonstrate that using the convolutional PML (CPML) has significant advantages in terms of implementation in FDTD method and reducing computer resources than using uniaxial PML (UPML). An extensive number of numerical experiments has been performed and results have shown that CPML is more efficient in electromagnetic waves absorption. Numerical code is prepared, several problems are analyzed and relative error is calculated and 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.

Global existence and exponential growth for a viscoelastic wave equation with dynamic boundary conditions

KAUST Repository

Gerbi, Sté phane; Said-Houari, Belkacem

2013-01-01

The goal of this work is to study a model of the wave equation with dynamic boundary conditions and a viscoelastic term. First, applying the Faedo-Galerkin method combined with the fixed point theorem, we show the existence and uniqueness of a local in time solution. Second, we show that under some restrictions on the initial data, the solution continues to exist globally in time. On the other hand, if the interior source dominates the boundary damping, then the solution is unbounded and grows as an exponential function. In addition, in the absence of the strong damping, then the solution ceases to exist and blows up in finite time.

Global existence and exponential growth for a viscoelastic wave equation with dynamic boundary conditions

KAUST Repository

Gerbi, Stéphane

2013-01-15

The goal of this work is to study a model of the wave equation with dynamic boundary conditions and a viscoelastic term. First, applying the Faedo-Galerkin method combined with the fixed point theorem, we show the existence and uniqueness of a local in time solution. Second, we show that under some restrictions on the initial data, the solution continues to exist globally in time. On the other hand, if the interior source dominates the boundary damping, then the solution is unbounded and grows as an exponential function. In addition, in the absence of the strong damping, then the solution ceases to exist and blows up in finite time.

Determinantal Representation of the Time-Dependent Stationary Correlation Function for the Totally Asymmetric Simple Exclusion Model

Directory of Open Access Journals (Sweden)

Nikolay M. Bogoliubov

2009-04-01

Full Text Available The basic model of the non-equilibrium low dimensional physics the so-called totally asymmetric exclusion process is related to the 'crystalline limit' (q → ∞ of the SU_q(2 quantum algebra. Using the quantum inverse scattering method we obtain the exact expression for the time-dependent stationary correlation function of the totally asymmetric simple exclusion process on a one dimensional lattice with the periodic boundary conditions.

Bayesian techniques for fatigue life prediction and for inference in linear time dependent PDEs

KAUST Repository

Scavino, Marco

2016-01-08

In this talk we introduce first the main characteristics of a systematic statistical approach to model calibration, model selection and model ranking when stress-life data are drawn from a collection of records of fatigue experiments. Focusing on Bayesian prediction assessment, we consider fatigue-limit models and random fatigue-limit models under different a priori assumptions. In the second part of the talk, we present a hierarchical Bayesian technique for the inference of the coefficients of time dependent linear PDEs, under the assumption that noisy measurements are available in both the interior of a domain of interest and from boundary conditions. We present a computational technique based on the marginalization of the contribution of the boundary parameters and apply it to inverse heat conduction problems.

Double absorbing boundaries for finite-difference time-domain electromagnetics

Energy Technology Data Exchange (ETDEWEB)

LaGrone, John, E-mail: jlagrone@smu.edu; Hagstrom, Thomas, E-mail: thagstrom@smu.edu

2016-12-01

We describe the implementation of optimal local radiation boundary condition sequences for second order finite difference approximations to Maxwell's equations and the scalar wave equation using the double absorbing boundary formulation. Numerical experiments are presented which demonstrate that the design accuracy of the boundary conditions is achieved and, for comparable effort, exceeds that of a convolution perfectly matched layer with reasonably chosen parameters. An advantage of the proposed approach is that parameters can be chosen using an accurate a priori error bound.

Computational complexity of time-dependent density functional theory

International Nuclear Information System (INIS)

Whitfield, J D; Yung, M-H; Tempel, D G; Aspuru-Guzik, A; Boixo, S

2014-01-01

Time-dependent density functional theory (TDDFT) is rapidly emerging as a premier method for solving dynamical many-body problems in physics and chemistry. The mathematical foundations of TDDFT are established through the formal existence of a fictitious non-interacting system (known as the Kohn–Sham system), which can reproduce the one-electron reduced probability density of the actual system. We build upon these works and show that on the interior of the domain of existence, the Kohn–Sham system can be efficiently obtained given the time-dependent density. We introduce a V-representability parameter which diverges at the boundary of the existence domain and serves to quantify the numerical difficulty of constructing the Kohn-Sham potential. For bounded values of V-representability, we present a polynomial time quantum algorithm to generate the time-dependent Kohn–Sham potential with controllable error bounds. (paper)

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.

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.

Simulating transient dynamics of the time-dependent time fractional Fokker–Planck systems

Energy Technology Data Exchange (ETDEWEB)

Kang, Yan-Mei, E-mail: ymkang@mail.xjtu.edu.cn

2016-09-16

For a physically realistic type of time-dependent time fractional Fokker–Planck (FP) equation, derived as the continuous limit of the continuous time random walk with time-modulated Boltzmann jumping weight, a semi-analytic iteration scheme based on the truncated (generalized) Fourier series is presented to simulate the resultant transient dynamics when the external time modulation is a piece-wise constant signal. At first, the iteration scheme is demonstrated with a simple time-dependent time fractional FP equation on finite interval with two absorbing boundaries, and then it is generalized to the more general time-dependent Smoluchowski-type time fractional Fokker–Planck equation. The numerical examples verify the efficiency and accuracy of the iteration method, and some novel dynamical phenomena including polarized motion orientations and periodic response death are discussed. - Highlights: • An iteration method is proposed for the transient dynamics of time-dependent time fractional Fokker–Planck equations. • The method is based on Fourier Series solution and the multi-step transition probability formula. • With the time-modulated subdiffusion on finite interval as example, the polarized motion orientation is disclosed. • With the time-modulated subdiffusion within a confined potential as example, the death of dynamic response is observed.

Simulating transient dynamics of the time-dependent time fractional Fokker–Planck systems

International Nuclear Information System (INIS)

Kang, Yan-Mei

2016-01-01

For a physically realistic type of time-dependent time fractional Fokker–Planck (FP) equation, derived as the continuous limit of the continuous time random walk with time-modulated Boltzmann jumping weight, a semi-analytic iteration scheme based on the truncated (generalized) Fourier series is presented to simulate the resultant transient dynamics when the external time modulation is a piece-wise constant signal. At first, the iteration scheme is demonstrated with a simple time-dependent time fractional FP equation on finite interval with two absorbing boundaries, and then it is generalized to the more general time-dependent Smoluchowski-type time fractional Fokker–Planck equation. The numerical examples verify the efficiency and accuracy of the iteration method, and some novel dynamical phenomena including polarized motion orientations and periodic response death are discussed. - Highlights: • An iteration method is proposed for the transient dynamics of time-dependent time fractional Fokker–Planck equations. • The method is based on Fourier Series solution and the multi-step transition probability formula. • With the time-modulated subdiffusion on finite interval as example, the polarized motion orientation is disclosed. • With the time-modulated subdiffusion within a confined potential as example, the death of dynamic response is observed.

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

Implicit time-dependent finite different algorithm for quench simulation

International Nuclear Information System (INIS)

Koizumi, Norikiyo; Takahashi, Yoshikazu; Tsuji, Hiroshi

1994-12-01

A magnet in a fusion machine has many difficulties in its application because of requirement of a large operating current, high operating field and high breakdown voltage. A cable-in-conduit (CIC) conductor is the best candidate to overcome these difficulties. However, there remained uncertainty in a quench event in the cable-in-conduit conductor because of a difficulty to analyze a fluid dynamics equation. Several scientists, then, developed the numerical code for the quench simulation. However, most of them were based on an explicit time-dependent finite difference scheme. In this scheme, a discrete time increment is strictly restricted by CFL (Courant-Friedrichs-Lewy) condition. Therefore, long CPU time was consumed for the quench simulation. Authors, then, developed a new quench simulation code, POCHI1, which is based on an implicit time dependent scheme. In POCHI1, the fluid dynamics equation is linearlized according to a procedure applied by Beam and Warming and then, a tridiagonal system can be offered. Therefore, no iteration is necessary to solve the fluid dynamics equation. This leads great reduction of the CPU time. Also, POCHI1 can cope with non-linear boundary condition. In this study, comparison with experimental results was carried out. The normal zone propagation behavior was investigated in two samples of CIC conductors which had different hydraulic diameters. The measured and simulated normal zone propagation length showed relatively good agreement. However, the behavior of the normal voltage shows a little disagreement. These results indicate necessity to improve the treatment of the heat transfer coefficient in the turbulent flow region and the electric resistivity of the copper stabilizer in high temperature and high field region. (author)

Design of broadband time-domain impedance boundary conditions using the oscillatory-diffusive representation of acoustical models.

Science.gov (United States)

Monteghetti, Florian; Matignon, Denis; Piot, Estelle; Pascal, Lucas

2016-09-01

A methodology to design broadband time-domain impedance boundary conditions (TDIBCs) from the analysis of acoustical models is presented. The derived TDIBCs are recast exclusively as first-order differential equations, well-suited for high-order numerical simulations. Broadband approximations are yielded from an elementary linear least squares optimization that is, for most models, independent of the absorbing material geometry. This methodology relies on a mathematical technique referred to as the oscillatory-diffusive (or poles and cuts) representation, and is applied to a wide range of acoustical models, drawn from duct acoustics and outdoor sound propagation, which covers perforates, semi-infinite ground layers, as well as cavities filled with a porous medium. It is shown that each of these impedance models leads to a different TDIBC. Comparison with existing numerical models, such as multi-pole or extended Helmholtz resonator, provides insights into their suitability. Additionally, the broadly-applicable fractional polynomial impedance models are analyzed using fractional calculus.

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

First-passage time asymptotics over moving boundaries for random walk bridges

NARCIS (Netherlands)

Sloothaak, F.; Zwart, B.; Wachtel, V.

2017-01-01

We study the asymptotic tail probability of the first-passage time over a moving boundary for a random walk conditioned to return to zero, where the increments of the random walk have finite variance. Typically, the asymptotic tail behavior may be described through a regularly varying function with

Compressible turbulent channel flow with impedance boundary conditions

Science.gov (United States)

Scalo, Carlo; Bodart, Julien; Lele, Sanjiva K.

2015-03-01

We have performed large-eddy simulations of isothermal-wall compressible turbulent channel flow with linear acoustic impedance boundary conditions (IBCs) for the wall-normal velocity component and no-slip conditions for the tangential velocity components. Three bulk Mach numbers, Mb = 0.05, 0.2, 0.5, with a fixed bulk Reynolds number, Reb = 6900, have been investigated. For each Mb, nine different combinations of IBC settings were tested, in addition to a reference case with impermeable walls, resulting in a total of 30 simulations. The adopted numerical coupling strategy allows for a spatially and temporally consistent imposition of physically realizable IBCs in a fully explicit compressible Navier-Stokes solver. The IBCs are formulated in the time domain according to Fung and Ju ["Time-domain impedance boundary conditions for computational acoustics and aeroacoustics," Int. J. Comput. Fluid Dyn. 18(6), 503-511 (2004)]. The impedance adopted is a three-parameter damped Helmholtz oscillator with resonant angular frequency, ωr, tuned to the characteristic time scale of the large energy-containing eddies. The tuning condition, which reads ωr = 2πMb (normalized with the speed of sound and channel half-width), reduces the IBCs' free parameters to two: the damping ratio, ζ, and the resistance, R, which have been varied independently with values, ζ = 0.5, 0.7, 0.9, and R = 0.01, 0.10, 1.00, for each Mb. The application of the tuned IBCs results in a drag increase up to 300% for Mb = 0.5 and R = 0.01. It is shown that for tuned IBCs, the resistance, R, acts as the inverse of the wall-permeability and that varying the damping ratio, ζ, has a secondary effect on the flow response. Typical buffer-layer turbulent structures are completely suppressed by the application of tuned IBCs. A new resonance buffer layer is established characterized by large spanwise-coherent Kelvin-Helmholtz rollers, with a well-defined streamwise wavelength λx, traveling downstream with

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)

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)

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)

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

DCXE, Time-Dependent Xe Diffusion in Non-Multiplying Slab

International Nuclear Information System (INIS)

Horikami, K.; Kawasaki, S.

1969-01-01

1 - Nature of physical problem solved: Programme DCXE was designed for the analysis of the diffusion phenomena of xenon in non-multiplying media (slab geometry) and solves the Fick's second diffusion equation with boundary conditions. Initial distribution of xenon in the media at time 0 and 0 value of distribution at both ends of media at any positive time. 2 - Method of solution: Difference approximation is used to solve Fick's second diffusion equation. Strict stability condition is chosen between time and spatial intervals

DCXE, Time-Dependent Xe Diffusion in Non-Multiplying Slab

Energy Technology Data Exchange (ETDEWEB)

Horikami, K; Kawasaki, S [Japan Atomic Energy Research Institute, Tokai, Ibaraki (Japan)

1969-01-01

1 - Nature of physical problem solved: Programme DCXE was designed for the analysis of the diffusion phenomena of xenon in non-multiplying media (slab geometry) and solves the Fick's second diffusion equation with boundary conditions. Initial distribution of xenon in the media at time 0 and 0 value of distribution at both ends of media at any positive time. 2 - Method of solution: Difference approximation is used to solve Fick's second diffusion equation. Strict stability condition is chosen between time and spatial intervals.

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

An MHD simulation model of time-dependent global solar corona with temporally varying solar-surface magnetic field maps

Science.gov (United States)

Hayashi, K.

2013-11-01

We present a model of a time-dependent three-dimensional magnetohydrodynamics simulation of the sub-Alfvenic solar corona and super-Alfvenic solar wind with temporally varying solar-surface boundary magnetic field data. To (i) accommodate observational data with a somewhat arbitrarily evolving solar photospheric magnetic field as the boundary value and (ii) keep the divergence-free condition, we developed a boundary model, here named Confined Differential Potential Field model, that calculates the horizontal components of the magnetic field, from changes in the vertical component, as a potential field confined in a thin shell. The projected normal characteristic method robustly simulates the solar corona and solar wind, in response to the temporal variation of the boundary Br. We conduct test MHD simulations for two periods, from Carrington Rotation number 2009 to 2010 and from Carrington Rotation 2074 to 2075 at solar maximum and minimum of Cycle 23, respectively. We obtained several coronal features that a fixed boundary condition cannot yield, such as twisted magnetic field lines at the lower corona and the transition from an open-field coronal hole to a closed-field streamer. We also obtained slight improvements of the interplanetary magnetic field, including the latitudinal component, at Earth.

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)

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.

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.

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

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

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

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.

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)

Uncertainty Propagation Analysis for the Monte Carlo Time-Dependent Simulations

International Nuclear Information System (INIS)

Shaukata, Nadeem; Shim, Hyung Jin

2015-01-01

In this paper, a conventional method to control the neutron population for super-critical systems is implemented. Instead of considering the cycles, the simulation is divided in time intervals. At the end of each time interval, neutron population control is applied on the banked neutrons. Randomly selected neutrons are discarded, until the size of neutron population matches the initial neutron histories at the beginning of time simulation. A time-dependent simulation mode has also been implemented in the development version of SERPENT 2 Monte Carlo code. In this mode, sequential population control mechanism has been proposed for modeling of prompt super-critical systems. A Monte Carlo method has been properly used in TART code for dynamic criticality calculations. For super-critical systems, the neutron population is allowed to grow over a period of time. The neutron population is uniformly combed to return it to the neutron population started with at the beginning of time boundary. In this study, conventional time-dependent Monte Carlo (TDMC) algorithm is implemented. There is an exponential growth of neutron population in estimation of neutron density tally for super-critical systems and the number of neutrons being tracked exceed the memory of the computer. In order to control this exponential growth at the end of each time boundary, a conventional time cut-off controlling population strategy is included in TDMC. A scale factor is introduced to tally the desired neutron density at the end of each time boundary. The main purpose of this paper is the quantification of uncertainty propagation in neutron densities at the end of each time boundary for super-critical systems. This uncertainty is caused by the uncertainty resulting from the introduction of scale factor. The effectiveness of TDMC is examined for one-group infinite homogeneous problem (the rod model) and two-group infinite homogeneous problem. The desired neutron density is tallied by the introduction of

Uncertainty Propagation Analysis for the Monte Carlo Time-Dependent Simulations

Energy Technology Data Exchange (ETDEWEB)

Shaukata, Nadeem; Shim, Hyung Jin [Seoul National University, Seoul (Korea, Republic of)

2015-10-15

In this paper, a conventional method to control the neutron population for super-critical systems is implemented. Instead of considering the cycles, the simulation is divided in time intervals. At the end of each time interval, neutron population control is applied on the banked neutrons. Randomly selected neutrons are discarded, until the size of neutron population matches the initial neutron histories at the beginning of time simulation. A time-dependent simulation mode has also been implemented in the development version of SERPENT 2 Monte Carlo code. In this mode, sequential population control mechanism has been proposed for modeling of prompt super-critical systems. A Monte Carlo method has been properly used in TART code for dynamic criticality calculations. For super-critical systems, the neutron population is allowed to grow over a period of time. The neutron population is uniformly combed to return it to the neutron population started with at the beginning of time boundary. In this study, conventional time-dependent Monte Carlo (TDMC) algorithm is implemented. There is an exponential growth of neutron population in estimation of neutron density tally for super-critical systems and the number of neutrons being tracked exceed the memory of the computer. In order to control this exponential growth at the end of each time boundary, a conventional time cut-off controlling population strategy is included in TDMC. A scale factor is introduced to tally the desired neutron density at the end of each time boundary. The main purpose of this paper is the quantification of uncertainty propagation in neutron densities at the end of each time boundary for super-critical systems. This uncertainty is caused by the uncertainty resulting from the introduction of scale factor. The effectiveness of TDMC is examined for one-group infinite homogeneous problem (the rod model) and two-group infinite homogeneous problem. The desired neutron density is tallied by the introduction of

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

Bifurcation of positive solutions to scalar reaction-diffusion equations with nonlinear boundary condition

Science.gov (United States)

Liu, Ping; Shi, Junping

2018-01-01

The bifurcation of non-trivial steady state solutions of a scalar reaction-diffusion equation with nonlinear boundary conditions is considered using several new abstract bifurcation theorems. The existence and stability of positive steady state solutions are proved using a unified approach. The general results are applied to a Laplace equation with nonlinear boundary condition and bistable nonlinearity, and an elliptic equation with superlinear nonlinearity and sublinear boundary conditions.

Three Dimensional Energy Transmitting Boundary in the Time Domain

Directory of Open Access Journals (Sweden)

Naohiro eNakamura

2015-11-01

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

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.

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

Science.gov (United States)

Reimer, Ashton S.; Cheviakov, Alexei F.

2013-03-01

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

Time-dependent angularly averaged inverse transport

International Nuclear Information System (INIS)

Bal, Guillaume; Jollivet, Alexandre

2009-01-01

This paper concerns the reconstruction of the absorption and scattering parameters in a time-dependent linear transport equation from knowledge of angularly averaged measurements performed at the boundary of a domain of interest. Such measurement settings find applications in medical and geophysical imaging. We show that the absorption coefficient and the spatial component of the scattering coefficient are uniquely determined by such measurements. We obtain stability results on the reconstruction of the absorption and scattering parameters with respect to the measured albedo operator. The stability results are obtained by a precise decomposition of the measurements into components with different singular behavior in the time domain

Quasilocal conservation laws in XXZ spin-1/2 chains: Open, periodic and twisted boundary conditions

Directory of Open Access Journals (Sweden)

Tomaž Prosen

2014-09-01

Full Text Available A continuous family of quasilocal exact conservation laws is constructed in the anisotropic Heisenberg (XXZ spin-1/2 chain for periodic (or twisted boundary conditions and for a set of commensurate anisotropies densely covering the entire easy plane interaction regime. All local conserved operators follow from the standard (Hermitian transfer operator in fundamental representation (with auxiliary spin s=1/2, and are all even with respect to a spin flip operation. However, the quasilocal family is generated by differentiation of a non-Hermitian highest weight transfer operator with respect to a complex auxiliary spin representation parameter s and includes also operators of odd parity. For a finite chain with open boundaries the time derivatives of quasilocal operators are not strictly vanishing but result in operators localized near the boundaries of the chain. We show that a simple modification of the non-Hermitian transfer operator results in exactly conserved, but still quasilocal operators for periodic or generally twisted boundary conditions. As an application, we demonstrate that implementing the new exactly conserved operator family for estimating the high-temperature spin Drude weight results, in the thermodynamic limit, in exactly the same lower bound as for almost conserved family and open boundaries. Under the assumption that the bound is saturating (suggested by agreement with previous thermodynamic Bethe ansatz calculations we propose a simple explicit construction of infinite time averages of local operators such as the spin current.

FULL GPU Implementation of Lattice-Boltzmann Methods with Immersed Boundary Conditions for Fast Fluid Simulations

Directory of Open Access Journals (Sweden)

G Boroni

2017-03-01

Full Text Available Lattice Boltzmann Method (LBM has shown great potential in fluid simulations, but performance issues and difficulties to manage complex boundary conditions have hindered a wider application. The upcoming of Graphic Processing Units (GPU Computing offered a possible solution for the performance issue, and methods like the Immersed Boundary (IB algorithm proved to be a flexible solution to boundaries. Unfortunately, the implicit IB algorithm makes the LBM implementation in GPU a non-trivial task. This work presents a fully parallel GPU implementation of LBM in combination with IB. The fluid-boundary interaction is implemented via GPU kernels, using execution configurations and data structures specifically designed to accelerate each code execution. Simulations were validated against experimental and analytical data showing good agreement and improving the computational time. Substantial reductions of calculation rates were achieved, lowering down the required time to execute the same model in a CPU to about two magnitude orders.

Time-dependent chemo-electro-mechanical behavior of hydrogel-based structures

Science.gov (United States)

Leichsenring, Peter; Wallmersperger, Thomas

2018-03-01

Charged hydrogels are ionic polymer gels and belong to the class of smart materials. These gels are multiphasic materials which consist of a solid phase, a fluid phase and an ionic phase. Due to the presence of bound charges these materials are stimuli-responsive to electrical or chemical loads. The application of electrical or chemical stimuli as well as mechanical loads lead to a viscoelastic response. On the macroscopic scale, the response is governed by a local reversible release or absorption of water which, in turn, leads to a local decrease or increase of mass and a respective volume change. Furthermore, the chemo-electro-mechanical equilibrium of a hydrogel depends on the chemical composition of the gel and the surrounding solution bath. Due to the presence of bound charges in the hydrogel, this system can be understood as an osmotic cell where differences in the concentration of mobile ions in the gel and solution domain lead to an osmotic pressure difference. In the present work, a continuum-based numerical model is presented in order to describe the time-dependent swelling behavior of hydrogels. The numerical model is based on the Theory of Porous Media and captures the fluid-solid, fluid-ion and ion-ion interactions. As a direct consequence of the chemo-electro-mechanical equilibrium, the corresponding boundary conditions are defined following the equilibrium conditions. For the interaction of the hydrogel with surrounding mechanical structures, also respective jump condtions are formulated. Finaly, numerical results of the time-dependent behavior of a hydrogel-based chemo-sensor will be presented.

A study on the boundary condition for analysis of bio-heat equation according to light irradiation

Energy Technology Data Exchange (ETDEWEB)

Ko, Dong Guk; Bae, Sung Woo; Im, Ik Tae [Chunbuk Natinal University, Junju (Korea, Republic of)

2015-11-15

In this study, the temperature change in an imitational biological tissue, when its surface is irradiated with bio-light, was measured by experiments. Using the experimental data, an equation for temperature as a function of time was developed in order to use it as a boundary condition in numerical studies for the model. The temperature profile was measured along the depth for several wavelengths and distances of the light source from the tissue. It was found that the temperature of the tissue increased with increasing wavelength and irradiation time; however, the difference in the temperatures with red light and near infrared light was not large. The numerical analysis results obtained by using the developed equation as boundary condition show good agreement with the measured temperatures.

Temperature dependence of the dynamics of zone boundary phonons in ZnO:Li

Science.gov (United States)

Yadav, Harish Kumar; Sreenivas, K.; Gupta, Vinay; Katiyar, R. S.

2008-12-01

Investigations of zone boundary phonons in ZnO:Li system (Li concentration: 10%) and their dynamics with temperature are reported. Additional modes at 127, 157, and 194 cm-1 are observed and assigned to zone boundary phonons at critical point M in the Brillouin zone [J. M. Calleja and M. Cardona, Phys. Rev. B 16, 3753 (1977)] due to breakdown of crystal translational symmetry with Li incorporation in ZnO. Anharmonicity in peak frequency and linewidth of the zone boundary phonons in a temperature range from 100 to 1000 K is also analyzed taking into account the decay of zone boundary phonons into three- and four-phonon modes (cubic and quadratic anharmonicities). The anharmonic behavior of peak frequency is found to be feebly dependent on three-phonon decay process but thermal expansion of lattice together with four-phonon decay process appropriately defines the temperature dependence. Linewidths, however, follow the simple four-phonon decay mechanism. E2(low) mode, on the other hand, shows a linear temperature dependency and therefore follows a three-phonon decay channel. The calculated values of phonon lifetimes at 100 K for the 127, 157, 194 cm-1, and E2(low) modes are 8.23, 6.54, 5.32, and 11.39 ps. Decay of the zone boundary phonon modes compared to E2(low) mode reveals that dopant induced disorder has a strong temperature dependency.

Reward-timing-dependent bidirectional modulation of cortical microcircuits during optical single-neuron operant conditioning.

Science.gov (United States)

Hira, Riichiro; Ohkubo, Fuki; Masamizu, Yoshito; Ohkura, Masamichi; Nakai, Junichi; Okada, Takashi; Matsuzaki, Masanori

2014-11-24

Animals rapidly adapt to environmental change. To reveal how cortical microcircuits are rapidly reorganized when an animal recognizes novel reward contingency, we conduct two-photon calcium imaging of layer 2/3 motor cortex neurons in mice and simultaneously reinforce the activity of a single cortical neuron with water delivery. Here we show that when the target neuron is not relevant to a pre-trained forelimb movement, the mouse increases the target neuron activity and the number of rewards delivered during 15-min operant conditioning without changing forelimb movement behaviour. The reinforcement bidirectionally modulates the activity of subsets of non-target neurons, independent of distance from the target neuron. The bidirectional modulation depends on the relative timing between the reward delivery and the neuronal activity, and is recreated by pairing reward delivery and photoactivation of a subset of neurons. Reward-timing-dependent bidirectional modulation may be one of the fundamental processes in microcircuit reorganization for rapid adaptation.

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

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

Mixed problems for linear symmetric hyperbolic systems with characteristic boundary conditions

International Nuclear Information System (INIS)

Secchi, P.

1994-01-01

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

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

Second-order domain derivative of normal-dependent boundary integrals

KAUST Repository

Balzer, Jonathan

2010-03-17

Numerous reconstruction tasks in (optical) surface metrology allow for a variational formulation. The occurring boundary integrals may be interpreted as shape functions. The paper is concerned with the second-order analysis of such functions. Shape Hessians of boundary integrals are considered difficult to find analytically because they correspond to third-order derivatives of an, in a sense equivalent, domain integral. We complement previous results by considering cost functions depending explicitly on the surface normal. The correctness and practicability of our calculations are verified in the context of a Newton-type shape reconstruction method. © 2010 Birkhäuser / Springer Basel AG.

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

The influence of boundary conditions on resonant frequencies of cavities in 3-D FDTD algorithm using non-orthogonal co-ordinates

Energy Technology Data Exchange (ETDEWEB)

Zhao, L.; Tong, L.S. [Southeast Univ., Nanjing (China). Research Inst. of Electronics; Carter, R.G. [Lancaster Univ. (United Kingdom). Engineering Dept.

1994-09-01

The 3-dimensional finite-difference time-domain method in non-orthogonal co-ordinates (non-standard FDTD) is used to calculate the frequencies of resonators. The numerical boundary conditions of the method are presented. The Influences of boundary conditions and discrete meshes on the numerical accuracy are investigated. The authors present the nonstandard FDTD method using the boundary-orthogonal mesh and equivalent dielectric constant so that the error is reduced from 8.66% to 3.0% for the cylindrical cavity loaded by a dielectric button.

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.

Attractor of Beam Equation with Structural Damping under Nonlinear Boundary Conditions

Directory of Open Access Journals (Sweden)

Danxia Wang

2015-01-01

Full Text Available Simultaneously, considering the viscous effect of material, damping of medium, and rotational inertia, we study a kind of more general Kirchhoff-type extensible beam equation utt-uxxtt+uxxxx-σ(∫0l‍(ux2dxuxx-ϕ(∫0l‍(ux2dxuxxt=q(x, in  [0,L]×R+ with the structural damping and the rotational inertia term. Little attention is paid to the longtime behavior of the beam equation under nonlinear boundary conditions. In this paper, under nonlinear boundary conditions, we prove not only the existence and uniqueness of global solutions by prior estimates combined with some inequality skills, but also the existence of a global attractor by the existence of an absorbing set and asymptotic compactness of corresponding solution semigroup. In addition, the same results also can be proved under the other nonlinear boundary conditions.

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.

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.

Scale-invariant solutions to partial differential equations of fractional order with a moving boundary condition

International Nuclear Information System (INIS)

Li Xicheng; Xu Mingyu; Wang Shaowei

2008-01-01

In this paper, we give similarity solutions of partial differential equations of fractional order with a moving boundary condition. The solutions are given in terms of a generalized Wright function. The time-fractional Caputo derivative and two types of space-fractional derivatives are considered. The scale-invariant variable and the form of the solution of the moving boundary are obtained by the Lie group analysis. A comparison between the solutions corresponding to two types of fractional derivative is also given

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

A new technique for the determination of coronal magnetic fields: A fixed mesh solution to Laplace's equation using line-of-sight boundary conditions

International Nuclear Information System (INIS)

Adams, J.; Pneuman, G.W.

1976-01-01

A new method for computing potential magnetic field configurations in the solar atmosphere is described. A discrete approximation to Laplace's equation is solved in the domain R(Sun) 1 , 0 1 being an arbitrary radial distance from the solar center). The method utilizes the measured line-of-sight magnetic fields directly as the boundary condition at the solar surface and constrains the field to become radial at the outer boundary, R 1 . First the differential equation and boundary conditions are reduced to a set of two-dimensional equations in r, theta by Fourier transforming out the periodic phi dependence. Next each transformed boundary condition is converted to a Dirichlet surface condition. Then each two-dimensional equation with standard Dirichlet-Dirichlet boundary conditions is solved for the Fourier coefficient it determines. Finally, the solution of the original three dimensional equation is obtained through inverse Fourier transformation. The primary numerical tools in this technique are the use of a finite fast Fourier transform technique and also a generalized cyclic reduction algorithm developed at NCAR. Any extraneous monopole component present in the data can be removed if so desired. (Auth.)

Non-equilibrium concentration fluctuations in binary liquids with realistic boundary conditions.

Science.gov (United States)

Ortiz de Zárate, J M; Kirkpatrick, T R; Sengers, J V

2015-09-01

Because of the spatially long-ranged nature of spontaneous fluctuations in thermal non-equilibrium systems, they are affected by boundary conditions for the fluctuating hydrodynamic variables. In this paper we consider a liquid mixture between two rigid and impervious plates with a stationary concentration gradient resulting from a temperature gradient through the Soret effect. For liquid mixtures with large Lewis and Schmidt numbers, we are able to obtain explicit analytical expressions for the intensity of the non-equilibrium concentration fluctuations as a function of the frequency ω and the wave number q of the fluctuations. In addition we elucidate the spatial dependence of the intensity of the non-equilibrium fluctuations responsible for a non-equilibrium Casimir effect.

Charged particle motion in a time-dependent flux-driven ring: an exactly solvable model

International Nuclear Information System (INIS)

Luan, P-G; Tang, C-S

2007-01-01

We consider a charged particle driven by a time-dependent flux threading a quantum ring. The dynamics of the charged particle is investigated using a classical treatment, a Fourier expansion technique, a time-evolution method, and the Lewis-Riesenfeld approach. We have shown that, by properly managing the boundary conditions, a time-dependent wavefunction can be obtained using a general non-Hermitian time-dependent invariant, which is a specific linear combination of initial angular-momentum and azimuthal-angle operators. It is shown that the linear invariant eigenfunction can be realized as a Gaussian-type wavepacket with a peak moving along the classical angular trajectory, while the distribution of the wavepacket is determined by the ratio of the coefficient of the initial angle to that of the initial canonical angular momentum. From the topologically nontrivial nature as well as the classical trajectory and angular momentum, one can determine the dynamical motion of the wavepacket. It should be noted that the peak position is no longer an expectation value of the angle operator, and hence the Ehrenfest theorem is not directly applicable in such a topologically nontrivial system

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

Synchronization of Reaction-Diffusion Neural Networks With Dirichlet Boundary Conditions and Infinite Delays.

Science.gov (United States)

Sheng, Yin; Zhang, Hao; Zeng, Zhigang

2017-10-01

This paper is concerned with synchronization for a class of reaction-diffusion neural networks with Dirichlet boundary conditions and infinite discrete time-varying delays. By utilizing theories of partial differential equations, Green's formula, inequality techniques, and the concept of comparison, algebraic criteria are presented to guarantee master-slave synchronization of the underlying reaction-diffusion neural networks via a designed controller. Additionally, sufficient conditions on exponential synchronization of reaction-diffusion neural networks with finite time-varying delays are established. The proposed criteria herein enhance and generalize some published ones. Three numerical examples are presented to substantiate the validity and merits of the obtained theoretical results.

Transport sensitivity studies for SITE-94: Time-dependent site-scale modelling of future glacial impact

International Nuclear Information System (INIS)

King-Clayton, L.; Smith, Paul

1996-10-01

The report details the methodology and preliminary results from the modelling of radionuclide transport from a hypothetical repository based at the Aespoe site in Sweden. The work complements and utilizes the results from regional-scale, variable density flow modelling by Provost, in which the groundwater flow field is time dependent, reflecting the impact of climate evolution over the next 130,000 years. The climate evolution include development of permafrost conditions and ice sheet advance and retreat. The results indicate that temporal changes in flow conditions owing to future climate changes have a significant effect on the transport of radionuclides from a repository. In all cases modelled with time-dependent boundary conditions, the greatest radionuclide fluxes occur towards the end of the main glacial periods, and correspond to periods of high groundwater discharge at the margin of the modelled ice sheets. Fluxes to the biosphere may, for limited periods (2000 years or less), be three times higher than those from the near field. The study provides a quantitative way of illustrating the possible effects of future glaciations on radionuclide transport from the repository. Such effects are likely to be significant in any potential siting area predicted to be affected by future periods of ice cover. 8 refs, 22 tabs, 119 figs

Nonlinear fractional differential equations and inclusions of arbitrary order and multi-strip boundary conditions

Directory of Open Access Journals (Sweden)

Bashir Ahmad

2012-06-01

Full Text Available We study boundary value problems of nonlinear fractional differential equations and inclusions of order $q in (m-1, m]$, $m ge 2$ with multi-strip boundary conditions. Multi-strip boundary conditions may be regarded as the generalization of multi-point boundary conditions. Our problem is new in the sense that we consider a nonlocal strip condition of the form: $$ x(1=sum_{i=1}^{n-2}alpha_i int^{eta_i}_{zeta_i} x(sds, $$ which can be viewed as an extension of a multi-point nonlocal boundary condition: $$ x(1=sum_{i=1}^{n-2}alpha_i x(eta_i. $$ In fact, the strip condition corresponds to a continuous distribution of the values of the unknown function on arbitrary finite segments $(zeta_i,eta_i$ of the interval $[0,1]$ and the effect of these strips is accumulated at $x=1$. Such problems occur in the applied fields such as wave propagation and geophysics. Some new existence and uniqueness results are obtained by using a variety of fixed point theorems. Some illustrative examples are also discussed.

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.

Seismic response of three-dimensional topographies using a time-domain boundary element method

Science.gov (United States)

Janod, François; Coutant, Olivier

2000-08-01

We present a time-domain implementation for a boundary element method (BEM) to compute the diffraction of seismic waves by 3-D topographies overlying a homogeneous half-space. This implementation is chosen to overcome the memory limitations arising when solving the boundary conditions with a frequency-domain approach. This formulation is flexible because it allows one to make an adaptive use of the Green's function time translation properties: the boundary conditions solving scheme can be chosen as a trade-off between memory and cpu requirements. We explore here an explicit method of solution that requires little memory but a high cpu cost in order to run on a workstation computer. We obtain good results with four points per minimum wavelength discretization for various topographies and plane wave excitations. This implementation can be used for two different aims: the time-domain approach allows an easier implementation of the BEM in hybrid methods (e.g. coupling with finite differences), and it also allows one to run simple BEM models with reasonable computer requirements. In order to keep reasonable computation times, we do not introduce any interface and we only consider homogeneous models. Results are shown for different configurations: an explosion near a flat free surface, a plane wave vertically incident on a Gaussian hill and on a hemispherical cavity, and an explosion point below the surface of a Gaussian hill. Comparison is made with other numerical methods, such as finite difference methods (FDMs) and spectral elements.

A hybrid method combining the FDTD and a time domain boundary-integral equation marching-on-in-time algorithm

Directory of Open Access Journals (Sweden)

A. Becker

2003-01-01

Full Text Available In this paper a hybrid method combining the FDTD/FIT with a Time Domain Boundary-Integral Marching-on-in-Time Algorithm (TD-BIM is presented. Inhomogeneous regions are modelled with the FIT-method, an alternative formulation of the FDTD. Homogeneous regions (which is in the presented numerical example the open space are modelled using a TD-BIM with equivalent electric and magnetic currents flowing on the boundary between the inhomogeneous and the homogeneous regions. The regions are coupled by the tangential magnetic fields just outside the inhomogeneous regions. These fields are calculated by making use of a Mixed Potential Integral Formulation for the magnetic field. The latter consists of equivalent electric and magnetic currents on the boundary plane between the homogeneous and the inhomogeneous region. The magnetic currents result directly from the electric fields of the Yee lattice. Electric currents in the same plane are calculated by making use of the TD-BIM and using the electric field of the Yee lattice as boundary condition. The presented hybrid method only needs the interpolations inherent in FIT and no additional interpolation. A numerical result is compared to a calculation that models both regions with FDTD.

IMPROVEMENT OF ACCURACY OF RADIATIVE HEAT TRANSFER DIFFERENTIAL APPROXIMATION METHOD FOR MULTI DIMENSIONAL SYSTEMS BY MEANS OF AUTO-ADAPTABLE BOUNDARY CONDITIONS

Directory of Open Access Journals (Sweden)

K. V. Dobrego

2015-01-01

Full Text Available Differential approximation is derived from radiation transfer equation by averaging over the solid angle. It is one of the more effective methods for engineering calculations of radia- tive heat transfer in complex three-dimensional thermal power systems with selective and scattering media. The new method for improvement of accuracy of the differential approximation based on using of auto-adaptable boundary conditions is introduced in the paper. The  efficiency  of  the  named  method  is  proved  for  the  test  2D-systems.  Self-consistent auto-adaptable boundary conditions taking into consideration the nonorthogonal component of the incident to the boundary radiation flux are formulated. It is demonstrated that taking in- to consideration of the non- orthogonal incident flux in multi-dimensional systems, such as furnaces, boilers, combustion chambers improves the accuracy of the radiant flux simulations and to more extend in the zones adjacent to the edges of the chamber.Test simulations utilizing the differential approximation method with traditional boundary conditions, new self-consistent boundary conditions and “precise” discrete ordinates method were performed. The mean square errors of the resulting radiative fluxes calculated along the boundary of rectangular and triangular test areas were decreased 1.5–2 times by using auto- adaptable boundary conditions. Radiation flux gaps in the corner points of non-symmetric sys- tems are revealed by using auto-adaptable boundary conditions which can not be obtained by using the conventional boundary conditions.

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.

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

Large Scale Skill in Regional Climate Modeling and the Lateral Boundary Condition Scheme

Science.gov (United States)

Veljović, K.; Rajković, B.; Mesinger, F.

2009-04-01

Several points are made concerning the somewhat controversial issue of regional climate modeling: should a regional climate model (RCM) be expected to maintain the large scale skill of the driver global model that is supplying its lateral boundary condition (LBC)? Given that this is normally desired, is it able to do so without help via the fairly popular large scale nudging? Specifically, without such nudging, will the RCM kinetic energy necessarily decrease with time compared to that of the driver model or analysis data as suggested by a study using the Regional Atmospheric Modeling System (RAMS)? Finally, can the lateral boundary condition scheme make a difference: is the almost universally used but somewhat costly relaxation scheme necessary for a desirable RCM performance? Experiments are made to explore these questions running the Eta model in two versions differing in the lateral boundary scheme used. One of these schemes is the traditional relaxation scheme, and the other the Eta model scheme in which information is used at the outermost boundary only, and not all variables are prescribed at the outflow boundary. Forecast lateral boundary conditions are used, and results are verified against the analyses. Thus, skill of the two RCM forecasts can be and is compared not only against each other but also against that of the driver global forecast. A novel verification method is used in the manner of customary precipitation verification in that forecast spatial wind speed distribution is verified against analyses by calculating bias adjusted equitable threat scores and bias scores for wind speeds greater than chosen wind speed thresholds. In this way, focusing on a high wind speed value in the upper troposphere, verification of large scale features we suggest can be done in a manner that may be more physically meaningful than verifications via spectral decomposition that are a standard RCM verification method. The results we have at this point are somewhat

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)

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.

Boundary Correct Real-Time Soft Shadows

DEFF Research Database (Denmark)

Jacobsen, Bjarke; Christensen, Niels Jørgen; Larsen, Bent Dalgaard

2004-01-01

This paper describes a method to determine correct shadow boundaries from an area light source using umbra and penumbra volumes. The light source is approximated by a circular disk as this gives a fast way to extrude the volumes. The method also gives a crude estimate of the visibility of the are...... for implementation on most programmable hardware. Though some crude approximations are used in the visibility function, the method can be used to produce soft shadows with correct boundaries in real time....

Thermal properties of nuclear matter under the periodic boundary condition

International Nuclear Information System (INIS)

Otuka, Naohiko; Ohnishi, Akira

1999-01-01

We present the thermal properties of nuclear matter under the periodic boundary condition by the use of our hadronic nucleus-nucleus cascade model (HANDEL) which is developed to treat relativistic heavy-ion collisions from BNL-AGS to CERN-SPS. We first show some results of p-p scattering calculation in our new version which is improved in order to treat isospin ratio and multiplicity more accurately. We then display the results of calculation of nuclear matter with baryon density ρ b = 0.77 fm 3 at some energy densities. Time evolution of particle abundance and temperature are shown. (author)

Effect of boundary conditions on measured water retention behavior within soils

Science.gov (United States)

Galindo-torres, S.; Scheuermann, A.; Pedroso, D.; Li, L.

2013-12-01

conducted to examine how the fluids are distributed inside the porous medium. This distribution is quantified by the measurement of the interfacial area which behaves also differently between the two configurations. Hassanizadeh proposed an unique relation among saturation, suction and interfacial area, which has been validated experimentally [4]. However we found that such relation is not 'unique' and instead depends on the flow and boundary conditions. While future experimental tests on these results need to be carried out, the simulated SWCC behaviors raise serious questions about the current experimental set-up for measuring the soil water retention characteristics. References. 1. Serrano, S.E., Modeling infiltration with approximate solutions to Richard's equation. Journal of Hydrologic Engineering, 2004. 9(5): p. 421-432. 2. Galindo-Torres, S.A., et al., A Lattice Boltzmann model for studying transient effects during imbibition-drainage cycles in unsaturated soils. Computer Physics Communications, 2013. 184(4): p. 1086-1093. 3. Drake, S.S., D.M. O'Carroll, and J.I. Gerhard, Wettability contrasts between fresh and weathered diesel fuels. Journal of contaminant hydrology, 2012. 4. Culligan, K.A., et al., Interfacial area measurements for unsaturated flow through a porous medium. Water Resources Research, 2004. 40(12).

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

Modeling long period swell in Southern California: Practical boundary conditions from buoy observations and global wave model predictions

Science.gov (United States)

Crosby, S. C.; O'Reilly, W. C.; Guza, R. T.

2016-02-01

Accurate, unbiased, high-resolution (in space and time) nearshore wave predictions are needed to drive models of beach erosion, coastal flooding, and alongshore transport of sediment, biota and pollutants. On highly sheltered shorelines, wave predictions are sensitive to the directions of onshore propagating waves, and nearshore model prediction error is often dominated by uncertainty in offshore boundary conditions. Offshore islands and shoals, and coastline curvature, create complex sheltering patterns over the 250km span of southern California (SC) shoreline. Here, regional wave model skill in SC was compared for different offshore boundary conditions created using offshore buoy observations and global wave model hindcasts (National Oceanographic and Atmospheric Administration Wave Watch 3, WW3). Spectral ray-tracing methods were used to transform incident offshore swell (0.04-0.09Hz) energy at high directional resolution (1-deg). Model skill is assessed for predictions (wave height, direction, and alongshore radiation stress) at 16 nearshore buoy sites between 2000 and 2009. Model skill using buoy-derived boundary conditions is higher than with WW3-derived boundary conditions. Buoy-driven nearshore model results are similar with various assumptions about the true offshore directional distribution (maximum entropy, Bayesian direct, and 2nd derivative smoothness). Two methods combining offshore buoy observations with WW3 predictions in the offshore boundary condition did not improve nearshore skill above buoy-only methods. A case example at Oceanside harbor shows strong sensitivity of alongshore sediment transport predictions to different offshore boundary conditions. Despite this uncertainty in alongshore transport magnitude, alongshore gradients in transport (e.g. the location of model accretion and erosion zones) are determined by the local bathymetry, and are similar for all predictions.

Short-time Asymptotics of the Heat Kernel on Bounded Domain with Piecewise Smooth Boundary Conditions and Its Applications to an Ideal Gas

Institute of Scientific and Technical Information of China (English)

E.M.E. ZAYED

2004-01-01

The asymptotic expansion of the heat kernel Θ(t)(∞∑=(i=0))exp (-λi) where({λi}∞i=1) Are the eigen-values of negative Laplacian( -△n=-n∑k=1(θ/θxk)2)in Rn(n=2 or 3) is studied for short-time t for a general bounded domainθΩwith a smooth boundary θΩ.In this paper, we consider the case of a finite number of the Dirichlet conditions φ=0 on Γi (i = J +1,….,J)and the Neumann conditions and (θφ/θ vi) = 0 on Γi (i = J+1,…,k) and the Robin condition (θφ/θ vi+γi) θ=(I=k+1,… m) where γi are piecewise smooth positive impedancem(θφ=mUi=1Γi. )We construct the required asymptotics in the form of a power series over t. The senior coe.cients inthis series are speci.ed as functionals of the geometric shape of the domain Ω.This result is applied to calculatethe one-particle partition function of a "special ideal gas", i.e., the set of non-interacting particles set up in abox with Dirichlet, Neumann and Robin boundary conditions for the appropriate wave function. Calculationof the thermodynamic quantities for the ideal gas such as the internal energy, pressure and speci.c heat revealsthat these quantities alone are incapable of distinguishing between two di.erent shapes of the domain. Thisconclusion seems to be intuitively clear because it is based on a limited information given by a one-particlepartition function; nevertheless, its formal theoretical motivation is of some interest.

Time-dependent patterns in quasivertical cylindrical binary convection

Science.gov (United States)

Alonso, Arantxa; Mercader, Isabel; Batiste, Oriol

2018-02-01

This paper reports on numerical investigations of the effect of a slight inclination α on pattern formation in a shallow vertical cylindrical cell heated from below for binary mixtures with a positive value of the Soret coefficient. By using direct numerical simulation of the three-dimensional Boussinesq equations with Soret effect in cylindrical geometry, we show that a slight inclination of the cell in the range α ≈0.036 rad =2∘ strongly influences pattern selection. The large-scale shear flow (LSSF) induced by the small tilt of gravity overcomes the squarelike arrangements observed in noninclined cylinders in the Soret regime, stratifies the fluid along the direction of inclination, and produces an enhanced separation of the two components of the mixture. The competition between shear effects and horizontal and vertical buoyancy alters significantly the dynamics observed in noninclined convection. Additional unexpected time-dependent patterns coexist with the basic LSSF. We focus on an unsual periodic state recently discovered in an experiment, the so-called superhighway convection state (SHC), in which ascending and descending regions of fluid move in opposite directions. We provide numerical confirmation that Boussinesq Navier-Stokes equations with standard boundary conditions contain the essential ingredients that allow for the existence of such a state. Also, we obtain a persistent heteroclinic structure where regular oscillations between a SHC pattern and a state of nearly stationary longitudinal rolls take place. We characterize numerically these time-dependent patterns and investigate the dynamics around the threshold of convection.

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

Time-Dependent Heat Conduction Problems Solved by an Integral-Equation Approach

International Nuclear Information System (INIS)

Oberaigner, E.R.; Leindl, M.; Antretter, T.

2010-01-01

Full text: A classical task of mathematical physics is the formulation and solution of a time dependent thermoelastic problem. In this work we develop an algorithm for solving the time-dependent heat conduction equation c p ρ∂ t T-kT, ii =0 in an analytical, exact fashion for a two-component domain. By the Green's function approach the formal solution of the problem is obtained. As an intermediate result an integral-equation for the temperature history at the domain interface is formulated which can be solved analytically. This method is applied to a classical engineering problem, i.e. to a special case of a Stefan-Problem. The Green's function approach in conjunction with the integral-equation method is very useful in cases were strong discontinuities or jumps occur. The initial conditions and the system parameters of the investigated problem give rise to two jumps in the temperature field. Purely numerical solutions are obtained by using the FEM (finite element method) and the FDM (finite difference method) and compared with the analytical approach. At the domain boundary the analytical solution and the FEM-solution are in good agreement, but the FDM results show a signicant smearing effect. (author)

Trace expansions for mixed boundary problems

Energy Technology Data Exchange (ETDEWEB)

Seeley, Robert T

2002-01-01

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

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

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.

The D(D3)-anyon chain: integrable boundary conditions and excitation spectra

International Nuclear Information System (INIS)

Finch, Peter E; Frahm, Holger

2013-01-01

Chains of interacting non-Abelian anyons with local interactions invariant under the action of the Drinfeld double of the dihedral group D 3 are constructed. Formulated as a spin chain the Hamiltonians are generated from commuting transfer matrices of an integrable vertex model for periodic and braided as well as open boundaries. A different anyonic model with the same local Hamiltonian is obtained within the fusion path formulation. This model is shown to be related to an integrable fusion interaction round the face model. Bulk and surface properties of the anyon chain are computed from the Bethe equations for the spin chain. The low-energy effective theories and operator content of the models (in both the spin chain and fusion path formulation) are identified from analytical and numerical studies of the finite-size spectra. For all boundary conditions considered the continuum theory is found to be a product of two conformal field theories. Depending on the coupling constants the factors can be a Z 4 parafermion or a M (5,6) minimal model. (paper)

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

Joining the Great Plains in Space, Place, and Time: Questioning a Time Zone Boundary

Science.gov (United States)

Kuper, Rob

2011-01-01

Standard time zone boundaries are invisible in the landscape, yet they abruptly delineate a temporal difference of one hour between two large areas located relative to one another on Earth. In most cases, standard time zone boundaries follow political ones and define areas within which daylight saving time (DST)--the seasonal advancement of…

Evaluation of Time-Dependent Behavior of Soils

DEFF Research Database (Denmark)

Augustesen, Anders; Liingaard, Morten; Lade, Poul V.

2004-01-01

The time-dependent behavior of soils has been investigated extensively through one-dimensional and triaxial test conditions. Most of the observations in literature have focused on the determination of the time-dependent behavior of clayey soils, whereas the reported experimental studies of granular...... situation for soils. That is whether the time-dependent behavior can be characterized as isotach or nonisotach. It seems that the isotach behavior is adequate for describing the time effects in clays in most situations. But for sand, the isotach description is inadequate. Further, the phenomenon...

Kinetics of drug release from ointments: Role of transient-boundary layer.

Science.gov (United States)

Xu, Xiaoming; Al-Ghabeish, Manar; Krishnaiah, Yellela S R; Rahman, Ziyaur; Khan, Mansoor A

2015-10-15

In the current work, an in vitro release testing method suitable for ointment formulations was developed using acyclovir as a model drug. Release studies were carried out using enhancer cells on acyclovir ointments prepared with oleaginous, absorption, and water-soluble bases. Kinetics and mechanism of drug release was found to be highly dependent on the type of ointment bases. In oleaginous bases, drug release followed a unique logarithmic-time dependent profile; in both absorption and water-soluble bases, drug release exhibited linearity with respect to square root of time (Higuchi model) albeit differences in the overall release profile. To help understand the underlying cause of logarithmic-time dependency of drug release, a novel transient-boundary hypothesis was proposed, verified, and compared to Higuchi theory. Furthermore, impact of drug solubility (under various pH conditions) and temperature on drug release were assessed. Additionally, conditions under which deviations from logarithmic-time drug release kinetics occur were determined using in situ UV fiber-optics. Overall, the results suggest that for oleaginous ointments containing dispersed drug particles, kinetics and mechanism of drug release is controlled by expansion of transient boundary layer, and drug release increases linearly with respect to logarithmic time. Published by Elsevier B.V.