A Rotational Pressure-Correction Scheme for Incompressible Two-Phase Flows with Open Boundaries

Science.gov (United States)

Dong, S.; Wang, X.

2016-01-01

Two-phase outflows refer to situations where the interface formed between two immiscible incompressible fluids passes through open portions of the domain boundary. We present several new forms of open boundary conditions for two-phase outflow simulations within the phase field framework, as well as a rotational pressure correction based algorithm for numerically treating these open boundary conditions. Our algorithm gives rise to linear algebraic systems for the velocity and the pressure that involve only constant and time-independent coefficient matrices after discretization, despite the variable density and variable viscosity of the two-phase mixture. By comparing simulation results with theory and the experimental data, we show that the method produces physically accurate results. We also present numerical experiments to demonstrate the long-term stability of the method in situations where large density contrast, large viscosity contrast, and backflows occur at the two-phase open boundaries. PMID:27163909

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

Acoustic Emission for on-line reactor pressure boundary monitoring

International Nuclear Information System (INIS)

Hutton, P.H.; Kurtz, R.J.; Pappas, R.A.

1985-01-01

The program objective is to develop AE for continuous surveillance to assess flaw growth in reactor pressure boundaries. Technology in the laboratory is being evaluated on structures. Results have demonstrated basic feasibility of the program objective. AE monitoring a long term fatigue test of a pressure vessel demonstrated an instrument system, and the ability to detect unexpected as well as well as known fatigue cracks. Monitoring a nuclear reactor system shows that the coolant flow noise problem is manageable and AE can be detected under simulated operating conditions

Analysis of specific factors causing RCS pressure boundary cracking

International Nuclear Information System (INIS)

Song, Taek-Ho; Jeong, Il-Seok

2007-01-01

As nuclear power plants become aged, pressure boundary integrity has become so important issue in domestic and foreign nuclear industry that many related research projects are on-going. KEPRI is going to embark a new research project for managing and preventing these kinds of cracks in nuclear power plants (NPPs). Many nuclear power plants experienced pressure boundary stress corrosion cracking (SCC) and shut downed because of it. In USA, V.C. Summer plant experienced reactor coolant pipe SCC near reactor outlet nozzle and Davis Vesse plant experienced reactor head crack around penetration pipe which is used to control rod drive mechanism. In this paper, RCS pressure boundary cracking cases and corrosion potential have been studied to find out what are the specific factors that have affected crack initiations in the reactor coolant pressure boundaries

Molecular dynamics simulations of a fully hydrated dipalmitoyl phosphatidylcholine bilayer with different macroscopic boundary conditions and parameters

NARCIS (Netherlands)

Tieleman, D.P; Berendsen, H.J.C.

1996-01-01

We compared molecular dynamics simulations of a bilayer of 128 fully hydrated phospholipid (DPPC) molecules, using different parameters and macroscopic boundary conditions. The same system was studied under constant pressure, constant volume, and constant surface tension boundary conditions, with

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

Transient pressure and productivity analysis in carbonate geothermal reservoirs with changing external boundary flux

Directory of Open Access Journals (Sweden)

Wang Dongying

2017-01-01

Full Text Available In this paper, a triple-medium flow model for carbonate geothermal reservoirs with an exponential external boundary flux is established. The pressure solution under constant production conditions in Laplace space is solved. The geothermal wellbore pressure change considering wellbore storage and skin factor is obtained by Stehfest numerical inversion. The well test interpretation charts and Fetkovich production decline chart for carbonate geothermal reservoirs are proposed for the first time. The proposed Fetkovich production decline curves are applied to analyze the production decline behavior. The results indicate that in carbonate geothermal reservoirs with exponential external boundary flux, the pressure derivative curve contains a triple dip, which represents the interporosity flow between the vugs or matrix and fracture system and the invading flow of the external boundary flux. The interporosity flow of carbonate geothermal reservoirs and changing external boundary flux can both slow down the extent of production decline and the same variation tendency is observed in the Fetkovich production decline curve.

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.

DNS of transcritical turbulent boundary layers at supercritical pressures under abrupt variations in thermodynamic properties

Science.gov (United States)

Kawai, Soshi

2014-11-01

In this talk, we first propose a numerical strategy that is robust and high-order accurate for enabling to simulate transcritical flows at supercritical pressures under abrupt variations in thermodynamic properties due to the real fluid effects. The method is based on introducing artificial density diffusion in a physically-consistent manner in order to capture the steep variation of thermodynamic properties in transcritical conditions robustly, while solving a pressure evolution equation to achieve pressure equilibrium at the transcritical interfaces. We then discuss the direct numerical simulation (DNS) of transcritical heated turbulent boundary layers on a zero-pressure-gradient flat plate at supercritical pressures. To the best of my knowledge, the present DNS is the first DNS of zero-pressure-gradient flat-plate transcritical turbulent boundary layer. The turbulent kinetic budget indicates that the compressibility effects (especially, pressure-dilatation correlation) are not negligible at the transcritical conditions even if the flow is subsonic. The unique and interesting interactions between the real fluid effects and wall turbulence, and their turbulence statistics, which have never been seen in the ideal-fluid turbulent boundary layers, are also discussed. This work was supported in part by Japan Society for the Promotion of Science (JSPS) Grant-in-Aid for Young Scientists (A) KAKENHI 26709066 and the JAXA International Top Young Fellowship Program.

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.

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

Science.gov (United States)

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

2012-11-01

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

Pressure Fluctuations Induced by a Hypersonic Turbulent Boundary Layer

Science.gov (United States)

Duan, Lian; Choudhari, Meelan M.; Zhang, Chao

2016-01-01

Direct numerical simulations (DNS) are used to examine the pressure fluctuations generated by a spatially-developed Mach 5.86 turbulent boundary layer. The unsteady pressure field is analyzed at multiple wall-normal locations, including those at the wall, within the boundary layer (including inner layer, the log layer, and the outer layer), and in the free stream. The statistical and structural variations of pressure fluctuations as a function of wall-normal distance are highlighted. Computational predictions for mean velocity pro les and surface pressure spectrum are in good agreement with experimental measurements, providing a first ever comparison of this type at hypersonic Mach numbers. The simulation shows that the dominant frequency of boundary-layer-induced pressure fluctuations shifts to lower frequencies as the location of interest moves away from the wall. The pressure wave propagates with a speed nearly equal to the local mean velocity within the boundary layer (except in the immediate vicinity of the wall) while the propagation speed deviates from the Taylor's hypothesis in the free stream. Compared with the surface pressure fluctuations, which are primarily vortical, the acoustic pressure fluctuations in the free stream exhibit a significantly lower dominant frequency, a greater spatial extent, and a smaller bulk propagation speed. The freestream pressure structures are found to have similar Lagrangian time and spatial scales as the acoustic sources near the wall. As the Mach number increases, the freestream acoustic fluctuations exhibit increased radiation intensity, enhanced energy content at high frequencies, shallower orientation of wave fronts with respect to the flow direction, and larger propagation velocity.

A parametric study of adverse pressure gradient turbulent boundary layers

International Nuclear Information System (INIS)

Monty, J.P.; Harun, Z.; Marusic, I.

2011-01-01

There are many open questions regarding the behaviour of turbulent boundary layers subjected to pressure gradients and this is confounded by the large parameter space that may affect these flows. While there have been many valuable investigations conducted within this parameter space, there are still insufficient data to attempt to reduce this parameter space. Here, we consider a parametric study of adverse pressure gradient turbulent boundary layers where we restrict our attention to the pressure gradient parameter, β, the Reynolds number and the acceleration parameter, K. The statistics analyzed are limited to the streamwise fluctuating velocity. The data show that the mean velocity profile in strong pressure gradient boundary layers does not conform to the classical logarithmic law. Moreover, there appears to be no measurable logarithmic region in these cases. It is also found that the large-scale motions scaling with outer variables are energised by the pressure gradient. These increasingly strong large-scale motions are found to be the dominant contributor to the increase in turbulence intensity (scaled with friction velocity) with increasing pressure gradient across the boundary layer.

A survey of repair practices for nuclear power plant containment metallic pressure boundaries

Energy Technology Data Exchange (ETDEWEB)

