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.

An energy-based equilibrium contact angle boundary condition on jagged surfaces for phase-field methods.

Science.gov (United States)

Frank, Florian; Liu, Chen; Scanziani, Alessio; Alpak, Faruk O; Riviere, Beatrice

2018-08-01

We consider an energy-based boundary condition to impose an equilibrium wetting angle for the Cahn-Hilliard-Navier-Stokes phase-field model on voxel-set-type computational domains. These domains typically stem from μCT (micro computed tomography) imaging of porous rock and approximate a (on μm scale) smooth domain with a certain resolution. Planar surfaces that are perpendicular to the main axes are naturally approximated by a layer of voxels. However, planar surfaces in any other directions and curved surfaces yield a jagged/topologically rough surface approximation by voxels. For the standard Cahn-Hilliard formulation, where the contact angle between the diffuse interface and the domain boundary (fluid-solid interface/wall) is 90°, jagged surfaces have no impact on the contact angle. However, a prescribed contact angle smaller or larger than 90° on jagged voxel surfaces is amplified. As a remedy, we propose the introduction of surface energy correction factors for each fluid-solid voxel face that counterbalance the difference of the voxel-set surface area with the underlying smooth one. The discretization of the model equations is performed with the discontinuous Galerkin method. However, the presented semi-analytical approach of correcting the surface energy is equally applicable to other direct numerical methods such as finite elements, finite volumes, or finite differences, since the correction factors appear in the strong formulation of the model. Copyright © 2018 Elsevier Inc. All rights reserved.

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

The boundary condition for vertical velocity and its interdependence with surface gas exchange

Directory of Open Access Journals (Sweden)

A. S. Kowalski

2017-07-01

Full Text Available The law of conservation of linear momentum is applied to surface gas exchanges, employing scale analysis to diagnose the vertical velocity (w in the boundary layer. Net upward momentum in the surface layer is forced by evaporation (E and defines non-zero vertical motion, with a magnitude defined by the ratio of E to the air density, as w = E/ρ. This is true even right down at the surface where the boundary condition is w|0 = E/ρ|0 (where w|0 and ρ|0 represent the vertical velocity and density of air at the surface. This Stefan flow velocity implies upward transport of a non-diffusive nature that is a general feature of the troposphere but is of particular importance at the surface, where it assists molecular diffusion with upward gas migration (of H2O, for example but opposes that of downward-diffusing species like CO2 during daytime. The definition of flux–gradient relationships (eddy diffusivities requires rectification to exclude non-diffusive transport, which does not depend on scalar gradients. At the microscopic scale, the role of non-diffusive transport in the process of evaporation from inside a narrow tube – with vapour transport into an overlying, horizontal airstream – was described long ago in classical mechanics and is routinely accounted for by chemical engineers, but has been neglected by scientists studying stomatal conductance. Correctly accounting for non-diffusive transport through stomata, which can appreciably reduce net CO2 transport and marginally boost that of water vapour, should improve characterisations of ecosystem and plant functioning.

A Snow Density Dataset for Improving Surface Boundary Conditions in Greenland Ice Sheet Firn Modeling

DEFF Research Database (Denmark)

S. Fausto, Robert; E. Box, Jason; Vandecrux, Baptiste Robert Marcel

2018-01-01

The surface snow density of glaciers and ice sheets is of fundamental importance in converting volume to mass in both altimetry and surface mass balance studies, yet it is often poorly constrained. Site-specific surface snow densities are typically derived from empirical relations based...... on temperature and wind speed. These parameterizations commonly calculate the average density of the top meter of snow, thereby systematically overestimating snow density at the actual surface. Therefore, constraining surface snow density to the top 0.1 m can improve boundary conditions in high-resolution firn......-evolution modeling. We have compiled an extensive dataset of 200 point measurements of surface snow density from firn cores and snow pits on the Greenland ice sheet. We find that surface snow density within 0.1 m of the surface has an average value of 315 kg m−3 with a standard deviation of 44 kg m−3, and has...

Ground Boundary Conditions for Thermal Convection Over Horizontal Surfaces at High Rayleigh Numbers

Science.gov (United States)

Hanjalić, K.; Hrebtov, M.

2016-07-01

We present "wall functions" for treating the ground boundary conditions in the computation of thermal convection over horizontal surfaces at high Rayleigh numbers using coarse numerical grids. The functions are formulated for an algebraic-flux model closed by transport equations for the turbulence kinetic energy, its dissipation rate and scalar variance, but could also be applied to other turbulence models. The three-equation algebraic-flux model, solved in a T-RANS mode ("Transient" Reynolds-averaged Navier-Stokes, based on triple decomposition), was shown earlier to reproduce well a number of generic buoyancy-driven flows over heated surfaces, albeit by integrating equations up to the wall. Here we show that by using a set of wall functions satisfactory results are found for the ensemble-averaged properties even on a very coarse computational grid. This is illustrated by the computations of the time evolution of a penetrative mixed layer and Rayleigh-Bénard (open-ended, 4:4:1 domain) convection, using 10 × 10 × 100 and 10 × 10 × 20 grids, compared also with finer grids (e.g. 60 × 60 × 100), as well as with one-dimensional treatment using 1 × 1 × 100 and 1 × 1 × 20 nodes. The approach is deemed functional for simulations of a convective boundary layer and mesoscale atmospheric flows, and pollutant transport over realistic complex hilly terrain with heat islands, urban and natural canopies, for diurnal cycles, or subjected to other time and space variations in ground conditions and stratification.

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

Heat and mass transfer boundary conditions at the surface of a heated sessile droplet

Science.gov (United States)

Ljung, Anna-Lena; Lundström, T. Staffan

2017-12-01

This work numerically investigates how the boundary conditions of a heated sessile water droplet should be defined in order to include effects of both ambient and internal flow. Significance of water vapor, Marangoni convection, separate simulations of the external and internal flow, and influence of contact angle throughout drying is studied. The quasi-steady simulations are carried out with Computational Fluid Dynamics and conduction, natural convection and Marangoni convection are accounted for inside the droplet. For the studied conditions, a noticeable effect of buoyancy due to evaporation is observed. Hence, the inclusion of moisture increases the maximum velocities in the external flow. Marangoni convection will, in its turn, increase the velocity within the droplet with up to three orders of magnitude. Results furthermore show that the internal and ambient flow can be simulated separately for the conditions studied, and the accuracy is improved if the internal temperature gradient is low, e.g. if Marangoni convection is present. Simultaneous simulations of the domains are however preferred at high plate temperatures if both internal and external flows are dominated by buoyancy and natural convection. The importance of a spatially resolved heat and mass transfer boundary condition is, in its turn, increased if the internal velocity is small or if there is a large variation of the transfer coefficients at the surface. Finally, the results indicate that when the internal convective heat transport is small, a rather constant evaporation rate may be obtained throughout the drying at certain conditions.

Trapping of diffusing particles by striped cylindrical surfaces. Boundary homogenization approach

Science.gov (United States)

Dagdug, Leonardo; Berezhkovskii, Alexander M.; Skvortsov, Alexei T.

2015-01-01

We study trapping of diffusing particles by a cylindrical surface formed by rolling a flat surface, containing alternating absorbing and reflecting stripes, into a tube. For an arbitrary stripe orientation with respect to the tube axis, this problem is intractable analytically because it requires dealing with non-uniform boundary conditions. To bypass this difficulty, we use a boundary homogenization approach which replaces non-uniform boundary conditions on the tube wall by an effective uniform partially absorbing boundary condition with properly chosen effective trapping rate. We demonstrate that the exact solution for the effective trapping rate, known for a flat, striped surface, works very well when this surface is rolled into a cylindrical tube. This is shown for both internal and external problems, where the particles diffuse inside and outside the striped tube, at three orientations of the stripe direction with respect to the tube axis: (a) perpendicular to the axis, (b) parallel to the axis, and (c) at the angle of π/4 to the axis. PMID:26093574

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

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

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

Surface capillary currents: Rediscovery of fluid-structure interaction by forced evolving boundary theory

Science.gov (United States)

Wang, Chunbai; Mitra, Ambar K.

2016-01-01

Any boundary surface evolving in viscous fluid is driven with surface capillary currents. By step function defined for the fluid-structure interface, surface currents are found near a flat wall in a logarithmic form. The general flat-plate boundary layer is demonstrated through the interface kinematics. The dynamics analysis elucidates the relationship of the surface currents with the adhering region as well as the no-slip boundary condition. The wall skin friction coefficient, displacement thickness, and the logarithmic velocity-defect law of the smooth flat-plate boundary-layer flow are derived with the advent of the forced evolving boundary method. This fundamental theory has wide applications in applied science and engineering.

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

Directory of Open Access Journals (Sweden)

P.BalaAnki Reddy

2017-12-01

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

Influence of boundary conditions on the existence and stability of minimal surfaces of revolution made of soap films

Science.gov (United States)

Salkin, Louis; Schmit, Alexandre; Panizza, Pascal; Courbin, Laurent

2014-09-01

Because of surface tension, soap films seek the shape that minimizes their surface energy and thus their surface area. This mathematical postulate allows one to predict the existence and stability of simple minimal surfaces. After briefly recalling classical results obtained in the case of symmetric catenoids that span two circular rings with the same radius, we discuss the role of boundary conditions on such shapes, working with two rings having different radii. We then investigate the conditions of existence and stability of other shapes that include two portions of catenoids connected by a planar soap film and half-symmetric catenoids for which we introduce a method of observation. We report a variety of experimental results including metastability—an hysteretic evolution of the shape taken by a soap film—explained using simple physical arguments. Working by analogy with the theory of phase transitions, we conclude by discussing universal behaviors of the studied minimal surfaces in the vicinity of their existence thresholds.

A Snow Density Dataset for Improving Surface Boundary Conditions in Greenland Ice Sheet Firn Modeling

Directory of Open Access Journals (Sweden)

Robert S. Fausto

2018-05-01

Full Text Available The surface snow density of glaciers and ice sheets is of fundamental importance in converting volume to mass in both altimetry and surface mass balance studies, yet it is often poorly constrained. Site-specific surface snow densities are typically derived from empirical relations based on temperature and wind speed. These parameterizations commonly calculate the average density of the top meter of snow, thereby systematically overestimating snow density at the actual surface. Therefore, constraining surface snow density to the top 0.1 m can improve boundary conditions in high-resolution firn-evolution modeling. We have compiled an extensive dataset of 200 point measurements of surface snow density from firn cores and snow pits on the Greenland ice sheet. We find that surface snow density within 0.1 m of the surface has an average value of 315 kg m−3 with a standard deviation of 44 kg m−3, and has an insignificant annual air temperature dependency. We demonstrate that two widely-used surface snow density parameterizations dependent on temperature systematically overestimate surface snow density over the Greenland ice sheet by 17–19%, and that using a constant density of 315 kg m−3 may give superior results when applied in surface mass budget modeling.

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.

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.

Effects of Thermal Radiation and Chemical Reaction on MHD Free Convection Flow past a Flat Plate with Heat Source and Convective Surface Boundary Condition

OpenAIRE

E.Hemalatha; N. Bhaskar Reddy

2015-01-01

This paper analyzes the radiation and chemical reaction effects on MHD steady two-dimensional laminar viscous incompressible radiating boundary layer flow over a flat plate in the presence of internal heat generation and convective boundary condition. It is assumed that lower surface of the plate is in contact with a hot fluid while a stream of cold fluid flows steadily over the upper surface with a heat source that decays exponentially. The Rosseland approximation is used to desc...

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.

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

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.

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.

SurfCut: Free-Boundary Surface Extraction

KAUST Repository

Algarni, Marei Saeed Mohammed

2016-09-15

We present SurfCut, an algorithm for extracting a smooth simple surface with unknown boundary from a noisy 3D image and a seed point. In contrast to existing approaches that extract smooth simple surfaces with boundary, our method requires less user input, i.e., a seed point, rather than a 3D boundary curve. Our method is built on the novel observation that certain ridge curves of a front propagated using the Fast Marching algorithm are likely to lie on the surface. Using the framework of cubical complexes, we design a novel algorithm to robustly extract such ridge curves and form the surface of interest. Our algorithm automatically cuts these ridge curves to form the surface boundary, and then extracts the surface. Experiments show the robustness of our method to errors in the data, and that we achieve higher accuracy with lower computational cost than comparable methods. © Springer International Publishing AG 2016.

HIFU scattering by the ribs: constrained optimisation with a complex surface impedance boundary condition

Science.gov (United States)

Gélat, P.; ter Haar, G.; Saffari, N.

2014-04-01

High intensity focused ultrasound (HIFU) enables highly localised, non-invasive tissue ablation and its efficacy has been demonstrated in the treatment of a range of cancers, including those of the kidney, prostate and breast. HIFU offers the ability to treat deep-seated tumours locally, and potentially bears fewer side effects than more established treatment modalities such as resection, chemotherapy and ionising radiation. There remains however a number of significant challenges which currently hinder its widespread clinical application. One of these challenges is the need to transmit sufficient energy through the ribcage to ablate tissue at the required foci whilst minimising the formation of side lobes and sparing healthy tissue. Ribs both absorb and reflect ultrasound strongly. This sometimes results in overheating of bone and overlying tissue during treatment, leading to skin burns. Successful treatment of a patient with tumours in the upper abdomen therefore requires a thorough understanding of the way acoustic and thermal energy is deposited. Previously, a boundary element (BE) approach based on a Generalised Minimal Residual (GMRES) implementation of the Burton-Miller formulation was developed to predict the field of a multi-element HIFU array scattered by human ribs, the topology of which was obtained from CT scan data [1]. Dissipative mechanisms inside the propagating medium have since been implemented, together with a complex surface impedance condition at the surface of the ribs. A reformulation of the boundary element equations as a constrained optimisation problem was carried out to determine the complex surface velocities of a multi-element HIFU array which generated the acoustic pressure field that best fitted a required acoustic pressure distribution in a least-squares sense. This was done whilst ensuring that an acoustic dose rate parameter at the surface of the ribs was kept below a specified threshold. The methodology was tested at an

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.

Long-term climate change: the evolution of shield surface boundary conditions

International Nuclear Information System (INIS)

Peltier, W.R.

2007-01-01

The Earths surface during the Pleistocene epoch has been repeatedly subjected to glacial cycles that have markedly influenced both the landscape and surface boundary conditions that, in part, governed past evolution of deep-seated Shield groundwater flow domains. As part of the Deep Geologic Repository Technology Programme simulations of the last Laurentide glacial episode have been undertaken with the University of Toronto Glacial System Model (GSM). The purpose of these simulations is to yield constrained predictions of the magnitude and time rate of change of peri-glacial, glacial and boreal regimes that have perturbed Shield flow domains in the geologic past. A detailed model of long timescale climate change has been developed, which is able to make useful predictions of the process of continental glaciation and deglaciation that has occurred in the past due to the small changes in the effective intensity of the Sun at the location of the Earth caused by gravitational many body effects in Solar System evolution. Based upon the success of this model we are able to assert that we have demonstrated a basic understanding of why this process has continually recurred in the past on a timescale of approximately 100 000 years. Continuing work with the Glacial Systems Model and efforts to provide explicit linkage to numerical analyses of sub-surface hydrology are beginning to yield a new understanding of groundwater flow system evolution and response to glacial perturbations. In so doing this understanding is not only providing a reasoned basis on which to examine issues of geosphere stability as relevant to the safety of a hypothetical repository for used nuclear fuel in Shield terrain, but is also offering an improved basis for the integrated interpretation of multi-disciplinary geo-scientific data necessary for development of a descriptive geosphere model that is seen as fundamental to the repository Safety Case. (author)

Sewing constraints for conformal field theories on surfaces with boundaries

International Nuclear Information System (INIS)

Lewellen, D.C.

1992-01-01

In a conformal field theory, correlation functions on any Riemann surface are in principle unambiguously defined by sewing together three-point functions on the sphere, provided that the four-point functions on the sphere are crossing symmetric, and the one-point functions on the torus are modular covariant. In this work we extend Sonoda's proof of this result to conformal field theories defined on surfaces with boundaries. Four additional sewing constraints arise; three on the half-plane and one on the cylinder. These relate the various OPE coefficients in the theory (bulk, boundary, and bulk-boundary) to one another. In rational theories these relations can be expressed in terms of data arising solely within the bulk theory: The matrix S which implements modular transformations on the characters, and the matrices implementing duality transformations on the four-point conformal-block functions. As an example we solve these relations for the boundary and bulk-boundary structure constants in the Ising model with all possible conformally invariant boundary conditions. The role of the basic sewing constraints in the construction of open string theories is discussed. (orig.)

SurfCut: Free-Boundary Surface Extraction

KAUST Repository

Algarni, Marei Saeed Mohammed; Sundaramoorthi, Ganesh

2016-01-01

We present SurfCut, an algorithm for extracting a smooth simple surface with unknown boundary from a noisy 3D image and a seed point. In contrast to existing approaches that extract smooth simple surfaces with boundary, our method requires less user

Thickened boundary layer theory for air film drag reduction on a van body surface

Science.gov (United States)

Xie, Xiaopeng; Cao, Lifeng; Huang, Heng

2018-05-01

To elucidate drag reduction mechanism on a van body surface under air film condition, a thickened boundary layer theory was proposed and a frictional resistance calculation model of the van body surface was established. The frictional resistance on the van body surface was calculated with different parameters of air film thickness. In addition, the frictional resistance of the van body surface under the air film condition was analyzed by computational fluid dynamics (CFD) simulation and different air film states that influenced the friction resistance on the van body surface were discussed. As supported by the CFD simulation results, the thickened boundary layer theory may provide reference for practical application of air film drag reduction on a van body surface.

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

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.

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

Impact of groundwater capillary rises as lower boundary conditions for soil moisture in a land surface model

Science.gov (United States)

Vergnes, Jean-Pierre; Decharme, Bertrand; Habets, Florence

2014-05-01

Groundwater is a key component of the global hydrological cycle. It sustains base flow in humid climate while it receives seepage in arid region. Moreover, groundwater influences soil moisture through water capillary rise into the soil and potentially affects the energy and water budget between the land surface and the atmosphere. Despite its importance, most global climate models do not account for groundwater and their possible interaction with both the surface hydrology and the overlying atmosphere. This study assesses the impact of capillary rise from shallow groundwater on the simulated water budget over France. The groundwater scheme implemented in the Total Runoff Integrated Pathways (TRIP) river routing model in a previous study is coupled with the Interaction between Soil Biosphere Atmosphere (ISBA) land surface model. In this coupling, the simulated water table depth acts as the lower boundary condition for the soil moisture diffusivity equation. An original parameterization accounting for the subgrid elevation inside each grid cell is proposed in order to compute this fully-coupled soil lower boundary condition. Simulations are performed at high (1/12°) and low (0.5°) resolutions and evaluated over the 1989-2009 period. Compared to a free-drain experiment, upward capillary fluxes at the bottom of soil increase the mean annual evapotranspiration simulated over the aquifer domain by 3.12 % and 1.54 % at fine and low resolutions respectively. This process logically induces a decrease of the simulated recharge from ISBA to the aquifers and contributes to enhance the soil moisture memory. The simulated water table depths are then lowered, which induces a slight decrease of the simulated mean annual river discharges. However, the fully-coupled simulations compare well with river discharge and water table depth observations which confirms the relevance of the coupling formalism.

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

Combined conduction and radiation in a two-layer planar medium with flux boundary condition

International Nuclear Information System (INIS)

Ho, C.H.; Ozisik, M.N.

1987-01-01

The interaction of conduction and radiation is investigated under both transient and steady-state conditions for an absorbing, emitting, and isotropically scattering two-layer slab having opaque coverings at both boundaries. The slab is subjected to an externally applied constant heat flux at one boundary surface and dissipates heat by radiation into external ambients from both boundary surfaces. An analytic approach is applied to solve the radiation part of the problem, and a finite-difference scheme is used to solve the conduction part. The effects of the conduction-to-radiation parameter, the single scattering albedo, the optical thickness, and the surface emissivity on the temperature distribution are examined

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

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

Boundary layer transition observations on a body of revolution with surface heating and cooling in water

Science.gov (United States)

Arakeri, V. H.

1980-04-01

Boundary layer flow visualization in water with surface heat transfer was carried out on a body of revolution which had the predicted possibility of laminar separation under isothermal conditions. Flow visualization was by in-line holographic technique. Boundary layer stabilization, including elimination of laminar separation, was observed to take place on surface heating. Conversely, boundary layer destabilization was observed on surface cooling. These findings are consistent with the theoretical predictions of Wazzan et al. (1970).

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

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

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.

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

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.

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)

Evolution of the lower planetary boundary layer over strongly contrasting surfaces

International Nuclear Information System (INIS)

Coulter, R.L.; Gao, W.; Martin, T.J.; Shannon, J.D.; Doran, J.C.; Hubbe, J.M.; Shaw, W.M.

1992-01-01

In a multilaboratory field study held near Boardman in northeastern Oregon in June 1991, various properties of the surface and lower atmospheric boundary layer over heavily irrigated cropland and adjacent desert steppe were investigated in the initial campaign of the Atmospheric Radiation Measurement (ARM) program. The locale was selected because its disparate characteristics over various spatial scales stress the ability of general circulation models (GCMS) to describe lower boundary conditions, particularly across the discontinuity between desert (in which turbulent flux of heat must be primarily as sensible heat) and large irrigated tracts (in which turbulent flux of latent heat should be the larger term). This campaign of ARM seeks to increase knowledge in three critical areas: (1) determination of the relationships between surface heat fluxes measured over multiple scales and the controlling surface parameters within each scale, (2) integration of local and nearly local heat flux estimates to produce estimates appropriate for GCM grid cells of 100-200 km horizontal dimension, and (3) characterization of the growth and development of the atmospheric boundary layer near transitions between surfaces with strongly contrasting moisture availabilities

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.

Change of Surface Roughness and Planetary Boundary Layer

DEFF Research Database (Denmark)

Jensen, Niels Otto

1978-01-01

The ratio between upstream and far downstream surface friction velocities relative to a change in surface roughness is given on the basis of results from surface Rossby number similarity theory. By simple theories for the internal boundary layer, which are found to compare quite well with recent...... numerical results from higher-order closure models, it is found that, even at a downwind distance such that the internal boundary layer has grown to the full height of the planetary boundary layers, the surface stress still considerably exceeds the equilibrium value...

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

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.

PC-version of RAM6-code for calculation of parameters of the effective logarithmic boundary condition at the absorbent rod surface in reactor

International Nuclear Information System (INIS)

Le Van Ngoc; Ngo Dang Nhan

1990-01-01

The RAM-6 code for calculation of parameters of the effective logarithmic boundary condition at the absorbent rod surface in reactor is suitably modofied to work on IBM PC, the instructions for its usage are presented and capabilities of the personal cpmputer oriented RAM-6 code are demonstrated. (author). 4 refs, 5 tabs, 2 figs

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

Energy Technology Data Exchange (ETDEWEB)

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

2002-06-01

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

Time-resolved PIV measurements of the atmospheric boundary layer over wind-driven surface waves

Science.gov (United States)

Markfort, Corey; Stegmeir, Matt

2017-11-01

Complex interactions at the air-water interface result in two-way coupling between wind-driven surface waves and the atmospheric boundary layer (ABL). Turbulence generated at the surface plays an important role in aquatic ecology and biogeochemistry, exchange of gases such as oxygen and carbon dioxide, and it is important for the transfer of energy and controlling evaporation. Energy transferred from the ABL promotes the generation and maintenance of waves. A fraction of the energy is transferred to the surface mixed layer through the generation of turbulence. Energy is also transferred back to the ABL by waves. There is a need to quantify the details of the coupled boundary layers of the air-water system to better understand how turbulence plays a role in the interactions. We employ time-resolved PIV to measure the detailed structure of the air and water boundary layers under varying wind and wave conditions in the newly developed IIHR Boundary-Layer Wind-Wave Tunnel. The facility combines a 30-m long recirculating water channel with an open-return boundary layer wind tunnel. A thick turbulent boundary layer is developed in the 1 m high air channel, over the water surface, allowing for the study of boundary layer turbulence interacting with a wind-driven wave field.

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

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

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

International Nuclear Information System (INIS)

Chang, C.T.

1982-02-01

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

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

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.

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

System and method for free-boundary surface extraction

KAUST Repository

Algarni, Marei

2017-10-26

A method of extracting surfaces in three-dimensional data includes receiving as inputs three-dimensional data and a seed point p located on a surface to be extracted. The method further includes propagating a front outwardly from the seed point p and extracting a plurality of ridge curves based on the propagated front. A surface boundary is detected based on a comparison of distances between adjacent ridge curves and the desired surface is extracted based on the detected surface boundary.

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.

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

Energy Technology Data Exchange (ETDEWEB)

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

2016-10-15

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

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

Modeling of Hydrophobic Surfaces by the Stokes Problem With the Stick–Slip Boundary Conditions

Czech Academy of Sciences Publication Activity Database

Kučera, R.; Šátek, V.; Haslinger, Jaroslav; Fialová, S.; Pochylý, F.

2017-01-01

Roč. 139, č. 1 (2017), č. článku 011202. ISSN 0098-2202 Institutional support: RVO:68145535 Keywords : algebra * boundary conditions * hydrophobicity * Lagrange multipliers * Navier Stokes equations Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 1.437, year: 2016 http://fluidsengineering.asmedigitalcollection.asme.org/article.aspx?articleid=2536532

The surface roughness and planetary boundary layer

Science.gov (United States)

Telford, James W.

1980-03-01

Applications of the entrainment process to layers at the boundary, which meet the self similarity requirements of the logarithmic profile, have been studied. By accepting that turbulence has dominating scales related in scale length to the height above the surface, a layer structure is postulated wherein exchange is rapid enough to keep the layers internally uniform. The diffusion rate is then controlled by entrainment between layers. It has been shown that theoretical relationships derived on the basis of using a single layer of this type give quantitatively correct factors relating the turbulence, wind and shear stress for very rough surface conditions. For less rough surfaces, the surface boundary layer can be divided into several layers interacting by entrainment across each interface. This analysis leads to the following quantitatively correct formula compared to published measurements. 1 24_2004_Article_BF00877766_TeX2GIFE1.gif {σ _w }/{u^* } = ( {2/{9Aa}} )^{{1/4}} ( {1 - 3^{{1/2}{ a/k{d_n }/z{σ _w }/{u^* }z/L} )^{{1/4}} = 1.28(1 - 0.945({{σ _w }/{u^* }}}) {{z/L}})^{{1/4 where u^* = ( {{tau/ρ}}^{{1/2}}, σ w is the standard deviation of the vertical velocity, z is the height and L is the Obukhov scale lenght. The constants a, A, k and d n are the entrainment constant, the turbulence decay constant, Von Karman's constant, and the layer depth derived from the theory. Of these, a and A, are universal constants and not empirically determined for the boundary layer. Thus the turbulence needed for the plume model of convection, which resides above these layers and reaches to the inversion, is determined by the shear stress and the heat flux in the surface layers. This model applies to convection in cool air over a warm sea. The whole field is now determined except for the temperature of the air relative to the water, and the wind, which need a further parameter describing sea surface roughness. As a first stop to describing a surface where roughness elements

Slarti: A boundary condition editor for a coupled climate model

Science.gov (United States)

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

2006-12-01

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

Boundary conditions 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

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

KAUST Repository

Li, Ping; Jiang, Li; Bagci, Hakan

2015-01-01

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

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

KAUST Repository

Li, Ping

2015-04-24

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

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)

Charge and current transport in open field lines turbulence: Influence of plasma-surface boundary conditions

Energy Technology Data Exchange (ETDEWEB)

Futtersack, R., E-mail: romain.futtersack@cea.fr [CEA, IRFM, F-13108 Saint-Paul-lez-Durance (France); Universite Paul Sabatier Toulouse, LAPLACE, 118 Route de Narbonne, F-31062 Toulouse Cedex 9 (France); Tamain, P. [CEA, IRFM, F-13108 Saint-Paul-lez-Durance (France); Hagelaar, G. [Universite Paul Sabatier Toulouse, LAPLACE, 118 Route de Narbonne, F-31062 Toulouse Cedex 9 (France); Ghendrih, Ph.; Simonin, A. [CEA, IRFM, F-13108 Saint-Paul-lez-Durance (France)

2013-07-15

We investigate the impact of both parallel and transverse boundary conditions on the current and charge transport in open field line systems using the TOKAM2D code, which solves a minimal model for interchange turbulence. Various limit test cases are discussed and analyzed. In the parallel direction, the sheath conductivity is found to play an essential role in the stabilization of large-scale potential structures, leading to the formation of transport channel or transport barrier respectively for an insulating end wall or a wall with an enhanced sheath conductivity. On another hand, the addition of transverse boundary conditions intrinsically changes the transport characteristics, influencing both radial profiles and probability density functions. It underlines that in some cases a detailed description of the plasma-wall interaction process is required to get a proper description of the current loop pattern that determines electrostatic turbulent transport.

Boundary conditions on the plasma emitter surface in the presence of a particle counter flow: I. Ion emitter

Energy Technology Data Exchange (ETDEWEB)

Astrelin, V. T., E-mail: V.T.Astrelin@inp.nsk.su; Kotelnikov, I. A. [Russian Academy of Sciences, Budker Institute of Nuclear Physics, Siberian Branch (Russian Federation)

2017-02-15

Emission of positively charged ions from a plasma emitter irradiated by a counterpropagating electron beam is studied theoretically. A bipolar diode with a plasma emitter in which the ion temperature is lower than the electron temperature and the counter electron flow is extracted from the ion collector is calculated in the one-dimensional model. An analog of Bohm’s criterion for ion emission in the presence of a counterpropagating electron beam is derived. The limiting density of the counterpropagating beam in a bipolar diode operating in the space-charge-limited-emission regime is calculated. The full set of boundary conditions on the plasma emitter surface that are required for operation of the high-current optics module in numerical codes used to simulate charged particle sources is formulated.

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.

Hydrodynamic boundary conditions for one-component liquid-gas flows on non-isothermal solid substrates

KAUST Repository

Xu, Xinpeng

2012-01-01

Recently, liquid-gas flows related to droplets, bubbles, and thin films on solid surfaces with thermal and wettability gradients have attracted widespread attention because of the many physical processes involved and their promising potential applications in biology, chemistry, and industry. Various new physical effects have been discovered at fluid-solid interfaces by experiments and molecular dynamics simulations, e.g., fluid velocity slip, temperature slip (Kapitza resistance), mechanical-thermal cross coupling, etc. There have been various models and theories proposed to explain these experimental and numerical observations. However, to the best of our knowledge,a continuum hydrodynamic model capable of predicting the temperature and velocity profiles of liquid-gas flows on non-isothermal, heterogeneous solid substrates is still absent. The purpose of this work is to construct a continuum model for simulating the liquid-gas flows on solid surfaces that are flat and rigid, and may involve wettability gradients and thermal gradients. This model is able to describe fluid velocity slip, temperature slip, and mechanical-thermal coupling that may occur at fluid-solid interfaces. For this purpose, we first employ the diffuse interface modeling to formulate the hydrodynamic equations for one-component liquid-gas flows in the bulk region. This reproduces the dynamic van der Waals theory of Onuki [Phys. Rev. Lett., 94: 054501, 2005]. We then extendWaldmann\\'s method [Z. Naturforsch. A, 22: 1269-1280, 1967] to formulate the boundary conditions at the fluid-solid interface that match the hydrodynamic equations in the bulk. The effects of the solid surface curvature are also briefly discussed in the appendix. The guiding principles of our model derivation are the conservation laws and the positive definiteness of entropy production together with the Onsager reciprocal relation. The derived model is self-consistent in the sense that the boundary conditions are

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.

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.

Combined measurement of surface, grain boundary and lattice diffusion coefficients on olivine bi-crystals

Science.gov (United States)

Marquardt, Katharina; Dohmen, Ralf; Wagner, Johannes

2014-05-01

Diffusion along interface and grain boundaries provides an efficient pathway and may control chemical transport in rocks as well as their mechanical strength. Besides the significant relevance of these diffusion processes for various geologic processes, experimental data are still very limited (e.g., Dohmen & Milke, 2010). Most of these data were measured using polycrystalline materials and the formalism of LeClaire (1951) to fit integrated concentration depth profiles. To correctly apply this formalism, certain boundary conditions of the diffusion problem need to be fulfilled, e.g., surface diffusion is ignored, and furthermore the lattice diffusion coefficient has to be known from other studies or is an additional fitting parameter, which produces some ambiguity in the derived grain boundary diffusion coefficients. We developed an experimental setup where we can measure the lattice and grain boundary diffusion coefficients simultaneously but independent and demonstrate the relevance of surface diffusion for typical grain boundary diffusion experiments. We performed Mg2SiO4 bicrystal diffusion experiments, where a single grain boundary is covered by a thin-film of pure Ni2SiO4 acting as diffusant source, produced by pulsed laser deposition. The investigated grain boundary is a 60° (011)/[100]. This specific grain boundary configuration was modeled using molecular dynamics for comparison with the experimental observations in the transmission electron microscope (TEM). Both, experiment and model are in good agreement regarding the misorientation, whereas there are still some disagreements regarding the strain fields along the grain boundary that are of outmost importance for the strengths of the material. The subsequent diffusion experiments were carried out in the temperature range between 800° and 1450° C. The inter diffusion profiles were measured using the TEMs energy dispersive x-ray spectrometer standardized using the Cliff-Lorimer equation and EMPA

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.

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

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

Selecting boundary conditions in physiological strain analysis of the femur: Balanced loads, inertia relief method and follower load.

Science.gov (United States)

Heyland, Mark; Trepczynski, Adam; Duda, Georg N; Zehn, Manfred; Schaser, Klaus-Dieter; Märdian, Sven

2015-12-01

Selection of boundary constraints may influence amount and distribution of loads. The purpose of this study is to analyze the potential of inertia relief and follower load to maintain the effects of musculoskeletal loads even under large deflections in patient specific finite element models of intact or fractured bone compared to empiric boundary constraints which have been shown to lead to physiological displacements and surface strains. The goal is to elucidate the use of boundary conditions in strain analyses of bones. Finite element models of the intact femur and a model of clinically relevant fracture stabilization by locking plate fixation were analyzed with normal walking loading conditions for different boundary conditions, specifically re-balanced loading, inertia relief and follower load. Peak principal cortex surface strains for different boundary conditions are consistent (maximum deviation 13.7%) except for inertia relief without force balancing (maximum deviation 108.4%). Influence of follower load on displacements increases with higher deflection in fracture model (from 3% to 7% for force balanced model). For load balanced models, follower load had only minor influence, though the effect increases strongly with higher deflection. Conventional constraints of fixed nodes in space should be carefully reconsidered because their type and position are challenging to justify and for their potential to introduce relevant non-physiological reaction forces. Inertia relief provides an alternative method which yields physiological strain results. Copyright © 2015 IPEM. Published by Elsevier Ltd. All rights reserved.

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

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.

Effective Stress Law in Unconventional Reservoirs under Different Boundary Conditions

Science.gov (United States)

Saurabh, S.; Harpalani, S.

2017-12-01

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

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.

Effect of inlet conditions for numerical modelling of the urban boundary layer

Science.gov (United States)

Gnatowska, Renata

2018-01-01

The paper presents the numerical results obtained with the use of the ANSYS FLUENT commercial code for analysing the flow structure around two rectangular inline surface-mounted bluff bodies immersed in a boundary layer. The effects of the inflow boundary layer for the accuracy of the numerical modelling of the flow field around a simple system of objects are described. The analysis was performed for two concepts. In the former case, the inlet velocity profile was defined using the power law, whereas the kinetic and dissipation energy was defined from the equations according to Richards and Hoxey [1]. In the latter case, the inlet conditions were calculated for the flow over the rough area composed of the rectangular components.

Streamline Patterns and their Bifurcations near a wall with Navier slip Boundary Conditions

DEFF Research Database (Denmark)

Tophøj, Laust; Møller, Søren; Brøns, Morten

2006-01-01

We consider the two-dimensional topology of streamlines near a surface where the Navier slip boundary condition applies. Using transformations to bring the streamfunction in a simple normal form, we obtain bifurcation diagrams of streamline patterns under variation of one or two external parameters....... Topologically, these are identical with the ones previously found for no-slip surfaces. We use the theory to analyze the Stokes flow inside a circle, and show how it can be used to predict new bifurcation phenomena. ©2006 American Institute of Physics...

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.

Automated Boundary Conditions for Wind Tunnel Simulations

Science.gov (United States)

Carlson, Jan-Renee

2018-01-01

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

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

International Nuclear Information System (INIS)

Adams, J.; Pneuman, G.W.

1976-01-01

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

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

Interaction between particles and grain boundaries under conditions of cooperative migration

International Nuclear Information System (INIS)

Marvina, L.A.; Marvin, V.B.

1996-01-01

The analysis of particle grain boundary interaction is performed for dispersion hardened alloys when cooperative migration takes place. It is shown that in a general case the particle experiences a Zener force and a force of grain boundary surface tension due to boundary bending between particles. Approximate numerical estimates are made for the force acting on a particle in dispersion hardened alloy Ni-HfO 2 . It is noted that during cooperative migration of particle and grain boundary the velocity of the particle is directed along the resulting force. The latter equals the sum of surface tension and Zener forces. 6 refs., 2 figs

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.

Influence of model boundary conditions on blood flow patterns in a patient specific stenotic right coronary artery.

Science.gov (United States)

Liu, Biyue; Zheng, Jie; Bach, Richard; Tang, Dalin

2015-01-01

inlet and the outlet did not have significant effect on the spatial distribution patterns of the PD, the WPG and the WSS on the lumen surface, regarding the locations of the maximum and the minimum of each quantity. The observations from this study indicated that the ways how pressure and velocity boundary conditions are imposed in computational models have considerable impact on flow velocity and shear stress predictions. Accuracy of in vivo measurements of blood pressure and velocity is of great importance for reliable model predictions.

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.

Recycling inflow method for simulations of spatially evolving turbulent boundary layers over rough surfaces

Science.gov (United States)

Yang, Xiang I. A.; Meneveau, Charles

2016-01-01

The technique by Lund et al. to generate turbulent inflow for simulations of developing boundary layers over smooth flat plates is extended to the case of surfaces with roughness elements. In the Lund et al. method, turbulent velocities on a sampling plane are rescaled and recycled back to the inlet as inflow boundary condition. To rescale mean and fluctuating velocities, appropriate length scales need be identified and for smooth surfaces, the viscous scale lν = ν/uτ (where ν is the kinematic viscosity and uτ is the friction velocity) is employed for the inner layer. Different from smooth surfaces, in rough wall boundary layers the length scale of the inner layer, i.e. the roughness sub-layer scale ld, must be determined by the geometric details of the surface roughness elements and the flow around them. In the proposed approach, it is determined by diagnosing dispersive stresses that quantify the spatial inhomogeneity caused by the roughness elements in the flow. The scale ld is used for rescaling in the inner layer, and the boundary layer thickness δ is used in the outer region. Both parts are then combined for recycling using a blending function. Unlike the blending function proposed by Lund et al. which transitions from the inner layer to the outer layer at approximately 0.2δ, here the location of blending is shifted upwards to enable simulations of very rough surfaces in which the roughness length may exceed the height of 0.2δ assumed in the traditional method. The extended rescaling-recycling method is tested in large eddy simulation of flow over surfaces with various types of roughness element shapes.

West-WRF Sensitivity to Sea Surface Temperature Boundary Condition in California Precipitation Forecasts of AR Related Events

Science.gov (United States)

Zhang, X.; Cornuelle, B. D.; Martin, A.; Weihs, R. R.; Ralph, M.

2017-12-01

We evaluated the merit in coastal precipitation forecasts by inclusion of high resolution sea surface temperature (SST) from blended satellite and in situ observations as a boundary condition (BC) to the Weather Research and Forecast (WRF) mesoscale model through simple perturbation tests. Our sensitivity analyses shows that the limited improvement of watershed scale precipitation forecast is credible. When only SST BC is changed, there is an uncertainty introduced because of artificial model state equilibrium and the nonlinear nature of the WRF model system. With the change of SST on the order of a fraction of a degree centigrade, we found that the part of random perturbation forecast response is saturated after 48 hours when it reaches to the order magnitude of the linear response. It is important to update the SST at a shorter time period, so that the independent excited nonlinear modes can cancel each other. The uncertainty in our SST configuration is quantitatively equivalent to adding to a spatially uncorrelated Guasian noise of zero mean and 0.05 degree of standard deviation to the SST. At this random noise perturbation magnitude, the ensemble average behaves well within a convergent range. It is also found that the sensitivity of forecast changes in response to SST changes. This is measured by the ratio of the spatial variability of mean of the ensemble perturbations over the spatial variability of the corresponding forecast. The ratio is about 10% for surface latent heat flux, 5 % for IWV, and less than 1% for surface pressure.

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

Science.gov (United States)

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

2017-08-01

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

Classically integrable boundary conditions for affine Toda field theories

International Nuclear Information System (INIS)

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

1995-01-01

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

Facilitating conditions for boundary-spanning behavior in governance networks

OpenAIRE

Meerkerk, Ingmar; Edelenbos, Jurian

2017-01-01

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

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

Directory of Open Access Journals (Sweden)

Yibao Li

2017-01-01

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

MHD Heat and Mass Transfer of Chemical Reaction Fluid Flow over a Moving Vertical Plate in Presence of Heat Source with Convective Surface Boundary Condition

Directory of Open Access Journals (Sweden)

B. R. Rout

2013-01-01

Full Text Available This paper aims to investigate the influence of chemical reaction and the combined effects of internal heat generation and a convective boundary condition on the laminar boundary layer MHD heat and mass transfer flow over a moving vertical flat plate. The lower surface of the plate is in contact with a hot fluid while the stream of cold fluid flows over the upper surface with heat source and chemical reaction. The basic equations governing the flow, heat transfer, and concentration are reduced to a set of ordinary differential equations by using appropriate transformation for variables and solved numerically by Runge-Kutta fourth-order integration scheme in association with shooting method. The effects of physical parameters on the velocity, temperature, and concentration profiles are illustrated graphically. A table recording the values of skin friction, heat transfer, and mass transfer at the plate is also presented. The discussion focuses on the physical interpretation of the results as well as their comparison with previous studies which shows good agreement as a special case of the problem.

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)

Influence of Spanwise Boundary Conditions on Slat Noise Simulations

Science.gov (United States)

Lockard, David P.; Choudhari, Meelan M.; Buning, Pieter G.

2015-01-01

The slat noise from the 30P/30N high-lift system is being investigated through computational fluid dynamics simulations with the OVERFLOW code in conjunction with a Ffowcs Williams-Hawkings acoustics solver. In the present study, two different spanwise grids are being used to investigate the effect of the spanwise extent and periodicity on the near-field unsteady structures and radiated noise. The baseline grid with periodic boundary conditions has a short span equal to 1/9th of the stowed chord, whereas the other, longer span grid adds stretched grids on both sides of the core, baseline grid to allow inviscid surface boundary conditions at both ends. The results indicate that the near-field mean statistics obtained using the two grids are similar to each other, as are the directivity and spectral shapes of the radiated noise. However, periodicity forces all acoustic waves with less than one wavelength across the span to be two-dimensional, without any variation in the span. The spanwise coherence of the acoustic waves is what is needed to make estimates of the noise that would be radiated from realistic span lengths. Simulations with periodic conditions need spans of at least six slat chords to allow spanwise variation in the low-frequencies associated with the peak of broadband slat noise. Even then, the full influence of the periodicity is unclear, so employing grids with a fine, central region and highly stretched meshes that go to slip walls may be a more efficient means of capturing the spanwise decorrelation of low-frequency acoustic phenomena.

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.

Riemann surfaces with boundaries and string theory

International Nuclear Information System (INIS)

Morozov, A.Yu.; Roslyj, A.A.

1989-01-01

A consideration of the cutting and joining operations for Riemann surfaces permits one to express the functional integral on a Riemann surface in terms of integrals over its pieces which are suarfaces with boundaries. This yields an expression for the determinant of the Laplacian on a Riemann surface in terms of Krichever maps for its pieces. Possible applications of the methods proposed to a study of the string perturbation theory in terms of an universal moduli space are mentioned

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.

Stepping towards new parameterizations for non-canonical atmospheric surface-layer conditions

Science.gov (United States)

Calaf, M.; Margairaz, F.; Pardyjak, E.

2017-12-01

Representing land-atmosphere exchange processes as a lower boundary condition remains a challenge. This is partially a result of the fact that land-surface heterogeneity exists at all spatial scales and its variability does not "average" out with decreasing scales. Such variability need not rapidly blend away from the boundary thereby impacting the near-surface region of the atmosphere. Traditionally, momentum and energy fluxes linking the land surface to the flow in NWP models have been parameterized using atmospheric surface layer (ASL) similarity theory. There is ample evidence that such representation is acceptable for stationary and planar-homogeneous flows in the absence of subsidence. However, heterogeneity remains a ubiquitous feature eliciting appreciable deviations when using ASL similarity theory, especially in scalars such moisture and air temperature whose blending is less efficient when compared to momentum. The focus of this project is to quantify the effect of surface thermal heterogeneity with scales Ο(1/10) the height of the atmospheric boundary layer and characterized by uniform roughness. Such near-canonical cases describe inhomogeneous scalar transport in an otherwise planar homogeneous flow when thermal stratification is weak or absent. In this work we present a large-eddy simulation study that characterizes the effect of surface thermal heterogeneities on the atmospheric flow using the concept of dispersive fluxes. Results illustrate a regime in which the flow is mostly driven by the surface thermal heterogeneities, in which the contribution of the dispersive fluxes can account for up to 40% of the total sensible heat flux. Results also illustrate an alternative regime in which the effect of the surface thermal heterogeneities is quickly blended, and the dispersive fluxes provide instead a quantification of the flow spatial heterogeneities produced by coherent turbulent structures result of the surface shear stress. A threshold flow

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.

Areal-averaged trace gas emission rates from long-range open-path measurements in stable boundary layer conditions

Directory of Open Access Journals (Sweden)

K. Schäfer

2012-07-01

Full Text Available Measurements of land-surface emission rates of greenhouse and other gases at large spatial scales (10 000 m² are needed to assess the spatial distribution of emissions. This can be readily done using spatial-integrating micro-meteorological methods like flux-gradient methods which were evaluated for determining land-surface emission rates of trace gases under stable boundary layers. Non-intrusive path-integrating measurements are utilized. Successful application of a flux-gradient method requires confidence in the gradients of trace gas concentration and wind, and in the applicability of boundary-layer turbulence theory; consequently the procedures to qualify measurements that can be used to determine the flux is critical. While there is relatively high confidence in flux measurements made under unstable atmospheres with mean winds greater than 1 m s⁻¹, there is greater uncertainty in flux measurements made under free convective or stable conditions. The study of N₂O emissions of flat grassland and NH₃ emissions from a cattle lagoon involves quality-assured determinations of fluxes under low wind, stable or night-time atmospheric conditions when the continuous "steady-state" turbulence of the surface boundary layer breaks down and the layer has intermittent turbulence. Results indicate that following the Monin-Obukhov similarity theory (MOST flux-gradient methods that assume a log-linear profile of the wind speed and concentration gradient incorrectly determine vertical profiles and thus flux in the stable boundary layer. An alternative approach is considered on the basis of turbulent diffusivity, i.e. the measured friction velocity as well as height gradients of horizontal wind speeds and concentrations without MOST correction for stability. It is shown that this is the most accurate of the flux-gradient methods under stable conditions.

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.

Discrete Boltzmann Method with Maxwell-Type Boundary Condition for Slip Flow

Science.gov (United States)

Zhang, Yu-Dong; Xu, Ai-Guo; Zhang, Guang-Cai; Chen, Zhi-Hua

2018-01-01

The rarefied effect of gas flow in microchannel is significant and cannot be well described by traditional hydrodynamic models. It has been known that discrete Boltzmann model (DBM) has the potential to investigate flows in a relatively wider range of Knudsen number because of its intrinsic kinetic nature inherited from Boltzmann equation. It is crucial to have a proper kinetic boundary condition for DBM to capture the velocity slip and the flow characteristics in the Knudsen layer. In this paper, we present a DBM combined with Maxwell-type boundary condition model for slip flow. The tangential momentum accommodation coefficient is introduced to implement a gas-surface interaction model. Both the velocity slip and the Knudsen layer under various Knudsen numbers and accommodation coefficients can be well described. Two kinds of slip flows, including Couette flow and Poiseuille flow, are simulated to verify the model. To dynamically compare results from different models, the relation between the definition of Knudsen number in hard sphere model and that in BGK model is clarified. Support of National Natural Science Foundation of China under Grant Nos. 11475028, 11772064, and 11502117 Science Challenge Project under Grant Nos. JCKY2016212A501 and TZ2016002

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.

Space-like surfaces with free boundary in the Lorentz-Minkowski space

International Nuclear Information System (INIS)

López, R; Pyo, J

2012-01-01

We investigate a variational problem in the Lorentz-Minkowski space L 3 whose critical points are space-like surfaces with a constant mean curvature and making a constant contact angle with a given support surface along its common boundary. We show that if the support surface is a pseudosphere, then the surface is a planar disc or a hyperbolic cap. We also study the problem of space-like hypersurfaces with free boundary in the higher dimensional Lorentz-Minkowski space L n+1 . (paper)

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.

Air Quality and Meteorological Boundary Conditions during the MCMA-2003 Field Campaign

Science.gov (United States)

Sosa, G.; Arriaga, J.; Vega, E.; Magaña, V.; Caetano, E.; de Foy, B.; Molina, L. T.; Molina, M. J.; Ramos, R.; Retama, A.; Zaragoza, J.; Martínez, A. P.; Márquez, C.; Cárdenas, B.; Lamb, B.; Velasco, E.; Allwine, E.; Pressley, S.; Westberg, H.; Reyes, R.

2004-12-01

A comprehensive field campaign to characterize photochemical smog in the Mexico City Metropolitan Area (MCMA) was conducted during April 2003. An important number of equipment was deployed all around the urban core and its surroundings to measure gas and particles composition from the various sources and receptor sites. In addition to air quality measurements, meteorology variables were also taken by regular weather meteorological stations, tethered balloons, radiosondes, sodars and lidars. One important issue with regard to the field campaign was the characterization of the boundary conditions in order to feed meteorological and air quality models. Four boundary sites were selected to measure continuously criteria pollutants, VOC and meteorological variables at surface level. Vertical meteorological profiles were measured at three other sites : radiosondes in Tacubaya site were launched every six hours daily; tethered balloons were launched at CENICA and FES-Cuautitlan sites according to the weather conditions, and one sodar was deployed at UNAM site in the south of the city. Additionally to these measurements, two fixed meteorological monitoring networks deployed along the city were available to complement these measurements. In general, we observed that transport of pollutants from the city to the boundary sites changes every day, according to the coupling between synoptic and local winds. This effect were less important at elevated sites such as Cerro de la Catedral and ININ, where synoptic wind were more dominant during the field campaign. Also, local sources nearby boundary sites hide the influence of pollution coming from the city some days, particularly at the La Reforma site.

Thermocouple Rakes for Measuring Boundary Layer Flows Extremely Close to Surface

Science.gov (United States)

Hwang, Danny P.; Fralick, Gustave C.; Martin, Lisa C.; Blaha, Charles A.

2001-01-01

Of vital interest to aerodynamic researchers is precise knowledge of the flow velocity profile next to the surface. This information is needed for turbulence model development and the calculation of viscous shear force. Though many instruments can determine the flow velocity profile near the surface, none of them can make measurements closer than approximately 0.01 in. from the surface. The thermocouple boundary-layer rake can measure much closer to the surface than conventional instruments can, such as a total pressure boundary layer rake, hot wire, or hot film. By embedding the sensors (thermocouples) in the region where the velocity is equivalent to the velocity ahead of a constant thickness strut, the boundary-layer flow profile can be obtained. The present device fabricated at the NASA Glenn Research Center microsystem clean room has a heater made of platinum and thermocouples made of platinum and gold. Equal numbers of thermocouples are placed both upstream and downstream of the heater, so that the voltage generated by each pair at the same distance from the surface is indicative of the difference in temperature between the upstream and downstream thermocouple locations. This voltage differential is a function of the flow velocity, and like the conventional total pressure rake, it can provide the velocity profile. In order to measure flow extremely close to the surface, the strut is made of fused quartz with extremely low heat conductivity. A large size thermocouple boundary layer rake is shown in the following photo. The latest medium size sensors already provide smooth velocity profiles well into the boundary layer, as close as 0.0025 in. from the surface. This is about 4 times closer to the surface than the previously used total pressure rakes. This device also has the advantage of providing the flow profile of separated flow and also it is possible to measure simultaneous turbulence levels within the boundary layer.

Global smooth solutions of 3-D null-form wave equations in exterior domains with Neumann boundary conditions

Science.gov (United States)

Jun, Li; Huicheng, Yin

2018-05-01

The paper is devoted to investigating long time behavior of smooth small data solutions to 3-D quasilinear wave equations outside of compact convex obstacles with Neumann boundary conditions. Concretely speaking, when the surface of a 3-D compact convex obstacle is smooth and the quasilinear wave equation fulfills the null condition, we prove that the smooth small data solution exists globally provided that the Neumann boundary condition on the exterior domain is given. One of the main ingredients in the current paper is the establishment of local energy decay estimates of the solution itself. As an application of the main result, the global stability to 3-D static compressible Chaplygin gases in exterior domain is shown under the initial irrotational perturbation with small amplitude.

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

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

International Nuclear Information System (INIS)

Lu Jia; Zhou Huaichun

2016-01-01

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

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.

Convective and global stability analysis of a Mach 5.8 boundary layer grazing a compliant surface

Science.gov (United States)

Dettenrieder, Fabian; Bodony, Daniel

2016-11-01

Boundary layer transition on high-speed vehicles is expected to be affected by unsteady surface compliance. The stability properties of a Mach 5.8 zero-pressure-gradient laminar boundary layer grazing a nominally-flat thermo-mechanically compliant panel is considered. The linearized compressible Navier-Stokes equations describe small amplitude disturbances in the fluid while the panel deformations are described by the Kirchhoff-Love plate equation and its thermal state by the transient heat equation. Compatibility conditions that couple disturbances in the fluid to those in the solid yield simple algebraic and robin boundary conditions for the velocity and thermal states, respectively. A local convective stability analysis shows that the panel can modify both the first and second Mack modes when, for metallic-like panels, the panel thickness exceeds the lengthscale δ99 Rex- 0 . 5 . A global stability analysis, which permits finite panel lengths with clamped-clamped boundary conditions, shows a rich eigenvalue spectrum with several branches. Unstable modes are found with streamwise-growing panel deformations leading to Mach wave-type radiation. Stable global modes are also found and have distinctly different panel modes but similar radiation patterns. Air Force Office of Scientific Research.

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

Science.gov (United States)

Colomo, F.; Pronko, A. G.

2005-05-01

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

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

Science.gov (United States)

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

2018-04-01

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

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

Energy Technology Data Exchange (ETDEWEB)

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

2015-11-15

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

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

Body surface adaptations to boundary-layer dynamics

NARCIS (Netherlands)

Videler, J.J.

1995-01-01

Evolutionary processes have adapted nektonic animals to interact efficiently with the water that surrounds them. Not all these adaptations serve the same purpose. This paper concentrates on reduction of drag due to friction in the boundary layer close to the body surface. Mucus, compliant skins,

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

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)

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)

A boundary integral equation method using auxiliary interior surface approach for acoustic radiation and scattering in two dimensions.

Science.gov (United States)

Yang, S A

2002-10-01

This paper presents an effective solution method for predicting acoustic radiation and scattering fields in two dimensions. The difficulty of the fictitious characteristic frequency is overcome by incorporating an auxiliary interior surface that satisfies certain boundary condition into the body surface. This process gives rise to a set of uniquely solvable boundary integral equations. Distributing monopoles with unknown strengths over the body and interior surfaces yields the simple source formulation. The modified boundary integral equations are further transformed to ordinary ones that contain nonsingular kernels only. This implementation allows direct application of standard quadrature formulas over the entire integration domain; that is, the collocation points are exactly the positions at which the integration points are located. Selecting the interior surface is an easy task. Moreover, only a few corresponding interior nodal points are sufficient for the computation. Numerical calculations consist of the acoustic radiation and scattering by acoustically hard elliptic and rectangular cylinders. Comparisons with analytical solutions are made. Numerical results demonstrate the efficiency and accuracy of the current solution method.

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.

Laminar boundary layer response to rotation of a finite diameter surface patch

International Nuclear Information System (INIS)

Klewicki, J.C.; Hill, R.B.

2003-01-01

The responses of the flat plate laminar boundary layer to perturbations generated by rotating a finite patch of the bounding surface are explored experimentally. The size of the surface patch was of the same order as the boundary layer thickness. The displacement thickness Reynolds number range of the boundary layers explored was 72-527. The rotation rates of the surface patch ranged from 2.14 to 62.8 s-1. Qualitative flow visualizations and quantitative molecular tagging velocimetry measurements revealed that rotation of a finite surface patch generates an asymmetric loop-like vortex. Significant features of this vortex include that, (i) the sign of the vorticity in the vortex head is opposite that of the boundary layer vorticity regardless of the sign of the input rotation, (ii) one leg of the vortex exhibits motion akin to solid body rotation while the other leg is best characterized as a spanwise shear layer, (iii) the vortex leg exhibiting near solid body rotation lifts more rapidly from the surface than the leg more like a shear layer, and (iv) the vortex leg exhibiting near solid body rotation always occurs on the side of the surface patch experiencing downstream motion. These asymmetries switch sides depending on the sign of the input rotation. The present results are interpreted and discussed relative to analytical solutions for infinite geometries. By way of analogy, plausible connections are drawn between the present results and the influences of wall normal vortices in turbulent boundary layer flows

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

DEFF Research Database (Denmark)

Escolano-Carrasco, José; Jacobsen, Finn

2007-01-01

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

Numerical model simulations of boundary-layer dynamics during winter conditions

DEFF Research Database (Denmark)

Melas, D.; Persson, T.; Bruin, H. de

2001-01-01

A mesoscale numerical model, incorporating a land-surface scheme based on Deardorffs' approach, is used to study the diurnal variation of the boundary layer structure and surface fluxes during four consecutive days with air temperatures well below zero, snow covered ground and changing synoptic f...

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.

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

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

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)

Drag reduction using wrinkled surfaces in high Reynolds number laminar boundary layer flows

Science.gov (United States)

Raayai-Ardakani, Shabnam; McKinley, Gareth H.

2017-09-01

Inspired by the design of the ribbed structure of shark skin, passive drag reduction methods using stream-wise riblet surfaces have previously been developed and tested over a wide range of flow conditions. Such textures aligned in the flow direction have been shown to be able to reduce skin friction drag by 4%-8%. Here, we explore the effects of periodic sinusoidal riblet surfaces aligned in the flow direction (also known as a "wrinkled" texture) on the evolution of a laminar boundary layer flow. Using numerical analysis with the open source Computational Fluid Dynamics solver OpenFOAM, boundary layer flow over sinusoidal wrinkled plates with a range of wavelength to plate length ratios ( λ / L ), aspect ratios ( 2 A / λ ), and inlet velocities are examined. It is shown that in the laminar boundary layer regime, the riblets are able to retard the viscous flow inside the grooves creating a cushion of stagnant fluid that the high-speed fluid above can partially slide over, thus reducing the shear stress inside the grooves and the total integrated viscous drag force on the plate. Additionally, we explore how the boundary layer thickness, local average shear stress distribution, and total drag force on the wrinkled plate vary with the aspect ratio of the riblets as well as the length of the plate. We show that riblets with an aspect ratio of close to unity lead to the highest reduction in the total drag, and that because of the interplay between the local stress distribution on the plate and stream-wise evolution of the boundary layer the plate has to exceed a critical length to give a net decrease in the total drag force.

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.

Chevron defect at the intersection of grain boundaries with free surfaces in Au

International Nuclear Information System (INIS)

Radetic, T.; Lancon, F.; Dahmen, U.

2002-01-01

We have identified a new defect at the intersection between grain boundaries and surfaces in Au using atomic resolution transmission electron microscopy. At the junction line of 90 deg. tilt grain boundaries of (110)-(001) orientation with the free surface, a small segment of the grain boundary, about 1 nm in length, dissociates into a triangular region with a chevronlike stacking disorder and a distorted hcp structure. The structure and stability of these defects are confirmed by atomistic simulations, and we point out the relationship with the one-dimensional incommensurate structure of the grain boundary

Modelling the artic stable boundary layer and its coupling to the surface

NARCIS (Netherlands)

Steeneveld, G.J.; Wiel, van de B.J.H.; Holtslag, A.A.M.

2006-01-01

The impact of coupling the atmosphere to the surface energy balance is examined for the stable boundary layer, as an extension of the first GABLS (GEWEX Atmospheric Boundary-Layer Study) one-dimensional model intercomparison. This coupling is of major importance for the stable boundary-layer

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.

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”

Nonlinear interaction of the surface waves at a plasma boundary

International Nuclear Information System (INIS)

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

1976-01-01

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

Euler–Lagrange simulation of gas–solid pipe flow with smooth and rough wall boundary conditions

DEFF Research Database (Denmark)

Mandø, Matthias; Yin, Chungen

2012-01-01

Numerical simulation of upward turbulent particle-laden pipe flow is performed with the intention to reveal the influence of surface roughness on the velocity statistics of the particle phase. A rough wall collision model, which models the surface as being sinusoidal, is proposed to account...... for the wall boundary condition ranging for smooth surfaces to very rough surfaces. This model accounts for the entire range of possible surface roughness found in pipes and industrial pneumatic equipment from smooth plastic pipes over machined steel pipes to cast iron surfaces. The model is based...... on a geometric interpretation of the wall collision process where the particle restitution coefficient is based on the data presented by Sommerfeld and Huber [1]. Simulations are performed using the Eulerian–Lagrangian methodology for the dilute one-way coupling regime. Results are reported for 3 different sizes...

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

Boundary surface and microstructure analysis of ceramic materials

International Nuclear Information System (INIS)

Woltersdorf, J.; Pippel, E.

1992-01-01

The article introduces the many possibilities of high voltage (HVEM) and high resolution electron microscopy (HREM) for boundary surface and microstructure analysis of ceramic materials. The investigations are limited to ceramic long fibre composites and a ceramic fibre/glass matrix system. (DG) [de

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.

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

Fictitious domain methods for elliptic problems with general boundary conditions with an application to the numerical simulation of two phase flows

International Nuclear Information System (INIS)

Ramiere, I.

2006-09-01

This work is dedicated to the introduction of two original fictitious domain methods for the resolution of elliptic problems (mainly convection-diffusion problems) with general and eventually mixed boundary conditions: Dirichlet, Robin or Neumann. The originality lies in the approximation of the immersed boundary by an approximate interface derived from the fictitious domain Cartesian mesh, which is generally not boundary-fitted to the physical domain. The same generic numerical scheme is used to impose the embedded boundary conditions. Hence, these methods require neither a surface mesh of the immersed boundary nor the local modification of the numerical scheme. We study two modelling of the immersed boundary. In the first one, called spread interface, the approximate immersed boundary is the union of the cells crossed by the physical immersed boundary. In the second one, called thin interface, the approximate immersed boundary lies on sides of mesh cells. Additional algebraic transmission conditions linking both flux and solution jumps through the thin approximate interface are introduced. The fictitious problem to solve as well as the treatment of the embedded boundary conditions are detailed for the two methods. A Q1 finite element scheme is implemented for the numerical validation of the spread interface approach while a new cell-centered finite volume scheme is derived for the thin interface approach with immersed jumps. Each method is then combined to multilevel local mesh refinement algorithms (with solution or flux residual) to increase the precision of the solution in the vicinity of the immersed interface. A convergence analysis of a Q1 finite element method with non-boundary fitted meshes is also presented. This study proves the convergence rates of the present methods. Among the various industrial applications, the simulation on a model of heat exchanger in french nuclear power plants enables us to appreciate the performances of the fictitious domain

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

Compressible turbulent channel flow with impedance boundary conditions

Science.gov (United States)

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

2015-03-01

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

Boundary-Layer Separation Control under Low-Pressure Turbine Airfoil Conditions using Glow-Discharge Plasma Actuators

Science.gov (United States)

Hultgren, Lennart S.; Ashpis, David E.

2003-01-01

Modem low-pressure turbines, in general, utilize highly loaded airfoils in an effort to improve efficiency and to lower the number of airfoils needed. Typically, the airfoil boundary layers are turbulent and fully attached at takeoff conditions, whereas a substantial fraction of the boundary layers on the airfoils may be transitional at cruise conditions due to the change of density with altitude. The strong adverse pressure gradients on the suction side of these airfoils can lead to boundary-layer separation at the latter low Reynolds number conditions. Large separation bubbles, particularly those which fail to reattach, cause a significant degradation of engine efficiency. A component efficiency drop of the order 2% may occur between takeoff and cruise conditions for large commercial transport engines and could be as large as 7% for smaller engines at higher altitude. An efficient means of of separation elimination/reduction is, therefore, crucial to improved turbine design. Because the large change in the Reynolds number from takeoff to cruise leads to a distinct change in the airfoil flow physics, a separation control strategy intended for cruise conditions will need to be carefully constructed so as to incur minimum impact/penalty at takeoff. A complicating factor, but also a potential advantage in the quest for an efficient strategy, is the intricate interplay between separation and transition for the situation at hand. Volino gives a comprehensive discussion of several recent studies on transition and separation under low-pressure-turbine conditions, among them one in the present facility. Transition may begin before or after separation, depending on the Reynolds number and other flow conditions. If the transition occurs early in the boundary layer then separation may be reduced or completely eliminated. Transition in the shear layer of a separation bubble can lead to rapid reattachment. This suggests using control mechanisms to trigger and enhance early

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

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.

Electron scattering at surfaces and grain boundaries in thin Au films

International Nuclear Information System (INIS)

Henriquez, Ricardo; Flores, Marcos; Moraga, Luis; Kremer, German; González-Fuentes, Claudio; Munoz, Raul C.

2013-01-01

The electron scattering at surfaces and grain boundaries is investigated using polycrystalline Au films deposited onto mica substrates. We vary the three length scales associated with: (i) electron scattering in the bulk, that at temperature T is characterized by the electronic mean free path in the bulk ℓ 0 (T); (ii) electron-surface scattering, that is characterized by the film thickness t; (iii) electron-grain boundary scattering, that is characterized by the mean grain diameter D. We varied independently the film thickness from approximately 50 nm to about 100 nm, and the typical grain size making up the samples from 12 nm to 160 nm. We also varied the scale of length associated with electron scattering in the bulk by measuring the resistivity of each specimen at temperatures T, 4 K 0 (T) by approximately 2 orders of magnitude. Detailed measurements of the grain size distribution as well as surface roughness of each sample were performed with a Scanning Tunnelling Microscope (STM). We compare, for the first time, theoretical predictions with resistivity data employing the two theories available that incorporate the effect of both electron-surface as well as electron-grain boundary scattering acting simultaneously: the theory of A.F. Mayadas and M. Shatzkes, Phys. Rev. 1 1382 (1970) (MS), and that of G. Palasantzas, Phys. Rev. B 58 9685 (1998). We eliminate adjustable parameters from the resistivity data analysis, by using as input the grain size distribution as well as the surface roughness measured with the STM on each sample. The outcome is that both theories provide a fair representation of both the temperature as well as the thickness dependence of the resistivity data, but yet there are marked differences between the resistivity predicted by these theories. In the case of the MS theory, when the average grain diameter D is significantly smaller than ℓ 0 (300) = 37 nm, the electron mean free path in the bulk at 300 K, the effect of electron

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

Science.gov (United States)

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

2005-05-01

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

Homogenized boundary conditions and resonance effects in Faraday cages

OpenAIRE

Hewett, DP; Hewitt, IJ

2016-01-01

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

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

Science.gov (United States)

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

2018-05-01

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

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

KAUST Repository

Reis, Tim

2012-01-01

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

Air Entrainment and Surface Ripples in a Turbulent Ship Hull Boundary Layer

Science.gov (United States)

Masnadi, Naeem; Erinin, Martin; Duncan, James H.

2017-11-01

The air entrainment and free-surface fluctuations caused by the interaction of a free surface and the turbulent boundary layer of a vertical surface-piercing plate is studied experimentally. In this experiment, a meter-wide stainless steel belt travels horizontally in a loop around two rollers with vertically oriented axes. This belt device is mounted inside a large water tank with the water level set just below the top edge of the belt. The belt, rollers, and supporting frame are contained within a sheet metal box to keep the device dry except for one 6-meter-long straight test section. The belt is accelerated suddenly from rest until reaching constant speed in order to create a temporally evolving boundary layer analogous to the spatially evolving boundary layer that would exist along a surface-piercing towed flat plate. Surface ripples are measured using a cinematic laser-induced fluorescence technique with the laser sheet oriented parallel or normal to the belt surface. Air entrainment events and bubble motions are recorded from underneath the water surface using a stereo imaging system. Measurements of small bubbles, that tend to stay submerged for a longer time, are planned via a high-speed digital in-line holographic system. The support of the Office of Naval Research is gratefully acknowledged.

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)

Influence of Last Glacial Maximum boundary conditions on the global water isotope distribution in an atmospheric general circulation model

Directory of Open Access Journals (Sweden)

T. Tharammal

2013-03-01

Full Text Available To understand the validity of δ18O proxy records as indicators of past temperature change, a series of experiments was conducted using an atmospheric general circulation model fitted with water isotope tracers (Community Atmosphere Model version 3.0, IsoCAM. A pre-industrial simulation was performed as the control experiment, as well as a simulation with all the boundary conditions set to Last Glacial Maximum (LGM values. Results from the pre-industrial and LGM simulations were compared to experiments in which the influence of individual boundary conditions (greenhouse gases, ice sheet albedo and topography, sea surface temperature (SST, and orbital parameters were changed each at a time to assess their individual impact. The experiments were designed in order to analyze the spatial variations of the oxygen isotopic composition of precipitation (δ18Oprecip in response to individual climate factors. The change in topography (due to the change in land ice cover played a significant role in reducing the surface temperature and δ18Oprecip over North America. Exposed shelf areas and the ice sheet albedo reduced the Northern Hemisphere surface temperature and δ18Oprecip further. A global mean cooling of 4.1 °C was simulated with combined LGM boundary conditions compared to the control simulation, which was in agreement with previous experiments using the fully coupled Community Climate System Model (CCSM3. Large reductions in δ18Oprecip over the LGM ice sheets were strongly linked to the temperature decrease over them. The SST and ice sheet topography changes were responsible for most of the changes in the climate and hence the δ18Oprecip distribution among the simulations.

Improved low frequency room responses by considering finiteness of room boundary surfaces

DEFF Research Database (Denmark)

Jeong, Cheol-Ho

2013-01-01

surface impedance values that are assigned to all the boundary surfaces, the suggested reflection coefficient is found to improve low frequency responses compared to the infinite panel theory; larger improvements are found for a more disproportionate room, more absorptive surfaces, and surfaces having...

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 D(D3)-anyon chain: integrable boundary conditions and excitation spectra

International Nuclear Information System (INIS)

Finch, Peter E; Frahm, Holger

2013-01-01

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

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

3D automatic segmentation method for retinal optical coherence tomography volume data using boundary surface enhancement

Directory of Open Access Journals (Sweden)

Yankui Sun

2016-03-01

Full Text Available With the introduction of spectral-domain optical coherence tomography (SD-OCT, much larger image datasets are routinely acquired compared to what was possible using the previous generation of time-domain OCT. Thus, there is a critical need for the development of three-dimensional (3D segmentation methods for processing these data. We present here a novel 3D automatic segmentation method for retinal OCT volume data. Briefly, to segment a boundary surface, two OCT volume datasets are obtained by using a 3D smoothing filter and a 3D differential filter. Their linear combination is then calculated to generate new volume data with an enhanced boundary surface, where pixel intensity, boundary position information, and intensity changes on both sides of the boundary surface are used simultaneously. Next, preliminary discrete boundary points are detected from the A-Scans of the volume data. Finally, surface smoothness constraints and a dynamic threshold are applied to obtain a smoothed boundary surface by correcting a small number of error points. Our method can extract retinal layer boundary surfaces sequentially with a decreasing search region of volume data. We performed automatic segmentation on eight human OCT volume datasets acquired from a commercial Spectralis OCT system, where each volume of datasets contains 97 OCT B-Scan images with a resolution of 496×512 (each B-Scan comprising 512 A-Scans containing 496 pixels; experimental results show that this method can accurately segment seven layer boundary surfaces in normal as well as some abnormal eyes.

Surface spectra of Weyl semimetals through self-adjoint extensions

Science.gov (United States)

Seradjeh, Babak; Vennettilli, Michael

2018-02-01

We apply the method of self-adjoint extensions of Hermitian operators to the low-energy, continuum Hamiltonians of Weyl semimetals in bounded geometries and derive the spectrum of the surface states on the boundary. This allows for the full characterization of boundary conditions and the surface spectra on surfaces both normal to the Weyl node separation as well as parallel to it. We show that the boundary conditions for quadratic bulk dispersions are, in general, specified by a U (2 ) matrix relating the wave function and its derivatives normal to the surface. We give a general procedure to obtain the surface spectra from these boundary conditions and derive them in specific cases of bulk dispersion. We consider the role of global symmetries in the boundary conditions and their effect on the surface spectrum. We point out several interesting features of the surface spectra for different choices of boundary conditions, such as a Mexican-hat shaped dispersion on the surface normal to Weyl node separation. We find that the existence of bound states, Fermi arcs, and the shape of their dispersion, depend on the choice of boundary conditions. This illustrates the importance of the physics at and near the boundaries in the general statement of bulk-boundary correspondence.

A Case Study of Offshore Advection of Boundary Layer Rolls over a Stably Stratified Sea Surface

Directory of Open Access Journals (Sweden)

Nina Svensson

2017-01-01

Full Text Available Streaky structures of narrow (8-9 km high wind belts have been observed from SAR images above the Baltic Sea during stably stratified conditions with offshore winds from the southern parts of Sweden. Case studies using the WRF model and in situ aircraft observations indicate that the streaks originate from boundary layer rolls generated over the convective air above Swedish mainland, also supported by visual satellite images showing the typical signature cloud streets. The simulations indicate that the rolls are advected and maintained at least 30–80 km off the coast, in agreement with the streaks observed by the SAR images. During evening when the convective conditions over land diminish, the streaky structures over the sea are still seen in the horizontal wind field; however, the vertical component is close to zero. Thus advected feature from a land surface can affect the wind field considerably for long times and over large areas in coastal regions. Although boundary layer rolls are a well-studied feature, no previous study has presented results concerning their persistence during situations with advection to a strongly stratified boundary layer. Such conditions are commonly encountered during spring in coastal regions at high latitudes.

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.

Analysis of the susceptibility in a fluid system with Neumann – plus boundary conditions

Directory of Open Access Journals (Sweden)

Djondjorov Peter

2018-01-01

Full Text Available The behaviour of the local and total susceptibilities of a fluid system bounded by different surfaces is studied in the framework of the Ginsburg-Landau Ising type model. The case of a plain geometry, Neumann-infinity boundary conditions under variations of the temperature and an external ordering field is considered. Exact analytic expressions for the order parameter, local and total susceptibilities in such a system are presented. They are used to analyse the phase behaviour of fluids confined in regions close to the bulk critical point of the respective infinite system.

Electron scattering at surfaces and grain boundaries in thin Au films

Energy Technology Data Exchange (ETDEWEB)

Henriquez, Ricardo [Departamento de Física, Universidad Técnica Federico Santa María, Av. España 1680, Casilla 110-V, Valparaíso (Chile); Flores, Marcos; Moraga, Luis [Departamento de Física, Facultad de Ciencias Físicas y Matemáticas, Universidad de Chile, Blanco Encalada 2008, Casilla 487-3, Santiago 8370449 (Chile); Kremer, German [Bachillerato, Universidad de Chile, Las Palmeras 3425, Santiago 7800024 (Chile); González-Fuentes, Claudio [Departamento de Física, Universidad Técnica Federico Santa María, Av. España 1680, Casilla 110-V, Valparaíso (Chile); Munoz, Raul C., E-mail: ramunoz@ing.uchile.cl [Departamento de Física, Facultad de Ciencias Físicas y Matemáticas, Universidad de Chile, Blanco Encalada 2008, Casilla 487-3, Santiago 8370449 (Chile)

2013-05-15

The electron scattering at surfaces and grain boundaries is investigated using polycrystalline Au films deposited onto mica substrates. We vary the three length scales associated with: (i) electron scattering in the bulk, that at temperature T is characterized by the electronic mean free path in the bulk ℓ{sub 0}(T); (ii) electron-surface scattering, that is characterized by the film thickness t; (iii) electron-grain boundary scattering, that is characterized by the mean grain diameter D. We varied independently the film thickness from approximately 50 nm to about 100 nm, and the typical grain size making up the samples from 12 nm to 160 nm. We also varied the scale of length associated with electron scattering in the bulk by measuring the resistivity of each specimen at temperatures T, 4 K < T < 300 K. Cooling the samples to 4 K increases ℓ{sub 0}(T) by approximately 2 orders of magnitude. Detailed measurements of the grain size distribution as well as surface roughness of each sample were performed with a Scanning Tunnelling Microscope (STM). We compare, for the first time, theoretical predictions with resistivity data employing the two theories available that incorporate the effect of both electron-surface as well as electron-grain boundary scattering acting simultaneously: the theory of A.F. Mayadas and M. Shatzkes, Phys. Rev. 1 1382 (1970) (MS), and that of G. Palasantzas, Phys. Rev. B 58 9685 (1998). We eliminate adjustable parameters from the resistivity data analysis, by using as input the grain size distribution as well as the surface roughness measured with the STM on each sample. The outcome is that both theories provide a fair representation of both the temperature as well as the thickness dependence of the resistivity data, but yet there are marked differences between the resistivity predicted by these theories. In the case of the MS theory, when the average grain diameter D is significantly smaller than ℓ{sub 0}(300) = 37 nm, the electron mean

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

Nanoscale roughness effect on Maxwell-like boundary conditions for the Boltzmann equation

Energy Technology Data Exchange (ETDEWEB)

Brull, S., E-mail: Stephane.Brull@math.u-bordeaux.fr; Charrier, P., E-mail: Pierre.Charrier@math.u-bordeaux.fr; Mieussens, L., E-mail: Luc.Mieussens@math.u-bordeaux.fr [University of Bordeaux, CNRS, Bordeaux INP, IMB, UMR 5251, F-33400 Talence (France)

2016-08-15

It is well known that the roughness of the wall has an effect on microscale gas flows. This effect can be shown for large Knudsen numbers by using a numerical solution of the Boltzmann equation. However, when the wall is rough at a nanometric scale, it is necessary to use a very small mesh size which is much too expansive. An alternative approach is to incorporate the roughness effect in the scattering kernel of the boundary condition, such as the Maxwell-like kernel introduced by the authors in a previous paper. Here, we explain how this boundary condition can be implemented in a discrete velocity approximation of the Boltzmann equation. Moreover, the influence of the roughness is shown by computing the structure scattering pattern of mono-energetic beams of the incident gas molecules. The effect of the angle of incidence of these molecules, of their mass, and of the morphology of the wall is investigated and discussed in a simplified two-dimensional configuration. The effect of the azimuthal angle of the incident beams is shown for a three-dimensional configuration. Finally, the case of non-elastic scattering is considered. All these results suggest that our approach is a promising way to incorporate enough physics of gas-surface interaction, at a reasonable computing cost, to improve kinetic simulations of micro- and nano-flows.

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

Science.gov (United States)

Wang, Mengjie; Herdeiro, Carlos; Jing, Jiliang

2017-11-01

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

Quantum communication through a spin ring with twisted boundary conditions

International Nuclear Information System (INIS)

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

2005-01-01

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

The thermodynamic Casimir effect with symmetry-preserving and symmetry-breaking boundary conditions; Der thermodynamische Casimir-Effekt mit symmetrieerhaltenden und symmetriebrechenden Randbedingungen

Energy Technology Data Exchange (ETDEWEB)

Schmidt, Felix

2014-07-16

When macroscopic bodies are immersed in fluctuating media, long-range forces between these bodies may occur. The fluctuation's spectrum is modified resulting in a dependence of the system's energy on the separation between the objects, straightforwardly leading to the existence of a force between the bodies. This work is dedicated to the analysis of how boundary conditions affect the thermodynamic Casimir effect where thermal fluctuations near a critical point induce these forces. O(n) symmetric φ 4 theories in d-dimensional slab geometries of thickness L are considered. When symmetry breaking external fields are present as well, the generic boundary conditions of these theories read ∂{sub n}φ-c{sub j}φ=-h{sub j} where the coefficients c{sub j} are surface couplings, serving as linearly extrapolated penetration depths into the surfaces in Landau theory, and h{sub j} are surface fields. The influence of the surface couplings c{sub j} on the Casimir force is investigated by means of the renormalization-group-improved perturbation theory in d=4-ε dimensions to two-loop order at the bulk critical point. Special attention is paid to the case of critical enhancement of the surface interactions which results in the existence of a zero mode leading to a breakdown of the usual loop expansion of the free energy and implicating the emergence of non-integer powers of ε in the ε expansion. These perturbative methods are restricted to the disordered phase with T≥T{sub c,∞}, c{sub j}≥c{sub sp}, and h{sub j}=0. In order to extend the analysis to the whole temperature axis, the exactly treatable limit n → ∞ of the three-dimensional φ 4 model is investigated. A set of self-consistent equations for the free energy is derived that can be solved numerically exact. Considering Dirichlet boundary conditions and vanishing external fields, one finds a temperature dependence of the Casimir force that exhibits the qualitative features of the experimentally

A high-order boundary integral method for surface diffusions on elastically stressed axisymmetric rods

OpenAIRE

Li, Xiaofan; Nie, Qing

2009-01-01

Many applications in materials involve surface diffusion of elastically stressed solids. Study of singularity formation and long-time behavior of such solid surfaces requires accurate simulations in both space and time. Here we present a high-order boundary integral method for an elastically stressed solid with axi-symmetry due to surface diffusions. In this method, the boundary integrals for isotropic elasticity in axi-symmetric geometry are approximated through modified alternating quadratu...

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.

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.

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.

Groundwater infiltration, surface water inflow and sewerage exfiltration considering hydrodynamic conditions in sewer systems.

Science.gov (United States)

Karpf, Christian; Hoeft, Stefan; Scheffer, Claudia; Fuchs, Lothar; Krebs, Peter

2011-01-01

Sewer systems are closely interlinked with groundwater and surface water. Due to leaks and regular openings in the sewer system (e.g. combined sewer overflow structures with sometimes reverse pressure conditions), groundwater infiltration and surface water inflow as well as exfiltration of sewage take place and cannot be avoided. In the paper a new hydrodynamic sewer network modelling approach will be presented, which includes--besides precipitation--hydrographs of groundwater and surface water as essential boundary conditions. The concept of the modelling approach and the models to describe the infiltration, inflow and exfiltration fluxes are described. The model application to the sewerage system of the City of Dresden during a flood event with complex conditions shows that the processes of infiltration, exfiltration and surface water inflows can be described with a higher reliability and accuracy, showing that surface water inflow causes a pronounced system reaction. Further, according to the simulation results, a high sensitivity of exfiltration rates on the in-sewer water levels and a relatively low influence of the dynamic conditions on the infiltration rates were found.

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

OpenAIRE

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

2006-01-01

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

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

KAUST Repository

Leise, Tanya

2008-11-01

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

Matrix factorisations for rational boundary conditions by defect fusion

International Nuclear Information System (INIS)

Behr, Nicolas; Fredenhagen, Stefan

2015-01-01

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

Matrix factorisations for rational boundary conditions by defect fusion

Energy Technology Data Exchange (ETDEWEB)

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

2015-05-11

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

Wind-driven rain as a boundary condition for HAM simulations: analysis of simplified modelling approaches

DEFF Research Database (Denmark)

Janssen, Hans; Blocken, Bert; Roels, Staf

2007-01-01

While the numerical simulation of moisture transfer inside building components is currently undergoing standardisation, the modelling of the atmospheric boundary conditions has received far less attention. This article analyses the modelling of the wind-driven-rain load on building facades...... though: the full variability with the perpendicular wind speed and horizontal rain intensity should be preserved, where feasible, for improved estimations of the moisture transfer in building components. In the concluding section, it is moreover shown that the dependence of the surface moisture transfer...

Physical modeling of emergency emission in the atmosphere (experimental investigation of Lagrangian turbulence characteristics in the surface and boundary layer of the atmosphere)

International Nuclear Information System (INIS)

Garger, E.K.

2013-01-01

Results of diffusion experiments simulating emergency emission in the surface and boundary layers of the atmosphere are presented. Interpretation of measurements in the surface layer of the atmosphere had been conducted on the basis of the Lagrangian similarity hypothesis., Results of measurements in the boundary layer of the atmosphere are interpreted with use of the homogeneous turbulence theory. Regimes of turbulent diffusion from land and low sources of admixtures predicted by the Lagrangian similarity hypothesis for various conditions of thermal stratification in the surface layer of the atmosphere are experimentally confirmed. Universal empirical constants for these regimes are received that allows to use their in practice. Calculation diffusion parameters and concentrations of an admixture from various sources in the surface layer of the atmosphere by model is presented. Results of calculation on this model are compared to independent measurements of mass concentration of a admixture in horizontal and vertical planes. Results of simultaneous measurements Eulerian and Lagrangian turbulence characteristics for various diffusion times in the boundary layer of the atmosphere have allowed to estimate turbulence time scales in Lagrangian variables for conditions close to neutral thermal stratification. The monograph is intended for scientists and students engaged in the field of meteorology, physics of the atmosphere and pollution air control, services of radiation and ecological safety

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

Natural frequency analysis of fluid conveying pipeline with different boundary conditions

International Nuclear Information System (INIS)

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

2010-01-01

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

Boundary Conditions, 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...

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

Czech Academy of Sciences Publication Activity Database

Kružík, Martin; Zimmer, J.

2012-01-01

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

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

Science.gov (United States)

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

2012-11-01

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

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.

Analysis of the Slab Temperature, Thermal Stresses and Fractures Computed with the Implementation of Local and Average Boundary Conditions in the Secondary Cooling Zones

Directory of Open Access Journals (Sweden)

Hadała B.

2016-12-01

Full Text Available The numerical simulations of the temperature fields have been accomplished for slab casting made of a low carbon steel. The casting process of slab of 1500 mm in width and 225 mm in height has been modeled. Two types of boundary condition models of heat transfer have been employed in numerical simulations. The heat transfer coefficient in the first boundary condition model was calculated from the formula which takes into account the slab surface temperature and water flow rate in each secondary cooling zone. The second boundary condition model defines the heat transfer coefficient around each water spray nozzle. The temperature fields resulting from the average in zones water flow rate and from the nozzles arrangement have been compared. The thermal stresses and deformations resulted from such temperature field have given higher values of fracture criterion at slab corners.

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.

Prediction of the Stress-Strain Behavior of Open-Cell Aluminum Foam under Compressive Loading and the Effects of Various RVE Boundary Conditions

Science.gov (United States)

Hamidi Ghaleh Jigh, Behrang; Farsi, Mohammad Ali; Hosseini Toudeshky, Hossein

2018-04-01

The prediction of the mechanical behavior of metallic foams with realistic microstructure and the effects of various boundary conditions on the mechanical behavior is an important and challenging issue in modeling representative volume elements (RVEs). A numerical investigation is conducted to determine the effects of various boundary conditions and cell wall cross sections on the compressive mechanical properties of aluminum foam, including the stiffness, plateau stress and onset strain of densification. The open-cell AA6101-T6 aluminum foam Duocel is used in the analyses in this study. Geometrical characteristics including the cell size, foam relative density, and cross-sectional shape and thickness of the cell walls are extracted from images of the foam. Then, the obtained foam microstructure is analyzed as a 2D model. The ligaments are modeled as shear deformable beams with elastic-plastic material behavior. To prevent interpenetration of the nodes and walls inside the cells with large deformations, self-contact-type frictionless interaction is stipulated between the internal surfaces. Sensitivity analyses are performed using several boundary conditions and cells wall cross-sectional shapes. The predicted results from the finite element analyses are compared with the experimental results. Finally, the most appropriate boundary conditions, leading to more consistent results with the experimental data, are introduced.

Prediction of the Stress-Strain Behavior of Open-Cell Aluminum Foam under Compressive Loading and the Effects of Various RVE Boundary Conditions

Science.gov (United States)

Hamidi Ghaleh Jigh, Behrang; Farsi, Mohammad Ali; Hosseini Toudeshky, Hossein

2018-05-01

The prediction of the mechanical behavior of metallic foams with realistic microstructure and the effects of various boundary conditions on the mechanical behavior is an important and challenging issue in modeling representative volume elements (RVEs). A numerical investigation is conducted to determine the effects of various boundary conditions and cell wall cross sections on the compressive mechanical properties of aluminum foam, including the stiffness, plateau stress and onset strain of densification. The open-cell AA6101-T6 aluminum foam Duocel is used in the analyses in this study. Geometrical characteristics including the cell size, foam relative density, and cross-sectional shape and thickness of the cell walls are extracted from images of the foam. Then, the obtained foam microstructure is analyzed as a 2D model. The ligaments are modeled as shear deformable beams with elastic-plastic material behavior. To prevent interpenetration of the nodes and walls inside the cells with large deformations, self-contact-type frictionless interaction is stipulated between the internal surfaces. Sensitivity analyses are performed using several boundary conditions and cells wall cross-sectional shapes. The predicted results from the finite element analyses are compared with the experimental results. Finally, the most appropriate boundary conditions, leading to more consistent results with the experimental data, are introduced.

Effective Boundary Slip Induced by Surface Roughness and Their Coupled Effect on Convective Heat Transfer of Liquid Flow

Directory of Open Access Journals (Sweden)

Yunlu Pan

2018-05-01

Full Text Available As a significant interfacial property for micro/nano fluidic system, the effective boundary slip can be induced by the surface roughness. However, the effect of surface roughness on the effective slip is still not clear, both increased and decreased effective boundary slip were found with increased roughness. The present work develops a simplified model to study the effect of surface roughness on the effective boundary slip. In the created rough models, the reference position of the rough surfaces to determinate effective boundary slip was set based on ISO/ASME standard and the surface roughness parameters including Ra (arithmetical mean deviation of the assessed profile, Rsm (mean width of the assessed profile elements and shape of the texture varied to form different surface roughness. Then, the effective boundary slip of fluid flow through the rough surface was analyzed by using COMSOL 5.3. The results show that the effective boundary slip induced by surface roughness of fully wetted rough surface keeps negative and further decreases with increasing Ra or decreasing Rsm. Different shape of roughness texture also results in different effective slip. A simplified corrected method for the measured effective boundary slip was developed and proved to be efficient when the Rsm is no larger than 200 nm. Another important finding in the present work is that the convective heat transfer firstly increases followed by an unobvious change with increasing Ra, while the effective boundary slip keeps decreasing. It is believed that the increasing Ra enlarges the area of solid-liquid interface for convective heat transfer, however, when Ra is large enough, the decreasing roughness-induced effective boundary slip counteracts the enhancement effect of roughness itself on the convective heat transfer.

Spheroidal and toroidal configurations as sources of the Kerr metric. Pt. 1. A kinematical approach. [Boyer surface condition, singularity

Energy Technology Data Exchange (ETDEWEB)

de Felice, F; Nobili, L [Padua Univ. (Italy). Istituto di Fisica; Calvani, M [Padua univ. (Italy). Istituto di Astronomia

1976-03-01

The existence of extended Kerr metric sources of perfect fluid is taken as work-hypothesis to investigate the structure of the boundaries which derive from Boyer's surface condition. We find closed spheroidal configurations which hide an internal cavity as well as toroidal configurations; however, both the boundary of the internal cavity and the tori touch the ring singularity. We judge this feature non-physical and conclude that Boyer's condition is not sufficient to completely define a well behaved physical source.

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

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

Science.gov (United States)

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

2015-07-01

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

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.

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

Directory of Open Access Journals (Sweden)

Tanmoy Chakraborty

2018-03-01

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

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.

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.

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.

Free surface simulation of a two-layer fluid by boundary element method

Directory of Open Access Journals (Sweden)

Weoncheol Koo

2010-09-01

Full Text Available A two-layer fluid with free surface is simulated in the time domain by a two-dimensional potential-based Numerical Wave Tank (NWT. The developed NWT is based on the boundary element method and a leap-frog time integration scheme. A whole domain scheme including interaction terms between two layers is applied to solve the boundary integral equation. The time histories of surface elevations on both fluid layers in the respective wave modes are verified with analytic results. The amplitude ratios of upper to lower elevation for various density ratios and water depths are also compared.

Stretched flow of Carreau nanofluid with convective boundary ...

Indian Academy of Sciences (India)

journal of. January 2016 physics pp. 3–17. Stretched flow of Carreau nanofluid with ... fluid over a flat plate subjected to convective surface condition. ... the steady laminar boundary layer flow over a permeable plate with a convective boundary.

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.

Analytical Solution of Heat Conduction for Hollow Cylinders with Time-Dependent Boundary Condition and Time-Dependent Heat Transfer Coefficient

Directory of Open Access Journals (Sweden)

Te-Wen Tu

2015-01-01

Full Text Available An analytical solution for the heat transfer in hollow cylinders with time-dependent boundary condition and time-dependent heat transfer coefficient at different surfaces is developed for the first time. The methodology is an extension of the shifting function method. By dividing the Biot function into a constant plus a function and introducing two specially chosen shifting functions, the system is transformed into a partial differential equation with homogenous boundary conditions only. The transformed system is thus solved by series expansion theorem. Limiting cases of the solution are studied and numerical results are compared with those in the literature. The convergence rate of the present solution is fast and the analytical solution is simple and accurate. Also, the influence of physical parameters on the temperature distribution of a hollow cylinder along the radial direction is investigated.

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.

Solvability conditions for dendritic growth in the boundary-layer model with capillary anisotropy

Science.gov (United States)

Langer, J. S.; Hong, D. C.

1986-01-01

This paper is concerned primarily with the development of an analytic approach to the theory of steady-state velocity selection in the boundary-layer model of dendritic solidification. The two-dimensional version of this model with a fourfold crystalline anisotropy alpha in the surface tension is considered. By extending a WKB method introduced in an earlier paper, the alpha dependence of the selected growth rate is determined in the limit of small alpha; and this rate is studied for large alphas in the limit in which the dimensionless undercooling approaches unity. Portions of the paper are devoted to a reinterpretation of the mathematical structure of the solvability condition in problems of this kind.

Effects of microstructures and creep conditions on the fractal dimension of grain boundary fracture in high-temperature creep of heat-resistant alloys

International Nuclear Information System (INIS)

Tanaka, Manabu

1993-01-01

The effects of microstructural aspects, such as grain size and grain boundary configuration, and creep conditions on the fractal dimension of the grain boundary fracture were examined using several heat-resistant alloys, principally in an analysis scale range between one grain boundary length and specimen size. Grain boundary fracture surface profiles in the heat-resistant alloys exhibited a fractal nature in the scale range between one grain boundary length and specimen size as well as in the scale range below one grain boundary length. The fractal dimension of the grain boundary fracture slightly increased with decreasing grain size and was generally a little larger in the specimens with serrated grain boundaries than in those with straight grain boundaries. The fractal dimension of the grain boundary and the number of grain boundary microcracks which affected the grain boundary fracture patterns were a little larger in the specimen with the smaller grain size, and were also larger in the specimen with serrated grain boundaries. The fractal dimension of the grain boundary fracture increased with decreasing creep stress in the temperature range from 973 to 1422 K in these alloys, since more grain boundary microcracks existed in the specimens ruptured under the lower stresses at the higher temperatures. (orig.) [de

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.

Impacts of synoptic condition and planetary boundary layer structure on the trans-boundary aerosol transport from Beijing-Tianjin-Hebei region to northeast China

Science.gov (United States)

Miao, Yucong; Guo, Jianping; Liu, Shuhua; Zhao, Chun; Li, Xiaolan; Zhang, Gen; Wei, Wei; Ma, Yanjun

2018-05-01

The northeastern China frequently experiences severe aerosol pollution in winter under unfavorable meteorological conditions. How and to what extent the meteorological factors affect the air quality there are not yet clearly understood. Thus, this study investigated the impacts of synoptic patterns on the aerosol transport and planetary boundary layer (PBL) structure in Shenyang from 1 to 3 December 2016, using surface observations, sounding measurements, satellite data, and three-dimensional simulations. Results showed that the aerosol pollution occurred in Shenyang was not only related to the local emissions, but also contributed by trans-boundary transport of aerosols from the Beiijng-Tianjin-Hebei (BTH) region. In the presence of the westerly and southwesterly synoptic winds, the aerosols emitted from BTH could be brought to Shenyang. From December 2 to 3, the aerosols emitted from BTH accounted for ∼20% of near-surface PM2.5 in Shenyang. In addition, the large-scale synoptic forcings could affect the vertical mixing of pollutants through modulating the PBL structure in Shenyang. The westerly and southwesterly synoptic winds not only brought the aerosols but also the warmer air masses from the southwest regions to Shenyang. The strong warm advections above PBL could enhance the already existing thermal inversion layers capping over PBL in Shenyang, leading to the suppressions of PBL. Both the trans-boundary transport of aerosols and the suppressions of PBL caused by the large-scale synoptic forcings should be partly responsible for the poor air quality in Shenyang, in addition to the high pollutant emissions. The present study revealed the physical mechanisms underlying the aerosol pollution in Shenyang, which has important implications for better forecasting and controlling the aerosols pollution.

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.

Numerical investigations of solute transport in bimodal porous media under dynamic boundary conditions

Science.gov (United States)

Cremer, Clemens; Neuweiler, Insa; Bechtold, Michel; Vanderborght, Jan

2016-04-01

behavior depends on the magnitude of the flow rates and hydraulic conductivity curves of the materials. Based on the unsaturated hydraulic conductivity at the intersection point of conductivity curves, we are able to define an estimate of flow rates at which the dynamic of the upper boundary condition significantly alters preferential flow paths through the system. If flow rates are low, with regard to the materials hydraulic conductivity at the intersection point, the influence of dynamic boundary conditions is small. If flow rates are in the range of the unsaturated hydraulic conductivity at intersection, solute is trapped in the fine material during upwards transport, which results in a more pronounced tailing. For flow rates exceeding the intersection conductivity, a redistribution at the soil surface can occur. References: Bechtold, M., S. Haber-Pohlmeier, J. Vanderborght, A. Pohlmeier, T.P.A. Ferré and H. Veerecken. 2011a. Near-surface solute redistribution during evaporation. Geophys. Res. Lett., 38, L17404, doi:10.1029/2011GL048147. Bechtold, M., J. Vanderborght, O. Ippisch and H. Vereecken. 2011b. Efficient random walk particle tracking algorithm for advective dispersive transport in media with discontinuous dispersion coefficients and water contents. Water Resour. Res., 47, W10526, doi: 10.1029/2010WR010267. Ippisch O., H.-J. Vogel and P. Bastian. 2006. Validity limits fort he van Genuchten-Mualem model and implications for parameter estimation and numerical simulation. Adv. Water Resour., 29, 1780-1789, doi: 10.1016/j.advwateres.2005.12.011. Lehmann, P. and D. Or. 2009. Evaporation and capillary coupling across vertical textural contrasts in porous media. Phys. Rev. E, 80, 046318, doi:10.1103/PhysRevE.80.046318.

Boundary layer for non-newtonian fluids on curved surfaces

International Nuclear Information System (INIS)

Stenger, N.

1981-04-01

By using the basic equation of fluid motion (conservation of mass and momentum) the boundary layer parameters for a Non-Newtonian, incompressible and laminar fluid flow, has been evaluated. As a test, the flat plate boundary layer is first analized and afterwards, a case with pressure gradient, allowing separation, is studied. In the case of curved surfaces, the problem is first developed in general and afterwards particularized to a circular cylinder. Finally suction and slip in the flow interface are examined. The power law model is used to represent the stress strain relationship in Non-Newtonian flow. By varying the fluid exponent one can then, have an idea of how the Non-Newtonian behavior of the flow influences the parameters of the boundary layer. Two equations, in an appropriate coordinate system have been obtained after an order of magnitude analysis of the terms in the equations of motion is performed. (Author) [pt

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.

Thalamocortical dynamics of the McCollough effect: boundary-surface alignment through perceptual learning.

Science.gov (United States)

Grossberg, Stephen; Hwang, Seungwoo; Mingolla, Ennio

2002-05-01

This article further develops the FACADE neural model of 3-D vision and figure-ground perception to quantitatively explain properties of the McCollough effect (ME). The model proposes that many ME data result from visual system mechanisms whose primary function is to adaptively align, through learning, boundary and surface representations that are positionally shifted due to the process of binocular fusion. For example, binocular boundary representations are shifted by binocular fusion relative to monocular surface representations, yet the boundaries must become positionally aligned with the surfaces to control binocular surface capture and filling-in. The model also includes perceptual reset mechanisms that use habituative transmitters in opponent processing circuits. Thus the model shows how ME data may arise from a combination of mechanisms that have a clear functional role in biological vision. Simulation results with a single set of parameters quantitatively fit data from 13 experiments that probe the nature of achromatic/chromatic and monocular/binocular interactions during induction of the ME. The model proposes how perceptual learning, opponent processing, and habituation at both monocular and binocular surface representations are involved, including early thalamocortical sites. In particular, it explains the anomalous ME utilizing these multiple processing sites. Alternative models of the ME are also summarized and compared with the present model.

Initializing a Mesoscale Boundary-Layer Model with Radiosonde Observations

Science.gov (United States)

Berri, Guillermo J.; Bertossa, Germán

2018-01-01

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

Surface influence upon vertical profiles in the nocturnal boundary layer

Science.gov (United States)

Garratt, J. R.

1983-05-01

Near-surface wind profiles in the nocturnal boundary layer, depth h, above relatively flat, tree-covered terrain are described in the context of the analysis of Garratt (1980) for the unstable atmospheric boundary layer. The observations at two sites imply a surface-based transition layer, of depth z *, within which the observed non-dimensional profiles Φ M 0 are a modified form of the inertial sub-layer relation Φ _M ( {{z L}} = ( {{{1 + 5_Z } L}} ) according to Φ _M^{{0}} ˜eq ( {{{1 + 5z} L}} )exp [ { - 0.7( {{{1 - z} z}_ * } )] , where z is height above the zero-plane displacement and L is the Monin-Obukhov length. At both sites the depth z * is significantly smaller than the appropriate neutral value ( z * N ) found from the previous analysis, as might be expected in the presence of a buoyant sink for turbulent kinetic energy.

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

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

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

Impact of Bay-Breeze Circulations on Surface Air Quality and Boundary Layer Export

Science.gov (United States)

Loughner, Christopher P.; Tzortziou, Maria; Follette-Cook, Melanie; Pickering, Kenneth E.; Goldberg, Daniel; Satam, Chinmay; Weinheimer, Andrew; Crawford, James H.; Knapp, David J.; Montzka, Denise D.;

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

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

International Nuclear Information System (INIS)

Ishak, Anuar; Nazar, Roslinda; Pop, Ioan

2008-01-01

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

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.

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.

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

Science.gov (United States)

Wouden, Alex; Cimbala, John; Lewis, Bryan

2014-11-01

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

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.

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

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.

LDV measurement of boundary layer on rotating blade surface in wind tunnel

Science.gov (United States)

Maeda, Takao; Kamada, Yasunari; Murata, Junsuke; Suzuki, Daiki; Kaga, Norimitsu; Kagisaki, Yosuke

2014-12-01

Wind turbines generate electricity due to extracting energy from the wind. The rotor aerodynamics strongly depends on the flow around blade. The surface flow on the rotating blade affects the sectional performance. The wind turbine surface flow has span-wise component due to span-wise change of airfoil section, chord length, twisted angle of blade and centrifugal force on the flow. These span-wise flow changes the boundary layer on the rotating blade and the sectional performance. Hence, the thorough understanding of blade surface flow is important to improve the rotor performance. For the purpose of clarification of the flow behaviour around the rotor blade, the velocity in the boundary layer on rotating blade surface of an experimental HAWT was measured in a wind tunnel. The velocity measurement on the blade surface was carried out by a laser Doppler velocimeter (LDV). As the results of the measurement, characteristics of surface flow are clarified. In optimum tip speed operation, the surface flow on leading edge and r/R=0.3 have large span-wise velocity which reaches 20% of sectional inflow velocity. The surface flow inboard have three dimensional flow patterns. On the other hand, the flow outboard is almost two dimensional in cross sectional plane.

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.

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

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

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.

The WELSONS experiment: overview and presentation of first results on the surface atmospheric boundary-layer in semiarid Spain

Directory of Open Access Journals (Sweden)

J.-P. Frangi

Full Text Available This study presents the preliminary results of the local energy budget and dynamic characteristics of the surface atmospheric boundary-layer (SBL during the WELSONS (wind erosion and losses of soil nutrients in semiarid Spain experiment. Some Mediterranean regions suffer land degradation by wind erosion as a consequence of their particular soil and climate conditions and inappropriate agricultural practice. In Spain, where land degradation by water erosion is well known, the lack of field studies to quantify soils losses by wind erosion resulted in the European Community organizing a scientific program for this specific issue. The European programme known as WELSONS was devoted to study the wind erosion process in central Aragon (NE Spain. This multidisciplinary experiment, which began in 1996 and finished in 1998, was carried out over an agricultural soil which was left fallow. Within the experimental field, two plots were delimited where two tillage treatments were applied, a mould-board ploughing (or conventional tillage denoted CT and chisel ploughing (reduced tillage denoted RT. This was to study on bare soil the influence of tillage method on surface conditions, saltation flux, vertical dust flux, erosion rates, dynamics characteristics such as friction velocity, roughness length, etc., and energy budget. The partitioning of the available energy, resulting from the dynamics of the SBL, are quite different over the two plots because of their own peculiar soil and surface properties. The first results show that the RT treatment seems to provide a wind erosion protection. Because of the long data recording time and particular phenomena (formation of a crust at the soil surface, very dry conditions, high wind speed for instance, these microclimatological data acquired during the WELSONS programmes may be helpful to test atmospheric boundary-layer models coupled with soil models.

Key words: Hydrology (desertification - Meterology and

The WELSONS experiment: overview and presentation of first results on the surface atmospheric boundary-layer in semiarid Spain

Directory of Open Access Journals (Sweden)

J.-P. Frangi

2000-03-01

Full Text Available This study presents the preliminary results of the local energy budget and dynamic characteristics of the surface atmospheric boundary-layer (SBL during the WELSONS (wind erosion and losses of soil nutrients in semiarid Spain experiment. Some Mediterranean regions suffer land degradation by wind erosion as a consequence of their particular soil and climate conditions and inappropriate agricultural practice. In Spain, where land degradation by water erosion is well known, the lack of field studies to quantify soils losses by wind erosion resulted in the European Community organizing a scientific program for this specific issue. The European programme known as WELSONS was devoted to study the wind erosion process in central Aragon (NE Spain. This multidisciplinary experiment, which began in 1996 and finished in 1998, was carried out over an agricultural soil which was left fallow. Within the experimental field, two plots were delimited where two tillage treatments were applied, a mould-board ploughing (or conventional tillage denoted CT and chisel ploughing (reduced tillage denoted RT. This was to study on bare soil the influence of tillage method on surface conditions, saltation flux, vertical dust flux, erosion rates, dynamics characteristics such as friction velocity, roughness length, etc., and energy budget. The partitioning of the available energy, resulting from the dynamics of the SBL, are quite different over the two plots because of their own peculiar soil and surface properties. The first results show that the RT treatment seems to provide a wind erosion protection. Because of the long data recording time and particular phenomena (formation of a crust at the soil surface, very dry conditions, high wind speed for instance, these microclimatological data acquired during the WELSONS programmes may be helpful to test atmospheric boundary-layer models coupled with soil models.Key words: Hydrology (desertification - Meterology and atmospheric

A three-dimensional, iterative mapping procedure for the implementation of an ionosphere-magnetosphere anisotropic Ohm's law boundary condition in global magnetohydrodynamic simulations

Directory of Open Access Journals (Sweden)

M. L. Goodman

1995-08-01

Full Text Available The mathematical formulation of an iterative procedure for the numerical implementation of an ionosphere-magnetosphere (IM anisotropic Ohm's law boundary condition is presented. The procedure may be used in global magnetohydrodynamic (MHD simulations of the magnetosphere. The basic form of the boundary condition is well known, but a well-defined, simple, explicit method for implementing it in an MHD code has not been presented previously. The boundary condition relates the ionospheric electric field to the magnetic field-aligned current density driven through the ionosphere by the magnetospheric convection electric field, which is orthogonal to the magnetic field B, and maps down into the ionosphere along equipotential magnetic field lines. The source of this electric field is the flow of the solar wind orthogonal to B. The electric field and current density in the ionosphere are connected through an anisotropic conductivity tensor which involves the Hall, Pedersen, and parallel conductivities. Only the height-integrated Hall and Pedersen conductivities (conductances appear in the final form of the boundary condition, and are assumed to be known functions of position on the spherical surface R=R1 representing the boundary between the ionosphere and magnetosphere. The implementation presented consists of an iterative mapping of the electrostatic potential ψ the gradient of which gives the electric field, and the field-aligned current density between the IM boundary at R=R1 and the inner boundary of an MHD code which is taken to be at R2>R1. Given the field-aligned current density on R=R2, as computed by the MHD simulation, it is mapped down to R=R1 where it is used to compute ψ by solving the equation that is the IM Ohm's law boundary condition. Then ψ is mapped out to R=R2, where it is used to update the electric field and the component of velocity perpendicular to B. The updated electric field and perpendicular velocity serve as new boundary

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

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

Tile Surface Thermocouple Measurement Challenges from the Orbiter Boundary Layer Transition Flight Experiment

Science.gov (United States)

Campbell, Charles H.; Berger, Karen; Anderson, Brian

2012-01-01

Hypersonic entry flight testing motivated by efforts seeking to characterize boundary layer transition on the Space Shuttle Orbiters have identified challenges in our ability to acquire high quality quantitative surface temperature measurements versus time. Five missions near the end of the Space Shuttle Program implemented a tile surface protuberance as a boundary layer trip together with tile surface thermocouples to capture temperature measurements during entry. Similar engineering implementations of these measurements on Discovery and Endeavor demonstrated unexpected measurement voltage response during the high heating portion of the entry trajectory. An assessment has been performed to characterize possible causes of the issues experienced during STS-119, STS-128, STS-131, STS-133 and STS-134 as well as similar issues encountered during other orbiter entries.

A Multigrid Algorithm for an Elliptic Problem with a Perturbed Boundary Condition

KAUST Repository

Bonito, Andrea; Pasciak, Joseph E.

2013-01-01

We discuss the preconditioning of systems coupling elliptic operators in Ω⊂Rd, d=2,3, with elliptic operators defined on hypersurfaces. These systems arise naturally when physical phenomena are affected by geometric boundary forces, such as the evolution of liquid drops subject to surface tension. The resulting operators are sums of interior and boundary terms weighted by parameters. We investigate the behavior of multigrid algorithms suited to this context and demonstrate numerical results which suggest uniform preconditioning bounds that are level and parameter independent.

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.

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)

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

Critical deflagration waves leading to detonation onset under different boundary conditions

International Nuclear Information System (INIS)

Lin Wei; Zhou Jin; Lin Zhi-Yong; Fan Xiao-Hua

2015-01-01

High-speed turbulent critical deflagration waves before detonation onset in H 2 –air mixture propagated into a square cross section channel, which was assembled of optional rigid rough, rigid smooth, or flexible walls. The corresponding propagation characteristic and the influence of the wall boundaries on the propagation were investigated via high-speed shadowgraph and a high-frequency pressure sampling system. As a comprehensive supplement to the different walls effect investigation, the effect of porous absorbing walls on the detonation propagation was also investigated via smoke foils and the high-frequency pressure sampling system. Results are as follows. In the critical deflagration stage, the leading shock and the closely following turbulent flame front travel at a speed of nearly half the CJ detonation velocity. In the preheated zone, a zonary flame arises from the overlapping part of the boundary layer and the pressure waves, and then merges into the mainstream flame. Among these wall boundary conditions, the rigid rough wall plays a most positive role in the formation of the zonary flame and thus accelerates the transition of the deflagration to detonation (DDT), which is due to the boost of the boundary layer growth and the pressure wave reflection. Even though the flexible wall is not conducive to the pressure wave reflection, it brings out a faster boundary layer growth, which plays a more significant role in the zonary flame formation. Additionally, the porous absorbing wall absorbs the transverse wave and yields detonation decay and velocity deficit. After the absorbing wall, below some low initial pressure conditions, no re-initiation occurs and the deflagration propagates in critical deflagration for a relatively long distance. (paper)

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

Surface condition effects on tritium permeation through the first wall of a water-cooled ceramic breeder blanket

Energy Technology Data Exchange (ETDEWEB)

Zhou, H.-S. [Institute of Plasma Physics, Chinese Academy of Sciences, P.O. Box 1126, Hefei (China); Xu, Y.-P.; Liu, H.-D. [Science Island Branch of Graduate School, University of Science and Technology of China, P.O. Box 1126, Hefei (China); Liu, F.; Li, X.-C.; Zhao, M.-Z.; Qi, Q.; Ding, F. [Institute of Plasma Physics, Chinese Academy of Sciences, P.O. Box 1126, Hefei (China); Luo, G.-N., E-mail: gnluo@ipp.ac.cn [Institute of Plasma Physics, Chinese Academy of Sciences, P.O. Box 1126, Hefei (China); Science Island Branch of Graduate School, University of Science and Technology of China, P.O. Box 1126, Hefei (China); Hefei Center for Physical Science and Technology, P.O. Box 1126, Hefei (China); Hefei Science Center of Chinese Academy of Science, P.O. Box 1126, Hefei (China)

2016-11-01

Highlights: • We investigate surface effects on T transport through the first wall. • We solve transport equations with various surface conditions. • The RAFMs walls w/and w/o W exhibit different T permeation behavior. • Diffusion in W has been found to be the rate-limiting step. - Abstract: Plasma-driven permeation of tritium (T) through the first wall of a water-cooled ceramic breeder (WCCB) blanket may raise safety and other issues. In the present work, surface effects on T transport through the first wall of a WCCB blanket have been investigated by theoretical calculation. Two types of wall structures, i.e., reduced activation ferritic/martensitic steels (RAFMs) walls with and without tungsten (W) armor, have been analyzed. Surface recombination is assumed to be the boundary condition for both the plasma-facing side and the coolant side. It has been found that surface conditions at both sides can affect T permeation flux and inventory. For the first wall using W as armor material, T permeation is not sensitive to the plasma-facing surface conditions. Contamination of the surfaces will lead to higher T inventory inside the first wall.

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.

Altered vocal fold kinematics in synthetic self-oscillating models that employ adipose tissue as a lateral boundary condition.

Science.gov (United States)

Saidi, Hiba; Erath, Byron D.

2015-11-01

The vocal folds play a major role in human communication by initiating voiced sound production. During voiced speech, the vocal folds are set into sustained vibrations. Synthetic self-oscillating vocal fold models are regularly employed to gain insight into flow-structure interactions governing the phonation process. Commonly, a fixed boundary condition is applied to the lateral, anterior, and posterior sides of the synthetic vocal fold models. However, physiological observations reveal the presence of adipose tissue on the lateral surface between the thyroid cartilage and the vocal folds. The goal of this study is to investigate the influence of including this substrate layer of adipose tissue on the dynamics of phonation. For a more realistic representation of the human vocal folds, synthetic multi-layer vocal fold models have been fabricated and tested while including a soft lateral layer representative of adipose tissue. Phonation parameters have been collected and are compared to those of the standard vocal fold models. Results show that vocal fold kinematics are affected by adding the adipose tissue layer as a new boundary condition.

Regional air-quality forecasting for the Pacific Northwest using MOPITT/TERRA assimilated carbon monoxide MOZART-4 forecasts as a near real-time boundary condition

Directory of Open Access Journals (Sweden)

F. L. Herron-Thorpe

2012-06-01

Full Text Available Results from a regional air quality forecast model, AIRPACT-3, were compared to AIRS carbon monoxide column densities for the spring of 2010 over the Pacific Northwest. AIRPACT-3 column densities showed high correlation (R > 0.9 but were significantly biased (~25% with consistent under-predictions for spring months when there is significant transport from Asia. The AIRPACT-3 CO bias relative to AIRS was eliminated by incorporating dynamic boundary conditions derived from NCAR's MOZART forecasts with assimilated MOPITT carbon monoxide. Changes in ozone-related boundary conditions derived from MOZART forecasts are also discussed and found to affect background levels by ± 10 ppb but not found to significantly affect peak ozone surface concentrations.

High Reynolds number rough wall turbulent boundary layer experiments using Braille surfaces

Science.gov (United States)

Harris, Michael; Monty, Jason; Nova, Todd; Allen, James; Chong, Min

2007-11-01

This paper details smooth, transitional and fully rough turbulent boundary layer experiments in the New Mexico State high Reynolds number rough wall wind tunnel. The initial surface tested was generated with a Braille printer and consisted of an uniform array of Braille points. The average point height being 0.5mm, the spacing between the points in the span was 0.5mm and the surface consisted of span wise rows separated by 4mm. The wavelength to peak ratio was 8:1. The boundary layer thickness at the measurement location was 190mm giving a large separation of roughness height to layer thickness. The maximum friction velocity was uτ=1.5m/s at Rex=3.8 x10^7. Results for the skin friction co-efficient show that this surface follows a Nikuradse type inflectional curve and that Townsends outer layer similarity hypothesis is valid for rough wall flows with a large separation of scales. Mean flow and turbulence statistics will be presented.

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.

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.

Spatial Manipulation of Heat Flow by Surface Boundaries at the Nanoscale

Science.gov (United States)

Malhotra, Abhinav; Maldovan, Martin

The precise manipulation of phonon transport properties is central to controlling thermal transport in semiconductor nanostructures. The physical understanding, prediction, and control of thermal phonon heat spectra and thermal conductivity accumulation functions - which establish the proportion of heat transported by phonons with different frequencies and mean-free-paths - has attracted significant attention in recent years. In this talk, we advance the possibilities of manipulating heat by spatially modulating thermal transport in nanostructures. We show that phonon scattering at interfaces impacts the most preferred physical pathway used by heat energy flow in thermal transport in nanostructures. The role of introducing boundaries with different surface conditions on resultant thermal flux is presented and methodologies to enhance these spatial modulations are discussed. This talk aims to advance the fundamental understanding on the nature of heat transport at nanoscale with potential applications in multiple research areas ranging from energy materials to optoelectronics.

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.

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.

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.

Elastic interaction between surface and spherical pore

International Nuclear Information System (INIS)

Ganeev, G.Z.; Kadyrzhanov, K.K.; Kislitsyn, S.B.; Turkebaev, T.Eh.

2000-01-01

The energy of elastic interaction of a gas-filled spherical cavity with a boundary of an elastic isotropic half-space is determined. The elastic field of a system of a spherical cavity - boundary is represented as an expansion in series of potential functions. The factors of expansions are determined by boundary conditions on a free surface of an elastic half-space and on a spherical surface of a cavity with pressure of gas P. Function of a Tresca-Miesesa on a surface of elastic surface is defined additionally with purpose creep condition determination caused by gas pressure in the cavity. (author)

Quantum effective action in spacetimes with branes and boundaries

International Nuclear Information System (INIS)

Barvinsky, A.O.; Nesterov, D.V.

2006-01-01

We construct quantum effective action in spacetime with branes/boundaries. This construction is based on the reduction of the underlying Neumann type boundary value problem for the propagator of the theory to that of the much more manageable Dirichlet problem. In its turn, this reduction follows from the recently suggested Neumann-Dirichlet duality which we extend beyond the tree-level approximation. In the one-loop approximation this duality suggests that the functional determinant of the differential operator subject to Neumann boundary conditions factorizes into the product of its Dirichlet counterpart and the functional determinant of a special operator on the brane--the inverse of the brane-to-brane propagator. As a byproduct of this relation we suggest a new method for surface terms of the heat kernel expansion. This method allows one to circumvent well-known difficulties in the heat kernel theory on manifolds with boundaries for a wide class of generalized Neumann boundary conditions. In particular, we easily recover several lowest-order surface terms in the case of Robin and oblique boundary onditions. We briefly discuss multiloop applications of the suggested Dirichlet reduction and the prospects of constructing the universal background-field method for systems with branes/boundaries, analogous to the Schwinger-DeWitt technique

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

Experimental and numerical study on inlet and outlet conditions of a bulb turbine with considering free surface

International Nuclear Information System (INIS)

Zhao, Y P; Liao, W L; Feng, H D; Ruan, H; Luo, X Q

2012-01-01

For a bulb turbine, it has a low head and a big runner diameter, and the free surface influences the flow at the inlet and outlet of the turbine, which bring many problems such as vibration, cracks and cavitation to the turbine. Therefore, it is difficult to get the precise internal flow characteristics through a numerical simulation with conventional ideal flow conditions. In this paper, both numerical and experimental methods are adopted to investigate the flow characteristics at the inlet and outlet of the bulb turbine with considering free surface. Firstly, experimental and numerical studies in a low head pressure pipeline are conducted, and the corresponding boundary condition according with reality is obtained through the comparison between the model test result and the CFD simulation result. Then, through an analysis of the velocity and pressure fields at the inlet of the bulb turbine at different heads, the flow characteristics and rules at the entrance of the bulb turbine have been revealed with considering free surface; Finally, the performance predictions for a bulb turbine have been conducted by using the obtained flow rules at the inlet as the boundary condition of a turbine, and the causes that lead to non-uniform forces on blades, cavitation and vibration have been illustrated in this paper, which also provide a theory basis for an accurate numerical simulation and optimization design of a bulb turbine.

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

Effect of residual stress and hardening on grain boundary sliding in welds of low-carbon stainless steels with surface machining

International Nuclear Information System (INIS)

Mori, Hiroaki; Mochizuki, Masahito; Nishimoto, Kazutoshi; Katsuyama, Jinya

2007-01-01

To clarify the effects of residual stress and hardening on intergranular stress corrosion cracking (IGSCC) behavior in welds of low-carbon austenitic stainless steels with surface machining, residual stress and hardness were evaluated by 3-dimentional thermo elastic-plastic analysis and grain boundary sliding behavior was examined using a constant strain rate tensile test. It was revealed that grain boundary sliding occurred in the material at 561K by the tensile test with the numerically simulated tensile residual stress due to welding and surface machining. In addition, it was clarified that the grain boundary energy is raised by the grain boundary sliding. On the basis of these results, it was concluded that the cause of IGSCC in the welds of low-carbon austenitic stainless steel with surface hardening is the increase in grain boundary energy due to grain boundary sliding accelerated by residual stress of multi pass welding and surface hardening. (author)

Effect of residual stress and hardening on grain boundary sliding in welds of low-carbon stainless steels with surface machining

International Nuclear Information System (INIS)

Mori, Hiroaki; Mochizuki, Masahito; Nishimoto, Kazutoshi; Katsuyama, Jinya

2008-01-01

To clarify the effects of residual stress and hardening on intergranular stress corrosion cracking (IGSCC) behavior in welds of low-carbon austenitic stainless steels with surface machining, residual stress and hardness were evaluated by 3-dimentional thermo elastic-plastic analysis and grain boundary sliding behavior was examined using a constant strain rate tensile test. It was revealed that grain boundary sliding occurred in the material at 561K by the tensile test with the numerically simulated tensile residual stress due to multi-pass welding and surface machining. In addition, it was clarified that the grain boundary energy is raised by the grain boundary sliding. On the basis of these results, it was concluded that the cause of IGSCC in the welds of low-carbon austenitic stainless steel with surface hardening is the increase in grain boundary energy due to grain boundary sliding induced by residual stress of multi pass welding and surface hardening. (author)