Oland, C.B.; Naus, D.J. [Oak Ridge National Lab., TN (United States)

1998-05-01

The Nuclear Regulatory Commission has initiated a program at the Oak Ridge National Laboratory to provide assistance in their assessment of the effects of potential degradation on the structural integrity and leaktightness of metal containment vessels and steel liners of concrete containments in nuclear power plants. One of the program objectives is to identify repair practices for restoring metallic containment pressure boundary components that have been damaged or degraded in service. This report presents issues associated with inservice condition assessments and continued service evaluations and identifies the rules and requirements for the repair and replacement of nonconforming containment pressure boundary components by welding or metal removal. Discussion topics include base and welding materials, welding procedure and performance qualifications, inspection techniques, testing methods, acceptance criteria, and documentation requirements necessary for making acceptable repairs and replacements so that the plant can be returned to a safe operating condition.

A survey of repair practices for nuclear power plant containment metallic pressure boundaries

International Nuclear Information System (INIS)

Oland, C.B.; Naus, D.J.

1998-05-01

The Nuclear Regulatory Commission has initiated a program at the Oak Ridge National Laboratory to provide assistance in their assessment of the effects of potential degradation on the structural integrity and leaktightness of metal containment vessels and steel liners of concrete containments in nuclear power plants. One of the program objectives is to identify repair practices for restoring metallic containment pressure boundary components that have been damaged or degraded in service. This report presents issues associated with inservice condition assessments and continued service evaluations and identifies the rules and requirements for the repair and replacement of nonconforming containment pressure boundary components by welding or metal removal. Discussion topics include base and welding materials, welding procedure and performance qualifications, inspection techniques, testing methods, acceptance criteria, and documentation requirements necessary for making acceptable repairs and replacements so that the plant can be returned to a safe operating condition

Crack growth behaviour of low-alloy steels for pressure boundary components under transient light water reactor operating conditions - CASTOC, Part II: WWER conditions

Energy Technology Data Exchange (ETDEWEB)

Ernestova, M.; Zamboch, M. [Nuclear Research Institute, NRI, Rez (Czech Republic); Devrient, B.; Roth, A. [Framatome ANP GmbH, Erlangen (Germany); Ehrnsten, U. [VTT Industrial Systems, Espoo (Finland); Foehl, J.; Weissenberg, T. [Staatliche Materialpruefungsanstalt, MPA, Stuttgart (Germany); Gomez-Briceno, D.; Lapena, J. [Centro de Investigaciones Energeticas Medioambientales y Tecnologicas, CIEMAT, Madrid (Spain); Ritter, S.; Seifert, H.P. [Paul Scherrer Institute, PSI, Villigen (Switzerland)

2004-07-01

One of the ageing phenomena of pressure boundary components of light water reactors (LWRs) is environmentally-assisted cracking (EAC). The project CASTOC (5. Framework Programme of the EU) was launched September 2000 with six European partners and terminated August 2003. It focused in particular on the EAC behaviour of low-alloy steels (LAS) and to some extent to weld metal, heat affected zone and the influence of an austenitic cladding. The main objective was directed to the clarification of crack growth behavior of LAS in high-temperature water due to EAC under constant load (steady-state power operation), to study the effect of transient conditions (during operation or start-up/shut-down of a plant) using their impact on time-based and cycle-based crack growth rates and to a more detailed understanding of the acting mechanisms. Autoclave tests were performed with Western and Russian type reactor pressure vessel steels under simulated boiling water reactor (BWR)/normal water chemistry (NWC) and pressurized water reactor (WWER) conditions. The investigations were performed with fracture mechanics specimens of different sizes and geometries. The applied loading comprised cyclic loads, static loads and load spectra where the static load was periodically interrupted by partial unloading. With regard to water chemistry, the oxygen content (WWER) and impurities of sulphate and chlorides (BWR) were varied beyond allowable limits for continuous operation. The current paper summarizes the most important crack growth results obtained under simulated WWER conditions. The influence of oxygen content and the effect of specimen size (C(T)25 versus C(T)50 specimens) on the crack growth rates are shown. The results are discussed in the context of the current crack growth rate curves in the corresponding nuclear codes. (authors)

Crack growth behaviour of low-alloy steels for pressure boundary components under transient light water reactor operating conditions - CASTOC, Part II: WWER conditions

International Nuclear Information System (INIS)

Ernestova, M.; Zamboch, M.; Devrient, B.; Roth, A.; Ehrnsten, U.; Foehl, J.; Weissenberg, T.; Gomez-Briceno, D.; Lapena, J.; Ritter, S.; Seifert, H.P.

2004-01-01

One of the ageing phenomena of pressure boundary components of light water reactors (LWRs) is environmentally-assisted cracking (EAC). The project CASTOC (5. Framework Programme of the EU) was launched September 2000 with six European partners and terminated August 2003. It focused in particular on the EAC behaviour of low-alloy steels (LAS) and to some extent to weld metal, heat affected zone and the influence of an austenitic cladding. The main objective was directed to the clarification of crack growth behavior of LAS in high-temperature water due to EAC under constant load (steady-state power operation), to study the effect of transient conditions (during operation or start-up/shut-down of a plant) using their impact on time-based and cycle-based crack growth rates and to a more detailed understanding of the acting mechanisms. Autoclave tests were performed with Western and Russian type reactor pressure vessel steels under simulated boiling water reactor (BWR)/normal water chemistry (NWC) and pressurized water reactor (WWER) conditions. The investigations were performed with fracture mechanics specimens of different sizes and geometries. The applied loading comprised cyclic loads, static loads and load spectra where the static load was periodically interrupted by partial unloading. With regard to water chemistry, the oxygen content (WWER) and impurities of sulphate and chlorides (BWR) were varied beyond allowable limits for continuous operation. The current paper summarizes the most important crack growth results obtained under simulated WWER conditions. The influence of oxygen content and the effect of specimen size (C(T)25 versus C(T)50 specimens) on the crack growth rates are shown. The results are discussed in the context of the current crack growth rate curves in the corresponding nuclear codes. (authors)

Crack growth behaviour of low-alloy steels for pressure boundary components under transient light water reactor operating conditions - CASTOC, Part I: BWR/NWC conditions

International Nuclear Information System (INIS)

Ritter, S.; Seifert, H.P.; Devrient, B.; Roth, A.; Ehrnsten, U.; Ernestova, M.; Zamboch, M.; Foehl, J.; Weissenberg, T.; Gomez-Briceno, D.; Lapena, J.

2004-01-01

One of the ageing phenomena of pressure boundary components of light water reactors (LWR) is environmentally-assisted cracking (EAC). The project CASTOC (5. Framework Programme of the EU) was launched September 2000 with six European partners and terminated August 2003. It was focused in particular on the EAC behaviour of low-alloy steels (LAS) and to some extent to weld metal, heat affected zone and the influence of an austenitic cladding. The main objective was directed to the clarification of EAC crack growth behaviour/mechanism of LAS in high-temperature water under steady-state power operation (constant load) and transient operating conditions (e.g., start-up/shut-down, transients in water chemistry and load). Autoclave tests were performed with Western and Russian type reactor pressure vessel steels under simulated boiling water reactor (BWR)/normal water chemistry (NWC) and pressurised water reactor (VVER) conditions. The investigations were performed with fracture mechanics specimens of different sizes and geometries. The applied loading comprised cyclic loads, static loads and load spectra where the static load was periodically interrupted by partial unloading. With regard to water chemistry, the oxygen content (VVER) and impurities of sulphate and chlorides (BWR) were varied beyond allowable limits for continuous operation. The current paper summarises the most important crack growth results obtained under simulated BWR/NWC conditions. The results are discussed in the context of the current crack growth rate curves in the corresponding nuclear codes. (authors)

Crack growth behaviour of low-alloy steels for pressure boundary components under transient light water reactor operating conditions - CASTOC, Part I: BWR/NWC conditions

Energy Technology Data Exchange (ETDEWEB)

Ritter, S.; Seifert, H.P. [Paul Scherrer Institute, PSI, Villigen (Switzerland); Devrient, B.; Roth, A. [Framatome ANP GmbH, Erlangen (Germany); Ehrnsten, U. [VTT Industrial Systems, Espoo (Finland); Ernestova, M.; Zamboch, M. [Nuclear Research Institute, NRI, Rez (Czech Republic); Foehl, J.; Weissenberg, T. [Staatliche Materialpruefungsanstalt, MPA, Stuttgart (Germany); Gomez-Briceno, D.; Lapena, J. [Centro de Investigaciones Energeticas Medioambientales y Tecnologicas, CIEMAT, Madrid (Spain)

2004-07-01

One of the ageing phenomena of pressure boundary components of light water reactors (LWR) is environmentally-assisted cracking (EAC). The project CASTOC (5. Framework Programme of the EU) was launched September 2000 with six European partners and terminated August 2003. It was focused in particular on the EAC behaviour of low-alloy steels (LAS) and to some extent to weld metal, heat affected zone and the influence of an austenitic cladding. The main objective was directed to the clarification of EAC crack growth behaviour/mechanism of LAS in high-temperature water under steady-state power operation (constant load) and transient operating conditions (e.g., start-up/shut-down, transients in water chemistry and load). Autoclave tests were performed with Western and Russian type reactor pressure vessel steels under simulated boiling water reactor (BWR)/normal water chemistry (NWC) and pressurised water reactor (VVER) conditions. The investigations were performed with fracture mechanics specimens of different sizes and geometries. The applied loading comprised cyclic loads, static loads and load spectra where the static load was periodically interrupted by partial unloading. With regard to water chemistry, the oxygen content (VVER) and impurities of sulphate and chlorides (BWR) were varied beyond allowable limits for continuous operation. The current paper summarises the most important crack growth results obtained under simulated BWR/NWC conditions. The results are discussed in the context of the current crack growth rate curves in the corresponding nuclear codes. (authors)

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

Dynamics and local boundary properties of the dawn-side magnetopause under conditions observed by Equator-S

Directory of Open Access Journals (Sweden)

M. W. Dunlop

Full Text Available Magnetic field measurements, taken by the magnetometer experiment (MAM on board the German Equator-S spacecraft, have been used to identify and categorise 131 crossings of the dawn-side magnetopause at low latitude, providing unusual, long duration coverage of the adjacent magnetospheric regions and near magnetosheath. The crossings occurred on 31 orbits, providing unbiased coverage over the full range of local magnetic shear from 06:00 to 10:40 LT. Apogee extent places the spacecraft in conditions associated with intermediate, rather than low, solar wind dynamic pressure, as it processes into the flank region. The apogee of the spacecraft remains close to the magnetopause for mean solar wind pressure. The occurrence of the magnetopause encounters are summarised and are found to compare well with predicted boundary location, where solar wind conditions are known. Most scale with solar wind pressure. Magnetopause shape is also documented and we find that the magnetopause orientation is consistently sunward of a model boundary and is not accounted for by IMF or local magnetic shear conditions. A number of well-established crossings, particularly those at high magnetic shear, or exhibiting unusually high-pressure states, were observed and have been analysed for their boundary characteristics and some details of their boundary and near magnetosheath properties are discussed. Of particular note are the occurrence of mirror-like signatures in the adjacent magnetosheath during a significant fraction of the encounters and a high number of multiple crossings over a long time period. The latter is facilitated by the spacecraft orbit which is designed to remain in the near magnetosheath for average solar wind pressure. For most encounters, a well-ordered, tangential (draped magnetosheath field is observed and there is little evidence of large deviations in local boundary orientations. Two passes corresponding to close conjunctions of the Geotail spacecraft

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.

Comparison of particle-wall interaction boundary conditions in the prediction of cyclone collection efficiency in computational fluid dynamics (CFD) modeling

International Nuclear Information System (INIS)

Valverde Ramirez, M.; Coury, J.R.; Goncalves, J.A.S.

2009-01-01

In recent years, many computational fluid dynamics (CFD) studies have appeared attempting to predict cyclone pressure drop and collection efficiency. While these studies have been able to predict pressure drop well, they have been only moderately successful in predicting collection efficiency. Part of the reason for this failure has been attributed to the relatively simple wall boundary conditions implemented in the commercially available CFD software, which are not capable of accurately describing the complex particle-wall interaction present in a cyclone. According, researches have proposed a number of different boundary conditions in order to improve the model performance. This work implemented the critical velocity boundary condition through a user defined function (UDF) in the Fluent software and compared its predictions both with experimental data and with the predictions obtained when using Fluent's built-in boundary conditions. Experimental data was obtained from eight laboratory scale cyclones with varying geometric ratios. The CFD simulations were made using the software Fluent 6.3.26. (author)

Effect of Reynolds Number and Periodic Unsteady Wake Flow Condition on Boundary Layer Development, Separation, and Intermittency Behavior Along the Suction Surface of a Low Pressure Turbine Blade

Science.gov (United States)

Schobeiri, M. T.; Ozturk, B.; Ashpis, David E.

2007-01-01

The paper experimentally studies the effects of periodic unsteady wake flow and different Reynolds numbers on boundary layer development, separation and re-attachment along the suction surface of a low pressure turbine blade. The experimental investigations were performed on a large scale, subsonic unsteady turbine cascade research facility at Turbomachinery Performance and Flow Research Laboratory (TPFL) of Texas A&M University. The experiments were carried out at Reynolds numbers of 110,000 and 150,000 (based on suction surface length and exit velocity). One steady and two different unsteady inlet flow conditions with the corresponding passing frequencies, wake velocities, and turbulence intensities were investigated. The reduced frequencies chosen cover the operating range of LP turbines. In addition to the unsteady boundary layer measurements, surface pressure measurements were performed. The inception, onset, and the extent of the separation bubble information collected from the pressure measurements were compared with the hot wire measurements. The results presented in ensemble-averaged, and the contour plot forms help to understand the physics of the separation phenomenon under periodic unsteady wake flow and different Reynolds number. It was found that the suction surface displayed a strong separation bubble for these three different reduced frequencies. For each condition, the locations defining the separation bubble were determined carefully analyzing and examining the pressure and mean velocity profile data. The location of the boundary layer separation was dependent of the Reynolds number. It is observed that starting point of the separation bubble and the re-attachment point move further downstream by increasing Reynolds number from 110,000 to 150,000. Also, the size of the separation bubble is smaller when compared to that for Re=110,000.

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

KAUST Repository

Reis, Tim; Dellar, Paul J.

2012-01-01

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

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

Energy Technology Data Exchange (ETDEWEB)

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

1996-12-01

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

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

International Nuclear Information System (INIS)

Staat, M.

1996-01-01

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

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.

Evaluation of nonequilibrium boundary conditions for hypersonic rarefied gas flows

Science.gov (United States)

Le, N. T. P.; Greenshields, Ch. J.; Reese, J. M.

2012-01-01

A new Computational Fluid Dynamics (CFD) solver for high-speed viscous §ows in the OpenFOAM code is validated against published experimental data and Direct Simulation Monte Carlo (DSMC) results. The laminar §at plate and circular cylinder cases are studied for Mach numbers, Ma, ranging from 6 to 12.7, and with argon and nitrogen as working gases. Simulation results for the laminar §at plate cases show that the combination of accommodation coefficient values σu = 0.7 and σT = 1.0 in the Maxwell/Smoluchowski conditions, and the coefficient values A1 = 1.5 and A2 = 1.0 in the second-order velocity slip condition, give best agreement with experimental data of surface pressure. The values σu = 0.7 and σT = 1.0 also give good agreement with DSMC data of surface pressure at the stagnation point in the circular cylinder case at Kn = 0.25. The Langmuir surface adsorption condition is also tested for the laminar §at plate case, but initial results were not as good as the Maxwell/Smoluchowski boundary conditions.

An optimal control method for fluid structure interaction systems via adjoint boundary pressure

Science.gov (United States)

Chirco, L.; Da Vià, R.; Manservisi, S.

2017-11-01

In recent year, in spite of the computational complexity, Fluid-structure interaction (FSI) problems have been widely studied due to their applicability in science and engineering. Fluid-structure interaction systems consist of one or more solid structures that deform by interacting with a surrounding fluid flow. FSI simulations evaluate the tensional state of the mechanical component and take into account the effects of the solid deformations on the motion of the interior fluids. The inverse FSI problem can be described as the achievement of a certain objective by changing some design parameters such as forces, boundary conditions and geometrical domain shapes. In this paper we would like to study the inverse FSI problem by using an optimal control approach. In particular we propose a pressure boundary optimal control method based on Lagrangian multipliers and adjoint variables. The objective is the minimization of a solid domain displacement matching functional obtained by finding the optimal pressure on the inlet boundary. The optimality system is derived from the first order necessary conditions by taking the Fréchet derivatives of the Lagrangian with respect to all the variables involved. The optimal solution is then obtained through a standard steepest descent algorithm applied to the optimality system. The approach presented in this work is general and could be used to assess other objective functionals and controls. In order to support the proposed approach we perform a few numerical tests where the fluid pressure on the domain inlet controls the displacement that occurs in a well defined region of the solid domain.

Integral Method of Boundary Characteristics: Neumann Condition

Science.gov (United States)

Kot, V. A.

2018-05-01

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

Solution of moving boundary problems with implicit boundary condition

International Nuclear Information System (INIS)

Moyano, E.A.

1990-01-01

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

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.

Grain boundary cavity growth under applied stress and internal pressure

International Nuclear Information System (INIS)

Mancuso, J.F.

1977-08-01

The growth of grain boundary cavities under applied stress and internal gas pressure was investigated. Methane gas filled cavities were produced by the C + 4H reversible CH4 reaction in the grain boundaries of type 270 nickel by hydrogen charging in an autoclave at 500 0 C with a hydrogen pressure of either 3.4 or 14.5 MPa. Intergranular fracture of nickel was achieved at a charging temperature of 300 0 C and 10.3 MPa hydrogen pressure. Cavities on the grain boundaries were observed in the scanning electron microscope after fracture. Photomicrographs of the cavities were produced in stereo pairs which were analyzed so as to correct for perspective distortion and also to determine the orientational dependence of cavity growth under an applied tensile stress

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

Wall-pressure fluctuations beneath a spatially evolving turbulent boundary layer

Science.gov (United States)

Mahesh, Krishnan; Kumar, Praveen

2016-11-01

Wall-pressure fluctuations beneath a turbulent boundary layer are important in applications dealing with structural deformation and acoustics. Simulations are performed for flat plate and axisymmetric, spatially evolving zero-pressure-gradient turbulent boundary layers at inflow Reynolds number of 1400 and 2200 based on momentum thickness. The simulations generate their own inflow using the recycle-rescale method. The results for mean velocity and second-order statistics show excellent agreement with the data available in literature. The spectral characteristics of wall-pressure fluctuations and their relation to flow structure will be discussed. This work is supported by ONR.

Parameterised Model of 2D Combustor Exit Flow Conditions for High-Pressure Turbine Simulations

Directory of Open Access Journals (Sweden)

Marius Schneider

2017-12-01

Full Text Available An algorithm is presented generating a complete set of inlet boundary conditions for Reynolds-averaged Navier–Stokes computational fluid dynamics (RANS CFD of high-pressure turbines to investigate their interaction with lean and rich burn combustors. The method shall contribute to understanding the sensitivities of turbine aerothermal performance in a systematic approach. The boundary conditions are based on a set of input parameters controlling velocity, temperature, and turbulence fields. All other quantities are derived from operating conditions and additional modelling assumptions. The algorithm is coupled with a CFD solver by applying the generated profiles as inlet boundary conditions. The successive steps to derive consistent flow profiles are described and results are validated against flow fields extracted from combustor CFD.

Temperature jump boundary conditions in radiation diffusion

International Nuclear Information System (INIS)

Alonso, C.T.

1976-12-01

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

Boundary condition histograms for modulated phases

International Nuclear Information System (INIS)

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

1997-11-01

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

Casimir pistons with general boundary conditions

Directory of Open Access Journals (Sweden)

Guglielmo Fucci

2015-02-01

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

A semi-analytical solution for viscothermal wave propagation in narrow gaps with arbitrary boundary conditions.

NARCIS (Netherlands)

Wijnant, Ysbrand H.; Spiering, R.M.E.J.; Blijderveen, M.; de Boer, Andries

2006-01-01

Previous research has shown that viscothermal wave propagation in narrow gaps can efficiently be described by means of the low reduced frequency model. For simple geometries and boundary conditions, analytical solutions are available. For example, Beltman [4] gives the acoustic pressure in the gap

Entropy Generation in Steady Laminar Boundary Layers with Pressure Gradients

Directory of Open Access Journals (Sweden)

Donald M. McEligot

2014-07-01

Full Text Available In an earlier paper in Entropy [1] we hypothesized that the entropy generation rate is the driving force for boundary layer transition from laminar to turbulent flow. Subsequently, with our colleagues we have examined the prediction of entropy generation during such transitions [2,3]. We found that reasonable predictions for engineering purposes could be obtained for flows with negligible streamwise pressure gradients by adapting the linear combination model of Emmons [4]. A question then arises—will the Emmons approach be useful for boundary layer transition with significant streamwise pressure gradients as by Nolan and Zaki [5]. In our implementation the intermittency is calculated by comparison to skin friction correlations for laminar and turbulent boundary layers and is then applied with comparable correlations for the energy dissipation coefficient (i.e., non-dimensional integral entropy generation rate. In the case of negligible pressure gradients the Blasius theory provides the necessary laminar correlations.

Peristaltic flow of Johnson-Segalman fluid in asymmetric channel with convective boundary conditions

Institute of Scientific and Technical Information of China (English)

H YASMIN; T HAYAT; A ALSAEDI; HH ALSULAMI

2014-01-01

This work is concerned with the peristaltic transport of the Johnson-Segalman fluid in an asymmetric channel with convective boundary conditions. The mathematical modeling is based upon the conservation laws of mass, linear momentum, and energy. The resulting equations are solved after long wavelength and low Reynolds number are used. The results for the axial pressure gradient, velocity, and temperature profiles are obtained for small Weissenberg number. The expressions of the pressure gra-dient, velocity, and temperature are analyzed for various embedded parameters. Pumping and trapping phenomena are also explored.

Formation of a Boundary-Free Dust Cluster in a Low-Pressure Gas-Discharge Plasma

International Nuclear Information System (INIS)

Usachev, A. D.; Zobnin, A. V.; Petrov, O. F.; Fortov, V. E.; Annaratone, B. M.; Thoma, M. H.; Hoefner, H.; Kretschmer, M.; Fink, M.; Morfill, G. E.

2009-01-01

An attraction between negatively charged micron-sized plastic particles was observed in the bulk of a low-pressure gas-discharge plasma under microgravity conditions. This attraction had led to the formation of a boundary-free dust cluster, containing one big central particle with a radius of about 6 μm and about 30 1 μm-sized particles situated on a sphere with a radius of 190 μm and with the big particle in the center. The stability of this boundary-free dust cluster was possible due to its confinement by the plasma flux on the central dust particle

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

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.

Integrable boundary conditions and modified Lax equations

International Nuclear Information System (INIS)

Avan, Jean; Doikou, Anastasia

2008-01-01

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

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

A domain-decomposition method to implement electrostatic free boundary conditions in the radial direction for electric discharges

Science.gov (United States)

Malagón-Romero, A.; Luque, A.

2018-04-01

At high pressure electric discharges typically grow as thin, elongated filaments. In a numerical simulation this large aspect ratio should ideally translate into a narrow, cylindrical computational domain that envelops the discharge as closely as possible. However, the development of the discharge is driven by electrostatic interactions and, if the computational domain is not wide enough, the boundary conditions imposed to the electrostatic potential on the external boundary have a strong effect on the discharge. Most numerical codes circumvent this problem by either using a wide computational domain or by calculating the boundary conditions by integrating the Green's function of an infinite domain. Here we describe an accurate and efficient method to impose free boundary conditions in the radial direction for an elongated electric discharge. To facilitate the use of our method we provide a sample implementation. Finally, we apply the method to solve Poisson's equation in cylindrical coordinates with free boundary conditions in both radial and longitudinal directions. This case is of particular interest for the initial stages of discharges in long gaps or natural discharges in the atmosphere, where it is not practical to extend the simulation volume to be bounded by two electrodes.

Electromagnetic stress at the boundary: Photon pressure or tension?

Science.gov (United States)

Wang, Shubo; Ng, Jack; Xiao, Meng; Chan, Che Ting

2016-03-01

It is well known that incident photons carrying momentum ℏk exert a positive photon pressure. But if light is impinging from a negative refractive medium in which ℏk is directed toward the source of radiation, should light exert a photon "tension" instead of a photon pressure? Using an ab initio method that takes the underlying microstructure of a material into account, we find that when an electromagnetic wave propagates from one material into another, the electromagnetic stress at the boundary is, in fact, indeterminate if only the macroscopic parameters are specified. Light can either pull or push the boundary, depending not only on the macroscopic parameters but also on the microscopic lattice structure of the polarizable units that constitute the medium. Within the context of an effective-medium approach, the lattice effect is attributed to electrostriction and magnetostriction, which can be accounted for by the Helmholtz stress tensor if we use the macroscopic fields to calculate the boundary optical stress.

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

Process-based investigation of cross-boundary environmental pressure from urban household consumption

International Nuclear Information System (INIS)

Yang, Dewei; Lin, Yanjie; Gao, Lijie; Sun, Yanwei; Wang, Run; Zhang, Guoqin

2013-01-01

Sustainability research at the city scale is increasingly focusing on urban household consumption in the context of global climate change. We use a complementary emergy accounting (EMA) and carbon footprint accounting (CFA) method to investigate the environmental pressure generated by household consumption in Xiamen, China. We distinguish between the resource extraction, consumption and disposal stages within an urban spatial conceptual framework, comprising the Urban Footprint Region (UFR) and Urban Sprawl Region (USR), and analyze five environmental footprint categories associated with cross-boundary household emergy and carbon flows. Cross-boundary activities, which link the USR with its UFR, contributed nearly 90% of total emergy and 70% of total GHG emissions in CFA. Transport fuel, building materials and food contribute most to environmental pressure in both EMA and CFA. The results indicate a significant cross-boundary resource burden and environmental footprint associated with household activities. The employed framework, method, and scope challenge the conventional spatial boundary of the urban system, and the results have important policy implications for urban sustainability and cross-boundary environmental management. - Highlights: ► We propose an urban spatial conceptual framework that includes USR and UFRs. ► A complementary EMA and CFA method is employed in urban household consumption system. ► Process-based cross-boundary environmental pressure of household consumption are evaluated. ► USR exerts pressure on its UFRs by extensive resource extraction and environmental emissions. ► We elucidate the USR–UFR environmental relationships and household energy policy

Einstein boundary conditions for the 3+1 Einstein equations

International Nuclear Information System (INIS)

Frittelli, Simonetta; Gomez, Roberto

2003-01-01

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

On the Boussinesq-Burgers equations driven by dynamic boundary conditions

Science.gov (United States)

Zhu, Neng; Liu, Zhengrong; Zhao, Kun

2018-02-01

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

Quadratic Functionals with General Boundary Conditions

International Nuclear Information System (INIS)

Dosla, Z.; Dosly, O.

1997-01-01

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

On filter boundary conditions in topology optimization

DEFF Research Database (Denmark)

Clausen, Anders; Andreassen, Erik

2017-01-01

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

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)

An Outflow Boundary Condition Model for Noninvasive Prediction of Fractional Flow Reserve in Diseased Coronary Arteries.

Science.gov (United States)

Fayssal, Iyad A; Moukalled, Fadl; Alam, Samir; Isma'eel, Hussain

2018-04-01

This paper reports on a new boundary condition formulation to model the total coronary myocardial flow and resistance characteristics of the myocardial vascular bed for any specific patient when considered for noninvasive diagnosis of ischemia. The developed boundary condition model gives an implicit representation of the downstream truncated coronary bed. Further, it is based on incorporating patient-specific physiological parameters that can be noninvasively extracted to account for blood flow demand to the myocardium at rest and hyperemic conditions. The model is coupled to a steady three-dimensional (3D) collocated pressure-based finite volume flow solver and used to characterize the "functional significance" of a patient diseased coronary artery segment without the need for predicting the hemodynamics of the entire arterial system. Predictions generated with this boundary condition provide a deep understanding of the inherent challenges behind noninvasive image-based diagnostic techniques when applied to human diseased coronary arteries. The overall numerical method and formulated boundary condition model are validated via two computational-based procedures and benchmarked with available measured data. The newly developed boundary condition is used via a designed computational methodology to (a) confirm the need for incorporating patient-specific physiological parameters when modeling the downstream coronary resistance, (b) explain the discrepancies presented in the literature between measured and computed fractional flow reserve (FFRCT), and (c) discuss the current limitations and future challenges in shifting to noninvasive assessment of ischemia.

Development of pressure boundaries leak detection technology for nuclear reactor

International Nuclear Information System (INIS)

Zhang Yao; Zhang Dafa; Chen Dengke; Zhang Liming

2008-01-01

The leak detection for the pressure boundaries is an important safeguard in nuclear reactor operation. In the paper, the status and the characters on the development of the pressure boundaries leak detection technology for the nuclear reactor were reviewed, especially, and the advance of the radiation leak detection technology and the acoustic emission leak detection technology were analyzed. The new advance trend of the leak detection technology was primarily explored. According to the analysis results, it is point out that the advancing target of the leak detection technology is to enhance its response speed, sensitivity, and reliability, and to provide effective information for operator and decision-maker. The realization of the global leak detection and the whole life cycle health monitoring for the nuclear boundaries is a significant advancing tendency of the leak detection technology. (authors)

Implicit Large-Eddy Simulations of Zero-Pressure Gradient, Turbulent Boundary Layer

Science.gov (United States)

Sekhar, Susheel; Mansour, Nagi N.

2015-01-01

A set of direct simulations of zero-pressure gradient, turbulent boundary layer flows are conducted using various span widths (62-630 wall units), to document their influence on the generated turbulence. The FDL3DI code that solves compressible Navier-Stokes equations using high-order compact-difference scheme and filter, with the standard recycling/rescaling method of turbulence generation, is used. Results are analyzed at two different Re values (500 and 1,400), and compared with spectral DNS data. They show that a minimum span width is required for the mere initiation of numerical turbulence. Narrower domains ((is) less than 100 w.u.) result in relaminarization. Wider spans ((is) greater than 600 w.u.) are required for the turbulent statistics to match reference DNS. The upper-wall boundary condition for this setup spawns marginal deviations in the mean velocity and Reynolds stress profiles, particularly in the buffer region.

Role of the vertical pressure gradient in wave boundary layers

DEFF Research Database (Denmark)

Jensen, Karsten Lindegård; Sumer, B. Mutlu; Vittori, Giovanna

2014-01-01

By direct numerical simulation (DNS) of the flow in an oscillatory boundary layer, it is possible to obtain the pressure field. From the latter, the vertical pressure gradient is determined. Turbulent spots are detected by a criterion involving the vertical pressure gradient. The vertical pressure...... gradient is also treated as any other turbulence quantity like velocity fluctuations and statistical properties of the vertical pressure gradient are calculated from the DNS data. The presence of a vertical pressure gradient in the near bed region has significant implications for sediment transport....

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.

Minimized Capillary End Effect During CO2 Displacement in 2-D Micromodel by Manipulating Capillary Pressure at the Outlet Boundary in Lattice Boltzmann Method

Science.gov (United States)

Kang, Dong Hun; Yun, Tae Sup

2018-02-01

We propose a new outflow boundary condition to minimize the capillary end effect for a pore-scale CO2 displacement simulation. The Rothman-Keller lattice Boltzmann method with multi-relaxation time is implemented to manipulate a nonflat wall and inflow-outflow boundaries with physically acceptable fluid properties in 2-D microfluidic chip domain. Introducing a mean capillary pressure acting at CO2-water interface to the nonwetting fluid at the outlet effectively prevents CO2 injection pressure from suddenly dropping upon CO2 breakthrough such that the continuous CO2 invasion and the increase of CO2 saturation are allowed. This phenomenon becomes most pronounced at capillary number of logCa = -5.5, while capillary fingering and massive displacement of CO2 prevail at low and high capillary numbers, respectively. Simulations with different domain length in homogeneous and heterogeneous domains reveal that capillary pressure and CO2 saturation near the inlet are reproducible compared with those with a proposed boundary condition. The residual CO2 saturation uniquely follows the increasing tendency with increasing capillary number, corroborated by experimental evidences. The determination of the mean capillary pressure and its sensitivity are also discussed. The proposed boundary condition is commonly applicable to other pore-scale simulations to accurately capture the spatial distribution of nonwetting fluid and corresponding displacement ratio.

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

CERN Document Server

Nier, Francis

2018-01-01

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

Adaptive boundary conditions for exterior flow problems

CERN Document Server

Boenisch, V; Wittwer, S

2003-01-01

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

Application of pressure-sensitive paint in shock-boundary layer interaction experiments

OpenAIRE

Seivwright, Douglas L.

1996-01-01

Approved for public release; distribution is unlimited A new type of pressure transducer, pressure-sensitive paint, was used to obtain pressure distributions associated with shock-boundary layer interaction. Based on the principle of photoluminescence and the process of oxygen quenching, pressure-sensitive paint provides a continous mapping of a pressure field over a surface of interest. The data measurement and acquisition system developed for use with the photoluminescence sensor was eva...

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.

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.

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.

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.

Single-phase liquid flow forced convection under a nearly uniform heat flux boundary condition in microchannels

KAUST Repository

Lee, Man

2012-02-22

A microchannel heat sink, integrated with pressure and temperature microsensors, is utilized to study single-phase liquid flow forced convection under a uniform heat flux boundary condition. Utilizing a waferbond-and-etch- back technology, the heat source, temperature and pressure sensors are encapsulated in a thin composite membrane capping the microchannels, thus allowing experimentally good control of the thermal boundary conditions. A three-dimensional physical model has been constructed to facilitate numerical simulations of the heat flux distribution. The results indicate that upstream the cold working fluid absorbs heat, while, within the current operating conditions, downstream the warmer working fluid releases heat. The Nusselt number is computed numerically and compared with experimental and analytical results. The wall Nusselt number in a microchannel can be estimated using classical analytical solutions only over a limited range of the Reynolds number, Re: both the top and bottom Nusselt numbers approach 4 for Re < 1, while the top and bottom Nusselt numbers approach 0 and 5.3, respectively, for Re > 100. The experimentally estimated Nusselt number for forced convection is highly sensitive to the location of the temperature measurements used in calculating the Nusselt number. © 2012 IOP Publishing Ltd.

Single-phase liquid flow forced convection under a nearly uniform heat flux boundary condition in microchannels

KAUST Repository

Lee, Man; Lee, Yi-Kuen; Zohar, Yitshak

2012-01-01

A microchannel heat sink, integrated with pressure and temperature microsensors, is utilized to study single-phase liquid flow forced convection under a uniform heat flux boundary condition. Utilizing a waferbond-and-etch- back technology, the heat source, temperature and pressure sensors are encapsulated in a thin composite membrane capping the microchannels, thus allowing experimentally good control of the thermal boundary conditions. A three-dimensional physical model has been constructed to facilitate numerical simulations of the heat flux distribution. The results indicate that upstream the cold working fluid absorbs heat, while, within the current operating conditions, downstream the warmer working fluid releases heat. The Nusselt number is computed numerically and compared with experimental and analytical results. The wall Nusselt number in a microchannel can be estimated using classical analytical solutions only over a limited range of the Reynolds number, Re: both the top and bottom Nusselt numbers approach 4 for Re < 1, while the top and bottom Nusselt numbers approach 0 and 5.3, respectively, for Re > 100. The experimentally estimated Nusselt number for forced convection is highly sensitive to the location of the temperature measurements used in calculating the Nusselt number. © 2012 IOP Publishing Ltd.

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.

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.

Entropy Stability and the No-Slip Wall Boundary Condition

KAUST Repository

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

2018-01-01

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

Entropy Stability and the No-Slip Wall Boundary Condition

KAUST Repository

Svärd, Magnus

2018-01-18

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

Exploring exotic states with twisted boundary conditions

International Nuclear Information System (INIS)

Agadjanov, Dimitri

2017-01-01

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

Exploring exotic states with twisted boundary conditions

Energy Technology Data Exchange (ETDEWEB)

Agadjanov, Dimitri

2017-09-11

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

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

Direct numerical simulation of thermally-stratified turbulent boundary layer subjected to adverse pressure gradient

International Nuclear Information System (INIS)

Hattori, Hirofumi; Kono, Amane; Houra, Tomoya

2016-01-01

Highlights: • We study various thermally-stratified turbulent boundary layers having adverse pressure gradient (APG) by means of DNS. • The detailed turbulent statistics and structures in various thermally-stratified turbulent boundary layers having APG are discussed. • It is found that the friction coefficient and Stanton number decrease along the streamwise direction due to the effects of stable thermal stratification and APG, but those again increase due to the APG effect in the case of weak stable thermal stratification. • In the case of strong stable stratification with or without APG, the flow separation is observed in the downstream region. - Abstract: The objective of this study is to investigate and observe turbulent heat transfer structures and statistics in thermally-stratified turbulent boundary layers subjected to a non-equilibrium adverse pressure gradient (APG) by means of direct numerical simulation (DNS). DNSs are carried out under conditions of neutral, stable and unstable thermal stratifications with a non-equilibrium APG, in which DNS results reveal heat transfer characteristics of thermally-stratified non-equilibrium APG turbulent boundary layers. In cases of thermally-stratified turbulent boundary layers affected by APG, heat transfer performances increase in comparison with a turbulent boundary layer with neutral thermal stratification and zero pressure gradient (ZPG). Especially, it is found that the friction coefficient and Stanton number decrease along the streamwise direction due to the effects of stable thermal stratification and APG, but those again increase due to the APG effect in the case of weak stable thermal stratification (WSBL). Thus, the analysis for both the friction coefficient and Stanton number in the case of WSBL with/without APG is conducted using the FIK identity in order to investigate contributions from the transport equations, in which it is found that both Reynolds-shear-stress and the mean convection terms

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.

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

A suitable boundary condition for bounded plasma simulation without sheath resolution

International Nuclear Information System (INIS)

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

1993-01-01

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

Nested Bethe Ansatz for Spin Ladder Model with Open Boundary Conditions

International Nuclear Information System (INIS)

Wu Junfang; Zhang Chunmin; Yue Ruihong; Li Runling

2005-01-01

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

Pressure-induced transition in the grain boundary of diamond

Science.gov (United States)

Chen, J.; Tang, L.; Ma, C.; Fan, D.; Yang, B.; Chu, Q.; Yang, W.

2017-12-01

Equation of state of diamond powder with different average grain sizes was investigated using in situ synchrotron x-ray diffraction and a diamond anvil cell (DAC). Comparison of compression curves was made for two samples with average grain size of 50nm and 100nm. The two specimens were pre-pressed into pellets and loaded in the sample pressure chamber of the DAC separately to minimized differences of possible systematic errors for the two samples. Neon gas was used as pressure medium and ruby spheres as pressure calibrant. Experiments were conducted at room temperature and high pressures up to 50 GPa. Fitting the compression data in the full pressure range into the third order Birch-Murnaghan equation of state yields bulk modulus (K) and its pressure derivative (K') of 392 GPa and 5.3 for 50nm sample and 398GPa and 4.5 for 100nm sample respectively. Using a simplified core-shell grain model, this result indicates that the grain boundary has an effective bulk modulus of 54 GPa. This value is similar to that observed for carbon nanotube[1] validating the recent theoretical diamond surface modeling[2]. Differential analysis of the compression cures demonstrates clear relative compressibility change at the pressure about 20 GPa. When fit the compression data below and above this pressure separately, the effect of grain size on bulk modulus reverses in the pressure range above 20 GPa. This observation indicates a possible transition of grain boundary structure, likely from sp2 hybridization at the surface[2] towards sp3like orbital structure which behaves alike the inner crystal. [1] Jie Tang, Lu-Chang Qin, Taizo Sasaki, Masako Yudasaka, Akiyuki Matsushita, and Sumio Iijima, Compressibility and Polygonization of Single-Walled Carbon Nanotubes under Hydrostatic Pressure, Physical Review Letters, 85(9), 1187-1198, 2000. [2] Shaohua Lu, Yanchao Wang, Hanyu Liu, Mao-sheng Miao, and Yanming Ma, Self-assembled ultrathin nanotubes on diamond (100) surface, Nature

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.

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.

Development of boundary layers

International Nuclear Information System (INIS)

Herbst, R.

1980-01-01

Boundary layers develop along the blade surfaces on both the pressure and the suction side in a non-stationary flow field. This is due to the fact that there is a strongly fluctuating flow on the downstream blade row, especially as a result of the wakes of the upstream blade row. The author investigates the formation of boundary layers under non-stationary flow conditions and tries to establish a model describing the non-stationary boundary layer. For this purpose, plate boundary layers are measured, at constant flow rates but different interferent frequency and variable pressure gradients. By introducing the sample technique, measurements of the non-stationary boundary layer become possible, and the flow rate fluctuation can be divided in its components, i.e. stochastic turbulence and periodical fluctuation. (GL) [de

Prevention of bolting degradation or failure in pressure boundary and support applications

International Nuclear Information System (INIS)

Merrick, E.A.; Rivers, A.; Bickford, J.; Marston, T.U.

1986-01-01

A discussion is presented of bolting degradation or failure experience in pressure boundary and component support applications in US commercial nuclear plants and the industry program to prevent failures in the future. The focus turns to steps which plant owners can take today to guard against pressure boundary bolt failure or degradation for existing plants or units being constructed. 'Tools' or products which the plant owner can expect from current industry programs which will be available in the near future to aid in understanding and improving bolting practices are described. (author)

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

International Nuclear Information System (INIS)

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

1982-01-01

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

On domain wall boundary conditions for the XXZ spin Hamiltonian

DEFF Research Database (Denmark)

Orlando, Domenico; Reffert, Susanne; Reshetikhin, Nicolai

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

Robust boundary treatment for open-channel flows in divergence-free incompressible SPH

Science.gov (United States)

Pahar, Gourabananda; Dhar, Anirban

2017-03-01

A robust Incompressible Smoothed Particle Hydrodynamics (ISPH) framework is developed to simulate specified inflow and outflow boundary conditions for open-channel flow. Being purely divergence-free, the framework offers smoothed and structured pressure distribution. An implicit treatment of Pressure Poison Equation and Dirichlet boundary condition is applied on free-surface to minimize error in velocity-divergence. Beyond inflow and outflow threshold, multiple layers of dummy particles are created according to specified boundary condition. Inflow boundary acts as a soluble wave-maker. Fluid particles beyond outflow threshold are removed and replaced with dummy particles with specified boundary velocity. The framework is validated against different cases of open channel flow with different boundary conditions. The model can efficiently capture flow evolution and vortex generation for random geometry and variable boundary conditions.

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

Directory of Open Access Journals (Sweden)

Nahed S. Hussein

2014-01-01

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

Simulations of QCD and QED with C* boundary conditions

Science.gov (United States)

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

2018-03-01

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

Critical effects of downstream boundary conditions on vortex breakdown

Science.gov (United States)

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

1992-01-01

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

Boundary conditions for plasma fluid models at the magnetic presheath entrance

International Nuclear Information System (INIS)

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

2012-01-01

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

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

International Nuclear Information System (INIS)

Jiang, Haiyan; Lu, Tiao; Cai, Wei

2014-01-01

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

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 symmetry of the boundary conditions of the volume potential

Science.gov (United States)

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

2017-09-01

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

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

International Nuclear Information System (INIS)

Jo, Chang Keun

1993-02-01

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

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

Science.gov (United States)

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

2017-08-01

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

Classically integrable boundary conditions for affine Toda field theories

International Nuclear Information System (INIS)

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

1995-01-01

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

Facilitating conditions for boundary-spanning behavior in governance networks

OpenAIRE

Meerkerk, Ingmar; Edelenbos, Jurian

2017-01-01

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

A Sharp-Interface Immersed Boundary Method with Improved Mass Conservation and Reduced Spurious Pressure Oscillations.

Science.gov (United States)

Seo, Jung Hee; Mittal, Rajat

2011-08-10

A method for reducing the spurious pressure oscillations observed when simulating moving boundary flow problems with sharp-interface immersed boundary methods (IBMs) is proposed. By first identifying the primary cause of these oscillations to be the violation of the geometric conservation law near the immersed boundary, we adopt a cut-cell based approach to strictly enforce geometric conservation. In order to limit the complexity associated with the cut-cell method, the cut-cell based discretization is limited only to the pressure Poisson and velocity correction equations in the fractional-step method and the small-cell problem tackled by introducing a virtual cell-merging technique. The method is shown to retain all the desirable properties of the original finite-difference based IBM while at the same time, reducing pressure oscillations for moving boundaries by roughly an order of magnitude.

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

Science.gov (United States)

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

2013-04-01

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

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

International Nuclear Information System (INIS)

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

2006-01-01

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

Scattering through a straight quantum waveguide with combined boundary conditions

Czech Academy of Sciences Publication Activity Database

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

2014-01-01

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

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

International Nuclear Information System (INIS)

Schuetz, G.

1992-11-01

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

Pressure estimation from single-snapshot tomographic PIV in a turbulent boundary layer

NARCIS (Netherlands)

Schneiders, J.F.G.; Pröbsting, S.; Dwight, R.P.; Van Oudheusden, B.W.; Scarano, F.

2016-01-01

A method is proposed to determine the instantaneous pressure field from a single tomographic PIV velocity snapshot and is applied to a flat-plate turbulent boundary layer. The main concept behind the single-snapshot pressure evaluation method is to approximate the flow acceleration using the

Consideration of loading conditions initiated by thermal transients in PWR pressure vessels

International Nuclear Information System (INIS)

Azodi; Glahn; Kersting; Schulz; Jansky.

1983-01-01

This report describes the present state of PWR-plants in the Federal Republic of Germany with respect to - the design of the primary pressure boundary - the analysis of thermal transients and resulting loads - the material conditions and neutron fluence - the requirements for protection against fast fracture. The experimental and analytical research and development programs are delineated together with some foreign R and D programs. It is shown that the parameters investigated (loading condition, crack shape and orientation etc.) cover a broad range. Extensive analytical investigations are emphasized. (orig./RW) [de

Gas migration in KBS-3 buffer bentonite. Sensitivity of test parameters to experimental boundary conditions

Energy Technology Data Exchange (ETDEWEB)

Harrington, J.F.; Horseman, S.T. [British Geological Survey, Nottingham (United Kingdom)

2003-01-01

In the current Swedish repository design concept, hydrogen gas can be generated inside a waste canister by anaerobic corrosion of the ferrous metal liner. If the gas generation rate exceeds the diffusion rate of gas molecules in the buffer porewater, gas will accumulate in the void-space of a canister until its pressure becomes large enough for it to enter the bentonite as a discrete gaseous phase. Three long tenn gas injection tests have been performed on cylinders of pre-compacted MX80 bentonite. Two of these tests were undertaken using a custom-designed constant volume and radial flow (CVRF) apparatus. Gas was injected at a centrally located porous filter installed in the clay before hydration. Arrangements were made for gas to flow to three independently monitored sink-filter arrays mounted around the specimen. Axial and radial total stresses and internal porewater pressures were continuously monitored. Breakthrough and peak gas pressures were substantially larger than the sum of the swelling pressure and the external porewater. The third test was performed. using an apparatus which radially constrains the specimen during gas flow. Observed sensitivity of the breakthrough and peak gas pressures to the test boundary conditions suggests that gas entry must be accompanied by dilation of the bentonite fabric. In other words, there is a tendency for the volume of the specimen to increase during this process. The experimental evidence is consistent with the flow of gas along a relatively small number of crack-like pathways which propagate through the clay as gas pressure increases. Gas entry and breakthrough under constant volume boundary conditions causes a substantial increase in the total stress and the internal porewater pressure. It is possible to determine the point at which gas enters the clay by monitoring changes in these parameters. Localisation of gas flow within multiple pathways results, in nonuniform discharge rates at the sinks. When gas injection

Gas migration in KBS-3 buffer bentonite. Sensitivity of test parameters to experimental boundary conditions

International Nuclear Information System (INIS)

Harrington, J.F.; Horseman, S.T.

2003-01-01

In the current Swedish repository design concept, hydrogen gas can be generated inside a waste canister by anaerobic corrosion of the ferrous metal liner. If the gas generation rate exceeds the diffusion rate of gas molecules in the buffer porewater, gas will accumulate in the void-space of a canister until its pressure becomes large enough for it to enter the bentonite as a discrete gaseous phase. Three long tenn gas injection tests have been performed on cylinders of pre-compacted MX80 bentonite. Two of these tests were undertaken using a custom-designed constant volume and radial flow (CVRF) apparatus. Gas was injected at a centrally located porous filter installed in the clay before hydration. Arrangements were made for gas to flow to three independently monitored sink-filter arrays mounted around the specimen. Axial and radial total stresses and internal porewater pressures were continuously monitored. Breakthrough and peak gas pressures were substantially larger than the sum of the swelling pressure and the external porewater. The third test was performed. using an apparatus which radially constrains the specimen during gas flow. Observed sensitivity of the breakthrough and peak gas pressures to the test boundary conditions suggests that gas entry must be accompanied by dilation of the bentonite fabric. In other words, there is a tendency for the volume of the specimen to increase during this process. The experimental evidence is consistent with the flow of gas along a relatively small number of crack-like pathways which propagate through the clay as gas pressure increases. Gas entry and breakthrough under constant volume boundary conditions causes a substantial increase in the total stress and the internal porewater pressure. It is possible to determine the point at which gas enters the clay by monitoring changes in these parameters. Localisation of gas flow within multiple pathways results, in nonuniform discharge rates at the sinks. When gas injection

Temperature field conduction solution by incomplete boundary condition

Energy Technology Data Exchange (ETDEWEB)

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

1977-01-01

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

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

International Nuclear Information System (INIS)

Frittelli, Simonetta; Gomez, Roberto

2004-01-01

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

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

Directory of Open Access Journals (Sweden)

Javier A. Dottori

2015-01-01

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

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

International Nuclear Information System (INIS)

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

2016-01-01

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

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

Science.gov (United States)

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

2016-04-01

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

Hypersonic Wind-Tunnel Measurements of Boundary-Layer Pressure Fluctuations

Science.gov (United States)

2009-08-01

Fluctuation Cone The Pressure-Fluctuation Cone was used for all wind-tunnel tests (Figure 3.7). The model is a 7◦ half-angle stainless - steel cone. It...analysis as a medium for fault detection: A review. Journal of Tribology , 130, January 2008. [80] L. M. Mack. Boundary layer linear stability theory. In

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.

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

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.

Temperature and stress distribution in pressure vessel by the boundary element method

International Nuclear Information System (INIS)

Alujevic, A.; Apostolovic, D.

1990-01-01

The aim of this paper is to demonstrate the applicability of boundary element method for the solution of temperatures and thermal stresses in the body of reactor pressure vessel of the NPP Krsko . In addition to the theory of boundary elements for thermo-elastic continua (2D, 3D) results are given of a numerically evaluated meridional cross-section. (author)

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.

Pressure tuning of the morphotropic phase boundary in piezoelectric lead zirconate titanate

International Nuclear Information System (INIS)

Rouquette, J.; Haines, J.; Bornand, V.; Pintard, M.; Papet, Ph.; Bousquet, C.; Konczewicz, L.; Gorelli, F. A.; Hull, S.

2004-01-01

Titanium-rich PZT solid solutions were studied under high pressure by neutron and x-ray diffraction, Raman spectroscopy and dielectric measurements. The results show that high pressure stabilizes the ferroelectric monoclinic phases, which are proposed to be responsible for the high piezoelectric properties characteristic of the morphotropic composition PbZr 0.52 Ti 0.48 O 3 . Pressure may thus be used to tune the morphotropic phase boundary in the composition-pressure plane to include a wide range of titanium-rich PZT compositions

Integrity of the reactor coolant boundary of the European pressurized water reactor (EPR)

Energy Technology Data Exchange (ETDEWEB)

Goetsch, D.; Bieniussa, K.; Schulz, H.; Jalouneix, J.

1997-04-01

This paper is an abstract of the work performed in the frame of the development of the IPSN/GRS approach in view of the EPR conceptual safety features. EPR is a pressurized water reactor which will be based on the experience gained by utilities and designers in France and in Germany. The reactor coolant boundary of a PWR includes the reactor pressure vessel (RPV), those parts of the steam generators (SGs) which contain primary coolant, the pressurizer (PSR), the reactor coolant pumps (RCPs), the main coolant lines (MCLs) with their branches as well as the other connecting pipes and all branching pipes including the second isolation valves. The present work covering the integrity of the reactor coolant boundary is mainly restricted to the integrity of the main coolant lines (MCLs) and reflects the design requirements for the main components of the reactor coolant boundary. In the following the conceptual aspects, i.e. design, manufacture, construction and operation, will be assessed. A main aspect is the definition of break postulates regarding overall safety implications.

Integrity of the reactor coolant boundary of the European pressurized water reactor (EPR)

International Nuclear Information System (INIS)

Goetsch, D.; Bieniussa, K.; Schulz, H.; Jalouneix, J.

1997-01-01

This paper is an abstract of the work performed in the frame of the development of the IPSN/GRS approach in view of the EPR conceptual safety features. EPR is a pressurized water reactor which will be based on the experience gained by utilities and designers in France and in Germany. The reactor coolant boundary of a PWR includes the reactor pressure vessel (RPV), those parts of the steam generators (SGs) which contain primary coolant, the pressurizer (PSR), the reactor coolant pumps (RCPs), the main coolant lines (MCLs) with their branches as well as the other connecting pipes and all branching pipes including the second isolation valves. The present work covering the integrity of the reactor coolant boundary is mainly restricted to the integrity of the main coolant lines (MCLs) and reflects the design requirements for the main components of the reactor coolant boundary. In the following the conceptual aspects, i.e. design, manufacture, construction and operation, will be assessed. A main aspect is the definition of break postulates regarding overall safety implications

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

Tuning of turbulent boundary layer anisotropy for improved surface pressure and trailing-edge noise modeling

DEFF Research Database (Denmark)

Bertagnolio, Franck; Fischer, Andreas; Zhu, Wei Jun

2014-01-01

The modeling of the surface pressure spectrum beneath a turbulent boundary layer is investigated, focusing on the case of airfoil flows and associated trailing edge noise prediction using the so-called TNO model. This type of flow is characterized by the presence of an adverse pressure gradient...... along the airfoil chord. It is shown that discrepancies between measurements and results from the TNO model increase as the pressure gradient increases. The original model is modified by introducing anisotropy in the definition of the turbulent vertical velocity spectrum across the boundary layer...

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.

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.

Analysis of grain boundaries, twin boundaries, and Te precipitates in CdZnTe grown by high-pressure Bridgeman method

International Nuclear Information System (INIS)

Heffelfinger, J.R.; Medlin, D.L.; James, R.B.

1998-03-01

Grain boundaries and twin boundaries in commercial Cd 1-x Zn x Te, which is prepared by a high pressure Bridgeman technique, have been investigated with transmission electron microscopy, scanning electron microscopy, infrared light microscopy and visible light microscopy. Boundaries inside these materials were found to be decorated with Te precipitates. The shape and local density of the precipitates were found to depend on the particular boundary. For precipitates that decorate grain boundaries, their microstructure was found to consist of a single, saucer shaped grain of hexagonal Te (space group P3 1 21). Analysis of a Te precipitate precipitates by selected area diffraction revealed the Te to be aligned with the surrounding Cd 1-x Zn x Te grains. This alignment was found to match the (111) Cd 1-x Z x Te planes with the (1 bar 101) planes of hexagonal Te. Crystallographic alignments between the Cd 1-x Zn x Te grains were also observed for a high angle grain boundary. The structure of the grain boundaries and the Te/Cd 1-x Zn x Te interface are discussed

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.

Transient integral boundary layer method to calculate the translesional pressure drop and the fractional flow reserve in myocardial bridges

Directory of Open Access Journals (Sweden)

Möhlenkamp Stefan

2006-06-01

Full Text Available Abstract Background The pressure drop – flow relations in myocardial bridges and the assessment of vascular heart disease via fractional flow reserve (FFR have motivated many researchers the last decades. The aim of this study is to simulate several clinical conditions present in myocardial bridges to determine the flow reserve and consequently the clinical relevance of the disease. From a fluid mechanical point of view the pathophysiological situation in myocardial bridges involves fluid flow in a time dependent flow geometry, caused by contracting cardiac muscles overlying an intramural segment of the coronary artery. These flows mostly involve flow separation and secondary motions, which are difficult to calculate and analyse. Methods Because a three dimensional simulation of the haemodynamic conditions in myocardial bridges in a network of coronary arteries is time-consuming, we present a boundary layer model for the calculation of the pressure drop and flow separation. The approach is based on the assumption that the flow can be sufficiently well described by the interaction of an inviscid core and a viscous boundary layer. Under the assumption that the idealised flow through a constriction is given by near-equilibrium velocity profiles of the Falkner-Skan-Cooke (FSC family, the evolution of the boundary layer is obtained by the simultaneous solution of the Falkner-Skan equation and the transient von-Kármán integral momentum equation. Results The model was used to investigate the relative importance of several physical parameters present in myocardial bridges. Results have been obtained for steady and unsteady flow through vessels with 0 – 85% diameter stenosis. We compare two clinical relevant cases of a myocardial bridge in the middle segment of the left anterior descending coronary artery (LAD. The pressure derived FFR of fixed and dynamic lesions has shown that the flow is less affected in the dynamic case, because the distal

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

Science.gov (United States)

Matzkin, A.

2018-03-01

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

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

Science.gov (United States)

Hejranfar, Kazem; Parseh, Kaveh

2017-09-01

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

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.

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

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.

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

Development of mapped stress-field boundary conditions based on a Hill-type muscle model.

Science.gov (United States)

Cardiff, P; Karač, A; FitzPatrick, D; Flavin, R; Ivanković, A

2014-09-01

Forces generated in the muscles and tendons actuate the movement of the skeleton. Accurate estimation and application of these musculotendon forces in a continuum model is not a trivial matter. Frequently, musculotendon attachments are approximated as point forces; however, accurate estimation of local mechanics requires a more realistic application of musculotendon forces. This paper describes the development of mapped Hill-type muscle models as boundary conditions for a finite volume model of the hip joint, where the calculated muscle fibres map continuously between attachment sites. The applied muscle forces are calculated using active Hill-type models, where input electromyography signals are determined from gait analysis. Realistic muscle attachment sites are determined directly from tomography images. The mapped muscle boundary conditions, implemented in a finite volume structural OpenFOAM (ESI-OpenCFD, Bracknell, UK) solver, are employed to simulate the mid-stance phase of gait using a patient-specific natural hip joint, and a comparison is performed with the standard point load muscle approach. It is concluded that physiological joint loading is not accurately represented by simplistic muscle point loading conditions; however, when contact pressures are of sole interest, simplifying assumptions with regard to muscular forces may be valid. Copyright © 2014 John Wiley & Sons, Ltd.

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

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

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.

Boundary layers affected by different pressure gradients investigated computationally by a zonal RANS-LES method

International Nuclear Information System (INIS)

Roidl, B.; Meinke, M.; Schröder, W.

2014-01-01

Highlights: • Reformulated synthetic turbulence generation method (RSTGM) is applied. • Zonal RANS-LES method is applied to boundary layers at pressure gradients. • Good agreement with the pure LES and other reference data is obtained. • The RSTGM is applicable to pressure gradient flows without modification. • RANS-to-LES boundary should be located where -1·10 6 6 is satisfied. -- Abstract: The reformulated synthetic turbulence generation (RSTG) method is used to compute by a fully coupled zonal RANS-LES approach turbulent non-zero-pressure gradient boundary layers. The quality of the RSTG method, which is based on the same shape functions and length scale distributions as in zero-pressure gradient flow, is discussed by comparing the zonal RANS-LES findings with pure LES, pure RANS, direct numerical simulation (DNS), and experimental data. For the favorable pressure gradient (FPG) simulation the RANS-to-LES transition occurs in the accelerated flow region and for the adverse pressure gradient (APG) case it is located in the decelerated flow region. The results of the time and spanwise averaged skin-friction distributions, velocity profiles, and Reynolds stress distributions of the zonal RANS-LES simulation show a satisfactory to good agreement with the pure LES, reference DNS, and experimental data. The quality of the findings shows that the rigorous formulation of the synthetic turbulence generation makes the RSTG method applicable without a priori knowledge of the flow properties but those determined by the RANS solution and without using additional control planes to regulate the shear stress budget to a wide range of Reynolds numbers and pressure gradients. The method is a promising approach to formulate embedded RANS-to-LES boundaries in flow regions where the Pohlhausen or acceleration parameter satisfies -1·10 -6 ⩽K⩽2·10 -6

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

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

Science.gov (United States)

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

2018-01-01

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

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

International Nuclear Information System (INIS)

Bhan, Jaemi; Kwon, Younghun

2007-01-01

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

A semi-implicit augmented IIM for Navier-Stokes equations with open, traction, or free boundary conditions.

Science.gov (United States)

Li, Zhilin; Xiao, Li; Cai, Qin; Zhao, Hongkai; Luo, Ray

2015-08-15

In this paper, a new Navier-Stokes solver based on a finite difference approximation is proposed to solve incompressible flows on irregular domains with open, traction, and free boundary conditions, which can be applied to simulations of fluid structure interaction, implicit solvent model for biomolecular applications and other free boundary or interface problems. For some problems of this type, the projection method and the augmented immersed interface method (IIM) do not work well or does not work at all. The proposed new Navier-Stokes solver is based on the local pressure boundary method, and a semi-implicit augmented IIM. A fast Poisson solver can be used in our algorithm which gives us the potential for developing fast overall solvers in the future. The time discretization is based on a second order multi-step method. Numerical tests with exact solutions are presented to validate the accuracy of the method. Application to fluid structure interaction between an incompressible fluid and a compressible gas bubble is also presented.

A semi-implicit augmented IIM for Navier–Stokes equations with open, traction, or free boundary conditions

Science.gov (United States)

Li, Zhilin; Xiao, Li; Cai, Qin; Zhao, Hongkai; Luo, Ray

2016-01-01

In this paper, a new Navier–Stokes solver based on a finite difference approximation is proposed to solve incompressible flows on irregular domains with open, traction, and free boundary conditions, which can be applied to simulations of fluid structure interaction, implicit solvent model for biomolecular applications and other free boundary or interface problems. For some problems of this type, the projection method and the augmented immersed interface method (IIM) do not work well or does not work at all. The proposed new Navier–Stokes solver is based on the local pressure boundary method, and a semi-implicit augmented IIM. A fast Poisson solver can be used in our algorithm which gives us the potential for developing fast overall solvers in the future. The time discretization is based on a second order multi-step method. Numerical tests with exact solutions are presented to validate the accuracy of the method. Application to fluid structure interaction between an incompressible fluid and a compressible gas bubble is also presented. PMID:27087702

New boundary conditions for 3D RF modelling

International Nuclear Information System (INIS)

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

1990-01-01

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

Antishocks in the ASEP with open boundaries conditioned on low current

International Nuclear Information System (INIS)

Belitsky, V; Schütz, G M

2013-01-01

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

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.

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

International Nuclear Information System (INIS)

Lindman, E.L.

1975-01-01

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

An Irrotational Flow Field That Approximates Flat Plate Boundary Conditions

OpenAIRE

Ruffa, Anthony A.

2004-01-01

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

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

DEFF Research Database (Denmark)

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

2014-01-01

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

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

Energy Technology Data Exchange (ETDEWEB)

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

2002-07-01

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

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

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.

Bruce Power's nuclear pressure boundary quality assurance program requirements, implementation and transition

International Nuclear Information System (INIS)

Krane, J.C.

2009-01-01

The development of a full scope nuclear pressure boundary quality assurance program in Canada requires extensive knowledge of the structure and detailed requirements of codes and standards published by the Canadian Standards Association (CSA) and American Society of Mechanical Engineers (ASME). Incorporation into company governance documents and implementation of these requirements while managing the transition to more recent revisions of these codes and standards represents a significant challenge for Bruce Power, Canada's largest independent nuclear operator. This paper explores the key developments and innovative changes that are used to ensure successful regulatory compliance and effective implementation of the Bruce Power Pressure Boundary Quality Assurance Program. Challenges and mitigating strategies to sustain this large compliance based program at Bruce Power's 8 unit nuclear power plant site will also be detailed. (author)

MDEP Technical Report TR-CSWG-03. Technical Report: fundamental attributes for the design and construction of reactor coolant pressure-boundary components

International Nuclear Information System (INIS)

2014-01-01

The primary, long-term goal of MDEP's CSWG is to achieve international harmonisation of codes and standards for pressure boundary components in nuclear power plants that are important to reactor safety. The key to achieving harmonisation is to understand the extent of similarities and differences amongst the pressure boundary codes and standards used in various countries. To assist the CSWG in its long-term goals, several standards development organisations (SDOs) from various countries performed a comparison of their pressure boundary codes and standards to identify the extent of similarities and differences in code requirements and the reasons for their differences. This CSWG document provides the fundamental attributes which have been developed for the codes and standards used in the design and construction of reactor coolant pressure boundary components in nuclear power plants. The fundamental attributes are the basic concepts to be considered in the design, materials, fabrication, installation, examination, testing and over-pressure protection requirements for pressure boundary components

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.

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

The Soil Water Characteristic Curve (SWCC) is a practical representation of the behavior of soil water by relating the suction (difference between the air and water pressures to the moisture content (water saturation). The SWCC is characterized by a hysteresis loop, which is thought to be unique in that any drainage-imbibition cycle lies within a main hysteresis loop limited by two different curves for drainage and imbibition. This 'uniqueness' is the main argument for considering the SWCC as a material-intrinsic feature that characterizes the pore structure and its interaction with fluids. Models have been developed with the SWCC as input data to describe the evolution of the water saturation and the suction within soils. One example of these models is the widely used Richard's equation [1]. In this work we present a series of numerical simulations to evaluate the 'unique' nature of the SWCC. The simulations involves the use of the Lattice Boltzmann Method (LBM) [2] within a regular soil, modelling the flow behavior of two immiscible fluids: wetting and non-wetting. The soil is packed within a cubic domain to resemble the experimental setups that are commonly used for measuring the SWCC[3]. The boundary conditions ensure that the non-wetting phase enters through one cubic face and the wetting phase enters trough the opposite phase, with no flow boundary conditions in the remaining 4 cubic faces. The SWCC known features are inspected including the presence of the common limit curves for different cycles involving varying limits for the suction. For this stage of simulations, the SWCC is indeed unique. Later, different boundary conditions are applied with the two fluids each injected from 3 opposing faces into the porous medium. The effect of this boundary condition change is a net flow direction, which is different from that in the previous case. A striking result is observed when both SWCC are compared and found to be noticeable different. Further analysis is

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”

Multiphysical Simulation of PT-CT Contact with Outer Boundary Condition

Energy Technology Data Exchange (ETDEWEB)

Chang, Se-Myong [Kunsan National Univ., Gunsan (Korea, Republic of); Kim, Hyoung Tae [Korea Atomic Energy Research Institute, Daejeon (Korea, Republic of)

2016-10-15

The present study is about preliminary calculation results for these ICSP activity works, where the COMSOL Multiphysics code is used to simulate plastic deformation of a pressure tube as a result of the interaction of stress and temperature. It is shown that the thermal stress model of COMSOL is compatible to simulate the multiple heat transfers (including the radiation heat transfer and heat conduction) and stress strain in the simplified 2-D problem. The benchmark test result for radiation heat transfer is in good agreement with the analytical solution for the concentric configuration of PT(pressure tube) and CT(calandria tube). In this paper, the authors did an open computation of these multi-physical phenomena by changing the outer boundary condition of CT according to the experimental result of ICSP. A series of simulation has been done based on the benchmark test proposed by IAEA/ICSP. The unsteady multi-physics was treated some numerical models with COMSOL. The comparison with CATHENA code is verified as a good agreement as we increase the accuracy of numerical method, Gaussian quadrature. The open computation for the validation of this numerical code is still on-going, and the temperature inside and outside the PT shows a very good agreement.

Multiphysical Simulation of PT-CT Contact with Outer Boundary Condition

International Nuclear Information System (INIS)

Chang, Se-Myong; Kim, Hyoung Tae

2016-01-01

The present study is about preliminary calculation results for these ICSP activity works, where the COMSOL Multiphysics code is used to simulate plastic deformation of a pressure tube as a result of the interaction of stress and temperature. It is shown that the thermal stress model of COMSOL is compatible to simulate the multiple heat transfers (including the radiation heat transfer and heat conduction) and stress strain in the simplified 2-D problem. The benchmark test result for radiation heat transfer is in good agreement with the analytical solution for the concentric configuration of PT(pressure tube) and CT(calandria tube). In this paper, the authors did an open computation of these multi-physical phenomena by changing the outer boundary condition of CT according to the experimental result of ICSP. A series of simulation has been done based on the benchmark test proposed by IAEA/ICSP. The unsteady multi-physics was treated some numerical models with COMSOL. The comparison with CATHENA code is verified as a good agreement as we increase the accuracy of numerical method, Gaussian quadrature. The open computation for the validation of this numerical code is still on-going, and the temperature inside and outside the PT shows a very good agreement

High-Reynolds-number turbulent-boundary-layer wall-pressure fluctuations with dilute polymer solutions

Science.gov (United States)

Elbing, Brian R.; Winkel, Eric S.; Ceccio, Steven L.; Perlin, Marc; Dowling, David R.

2010-08-01

Wall-pressure fluctuations were investigated within a high-Reynolds-number turbulent boundary layer (TBL) modified by the addition of dilute friction-drag-reducing polymer solutions. The experiment was conducted at the U.S. Navy's Large Cavitation Channel on a 12.9 m long flat-plate test model with the surface hydraulically smooth (k+<0.2) and achieving downstream-distance-based Reynolds numbers to 220×106. The polymer (polyethylene oxide) solution was injected into the TBL through a slot in the surface. The primary flow diagnostics were skin-friction drag balances and an array of flush-mounted dynamic pressure transducers 9.8 m from the model leading edge. Parameters varied included the free-stream speed (6.7, 13.4, and 20.2 m s-1) and the injection condition (polymer molecular weight, injection concentration, and volumetric injection flux). The behavior of the pressure spectra, convection velocity, and coherence, regardless of the injection condition, were determined primarily based on the level of drag reduction. Results were divided into two regimes dependent on the level of polymer drag reduction (PDR), nominally separated at a PDR of 40%. The low-PDR regime is characterized by decreasing mean-square pressure fluctuations and increasing convection velocity with increasing drag reduction. This shows that the decrease in the pressure spectra with increasing drag reduction is due in part to the moving of the turbulent structures from the wall. Conversely, with further increases in drag reduction, the high-PDR regime has negligible variation in the mean-squared pressure fluctuations and convection velocity. The convection velocity remains constant at approximately 10% above the baseline-flow convection velocity, which suggests that the turbulent structures no longer move farther from the wall with increasing drag reduction. In light of recent numerical work, the coherence results indicate that in the low-PDR regime, the turbulent structures are being elongated in

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

Science.gov (United States)

Yang, Chuan-Fu

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

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

DEFF Research Database (Denmark)

Svec, Oldrich; Skoček, Jan

2013-01-01

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

New Models for Velocity/Pressure-Gradient Correlations in Turbulent Boundary Layers

Science.gov (United States)

Poroseva, Svetlana; Murman, Scott

2014-11-01

To improve the performance of Reynolds-Averaged Navier-Stokes (RANS) turbulence models, one has to improve the accuracy of models for three physical processes: turbulent diffusion, interaction of turbulent pressure and velocity fluctuation fields, and dissipative processes. The accuracy of modeling the turbulent diffusion depends on the order of a statistical closure chosen as a basis for a RANS model. When the Gram-Charlier series expansions for the velocity correlations are used to close the set of RANS equations, no assumption on Gaussian turbulence is invoked and no unknown model coefficients are introduced into the modeled equations. In such a way, this closure procedure reduces the modeling uncertainty of fourth-order RANS (FORANS) closures. Experimental and direct numerical simulation data confirmed the validity of using the Gram-Charlier series expansions in various flows including boundary layers. We will address modeling the velocity/pressure-gradient correlations. New linear models will be introduced for the second- and higher-order correlations applicable to two-dimensional incompressible wall-bounded flows. Results of models' validation with DNS data in a channel flow and in a zero-pressure gradient boundary layer over a flat plate will be demonstrated. A part of the material is based upon work supported by NASA under award NNX12AJ61A.

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

Boundary-Layer Characteristics of Persistent Regional Haze Events and Heavy Haze Days in Eastern China

Directory of Open Access Journals (Sweden)

Peng Huaqing

2016-01-01

Full Text Available This paper analyzed the surface conditions and boundary-layer climate of regional haze events and heavy haze in southern Jiangsu Province in China. There are 5 types with the surface conditions which are equalized pressure (EQP, the advancing edge of a cold front (ACF, the base of high pressure (BOH, the backside of high pressure (BAH, the inverted trough of low pressure (INT, and saddle pressure (SAP with the haze days. At that time, 4 types are divided with the regional haze events and each of which has a different boundary-layer structure. During heavy haze, the surface mainly experiences EQP, ACF, BOH, BAH, and INT which also have different boundary-layer structures.

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

Science.gov (United States)

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

2017-12-01

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

Mixed problems for linear symmetric hyperbolic systems with characteristic boundary conditions

International Nuclear Information System (INIS)

Secchi, P.

1994-01-01

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

Numerical study of ambient pressure for laser-induced bubble near a rigid boundary

Science.gov (United States)

Li, BeiBei; Zhang, HongChao; Han, Bing; Lu, Jian

2012-07-01

The dynamics of the laser-induced bubble at different ambient pressures was numerically studied by Finite Volume Method (FVM). The velocity of the bubble wall, the liquid jet velocity at collapse, and the pressure of the water hammer while the liquid jet impacting onto the boundary are found to increase nonlinearly with increasing ambient pressure. The collapse time and the formation time of the liquid jet are found to decrease nonlinearly with increasing ambient pressure. The ratios of the jet formation time to the collapse time, and the displacement of the bubble center to the maximal radius while the jet formation stay invariant when ambient pressure changes. These ratios are independent of ambient pressure.

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.

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

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.

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

International Nuclear Information System (INIS)

Follin, S.

1999-06-01

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

Boundary 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

Experimental study on effects of inlet boundary layer thickness and boundary layer fence in a turbine cascade

International Nuclear Information System (INIS)

Jun, Y. M.; Chung, J. T.

2000-01-01

The working fluid from the combustor to the turbine stage of a gas turbine makes various boundary layer thickness. Since the inlet boundary layer thickness is one of the important factors that affect the turbine efficiency, It is necessary to investigate secondary flow and loss with various boundary layer thickness conditions. In the present study, the effect of various inlet boundary layer thickness on secondary flow and loss and the proper height of the boundary layer fences for various boundary layer thickness were investigated. Measurements of secondary flow velocity and total pressure loss within and downstream of the passage were taken under 5 boundary layer thickness conditions, 16, 36, 52, 69, 110mm. It was found that total pressure loss and secondary flow areas were increased with increase of thickness but they were maintained almost at the same position. At the following research about the boundary layer fences, 1/6, 1/3, 1/2 of each inlet boundary layer thickness and 12mm were used as the fence heights. As a result, it was observed that the proper height of the fences was generally constant since the passage vortex remained almost at the same position. Therefore once the geometry of a cascade is decided, the location of the passage vortex and the proper fence height are appeared to be determined at the same time. When the inlet boundary layer thickness is relatively small, the loss caused by the proper fence becomes bigger than end wall loss so that it dominates secondary loss. In these cases the proper fence height is decided not by the cascade geometry but by the inlet boundary layer thickness as previous investigations

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.

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

Directory of Open Access Journals (Sweden)

Kyncl Martin

2017-01-01

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

Homogenized boundary conditions and resonance effects in Faraday cages

OpenAIRE

Hewett, DP; Hewitt, IJ

2016-01-01

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

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

Science.gov (United States)

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

2018-05-01

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

On the Unsteadiness of a Transitional Shock Wave-Boundary Layer Interaction Using Fast-Response Pressure-Sensitive Paint

Science.gov (United States)

Lash, E. Lara; Schmisseur, John

2017-11-01

Pressure-sensitive paint has been used to evaluate the unsteady dynamics of transitional and turbulent shock wave-boundary layer interactions generated by a vertical cylinder on a flat plate in a Mach 2 freestream. The resulting shock structure consists of an inviscid bow shock that bifurcates into a separation shock and trailing shock. The primary features of interest are the separation shock and an upstream influence shock that is intermittently present in transitional boundary layer interactions, but not observed in turbulent interactions. The power spectral densities, frequency peaks, and normalized wall pressures are analyzed as the incoming boundary layer state changes from transitional to fully turbulent, comparing both centerline and outboard regions of the interaction. The present study compares the scales and frequencies of the dynamics of the separation shock structure in different boundary layer regimes. Synchronized high-speed Schlieren imaging provides quantitative statistical analyses as well as qualitative comparisons to the fast-response pressure sensitive paint measurements. Materials based on research supported by the U.S. Office of Naval Research under Award Number N00014-15-1-2269.

F-8 supercritical wing flight pressure, Boundary layer, and wake measurements and comparisons with wind tunnel data

Science.gov (United States)

Montoya, L. C.; Banner, R. D.

1977-01-01

Data for speeds from Mach 0.50 to Mach 0.99 are presented for configurations with and without fuselage area-rule additions, with and without leading-edge vortex generators, and with and without boundary-layer trips on the wing. The wing pressure coefficients are tabulated. Comparisons between the airplane and model data show that higher second velocity peaks occurred on the airplane wing than on the model wing. The differences were attributed to wind tunnel wall interference effects that caused too much rear camber to be designed into the wing. Optimum flow conditions on the outboard wing section occurred at Mach 0.98 at an angle of attack near 4 deg. The measured differences in section drag with and without boundary-layer trips on the wing suggested that a region of laminar flow existed on the outboard wing without trips.

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

International Nuclear Information System (INIS)

Nagai, K.; Hara, J.

1988-01-01

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

The neutron transport with general boundary conditions (II)

International Nuclear Information System (INIS)

Boulanouar, Mohamed

2012-01-01

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

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.

Acoustic emission and estimation of flaw significance in reactor pressure boundaries

International Nuclear Information System (INIS)

Hutton, P.H.; Kurtz, R.J.

1982-01-01

The work discussed is intended to establish the feasibility of using acoustic emission (AE) to detect and evaluate growing flaws in nuclear reactor pressure boundaries. Basic AE identification and interpretation methods have grown out of Phase 1. Phases 2 and 3 to test and demonstrate developed methodology on a vessel test and on a reactor are in progress

Direct measurement of methane hydrate composition along the hydrate equilibrium boundary

Science.gov (United States)

Circone, S.; Kirby, S.H.; Stern, L.A.

2005-01-01

The composition of methane hydrate, namely nW for CH 4??nWH2O, was directly measured along the hydrate equilibrium boundary under conditions of excess methane gas. Pressure and temperature conditions ranged from 1.9 to 9.7 MPa and 263 to 285 K. Within experimental error, there is no change in hydrate composition with increasing pressure along the equilibrium boundary, but nW may show a slight systematic decrease away from this boundary. A hydrate stoichiometry of n W = 5.81-6.10 H2O describes the entire range of measured values, with an average composition of CH4??5.99(??0.07) H2O along the equilibrium boundary. These results, consistent with previously measured values, are discussed with respect to the widely ranging values obtained by thermodynamic analysis. The relatively constant composition of methane hydrate over the geologically relevant pressure and temperature range investigated suggests that in situ methane hydrate compositions may be estimated with some confidence. ?? 2005 American Chemical Society.

The role of mechanical boundary conditions in the soft mode dynamics of PbTiO3.

Science.gov (United States)

McCash, Kevin; Mani, B K; Chang, C-M; Ponomareva, I

2014-10-29

The role of different mechanical boundary conditions in the soft mode dynamics of ferroelectric PbTiO3 is systematically investigated using first-principles-based simulations and analytical model. The change in the soft mode dynamics due to hydrostatic pressure, uniaxial and biaxial stresses and biaxial strains is studied in a wide temperature range. Our computations predict: (i) the existence of Curie-Weiss laws that relate the soft mode frequency to the stress or strain; (ii) a non-trivial temperature evolution of the associated Curie-Weiss constants; (iii) a qualitative difference between the soft mode response to stresses/strains and hydrostatic pressure. The latter finding implies that the Curie-Weiss pressure law commonly used for residual stress estimation may not apply for the cases of uniaxial and biaxial stresses and strains. On the other hand, our systematic study offers a way to eliminate this difficulty through the establishment of Curie-Weiss stress and strain laws. Implications of our predictions for some available experimental data are discussed.

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

Directory of Open Access Journals (Sweden)

Leendert van Maanen

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

The Hardy inequality with boundary or intermediate conditions

Czech Academy of Sciences Publication Activity Database

Kufner, Alois

2017-01-01

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

The Hardy inequality with boundary or intermediate conditions

Czech Academy of Sciences Publication Activity Database

Kufner, Alois

2017-01-01

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

Initial boundary-value problem for the spherically symmetric Einstein equations with fluids with tangential pressure.

Science.gov (United States)

Brito, Irene; Mena, Filipe C

2017-08-01

We prove that, for a given spherically symmetric fluid distribution with tangential pressure on an initial space-like hypersurface with a time-like boundary, there exists a unique, local in time solution to the Einstein equations in a neighbourhood of the boundary. As an application, we consider a particular elastic fluid interior matched to a vacuum exterior.

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.

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

Correlations for modeling transitional boundary layers under influences of freestream turbulence and pressure gradient

International Nuclear Information System (INIS)

Suluksna, Keerati; Dechaumphai, Pramote; Juntasaro, Ekachai

2009-01-01

This paper presents mathematical expressions for two significant parameters which control the onset location and length of transition in the γ-Re θ transition model of Menter et al. [Menter, F.R., Langtry, R.B., Volker, S., Huang, P.G., 2005. Transition modelling for general purpose CFD codes. In: ERCOFTAC International Symposium on Engineering Turbulence Modelling and Measurements]. The expressions are formulated and calibrated by means of numerical experiments for predicting transitional boundary layers under the influences of freestream turbulence and pressure gradient. It was also found that the correlation for transition momentum thickness Reynolds number needs only to be expressed in terms of local turbulence intensity, so that the more complex form that includes pressure gradient effects is unnecessary. Transitional boundary layers on a flat plate both with and without pressure gradients are employed to assess the performance of these two expressions for predicting the transition. The results show that the proposed expressions can work well with the model of Menter et al. (2005)

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.

The response of the Goddard general circulation model to sea ice boundary conditions

Science.gov (United States)

Herman, G.; Johnson, W. T.

1979-01-01

The effect of variation in the location of Arctic sea ice boundaries on the model's mean monthly climatology was examined. When sea ice boundaries were at their maximum extent the differences resulted in the January-February climatology. Sea level pressure was higher over the Barents Sea, in the Davis Strait, and in the Sea of Okhotsk. Pressure was lower by as much as 8 mb in the North Atlantic between Iceland and the British Isles, and in the Gulf of Alaska. Pressure rises in the eastern subtropical regions of the North Atlantic and North Pacific accompanied pressure falls in the Gulf of Alaska and Icelandic region. Geopotential heights at 500 mb were more than 100 gpm lower in the Bering Sea, and more than 120 gpm lower in the Icelandic region. Zonally averaged temperatures were cooler by 4 deg C below 3800 mb between 50 deg and 70 deg N with little change elsewhere. Zonally averaged geopotentials were lower by as much as 70 gpm in the mid-troposphere between 50/-70 deg N and zonal winds increased by as much as 3 m s in the mid-troposphere between 35/-50 deg N.

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

Pressure and tension waves from bubble collapse near a solid boundary: A numerical approach.

Science.gov (United States)

Lechner, Christiane; Koch, Max; Lauterborn, Werner; Mettin, Robert

2017-12-01

The acoustic waves being generated during the motion of a bubble in water near a solid boundary are calculated numerically. The open source package OpenFOAM is used for solving the Navier-Stokes equation and extended to include nonlinear acoustic wave effects via the Tait equation for water. A bubble model with a small amount of gas is chosen, the gas obeying an adiabatic law. A bubble starting from a small size with high internal pressure near a flat, solid boundary is studied. The sequence of events from bubble growth via axial microjet formation, jet impact, annular nanojet formation, torus-bubble collapse, and bubble rebound to second collapse is described. The different pressure and tension waves with their propagation properties are demonstrated.

Beta limitation of matter-antimatter boundary layers

International Nuclear Information System (INIS)

Lehnert, B.

1987-08-01

A model has earlier been proposed for a boundary layer which separates a cloud of matter from one of antimatter in a magnetized ambiplasma. In this model steady pressure equilibrium ceases to exist when a certain beta limit is exceeded. The latter is defined as the ratio between the ambiplasma and magnetic field pressures which balance each other in the boundary layer. Thus, at an increasing density, the high-energy particles created by annihilation within the layer are 'pumped up' to a pressure which cannot be balanced by a given magnetic field. The boundary layer then 'disrupts'. The critical beta limit thus obtained falls within the observed parameter ranges of galaxies and other large cosmical objects. Provided that the considered matter-antimatter balance holds true, this limit is thus expected to impose certain existence conditions on matter-antimatter boundary layers. Such a limitation may apply to certain cosmical objects and cosmological models. The maximum time scale for the corresponding disruption development has been estimated to be in the range from about 10 -4 to 10 2 seconds for boundary layers at ambiplasma particle densities in the range from 10 4 to 10 -2 m -3 , respectively. (author)

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.

Acoustic Radiation From a Mach 14 Turbulent Boundary Layer

Science.gov (United States)

Zhang, Chao; Duan, Lian; Choudhari, Meelan M.

2016-01-01

Direct numerical simulations (DNS) are used to examine the turbulence statistics and the radiation field generated by a high-speed turbulent boundary layer with a nominal freestream Mach number of 14 and wall temperature of 0:18 times the recovery temperature. The flow conditions fall within the range of nozzle exit conditions of the Arnold Engineering Development Center (AEDC) Hypervelocity Tunnel No. 9 facility. The streamwise domain size is approximately 200 times the boundary-layer thickness at the inlet, with a useful range of Reynolds number corresponding to Re 450 ?? 650. Consistent with previous studies of turbulent boundary layer at high Mach numbers, the weak compressibility hypothesis for turbulent boundary layers remains applicable under this flow condition and the computational results confirm the validity of both the van Driest transformation and Morkovin's scaling. The Reynolds analogy is valid at the surface; the RMS of fluctuations in the surface pressure, wall shear stress, and heat flux is 24%, 53%, and 67% of the surface mean, respectively. The magnitude and dominant frequency of pressure fluctuations are found to vary dramatically within the inner layer (z/delta 0.< or approx. 0.08 or z+ < or approx. 50). The peak of the pre-multiplied frequency spectrum of the pressure fluctuation is f(delta)/U(sub infinity) approx. 2.1 at the surface and shifts to a lower frequency of f(delta)/U(sub infinity) approx. 0.7 in the free stream where the pressure signal is predominantly acoustic. The dominant frequency of the pressure spectrum shows a significant dependence on the freestream Mach number both at the wall and in the free stream.

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

Directory of Open Access Journals (Sweden)

A Nemati

2013-10-01

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

Molecular Dynamics with Helical Periodic Boundary Conditions

Czech Academy of Sciences Publication Activity Database

Kessler, Jiří; Bouř, Petr

2014-01-01

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

Boundary conditions in conformal and integrable theories

CERN Document Server

Petkova, V B

2000-01-01

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

Aging of elastomers in CANDU pressure boundary service

International Nuclear Information System (INIS)

VanBerlo, C.; Leidner, J.

1987-09-01

This report describes the properties and aging of elastomers, and examines the performance of major elastomeric components in CANDU pressure boundary service. The components examined are vacuum building roof seals, pressure relief duct seals, airlock door seals, fuelling machine hoses, and cable penetrations. For each of these components, the design requirements, technical specifications and component testing procedures are compared with applicable standards. Information on actual and recommended monitoring and maintenance methods is presented. Operational and environmental stressors are identified. Component failure modes, causes and frequencies are described, as well as the remedial action taken. Many different elastomers are used in CANDU plants, for many different applications. Standards and manufacturers' recommendations are not consistent and may vary from one component to another. Accordingly, the monitoring, maintenance and replacement practices tend to vary from one application to another, and may also be different at different stations. Recommendations are given in this report for improved monitoring and maintenance, in an attempt to provide more consistency in approach. A summary of some experiences with elastomers from non-Canadian sources is contained in the last section. 125 refs

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

International Nuclear Information System (INIS)

Dyer, Ethan; Hinterbichler, Kurt

2009-01-01

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

Role of open boundary conditions on the hysteretic behaviour of one-dimensional spin crossover nanoparticles

Energy Technology Data Exchange (ETDEWEB)

Chiruta, Daniel [GEMaC, Université de Versailles Saint-Quentin-en-Yvelines, CNRS-UVSQ (UMR 8635), 78035 Versailles Cedex (France); Faculty of Electrical Engineering and Computer Science, Stefan cel Mare University, Suceava 720229 (Romania); Linares, Jorge, E-mail: jorge.linares@uvsq.fr, E-mail: miya@spin.phys.s.u-tokyo.ac.jp; Boukheddaden, Kamel [GEMaC, Université de Versailles Saint-Quentin-en-Yvelines, CNRS-UVSQ (UMR 8635), 78035 Versailles Cedex (France); Miyashita, Seiji, E-mail: jorge.linares@uvsq.fr, E-mail: miya@spin.phys.s.u-tokyo.ac.jp [Department of Physics, Graduate School of Science, University of Tokyo, Bunkyo-Ku, Tokyo (Japan); CREST, JST, 4-1-8 Honcho Kawaguchi, Saitama 332-0012 (Japan)

2014-05-21

In order to explain clearly the role of the open boundary conditions (OBCs) on phase transition in one dimensional system, we consider an Ising model with both short-range (J) and long-range (G) interactions, which has allowed us to study the cooperative nature of spin-crossover (SCO) materials at the nanometer scale. At this end, we developed a transfer-matrix method for one-dimensional (1D) SCO system with free boundary conditions, and we give numerical evidences for how the thermal spin transition curves vary as a function of the physical parameters (J, G) or an applied pressure. Moreover for OBCs case, we have derived the bulk, surface and finite-size contributions to the free energy and we have investigated the variation of these energies as function of J and system size. We have found that the surface free energy behaves like J〈σ〉{sup 2}, where 〈σ〉 is the average magnetization per site. Since the properties of the nanometric scale are dramatically influenced by the system's size (N), our analytical outcomes for the size dependence represent a step to achieve new characteristic of the future devices and also a way to find various novel properties which are absent in the bulk materials.

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

International Nuclear Information System (INIS)

Tessarotto, Massimo; Cremaschini, Claudio

2014-01-01

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

CFD simulation of pressure and discharge surge in Francis turbine at off-design conditions

International Nuclear Information System (INIS)

Chirkov, D; Avdyushenko, A; Panov, L; Bannikov, D; Cherny, S; Skorospelov, V; Pylev, I

2012-01-01

A hybrid 1D-3D CFD model is developed for the numerical simulation of pressure and discharge surge in hydraulic power plants. The most essential part – the turbine itself – is simulated directly using 3D unsteady equations of turbulent motion of fluid-vapor mixture, while the rest of the hydraulic system is simulated in frames of 1D hydro-acoustic model. Thus the model accounts for the main factors responsible for excitation and propagation of pressure and discharge waves in hydraulic power plant. Boundary conditions at penstock inlet and draft tube outlet are discussed in detail. Then simulations of dynamic behavior at part load and full load operating points are performed. It is shown that the numerical model is able to capture self-excited oscillations in full load conditions. The influence of penstock length and flow structure behind the runner are investigated. The presented approach seems to be a promising tool for prediction and investigation the dynamic behavior in hydraulic power plants.

On some boundary value problems in quantum statistical mechanics

International Nuclear Information System (INIS)

Angelescu, N.

1978-01-01

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

Crack growth behaviour of low alloy steels for pressure boundary components under transient light water reactor operating conditions (CASTOC)

International Nuclear Information System (INIS)

Foehl, J.; Weissenberg, T.; Gomez-Briceno, D.; Lapena, J.; Ernestova, M.; Zamboch, M.; Seifert, H.P.; Ritter, S.; Roth, A.; Devrient, B.; Ehrnsten, U.

2004-01-01

The CASTOC project addresses environmentally assisted cracking (EAC) phenomena in low alloy steels used for pressure boundary components in both Western type boiling water reactors (BWR) and Russian type pressurised water reactors (VVER). It comprises the four work packages (WP): inter-laboratory comparison test (WP1); EAC behaviour under static load (WP2), EAC behaviour under cyclic load and load transients (WP3); evaluation of the results with regard to their relevance for components in practice (WP4). The use of sophisticated test facilities and measurement techniques for the on-line detection of crack advances have provided a more detailed understanding of the mechanisms of environmentally assisted cracking and provided quantitative data of crack growth rates as a function of loading events and time, respectively. The effect of several major parameters controlling EAC was investigated with particular emphasis on the transferability of the results to components in service. The obtained crack growth rate data were reflected on literature data and on commonly applied prediction curves as presented in the appropriate Code. At relevant stress intensity factors it could be shown that immediate cessation of growing cracks occurs after changing from cyclic to static load in high purity oxygenated BWR water and oxygen-free VVER water corresponding to steady state operation conditions. Susceptibility to environmentally assisted cracking under static load was observed for a heat affected zone material in oxygenated high purity water and also in base materials during a chloride transient representing BWR water condition below Action Level 1 of the EPRI Water Chemistry Guidelines according to the lectrical conductivity of the water but in the range of Action Level 2 according to the content of chlorides. Time based crack growth was also observed in one Russian type base material in oxygenated VVER water and in one Western type base material in oxygenated high purity BWR

Local characteristics of the nocturnal boundary layer in response to external pressure forcing

NARCIS (Netherlands)

van der Linden, S.J.A.; Baas, P.; van Hooft, J.A.; van Hooijdonk, I.G.S.; Bosveld, F.C.; van de Wiel, B.J.H.

2017-01-01

Geostrophic wind speed data, derived from pressure observations, are used in combination with tower measurements to investigate the nocturnal stable boundary layer at Cabauw (The Netherlands). Since the geostrophic wind speed is not directly influenced by local nocturnal stability, it may be

Slarti: A boundary condition editor for a coupled climate model

Science.gov (United States)

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

2006-12-01

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

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

OpenAIRE

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

2006-01-01

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

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

KAUST Repository

Leise, Tanya

2008-11-01

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

The online sealing performance test of the primary circuit pressure boundary check valve in nuclear power plants

International Nuclear Information System (INIS)

Yang Yunfei; Huang Huimin

2013-01-01

The primary circuit pressure boundary check valves of 320 MW pressurized water reactor is a nuclear grade I key equipment. The sealing demand is very high, which is directly related to the internal leakage rate of the primary circuit system. After the welding check valve is repaired, the sealing performance is judged by color printing checks. If there is water or humid vapor in the pipe, it will affect the accuracy of the color printing checks. For the particularity of the online check valve tightness test, online detecting device is designed by the hydraulic pressure drop method in other nuclear power plants, but the method has some shortcomings and restrictions. In this paper, we design a reliable and portable test equipment by the low-pressure gas seal test flow measurement, which make accurate and quantitative judgment of sealing property after the pressure boundary check valves are repaired. (authors)

Matrix factorisations for rational boundary conditions by defect fusion

International Nuclear Information System (INIS)

Behr, Nicolas; Fredenhagen, Stefan

2015-01-01

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

Matrix factorisations for rational boundary conditions by defect fusion

Energy Technology Data Exchange (ETDEWEB)

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

2015-05-11

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

Comparison of mass transport using average and transient rainfall boundary conditions

International Nuclear Information System (INIS)

Duguid, J.O.; Reeves, M.

1976-01-01

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

Generating wind fluctuations for Large Eddy Simulation inflow boundary condition

International Nuclear Information System (INIS)

Bekele, S.A.; Hangan, H.

2004-01-01

Large Eddy Simulation (LES) studies of flows over bluff bodies immersed in a boundary layer wind environment require instantaneous wind characteristics. The influences of the wind environment on the building pressure distribution are a well-established fact in the experimental study of wind engineering. Measured wind data of full or model scale are available only at a limited number of points. A method of obtaining instantaneous wind data at all mesh points of the inlet boundary for LES computation is necessary. Herein previous and new wind inflow generation techniques are presented. The generated wind data is then applied to a LES computation of a channel flow. The characteristics of the generated wind fluctuations in comparison to the measured data and the properties of the flow field computed from these two wind data are discussed. (author)

Natural frequency analysis of fluid conveying pipeline with different boundary conditions

International Nuclear Information System (INIS)

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

2010-01-01

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

A new boundary control scheme for simultaneous achievement of H-mode and radiative cooling (SHC boundary)

International Nuclear Information System (INIS)

Ohyabu, N.

1995-05-01

We have proposed a new boundary control scheme (SHC boundary), which could allow simultaneous achievement of the H-mode type confinement improvement and radiative cooling with wide heat flux distribution. In our proposed configuration, a low m island layer sharply separates a plasma confining region from an open 'ergodic' boundary. The degree of openness in the ergodic boundary must be high enough to make the plasma pressure constant along the field line, which in turn separates low density plasma just outside the plasma confining region (the key external condition for achieving a good H-mode discharge) from very high density, cold radiative plasma near the wall (required for effective edge radiative cooling). Examples of such proposed SHC boundaries for Heliotron typed devices and tokamaks are presented. (author)

Boundary conditions for free surface inlet and outlet problems

KAUST Repository

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

2012-01-01

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

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

National Research Council Canada - National Science Library

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

2007-01-01

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

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)

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

Czech Academy of Sciences Publication Activity Database

Kružík, Martin; Zimmer, J.

2012-01-01

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

Effects of irradiation and thermal aging upon fatigue-crack growth behavior of reactor pressure boundary materials. [Neutrons

Energy Technology Data Exchange (ETDEWEB)

James, L. A.

1978-10-01

Two processes that have the potential to produce degradation in the properties of pressure boundary materials are neutron irradiation and long-time thermal aging. This paper uses linear-elastic fracture mechanics techniques to assess the effect of these two processes upon the fatigue-crack growth behavior of a number of alloys commonly employed in reactor pressure boundaries. The materials evaluated include ferritic steels, austenitic stainless steels, and nickel-base alloys typical of those employed in a number of reactor types including water-cooled, gas-cooled, and liquid-metal-cooled designs.

Effects of inlet boundary conditions, on the computed flow in the Turbine-99 draft tube, using OpenFOAM and CFX

Science.gov (United States)

Nilsson, H.; Cervantes, M. J.

2012-11-01

The flow in the Turbine-99 Kaplan draft tube was thoroughly investigated at three workshops (1999, 2001, 2005), which aimed at determining the state of the art of draft tube simulations. The flow is challenging due to the different flow phenomena appearing simultaneously such as unsteadiness, separation, swirl, turbulence, and a strong adverse pressure gradient. The geometry and the experimentally determined inlet boundary conditions were provided to the Turbine-99 workshop participants. At the final workshop, angular resolved inlet velocity boundary conditions were provided. The rotating non-axi-symmetry of the inlet flow due to the runner blades was thus included. The effect of the rotating angular resolution was however not fully investigated at that workshop. The first purpose of this work is to further investigate this effect. Several different inlet boundary conditions are applied - the angular resolved experimental data distributed at the Turbine-99 workshop, the angular resolved results of a runner simulation with interpolated values using different resolution in the tangential and radial directions, and an axi-symmetric variant of the same numerical data. The second purpose of this work is to compare the results from the OpenFOAM and CFX CFD codes, using as similar settings as possible. The present results suggest that the experimental angular inlet boundary conditions proposed to the workshop are not adequate to simulate accurately the flow in the T-99 draft tube. The reason for this is that the experimental phase-averaged data has some important differences compared to the previously measured time-averaged data. Using the interpolated data from the runner simulation as inlet boundary condition however gives good results as long as the resolution of that data is sufficient. It is shown that the difference between the results using the angular-resolved and the corresponding symmetric inlet data is very small, suggesting that the importance of the angular

Effects of inlet boundary conditions, on the computed flow in the Turbine-99 draft tube, using OpenFOAM and CFX

International Nuclear Information System (INIS)

Nilsson, H; Cervantes, M J

2012-01-01

The flow in the Turbine-99 Kaplan draft tube was thoroughly investigated at three workshops (1999, 2001, 2005), which aimed at determining the state of the art of draft tube simulations. The flow is challenging due to the different flow phenomena appearing simultaneously such as unsteadiness, separation, swirl, turbulence, and a strong adverse pressure gradient. The geometry and the experimentally determined inlet boundary conditions were provided to the Turbine-99 workshop participants. At the final workshop, angular resolved inlet velocity boundary conditions were provided. The rotating non-axi-symmetry of the inlet flow due to the runner blades was thus included. The effect of the rotating angular resolution was however not fully investigated at that workshop. The first purpose of this work is to further investigate this effect. Several different inlet boundary conditions are applied – the angular resolved experimental data distributed at the Turbine-99 workshop, the angular resolved results of a runner simulation with interpolated values using different resolution in the tangential and radial directions, and an axi-symmetric variant of the same numerical data. The second purpose of this work is to compare the results from the OpenFOAM and CFX CFD codes, using as similar settings as possible. The present results suggest that the experimental angular inlet boundary conditions proposed to the workshop are not adequate to simulate accurately the flow in the T-99 draft tube. The reason for this is that the experimental phase-averaged data has some important differences compared to the previously measured time-averaged data. Using the interpolated data from the runner simulation as inlet boundary condition however gives good results as long as the resolution of that data is sufficient. It is shown that the difference between the results using the angular-resolved and the corresponding symmetric inlet data is very small, suggesting that the importance of the angular

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.

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.

Boundary conditions estimation of scoop dosimeter for primary sorting during earthworks

International Nuclear Information System (INIS)

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

2005-01-01

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

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

DEFF Research Database (Denmark)

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

2008-01-01

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

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

Science.gov (United States)

Zhang, Shuhai; Oskay, Caglar

2015-04-01

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

A new technique for observationally derived boundary conditions for space weather

Science.gov (United States)

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

2018-04-01

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

Reduced Need of Lubricity Additives in Soybean Oil Blends Under Boundary Lubrication Conditions

Science.gov (United States)

Converging prices of vegetable oils and petroleum, along with increased focus on renewable resources, gave more momentum to vegetable oil lubricants. Boundary lubrication properties of four Extreme Pressure (EP) additive blends in conventional Soy Bean Oil (SBO) and Paraffinic Mineral Oil (PMO) of ...

Evolution of passive movement in advective environments: General boundary condition

Science.gov (United States)

Zhou, Peng; Zhao, Xiao-Qiang

2018-03-01

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

Two-dimensional properties of n-inversion layers in InSb grain boundaries under high hydrostatic pressure

International Nuclear Information System (INIS)

Kraak, W.; Herrmann, R.; Nachtwei, G.

1985-01-01

Magnetotransport properties of n-inversion layers in grain boundaries of p-InSb bicrystals are investigated under high hydrostatic pressure up to 10 3 MPa. A rapid decrease of the carrier concentration in the inversion layer is observed when hydrostatic pressure is applied. A simple model taking into account the pressure dependence of the energy band structure of pure InSb is proposed to describe this behaviour. (author)

The interior regularity of pressure associated with a weak solution to the Navier-Stokes equations with the Navier-type boundary conditions

Czech Academy of Sciences Publication Activity Database

Neustupa, Jiří; Al Baba, Hind

2018-01-01

Roč. 463, č. 1 (2018), s. 222-234 ISSN 0022-247X R&D Projects: GA ČR(CZ) GA17-01747S Institutional support: RVO:67985840 Keywords : Navier-Stokes equation * Navier-type boundary conditions * interior regularity Subject RIV: BA - General Mathematics OBOR OECD: Pure mathematics Impact factor: 1.064, year: 2016 https://www.sciencedirect.com/science/article/pii/S0022247X18302233?via%3Dihub

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

Directory of Open Access Journals (Sweden)

Weiye Huo

2015-02-01

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

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

Science.gov (United States)

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

2018-06-01

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

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.

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.

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

Science.gov (United States)

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

1979-01-01

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

Models for Predicting Boundary Conditions in L-Mode Tokamak Plasma

International Nuclear Information System (INIS)

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

2009-07-01

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

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

Science.gov (United States)

Mantile, Andrea; Posilicano, Andrea; Sini, Mourad

2016-07-01

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

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

Science.gov (United States)

Qian, Wei; Werner, Wendelin

2018-06-01

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

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

Science.gov (United States)

Verschaeve, Joris C G

2011-06-13

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

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

Energy Technology Data Exchange (ETDEWEB)

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

2004-02-06

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

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.

The interior regularity of pressure associated with a weak solution to the Navier-Stokes equations with the Navier-type boundary conditions

Czech Academy of Sciences Publication Activity Database

Neustupa, Jiří; Al Baba, Hind

2018-01-01

Roč. 463, č. 1 (2018), s. 222-234 ISSN 0022-247X R&D Projects: GA ČR(CZ) GA17-01747S Institutional support: RVO:67985840 Keywords : Navier-Stokes equation * Navier-type boundary conditions * interior regularity Subject RIV: BA - General Mathematics OBOR OECD: Pure mathematics Impact factor: 1.064, year: 2016 https://www. science direct.com/ science /article/pii/S0022247X18302233?via%3Dihub

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

OpenAIRE

Dipierro, Serena; Karakhanyan, Aram; Valdinoci, Enrico

2017-01-01

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

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

Energy Technology Data Exchange (ETDEWEB)

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

2012-08-23

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

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

International Nuclear Information System (INIS)

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

2012-01-01

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

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.

Confinement dynamics and boundary condition studies in the Reversed Field Pinch

International Nuclear Information System (INIS)

Schoenberg, K.F.; Ingraham, J.C.; Moses, R.W. Jr.

1988-01-01

The study of confinement dynamics, including investigation of the boundary conditions required for plasma sustainment, are central to the development of the Reversed Field Pinch (RFP) concept. Recently, several insights into confinement have emerged from a detailed investigation RFP electron and ion dynamics. These insights derive from the recognition that both magnetohydrodynamic (MHD) and electron kinetic effects play an important and coupled role in RFP stability, sustainment, and confinement. In this paper, we summarize the results of confinement studies on the ZT-40M experiment, and boundary condition studies on the Wisconsin non-circular RFP experiment. A brief description of the newly commissioned Madison Symmetric Torus (MST) is also presented. 28 refs., 3 figs

On flows of viscoelastic fluids under threshold-slip boundary conditions

Science.gov (United States)

Baranovskii, E. S.

2018-03-01

We investigate a boundary-value problem for the steady isothermal flow of an incompressible viscoelastic fluid of Oldroyd type in a 3D bounded domain with impermeable walls. We use the Fujita threshold-slip boundary condition. This condition states that the fluid can slip along a solid surface when the shear stresses reach a certain critical value; otherwise the slipping velocity is zero. Assuming that the flow domain is not rotationally symmetric, we prove an existence theorem for the corresponding slip problem in the framework of weak solutions. The proof uses methods for solving variational inequalities with pseudo-monotone operators and convex functionals, the method of introduction of auxiliary viscosity, as well as a passage-to-limit procedure based on energy estimates of approximate solutions, Korn’s inequality, and compactness arguments. Also, some properties and estimates of weak solutions are established.

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.

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

International Nuclear Information System (INIS)

Krieg, R.; Hailfinger, G.

1979-01-01

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

Analytical prediction on the pump-induced pulsating pressure in a reactor coolant pipe

International Nuclear Information System (INIS)

Lee, K.B.; Im, I.Y.; Lee, S.K.

1992-01-01

An analytical method is presented for predicting the amplitudes of pump-induced fluctuating pressures in a reactor coolant pipe using a linear transformation technique which reduces a homogeneous differential equation with non-homogeneous boundary conditions into a nonhomogeneous differential equation with homogeneous boundary conditions. At the end of the pipe, three types of boundary conditions are considered-open, closed and piston-spring supported. Numerical examples are given for a typical reactor. Comparisons of measured pressure amplitudes in the pipe with model prediction are shown to be in good agreement for the forcing frequencies. (author)

On a stochastic Burgers equation with Dirichlet boundary conditions

Directory of Open Access Journals (Sweden)

Ekaterina T. Kolkovska

2003-01-01

Full Text Available We consider the one-dimensional Burgers equation perturbed by a white noise term with Dirichlet boundary conditions and a non-Lipschitz coefficient. We obtain existence of a weak solution proving tightness for a sequence of polygonal approximations for the equation and solving a martingale problem for the weak limit.

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.

Influence of Boundary Conditions on the Simulation of a Diamond-Type Lattice Structure: A Preliminary Study

Directory of Open Access Journals (Sweden)

Patrick Terriault

2017-01-01

Full Text Available Emergent additive manufacturing processes allow the use of metallic porous structures in various industrial applications. Because these structures comprise a large number of ordered unit cells, their design using conventional modeling approaches, such as finite elements, becomes a real challenge. A homogenization technique, in which the lattice structure is simulated as a fully dense volume having equivalent material properties, can then be employed. To determine these equivalent material properties, numerical simulations can be performed on a single unit cell of the lattice structure. However, a critical aspect to consider is the boundary conditions applied to the external faces of the unit cell. In the literature, different types of boundary conditions are used, but a comparative study is definitely lacking. In this publication, a diamond-type unit cell is studied in compression by applying different boundary conditions. If the porous structure’s boundaries are free to deform, then the periodic boundary condition is found to be the most representative, but constraint equations must be introduced in the model. If, instead, the porous structure is inserted in a rigid enclosure, it is then better to use frictionless boundary conditions. These preliminary results remain to be validated for other types of unit cells loaded beyond the yield limit of the material.

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.

Instability waves and transition in adverse-pressure-gradient boundary layers

Science.gov (United States)

Bose, Rikhi; Zaki, Tamer A.; Durbin, Paul A.

2018-05-01

Transition to turbulence in incompressible adverse-pressure-gradient (APG) boundary layers is investigated by direct numerical simulations. Purely two-dimensional instability waves develop on the inflectional base velocity profile. When the boundary layer is perturbed by isotropic turbulence from the free stream, streamwise elongated streaks form and may interact with the instability waves. Subsequent mechanisms that trigger transition depend on the intensity of the free-stream disturbances. All evidence from the present simulations suggest that the growth rate of instability waves is sufficiently high to couple with the streaks. Under very low levels of free-stream turbulence (˜0.1 % ), transition onset is highly sensitive to the inlet disturbance spectrum and is accelerated if the spectrum contains frequency-wave-number combinations that are commensurate with the instability waves. Transition onset and completion in this regime is characterized by formation and breakdown of Λ vortices, but they are more sporadic than in natural transition. Beneath free-stream turbulence with higher intensity (1-2 % ), bypass transition mechanisms are dominant, but instability waves are still the most dominant disturbances in wall-normal and spanwise perturbation spectra. Most of the breakdowns were by disturbances with critical layers close to the wall, corresponding to inner modes. On the other hand, the propensity of an outer mode to occur increases with the free-stream turbulence level. Higher intensity free-stream disturbances induce strong streaks that favorably distort the boundary layer and suppress the growth of instability waves. But the upward displacement of high amplitude streaks brings them to the outer edge of the boundary layer and exposes them to ambient turbulence. Consequently, high-amplitude streaks exhibit an outer-mode secondary instability.

Generalized wall function and its application to compressible turbulent boundary layer over a flat plate

Science.gov (United States)

Liu, J.; Wu, S. P.

2017-04-01

Wall function boundary conditions including the effects of compressibility and heat transfer are improved for compressible turbulent boundary flows. Generalized wall function formulation at zero-pressure gradient is proposed based on coupled velocity and temperature profiles in the entire near-wall region. The parameters in the generalized wall function are well revised. The proposed boundary conditions are integrated into Navier-Stokes computational fluid dynamics code that includes the shear stress transport turbulence model. Numerical results are presented for a compressible boundary layer over a flat plate at zero-pressure gradient. Compared with experimental data, the computational results show that the generalized wall function reduces the first grid spacing in the directed normal to the wall and proves the feasibility and effectivity of the generalized wall function method.

Robust and general method for determining surface fluid flow boundary conditions in articular cartilage contact mechanics modeling.

Science.gov (United States)

Pawaskar, Sainath Shrikant; Fisher, John; Jin, Zhongmin

2010-03-01

Contact detection in cartilage contact mechanics is an important feature of any analytical or computational modeling investigation when the biphasic nature of cartilage and the corresponding tribology are taken into account. The fluid flow boundary conditions will change based on whether the surface is in contact or not, which will affect the interstitial fluid pressurization. This in turn will increase or decrease the load sustained by the fluid phase, with a direct effect on friction, wear, and lubrication. In laboratory experiments or clinical hemiarthroplasty, when a rigid indenter or metallic prosthesis is used to apply load to the cartilage, there will not be any fluid flow normal to the surface in the contact region due to the impermeable nature of the indenter/prosthesis. In the natural joint, on the other hand, where two cartilage surfaces interact, flow will depend on the pressure difference across the interface. Furthermore, in both these cases, the fluid would flow freely in non-contacting regions. However, it should be pointed out that the contact area is generally unknown in advance in both cases and can only be determined as part of the solution. In the present finite element study, a general and robust algorithm was proposed to decide nodes in contact on the cartilage surface and, accordingly, impose the fluid flow boundary conditions. The algorithm was first tested for a rigid indenter against cartilage model. The algorithm worked well for two-dimensional four-noded and eight-noded axisymmetric element models as well as three-dimensional models. It was then extended to include two cartilages in contact. The results were in excellent agreement with the previous studies reported in the literature.

Predictability of Japan/East Sea (JES) System to Uncertain Initial/Lateral Boundary Conditions and Surface Winds

National Research Council Canada - National Science Library

Fang, Chin-Lung

2003-01-01

.... Change in either initial or boundary condition leads to a variety of model solutions. It is necessary to specify realistic initial and boundary conditions to achieve better understanding and prediction of the ocean behavior...

Effects of microscopic boundary conditions on plastic deformations of small-sized single crystals

DEFF Research Database (Denmark)

Kuroda, Mitsutoshi; Tvergaard, Viggo

2009-01-01

The finite deformation version of the higher-order gradient crystal plasticity model proposed by the authors is applied to solve plane strain boundary value problems, in order to obtain an understanding of the effect of the higher-order boundary conditions. Numerical solutions are carried out...

Boundary Conditions and SGS Models for LES of Wall-Bounded Separated Flows: An Application to Engine-Like Geometries

Directory of Open Access Journals (Sweden)

Piscaglia F.

2013-11-01

Full Text Available The implementation and the combination of advanced boundary conditions and subgrid scale models for Large Eddy Simulations are presented. The goal is to perform reliable cold flow LES simulations in complex geometries, such as in the cylinders of internal combustion engines. The implementation of an inlet boundary condition for synthetic turbulence generation and of two subgrid scale models, the local Dynamic Smagorinsky and the Wall-Adapting Local Eddy-viscosity SGS model ( WALE is described. The WALE model is based on the square of the velocity gradient tensor and it accounts for the effects of both the strain and the rotation rate of the smallest resolved turbulent fluctuations and it recovers the proper y3 near-wall scaling for the eddy viscosity without requiring dynamic pressure; hence, it is supposed to be a very reliable model for ICE simulation. Model validation has been performed separately on two steady state flow benches: a backward facing step geometry and a simple IC engine geometry with one axed central valve. A discussion on the completeness of the LES simulation (i.e. LES simulation quality is given.

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)

Mode solutions for a Klein-Gordon field in anti-de Sitter spacetime with dynamical boundary conditions of Wentzell type

Science.gov (United States)

Dappiaggi, Claudio; Ferreira, Hugo R. C.; Juárez-Aubry, Benito A.

2018-04-01

We study a real, massive Klein-Gordon field in the Poincaré fundamental domain of the (d +1 )-dimensional anti-de Sitter (AdS) spacetime, subject to a particular choice of dynamical boundary conditions of generalized Wentzell type, whereby the boundary data solves a nonhomogeneous, boundary Klein-Gordon equation, with the source term fixed by the normal derivative of the scalar field at the boundary. This naturally defines a field in the conformal boundary of the Poincaré fundamental domain of AdS. We completely solve the equations for the bulk and boundary fields and investigate the existence of bound state solutions, motivated by the analogous problem with Robin boundary conditions, which are recovered as a limiting case. Finally, we argue that both Robin and generalized Wentzell boundary conditions are distinguished in the sense that they are invariant under the action of the isometry group of the AdS conformal boundary, a condition which ensures in addition that the total flux of energy across the boundary vanishes.

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.

Axisymmetric Tornado Simulations with a Semi-Slip Boundary

Directory of Open Access Journals (Sweden)

Brian H. Fiedler

2017-12-01

Full Text Available The structure of natural tornadoes and simulated analogs are sensitive to the lower boundary condition for friction. Three-dimensional numerical simulations of storms require a choice for turbulence parameterizations and resolution of wind near the lower boundary. This article explores some of the consequences of choices of a surface drag coefficient on the structure of a mature simulated tornado, using a conventional axisymmetric model. The surface drag parameterization is explored over the range of the semi-slip condition, including the extremes of no-slip and free-slip. A moderate semi-slip condition allows for an extreme pressure deficit, but without the unrealistic vortex breakdown of the no-slip condition.

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

NARCIS (Netherlands)

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

2005-01-01

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

Influence of the Outer Boundary Condition on models of AGB stars

Science.gov (United States)

Wagstaff, G.; Weiss, A.

2018-04-01

Current implementations of the stellar atmosphere typically derive boundary conditions for the interior model from either grey plane-parallel atmospheres or scaled solar atmospheres, neither of which can be considered to have appropriate underlying assumptions for the Thermally Pulsing Asymptotic Giant Branch (TP-AGB). This paper discusses the treatment and influence of the outer boundary condition within stellar evolution codes, and the resulting effects on the AGB evolution. The complex interaction of processes, such as the third dredge up and mass loss, governing the TP-AGB can be affected by varying the treatment of this boundary condition. Presented here are the results from altering the geometry, opacities and the implementation of a grid of MARCS/COMARCS model atmospheres in order to improve this treatment. Although there are changes in the TP-AGB evolution, observable quantities, such as the final core mass, are not significantly altered as a result of the change of atmospheric treatment. During the course of the investigation, a previously unseen phenomena in the AGB models was observed and further investigated. This is believed to be physical, although arising from specific conditions which make its presence unlikely. If it were present in stars, this phenomenon would increase the carbon-star lifetime above 10Myr and increase the final core mass by ˜0.1M⊙ in the narrow initial-mass range where it was observed (˜2 - 2.3M⊙).

Enhancement of fatigue crack growth rates in pressure boundary materials due to light-water-reactor environments

International Nuclear Information System (INIS)

VanDerSluys, W.A.; Emanuelson, R.H.

1988-01-01

The high level of reliability required of the primary-coolant pressure boundary in a nuclear reactor system leads to a continuing interest in the interaction among the coolant, pressure boundary materials, and service loadings. One area of concern involves the possible enhancement of the growth rate of fatigue cracks due to the coolant. Advances have occurred recently toward a better understanding of the variables influencing the material/environment interactions and methods of addressing this interaction. Sulfur now appears to be one of the principal agents responsible for the observed enhancement of the fatigue crack growth rates in light-water-reactor (LWR) environments. This paper presents the results of investigations on the effect of sulfur in the steel, bulk water environment, and at the crack tip

Experimental verification of electrostatic boundary conditions in gate-patterned quantum devices

Science.gov (United States)

Hou, H.; Chung, Y.; Rughoobur, G.; Hsiao, T. K.; Nasir, A.; Flewitt, A. J.; Griffiths, J. P.; Farrer, I.; Ritchie, D. A.; Ford, C. J. B.

2018-06-01

In a model of a gate-patterned quantum device, it is important to choose the correct electrostatic boundary conditions (BCs) in order to match experiment. In this study, we model gated-patterned devices in doped and undoped GaAs heterostructures for a variety of BCs. The best match is obtained for an unconstrained surface between the gates, with a dielectric region above it and a frozen layer of surface charge, together with a very deep back boundary. Experimentally, we find a  ∼0.2 V offset in pinch-off characteristics of 1D channels in a doped heterostructure before and after etching off a ZnO overlayer, as predicted by the model. Also, we observe a clear quantised current driven by a surface acoustic wave through a lateral induced n-i-n junction in an undoped heterostructure. In the model, the ability to pump electrons in this type of device is highly sensitive to the back BC. Using the improved boundary conditions, it is straightforward to model quantum devices quite accurately using standard software.

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

Environment sensitive cracking in light water reactor pressure boundary materials

International Nuclear Information System (INIS)

Haenninen, H.; Aho-Mantila, I.

1985-01-01

The purpose of the paper is to review the available methods and the most promising future possibilities of preventive maintenance to counteract the various forms of environment sensitive cracking of pressure boundary materials in light water reactors. Environment sensitive cracking is considered from the metallurgical, mechanical and environmental point of view. The main emphasis is on intergranular stress corrosion cracking of austenitic stainless steels and high strenght Ni-base alloys, as well as on corrosion fatigue of low alloy and stainless steels. Finally, some general ideas how to predict, reduce or eliminate environment sensitive cracking in service are presented

Simulation of microcirculatory hemodynamics: estimation of boundary condition using particle swarm optimization.

Science.gov (United States)

Pan, Qing; Wang, Ruofan; Reglin, Bettina; Fang, Luping; Pries, Axel R; Ning, Gangmin

2014-01-01

Estimation of the boundary condition is a critical problem in simulating hemodynamics in microvascular networks. This paper proposed a boundary estimation strategy based on a particle swarm optimization (PSO) algorithm, which aims to minimize the number of vessels with inverted flow direction in comparison to the experimental observation. The algorithm took boundary values as the particle swarm and updated the position of the particles iteratively to approach the optimization target. The method was tested in a real rat mesenteric network. With random initial boundary values, the method achieved a minimized 9 segments with an inverted flow direction in the network with 546 vessels. Compared with reported literature, the current work has the advantage of a better fit with experimental observations and is more suitable for the boundary estimation problem in pulsatile hemodynamic models due to the experiment-based optimization target selection.

Active Brownian particles near straight or curved walls: Pressure and boundary layers

Science.gov (United States)

Duzgun, Ayhan; Selinger, Jonathan V.

2018-03-01

Unlike equilibrium systems, active matter is not governed by the conventional laws of thermodynamics. Through a series of analytic calculations and Langevin dynamics simulations, we explore how systems cross over from equilibrium to active behavior as the activity is increased. In particular, we calculate the profiles of density and orientational order near straight or circular walls and show the characteristic width of the boundary layers. We find a simple relationship between the enhancements of density and pressure near a wall. Based on these results, we determine how the pressure depends on wall curvature and hence make approximate analytic predictions for the motion of curved tracers, as well as the rectification of active particles around small openings in confined geometries.

Universal parity effects in the entanglement entropy of XX chains with open boundary conditions

International Nuclear Information System (INIS)

Fagotti, Maurizio; Calabrese, Pasquale

2011-01-01

We consider the Rényi entanglement entropies in the one-dimensional XX spin-chains with open boundary conditions in the presence of a magnetic field. In the case of a semi-infinite system and a block starting from the boundary, we derive rigorously the asymptotic behavior for large block sizes on the basis of a recent mathematical theorem for the determinant of Toeplitz plus Hankel matrices. We conjecture a generalized Fisher–Hartwig form for the corrections to the asymptotic behavior of this determinant that allows the exact characterization of the corrections to the scaling at order o(l -1 ) for any n. By combining these results with conformal field theory arguments, we derive exact expressions also in finite chains with open boundary conditions and in the case when the block is detached from the boundary

Validation of Boundary Conditions for CFD Simulations on Ventilated Rooms

DEFF Research Database (Denmark)

Topp, Claus; Jensen, Rasmus Lund; Pedersen, D.N.

2001-01-01

The application of Computational Fluid Dynamics (CFD) for ventilation research and design of ventilation systems has increased during the recent years. This paper provides an investigation of direct description of boundary conditions for a complex inlet diffuser and a heated surface. A series...

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.

Bound states in waveguides with complex Robin boundary conditions

Czech Academy of Sciences Publication Activity Database

Novák, Radek

2016-01-01

Roč. 96, 3-4 (2016), s. 251-281 ISSN 0921-7134 R&D Projects: GA ČR(CZ) GA14-06818S Institutional support: RVO:61389005 Keywords : non-self-adjointness * waveguide * Robin boundary conditions * spectral analysis * essential spectrum * weak coupling * Birman-Schwinger principle * reality of the spectrum Subject RIV: BE - Theoretical Physics Impact factor: 0.933, year: 2016

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

NARCIS (Netherlands)

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

2005-01-01

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

Explicit formulation for natural frequencies of double-beam system with arbitrary boundary conditions

Energy Technology Data Exchange (ETDEWEB)

Mirzabeigy, Alborz; Madoliat, Reza [Iran University of Science and Technology, Narmak, Tehran (Iran, Islamic Republic of); Dabbagh, Vahid [University of Malaya, Kuala Lumpur (Malaysia)

2017-02-15

In this paper, free transverse vibration of two parallel beams connected through Winkler type elastic layer is investigated. Euler- Bernoulli beam hypothesis has been applied and it is assumed that boundary conditions of upper and lower beams are similar while arbitrary without any limitation even for non-ideal boundary conditions. Material properties and cross-section geometry of beams could be different from each other. The motion of the system is described by a homogeneous set of two partial differential equations, which is solved by using the classical Bernoulli-Fourier method. Explicit expressions are derived for the natural frequencies. In order to verify accuracy of results, the problem once again solved using modified Adomian decomposition method. Comparison between results indicates excellent accuracy of proposed formulation for any arbitrary boundary conditions. Derived explicit formulation is simplest method to determine natural frequencies of double-beam systems with high level of accuracy in comparison with other methods in literature.

Adapting a Fourier pseudospectral method to Dirichlet boundary conditions for Rayleigh-Bénard convection

Directory of Open Access Journals (Sweden)

I. C. Ramos

2015-10-01

Full Text Available We present the adaptation to non-free boundary conditions of a pseudospectral method based on the (complex Fourier transform. The method is applied to the numerical integration of the Oberbeck-Boussinesq equations in a Rayleigh-Bénard cell with no-slip boundary conditions for velocity and Dirichlet boundary conditions for temperature. We show the first results of a 2D numerical simulation of dry air convection at high Rayleigh number (. These results are the basis for the later study, by the same method, of wet convection in a solar still. Received: 20 Novembre 2014, Accepted: 15 September 2015; Edited by: C. A. Condat, G. J. Sibona; DOI:http://dx.doi.org/10.4279/PIP.070015 Cite as: I C Ramos, C B Briozzo, Papers in Physics 7, 070015 (2015

RECTC/RECTCF, 2. Order Elliptical Partial Differential Equation, Arbitrary Boundary Conditions

International Nuclear Information System (INIS)

Hackbusch, W.

1983-01-01

1 - Description of problem or function: A general linear elliptical second order partial differential equation on a rectangle with arbitrary boundary conditions is solved. 2 - Method of solution: Multi-grid iteration

The structure of a three-dimensional boundary layer subjected to streamwise-varying spanwise-homogeneous pressure gradient

International Nuclear Information System (INIS)

Bentaleb, Y.; Leschziner, M.A.

2013-01-01

Highlights: • We study a spatially-evolving three-dimensional boundary layer. • We impose a streamwise-varying spanwise-homogeneous pressure gradient. • A collateral flow is formed close to the wall, and this is investigated alongside the skewed upper part of the boundary layer. • A wide range of flow-physical properties have been studied. -- Abstract: A spatially-evolving three-dimensional boundary layer, subjected to a streamwise-varying spanwise-homogeneous pressure gradient, equivalent to a body force, is investigated by way of direct numerical simulation. The pressure gradient, prescribed to change its sign half-way along the boundary layer, provokes strong skewing of the velocity vector, with a layer of nearly collateral flow forming close to the wall up to the position of maximum spanwise velocity. A wide range of flow-physical properties have been studied, with particular emphasis on the near-wall layer, including second-moments, major budget contributions and wall-normal two-point correlations of velocity fluctuations and their angles, relative to wall-shear fluctuations. The results illustrate the complexity caused by skewing, including a damping in turbulent mixing and a significant lag between strains and stresses. The study has been undertaken in the context of efforts to develop and test novel hybrid LES–RANS schemes for non-equilibrium near-wall flows, with an emphasis on three-dimensional near-wall straining. Fundamental flow-physical issues aside, the data derived should be of particular relevance to a priori studies of second-moment RANS closure and the development and validation of RANS-type near-wall approximations implemented in LES schemes for high-Reynolds-number complex flows

Effect of magnetization boundary condition on cavity magnon polariton of YIG thin film

Science.gov (United States)

Jiang, H. H.; Xiao, Y.; Hu, C. M.; Guo, H.; Xia, K.

2018-06-01

Motivated by recent studies of cavity magnon polariton (CMP), we extended a previous theoretical work to generalize microwave transmission calculation with various magnetization boundary condition of YIG thin film embedded in cavity. It is found that numerical implementation given in this paper can be easily applied to other magnetization boundary condition and extended to magnetic multilayers. Numerical results show that ferromagnetic resonance mode of microwave transmission spectrum, which is absent in previous calculation, can be recovered by altering the pinning condition of surface spins. The demonstrated reliability of our theory opens attractive perspectives for studying CMP of thin film with complicated surface magnetization distribution and magnetic multilayers.

Global exponential stability and periodicity of reaction-diffusion delayed recurrent neural networks with Dirichlet boundary conditions

International Nuclear Information System (INIS)

Lu Junguo

2008-01-01

In this paper, the global exponential stability and periodicity for a class of reaction-diffusion delayed recurrent neural networks with Dirichlet boundary conditions are addressed by constructing suitable Lyapunov functionals and utilizing some inequality techniques. We first prove global exponential converge to 0 of the difference between any two solutions of the original reaction-diffusion delayed recurrent neural networks with Dirichlet boundary conditions, the existence and uniqueness of equilibrium is the direct results of this procedure. This approach is different from the usually used one where the existence, uniqueness of equilibrium and stability are proved in two separate steps. Furthermore, we prove periodicity of the reaction-diffusion delayed recurrent neural networks with Dirichlet boundary conditions. Sufficient conditions ensuring the global exponential stability and the existence of periodic oscillatory solutions for the reaction-diffusion delayed recurrent neural networks with Dirichlet boundary conditions are given. These conditions are easy to check and have important leading significance in the design and application of reaction-diffusion recurrent neural networks with delays. Finally, two numerical examples are given to show the effectiveness of the obtained results

On the physical solutions to the heat equation subjected to nonlinear boundary conditions

International Nuclear Information System (INIS)

Gama, R.M.S. da.

1990-01-01

This work consists of a discussion on the physical solutions to the steady-state heat transfer equation, when it is subjected to nonlinear boundary conditions. It will be presented a functional, whose minimum occurs for the (unique) physical solution to the condidered heat transfer problem, suitable for a large class of typical (nonlinear) boundary conditions (representing the radiative/convective loss from the body to the environment). It will be demonstrated that these problems admit-always one, and only one, physical solution (which represents the absolute temperature). (author)

Heat conduction in a plate-type fuel element with time-dependent boundary conditions

International Nuclear Information System (INIS)

Faya, A.J.G.; Maiorino, J.R.

1981-01-01

A method for the solution of boundary-value problems with variable boundary conditions is applied to solve a heat conduction problem in a plate-type fuel element with time dependent film coefficient. The numerical results show the feasibility of the method in the solution of this class of problems. (Author) [pt

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

Directory of Open Access Journals (Sweden)

Yanmei Sun

2012-01-01

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

Free-boundary Full-pressure Island Healing in a Stellarator: Coil-healing

International Nuclear Information System (INIS)

Hudson, S.R.; Reiman, A.; Strickler, D.; Brooks, A.; Monticello, D.A.; Hirshman, S.P.

2002-01-01

The lack of axisymmetry in stellarators guarantees that in general magnetic islands and chaotic magnetic field lines will exist. As particle transport is strongly tied to the magnetic field lines, magnetic islands and chaotic field lines result in poor plasma confinement. For stellarators to be feasible candidates for fusion power stations it is essential that, to a good approximation, the magnetic field lines lie on nested flux-surfaces, and the suppression of magnetic islands is a critical issue for stellarator coil design, particularly for small aspect ratio devices. A procedure for modifying stellarator coil designs to eliminate magnetic islands in free-boundary full-pressure magnetohydrodynamic equilibria is presented. Islands may be removed from coil-plasma free-boundary equilibria by making small changes to the coil geometry and also by variation of trim coil currents. A plasma and coil design relevant to the National Compact Stellarator Experiment is used to illustrate the technique

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)

Particles in a magnetic field and plasma analogies: doubly periodic boundary conditions

International Nuclear Information System (INIS)

Forrester, P J

2006-01-01

The N-particle free fermion state for quantum particles in the plane subject to a perpendicular magnetic field, and with doubly periodic boundary conditions, is written in a product form. The absolute value of this is used to formulate an exactly solvable one-component plasma model and further motivates the formulation of an exactly solvable two-species Coulomb gas. The large N expansion of the free energy of both these models exhibits the same O(1) term. On the basis of a relationship to the Gaussian free field, this term is predicted to be universal for conductive Coulomb systems in doubly periodic boundary conditions

Identification of complex model thermal boundary conditions based on exterior temperature measurement

International Nuclear Information System (INIS)

Lu Jianming; Ouyang Guangyao; Zhang Ping; Rong Bojun

2012-01-01

Combining the advantages of the finite element software in temperature field analyzing with the multivariate function optimization arithmetic, a feasibility method based on the exterior temperature was proposed to get the thermal boundary conditions, which was required in temperature field analyzing. The thermal boundary conditions can be obtained only by some temperature measurement values. Taking the identification of the convection heat transfer coefficient of a high power density diesel engine cylinder head as an example, the calculation result shows that when the temperature measurement error was less than 0.5℃, the maximum relative error was less than 2%. It is shown that the new method was feasible (authors)

Effect of disjoining pressure in a thin film equation with non-uniform forcing

KAUST Repository

MOULTON, D. E.

2013-08-02

We explore the effect of disjoining pressure on a thin film equation in the presence of a non-uniform body force, motivated by a model describing the reverse draining of a magnetic film. To this end, we use a combination of numerical investigations and analytical considerations. The disjoining pressure has a regularizing influence on the evolution of the system and appears to select a single steady-state solution for fixed height boundary conditions; this is in contrast with the existence of a continuum of locally attracting solutions that exist in the absence of disjoining pressure for the same boundary conditions. We numerically implement matched asymptotic expansions to construct equilibrium solutions and also investigate how they behave as the disjoining pressure is sent to zero. Finally, we consider the effect of the competition between forcing and disjoining pressure on the coarsening dynamics of the thin film for fixed contact angle boundary conditions. Copyright © Cambridge University Press 2013.

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

Grain boundary microchemistry and metallurgical characterization of Eurofer'97 after simulated service conditions

International Nuclear Information System (INIS)

Fernandez, P.; Garcia-Mazario, M.; Lancha, A.M.; Lapena, J.

2004-01-01

The aim of this paper is to describe the microstructural investigations, the mechanical properties (hardness, tensile and charpy) and the grain boundary microchemistry studied by Auger electron spectroscopy (AES), of the Eurofer'97 steel aged in the range of temperatures from 400 to 600 deg. C up to 10,000 h. After these thermal aging treatments the steel showed a high microstructural stability, and similar values of hardness, ultimate tensile strength and 0.2% proof stress regardless of the material condition. A slight DBTT increase was observed in the material aged at 600 deg. C for 5000 and 10,000 h. The Auger results showed chromium enrichment at grain boundaries in all material conditions. In addition, phosphorus was detected at the grain boundaries after the aging treatments at 500 deg. C

Liouville equation with boundary conditions derived from classical strings

International Nuclear Information System (INIS)

Marnelius, R.

1983-01-01

It is shown in terms of the classical string theory that a breaking of the Weyl invariance necessarily requires the Liouville equation for the variable phi=1n rho, where rho is the variable that appears in the conformal gauge gsub(α#betta#)=rhoetasub(α#betta#). Appropriate boundary conditions on phi for open and closed strings are then derived. (orig.)

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

Energy Technology Data Exchange (ETDEWEB)

Rian, Kjell Erik

2003-07-01

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

Analysis of Numerical Simulation Database for Pressure Fluctuations Induced by High-Speed Turbulent Boundary Layers

Science.gov (United States)

Duan, Lian; Choudhari, Meelan M.

2014-01-01

Direct numerical simulations (DNS) of Mach 6 turbulent boundary layer with nominal freestream Mach number of 6 and Reynolds number of Re(sub T) approximately 460 are conducted at two wall temperatures (Tw/Tr = 0.25, 0.76) to investigate the generated pressure fluctuations and their dependence on wall temperature. Simulations indicate that the influence of wall temperature on pressure fluctuations is largely limited to the near-wall region, with the characteristics of wall-pressure fluctuations showing a strong temperature dependence. Wall temperature has little influence on the propagation speed of the freestream pressure signal. The freestream radiation intensity compares well between wall-temperature cases when normalized by the local wall shear; the propagation speed of the freestream pressure signal and the orientation of the radiation wave front show little dependence on the wall temperature.

Study of Boundary Layer Convective Heat Transfer with Low Pressure Gradient Over a Flat Plate Via He's Homotopy Perturbation Method

International Nuclear Information System (INIS)

Fathizadeh, M.; Aroujalian, A.

2012-01-01

The boundary layer convective heat transfer equations with low pressure gradient over a flat plate are solved using Homotopy Perturbation Method, which is one of the semi-exact methods. The nonlinear equations of momentum and energy solved simultaneously via Homotopy Perturbation Method are in good agreement with results obtained from numerical methods. Using this method, a general equation in terms of Pr number and pressure gradient (λ) is derived which can be used to investigate velocity and temperature profiles in the boundary layer.

Reynolds stress structures in a self-similar adverse pressure gradient turbulent boundary layer at the verge of separation.

Science.gov (United States)

Atkinson, C.; Sekimoto, A.; Jiménez, J.; Soria, J.

2018-04-01

Mean Reynolds stress profiles and instantaneous Reynolds stress structures are investigated in a self-similar adverse pressure gradient turbulent boundary layer (APG-TBL) at the verge of separation using data from direct numerical simulations. The use of a self-similar APG-TBL provides a flow domain in which the flow gradually approaches a constant non-dimensional pressure gradient, resulting in a flow in which the relative contribution of each term in the governing equations is independent of streamwise position over a domain larger than two boundary layer thickness. This allows the flow structures to undergo a development that is less dependent on the upstream flow history when compared to more rapidly decelerated boundary layers. This APG-TBL maintains an almost constant shape factor of H = 2.3 to 2.35 over a momentum thickness based Reynolds number range of Re δ 2 = 8420 to 12400. In the APG-TBL the production of turbulent kinetic energy is still mostly due to the correlation of streamwise and wall-normal fluctuations, 〈uv〉, however the contribution form the other components of the Reynolds stress tensor are no longer negligible. Statistical properties associated with the scale and location of sweeps and ejections in this APG-TBL are compared with those of a zero pressure gradient turbulent boundary layer developing from the same inlet profile, resulting in momentum thickness based range of Re δ 2 = 3400 to 3770. In the APG-TBL the peak in both the mean Reynolds stress and the production of turbulent kinetic energy move from the near wall region out to a point consistent with the displacement thickness height. This is associated with a narrower distribution of the Reynolds stress and a 1.6 times higher relative number of wall-detached negative uv structures. These structures occupy 5 times less of the boundary layer volume and show a similar reduction in their streamwise extent with respect to the boundary layer thickness. A significantly lower percentage

Pressure cylinders under fire condition

Directory of Open Access Journals (Sweden)

Jan Hora

2016-03-01

Full Text Available The presence of pressure cylinders under fire conditions significantly increases the risk rate for the intervening persons. It is considerably problematic to predict the pressure cylinders behaviour during heat exposition, its destruction progress and possible following explosion of the produced air–gas mixture because pressure cylinders and its environment generate a highly complicated dynamic system during an uncontrolled destruction. The large scale tests carried out by the Pilsen Fire and Rescue Department and the Rapid Response Unit of the Czech Republic Police in October 2012 and in May 2014 in the Military area Brdy and in the area of the former Lachema factory in Kaznějov had several objectives, namely, to record, qualify and quantify some of the aspects of an uncontrolled heat destruction procedure of an exposed pressure cylinder in an enclosed space and to qualify and describe the process of a controlled destruction of a pressure cylinder by shooting through it including basic tactical concepts. The article describes the experiments that were carried out.

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

NARCIS (Netherlands)

JUFFER, AH; BERENDSEN, HJC

1993-01-01

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

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.

Explicitly represented polygon wall boundary model for the explicit MPS method

Science.gov (United States)

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

2015-05-01

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

Interaction between a normal shock wave and a turbulent boundary layer at high transonic speeds. I - Pressure distribution

Science.gov (United States)

Messiter, A. F.

1980-01-01

Asymptotic solutions are derived for the pressure distribution in the interaction of a weak normal shock wave with a turbulent boundary layer. The undisturbed boundary layer is characterized by the law of the wall and the law of the wake for compressible flow. In the limiting case considered, for 'high' transonic speeds, the sonic line is very close to the wall. Comparisons with experiment are shown, with corrections included for the effect of longitudinal wall curvature and for the boundary-layer displacement effect in a circular pipe.

A preliminary investigation of boundary-layer transition along a flat plate with adverse pressure gradient

Science.gov (United States)

Von Doenhoff, Albert E

1938-01-01

Boundary-layer surveys were made throughout the transition region along a smooth flat plate placed in an airstream of practically zero turbulence and with an adverse pressure gradient. The boundary-layer Reynolds number at the laminar separation point was varied from 1,800 to 2,600. The test data, when considered in the light of certain theoretical deductions, indicated that transition probably began with separation of the laminar boundary layer. The extent of the transition region, defined as the distance from a calculated laminar separation point to the position of the first fully developed turbulent boundary-layer profile, could be expressed as a constant Reynolds number run of approximately 70,000. Some speculations are presented concerning the application of the foregoing concepts, after certain assumptions have been made, to the problem of the connection between transition on the upper surface of an airfoil at high angles of attack and the maximum lift.

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.

Asymptotic stability and blow up for a semilinear damped wave equation with dynamic boundary conditions

KAUST Repository

Gerbi, Stéphane

2011-12-01

In this paper we consider a multi-dimensional wave equation with dynamic boundary conditions, related to the KelvinVoigt damping. Global existence and asymptotic stability of solutions starting in a stable set are proved. Blow up for solutions of the problem with linear dynamic boundary conditions with initial data in the unstable set is also obtained. © 2011 Elsevier Ltd. All rights reserved.

Asymptotic stability and blow up for a semilinear damped wave equation with dynamic boundary conditions

KAUST Repository

Gerbi, Sté phane; Said-Houari, Belkacem

2011-01-01

In this paper we consider a multi-dimensional wave equation with dynamic boundary conditions, related to the KelvinVoigt damping. Global existence and asymptotic stability of solutions starting in a stable set are proved. Blow up for solutions of the problem with linear dynamic boundary conditions with initial data in the unstable set is also obtained. © 2011 Elsevier Ltd. All rights reserved.

Nonlinear $q$-fractional differential equations with nonlocal and sub-strip type boundary conditions

Directory of Open Access Journals (Sweden)

Bashir Ahmad

2014-06-01

Full Text Available This paper is concerned with new boundary value problems of nonlinear $q$-fractional differential equations with nonlocal and sub-strip type boundary conditions. Our results are new in the present setting and rely on the contraction mapping principle and a fixed point theorem due to O'Regan. Some illustrative examples are also presented.

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.

The effect of spacer grid critical component on pressure drop under both single and two phase flow conditions

International Nuclear Information System (INIS)

Han, B.; Yang, B.W.; Zhang, H.; Mao, H.; Zha, Y.

2016-01-01

As pressure drop is one of the most critical thermal hydraulic parameters for spacer grids the accurate estimation of it is the key to the design and development of spacer grids. Most of the available correlations for pressure drop do not contain any real geometrical parameters that characterize the grid effect. The main functions for spacer grid are structural support and flow mixing. Once the boundary sublayer near the rod bundle is disturbed, the liquid forms swirls or flow separation that affect pressure drop. However, under two phase flow conditions, due to the existence of steam bubble, the complexity for spacer grid are multiplied and pressure drop calculation becomes much more challenging. The influence of the dimple location, distance of mixing vane to the nearest strip, and the effect of inter-subchannel mixing among neighboring subchannels on pressure drop and downstream flow fields are analyzed in this paper. Based on this study, more detailed space grid geometry parameters are recommended for adding into the correlation when predicting pressure drop.

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.

Influence of Contact Angle Boundary Condition on CFD Simulation of T-Junction

Science.gov (United States)

Arias, S.; Montlaur, A.

2018-03-01

In this work, we study the influence of the contact angle boundary condition on 3D CFD simulations of the bubble generation process occurring in a capillary T-junction. Numerical simulations have been performed with the commercial Computational Fluid Dynamics solver ANSYS Fluent v15.0.7. Experimental results serve as a reference to validate numerical results for four independent parameters: the bubble generation frequency, volume, velocity and length. CFD simulations accurately reproduce experimental results both from qualitative and quantitative points of view. Numerical results are very sensitive to the gas-liquid-wall contact angle boundary conditions, confirming that this is a fundamental parameter to obtain accurate CFD results for simulations of this kind of problems.

Physically-consistent wall boundary conditions for the k-ω turbulence model

DEFF Research Database (Denmark)

Fuhrman, David R.; Dixen, Martin; Jacobsen, Niels Gjøl

2010-01-01

A model solving Reynolds-averaged Navier–Stokes equations, coupled with k-v turbulence closure, is used to simulate steady channel flow on both hydraulically smooth and rough beds. Novel experimental data are used as model validation, with k measured directly from all three components of the fluc......A model solving Reynolds-averaged Navier–Stokes equations, coupled with k-v turbulence closure, is used to simulate steady channel flow on both hydraulically smooth and rough beds. Novel experimental data are used as model validation, with k measured directly from all three components...... of the fluctuating velocity signal. Both conventional k = 0 and dk/dy = 0 wall boundary conditions are considered. Results indicate that either condition can provide accurate solutions, for the bulk of the flow, over both smooth and rough beds. It is argued that the zero-gradient condition is more consistent...... with the near wall physics, however, as it allows direct integration through a viscous sublayer near smooth walls, while avoiding a viscous sublayer near rough walls. This is in contrast to the conventional k = 0 wall boundary condition, which forces resolution of a viscous sublayer in all circumstances...

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.

Derivation of Zagarola-Smits scaling in zero-pressure-gradient turbulent boundary layers

Science.gov (United States)

Wei, Tie; Maciel, Yvan

2018-01-01

This Rapid Communication derives the Zagarola-Smits scaling directly from the governing equations for zero-pressure-gradient turbulent boundary layers (ZPG TBLs). It has long been observed that the scaling of the mean streamwise velocity in turbulent boundary layer flows differs in the near surface region and in the outer layer. In the inner region of small-velocity-defect boundary layers, it is generally accepted that the proper velocity scale is the friction velocity, uτ, and the proper length scale is the viscous length scale, ν /uτ . In the outer region, the most generally used length scale is the boundary layer thickness, δ . However, there is no consensus on velocity scales in the outer layer. Zagarola and Smits [ASME Paper No. FEDSM98-4950 (1998)] proposed a velocity scale, U ZS=(δ1/δ ) U∞ , where δ1 is the displacement thickness and U∞ is the freestream velocity. However, there are some concerns about Zagarola-Smits scaling due to the lack of a theoretical base. In this paper, the Zagarola-Smits scaling is derived directly from a combination of integral, similarity, and order-of-magnitude analysis of the mean continuity equation. The analysis also reveals that V∞, the mean wall-normal velocity at the edge of the boundary layer, is a proper scale for the mean wall-normal velocity V . Extending the analysis to the streamwise mean momentum equation, we find that the Reynolds shear stress in ZPG TBLs scales as U∞V∞ in the outer region. This paper also provides a detailed analysis of the mass and mean momentum balance in the outer region of ZPG TBLs.

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

Hawking radiation, effective actions and covariant boundary conditions

International Nuclear Information System (INIS)

Banerjee, Rabin; Kulkarni, Shailesh

2008-01-01

From an appropriate expression for the effective action, the Hawking radiation from charged black holes is derived, using only covariant boundary conditions at the event horizon. The connection of our approach with the Unruh vacuum and the recent analysis [S.P. Robinson, F. Wilczek, Phys. Rev. Lett. 95 (2005) 011303, (gr-qc/0502074); S. Iso, H. Umetsu, F. Wilczek, Phys. Rev. Lett. 96 (2006) 151302, (hep-th/0602146); R. Banerjee, S. Kulkarni, (arXiv: 0707.2449 [hep-th])] of Hawking radiation using anomalies is established

The role of lateral boundary conditions in simulations of mineral aerosols by a regional climate model of Southwest Asia

Energy Technology Data Exchange (ETDEWEB)

Marcella, Marc Pace [Ralph M. Parsons Laboratory, Department of Civil and Environmental Engineering, Cambridge, MA (United States); Massachusetts Institute of Technology, Cambridge, Massachusetts (United States); Eltahir, Elfatih A.B. [Massachusetts Institute of Technology, Cambridge, Massachusetts (United States)

2012-01-15

The importance of specifying realistic lateral boundary conditions in the regional modeling of mineral aerosols has not been examined previously. This study examines the impact of assigning values for mineral aerosol (dust) concentrations at the lateral boundaries of Regional Climate Model version 3 (RegCM3) and its aerosol model over Southwest Asia. Currently, the dust emission module of RegCM3 operates over the interior of the domain, allowing dust to be transported to the boundaries, but neglecting any dust emitted at these points or from outside the domain. To account for possible dust occurring at, or entering from the boundaries, mixing ratios of dust concentrations from a larger domain RegCM3 simulation are specified at the boundaries of a smaller domain over Southwest Asia. The lateral boundary conditions are monthly averaged concentration values ({mu}g of dust per kg of dry air) resolved in the vertical for all four dust bin sizes within RegCM3's aerosol model. RegCM3 simulations with the aerosol/dust model including lateral boundary conditions for dust are performed for a five year period and compared to model simulations without prescribed dust concentrations at the boundaries. Results indicate that specifying boundary conditions has a significant impact on dust loading across the entire domain over Southwest Asia. More specifically, a nearly 30% increase in aerosol optical depth occurs during the summer months from specifying realistic dust boundary conditions, bringing model results closer to observations such as MISR. In addition, smaller dust particles at the boundaries have a more important impact than large particles in affecting the dust loading within the interior of this domain. Moreover, increases in aerosol optical depth and dust concentrations within the interior domain are not entirely caused by inflow from the boundaries; results indicate that an increase in the gradient of concentration at the boundaries causes an increase of

Thermal boundary conditions for electrons in a weakly ionized gas near a catalytic wall

International Nuclear Information System (INIS)

Chekmarev, I.

1981-01-01

A technique of matched asymptotic expansions is used to examine the derivation of hydrodynamic transport equations for the external region of a weakly ionized multitemperature gas near an absorbing and conducting wall. An approximate moment solution is constructed for the Knudsen boundary layer. The conditions for the matching of the external and internal expansions lead to a new form of the hydrodynamic boundary conditions, from which the singular behavior of the energy equation for electrons near the wall has been eliminated

Recharge and Lateral Groundwater Flow Boundary Conditions for the Saturated Zone Site-Scale Flow and Transport Model

Energy Technology Data Exchange (ETDEWEB)

B. Arnold; T. Corbet

2001-12-18

The purpose of the flow boundary conditions analysis is to provide specified-flux boundary conditions for the saturated zone (SZ) site-scale flow and transport model. This analysis is designed to use existing modeling and analysis results as the basis for estimated groundwater flow rates into the SZ site-scale model domain, both as recharge at the upper (water table) boundary and as underflow at the lateral boundaries. The objective is to provide consistency at the boundaries between the SZ site-scale flow model and other groundwater flow models. The scope of this analysis includes extraction of the volumetric groundwater flow rates simulated by the SZ regional-scale flow model to occur at the lateral boundaries of the SZ site-scale flow model and the internal qualification of the regional-scale model for use in this analysis model report (AMR). In addition, the scope includes compilation of information on the recharge boundary condition taken from three sources: (1) distributed recharge as taken from the SZ regional-scale flow model, (2) recharge below the area of the unsaturated zone (UZ) site-scale flow model, and (3) focused recharge along the Fortymile Wash channel.

A free boundary problem describing the saturated-unsaturated flow in a porous medium

Directory of Open Access Journals (Sweden)

Gabriela Marinoschi

2004-01-01

Full Text Available This paper presents a functional approach to a nonlinear model describing the complete physical process of water infiltration into an unsaturated soil, including the saturation occurrence and the advance of the wetting front. The model introduced in this paper involves a multivalued operator covering the simultaneous saturated and unsaturated flow behaviors and enhances the study of the displacement of the free boundary between these two flow regimes. The model resides in Richards' equation written in pressure form with an initial condition and boundary conditions which in this work express the inflow due to the rain on the soil surface on the one hand, and characterize a certain permeability corresponding to the underground boundary, on the other hand. Existence, uniqueness, and regularity results for the transformed model in diffusive form, that is, for the moisture of the soil, and the existence of the weak solution for the pressure form are proved in the 3D case. The main part of the paper focuses on the existence of the free boundary between the saturated and unsaturated parts of the soil, and this is proved, in the 1D case, for certain stronger assumptions on the initial data and boundary conditions.

Recursive solution for dynamic response of one-dimensional structures with time-dependent boundary conditions

Energy Technology Data Exchange (ETDEWEB)

Abadi, Mohammad Tahaye [Aerospace Research Institute, Tehran (Iran, Islamic Republic of)

2015-10-15

A recursive solution method is derived for the transient response of one-dimensional structures subjected to the general form of time dependent boundary conditions. Unlike previous solution methods that assumed the separation of variables, the present method involves formulating and solving the dynamic problems using the summation of two single-argument functions satisfying the motion equation. Based on boundary and initial conditions, a recursive procedure is derived to determine the single-argument functions. Such a procedure is applied to the general form of boundary conditions, and an analytical solution is derived by solving the recursive equation. The present solution method is implemented for base excitation problems, and the results are compared with those of the previous analytical solution and the Finite element (FE) analysis. The FE results converge to the present analytical solution, although considerable error is found in predicting a solution method on the basis of the separation of variables. The present analytical solution predicts the transient response for wave propagation problems in broadband excitation frequencies.

Recursive solution for dynamic response of one-dimensional structures with time-dependent boundary conditions

International Nuclear Information System (INIS)

Abadi, Mohammad Tahaye

2015-01-01

A recursive solution method is derived for the transient response of one-dimensional structures subjected to the general form of time dependent boundary conditions. Unlike previous solution methods that assumed the separation of variables, the present method involves formulating and solving the dynamic problems using the summation of two single-argument functions satisfying the motion equation. Based on boundary and initial conditions, a recursive procedure is derived to determine the single-argument functions. Such a procedure is applied to the general form of boundary conditions, and an analytical solution is derived by solving the recursive equation. The present solution method is implemented for base excitation problems, and the results are compared with those of the previous analytical solution and the Finite element (FE) analysis. The FE results converge to the present analytical solution, although considerable error is found in predicting a solution method on the basis of the separation of variables. The present analytical solution predicts the transient response for wave propagation problems in broadband excitation frequencies.

Optimum heat power cycles for specified boundary conditions

International Nuclear Information System (INIS)

Ibrahim, O.M.; Klein, S.A.; Mitchell, J.W.

1991-01-01

In this paper optimization of the power output of Carnot and closed Brayton cycles is considered for both finite and infinite thermal capacitance rates of the external fluid streams. The method of Lagrange multipliers is used to solve for working fluid temperatures that yield maximum power. Analytical expressions for the maximum power and the cycle efficiency at maximum power are obtained. A comparison of the maximum power from the two cycles for the same boundary conditions, i.e., the same heat source/sink inlet temperatures, thermal capacitance rates, and heat exchanger conductances, shows that the Brayton cycle can produce more power than the Carnot cycle. This comparison illustrates that cycles exist that can produce more power than the Carnot cycle. The optimum heat power cycle, which will provide the upper limit of power obtained from any thermodynamic cycle for specified boundary conditions and heat exchanger conductances is considered. The optimum heat power cycle is identified by optimizing the sum of the power output from a sequence of Carnot cycles. The shape of the optimum heat power cycle, the power output, and corresponding efficiency are presented. The efficiency at maximum power of all cycles investigated in this study is found to be equal to (or well approximated by) η = 1 - sq. root T L.in /φT H.in where φ is a factor relating the entropy changes during heat rejection and heat addition

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

Science.gov (United States)

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

2017-12-01

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

Sealing of rotary drums for operation under pressurized conditions

International Nuclear Information System (INIS)

Shirvani, M.; Khanof, M. H.; Yousefi, M. R.; Sadighi, S.

2006-01-01

In practice, rotary drums are always designed for operation under vacuum conditions. In this paper, a novel technique is proposed for sealing the rotary drums under pressurized conditions. The proposed system is based on applying a secondary pressurized volume around the leaking gap of the drum. By controlling the pressure of this volume above the pressure of the drum, it will be possible to prevent from any leakage of gases to the ambient. The objective of a controller in this system is that the pressure of secondary volume be kept above the pressure of the drum in spite of the disturbances which may be exerted on the system by the wind outside the drum. The control system is also required to trace the variations in the drum pressure with the least fluctuations in the pressure difference among the drum and the volume

Environment sensitive cracking in pressure boundary materials of light water reactors

International Nuclear Information System (INIS)

Hanninen, H.; Aho-Mantila, I.; Torronen, K.

1987-08-01

A review of the various forms of environment sensitive cracking in pressure boundary materials of light water reactors is presented. The available methods and the most promising future possibilities of preventive maintenance to counteract the environmental degradation are evaluated. Environment sensitive cracking is considered from the metallurgical, mechanical and environmental point of view. The main emphasis is on intergranular stress corrosion cracking of austenitic stainless steels and high strength Ni-base alloys as well as on corrosion fatigue of low alloy and stainless steels. Additionally, some general ideas on how to predict, reduce, monitor or eliminate environment sensitive cracking in service are presented

Toroidal current asymmetry and boundary conditions in disruptions

Science.gov (United States)

Strauss, Henry

2014-10-01

It was discovered on JET that disruptions were accompanied by toroidal asymmetry of the plasma current. The toroidal current asymmetry ΔIϕ is proportional to the vertical current moment ΔMIZ , with positive sign for an upward vertical displacement event (VDE) and negative sign for a downward VDE. It was claimed that this could only be explained by Hiro current. It is shown that instead it is essentially a kinematic effect produced by the VDE displacement of a 3D magnetic perturbation. This is verified by M3D simulations. The simulation results do not require penetration of plasma into the boundary, as in the Hiro current model. It is shown that the normal velocity perpendicular to the magnetic field vanishes at the wall, in the small Larmor radius limit of electromagnetic sheath boundary conditions. Plasma is absorbed into the wall only via the parallel velocity, which is small, penetrates only an infinitesimal distance into the wall, and does not affect forces exerted by the plasma on the wall. Supported by USDOE and ITER.

Effects of boundary conditions on temperature and density in an EXTRAP Z-pinch

International Nuclear Information System (INIS)

Drake, J.R.; Karlsson, P.

1985-08-01

Using the fluid equations, we examine transport in an Extrap configuration by carrying out calculations incorporating model profiles for the density and temperature. The goal of this analysis is to examine the scaling of the pinch equilibrium plasma density, temperature and radius with parameters that are characteristic for Extrap Z-pinches. These parameters include the discharge current, the neutral hydrogen filling density, an oxygen impurity fractional concentration and the condition at the pinch boundary. An Extrap Z-pinch is a pinch discharge where the current channel has a characteristic non-circular cross-section achieved by bounding the discharge by a magnetic separatrix produced when a vacuum octupole magnetic field, generated by currents in external conductors, combines with the self-magnetic field produced by the discharge current. The pinch boundary is changed from a plasma-vacuum boundary to an interface between a high-beta pinch plasma and a low-beta plasma contained in the vacuum magnetic field. The energy that is lost from the pinch region sustains this boundary layer. The introduction of a separatrix boundary around the pinch with four X-point nulls deteriorates the containment of the pinch somewhat. However the presence of the warm, low-beta plasma scrape-off layer, which provides a boundary condition on the pinch, tends to counteract the negative effects of the poorer confinement. Thus the equilibrium parameters that characterize the pinch may not be severely deteriorated by the introduction of the separatrix when the entire configuration, including the scrape-off layer, is considered. (author)

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

Science.gov (United States)

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

2017-10-01

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

Edge states and conformal boundary conditions in super spin chains and super sigma models

International Nuclear Information System (INIS)

Bondesan, Roberto; Jacobsen, Jesper L.; Saleur, Hubert

2011-01-01

The sigma models on projective superspaces CP N+M-1|N with topological angle θ=πmod2π flow to non-unitary, logarithmic conformal field theories in the low-energy limit. In this paper, we determine the exact spectrum of these theories for all open boundary conditions preserving the full global symmetry of the model, generalizing recent work on the particular case M=0 [C. Candu et al., JHEP 1002 (2010) 015]. In the sigma model setting, these boundary conditions are associated with complex line bundles, and are labelled by an integer, related with the exact value of θ. Our approach relies on a spin chain regularization, where the boundary conditions now correspond to the introduction of additional edge states. The exact values of the exponents then follow from a lengthy algebraic analysis, a reformulation of the spin chain in terms of crossing and non-crossing loops (represented as a certain subalgebra of the Brauer algebra), and earlier results on the so-called one- and two-boundary Temperley-Lieb algebras (also known as blob algebras). A remarkable result is that the exponents, in general, turn out to be irrational. The case M=1 has direct applications to the spin quantum Hall effect, which will be discussed in a sequel.