Growth rate degeneracies in kinematic dynamos

NASA Astrophysics Data System (ADS)

Favier, B.; Proctor, M. R. E.

2013-09-01

We consider the classical problem of kinematic dynamo action in simple steady flows. Due to the adjointness of the induction operator, we show that the growth rate of the dynamo will be exactly the same for two types of magnetic boundary conditions: the magnetic field can be normal (infinite magnetic permeability, also called pseudovacuum) or tangent (perfect electrical conductor) to the boundaries of the domain. These boundary conditions correspond to well-defined physical limits often used in numerical models and relevant to laboratory experiments. The only constraint is for the velocity field u to be reversible, meaning there exists a transformation changing u into -u. We illustrate this surprising property using S2T2 type of flows in spherical geometry inspired by [Dudley and James, Proc. R. Soc. London A1364-502110.1098/rspa.1989.0112 425, 407 (1989)]. Using both types of boundary conditions, it is shown that the growth rates of the dynamos are identical, although the corresponding magnetic eigenmodes are drastically different.

Repulsive Casimir force in Bose–Einstein Condensate

NASA Astrophysics Data System (ADS)

Mehedi Faruk, Mir; Biswas, Shovon

2018-04-01

We study the Casimir effect for a three dimensional system of ideal free massive Bose gas in a slab geometry with Zaremba and anti-periodic boundary conditions. It is found that for these type of boundary conditions the resulting Casimir force is repulsive in nature, in contrast with usual periodic, Dirichlet or Neumann boundary condition where the Casimir force is attractive (Martin and Zagrebnov 2006 Europhys. Lett. 73 15). Casimir forces in these boundary conditions also maintain a power law decay function below condensation temperature and exponential decay function above the condensation temperature albeit with a positive sign, identifying the repulsive nature of the force.

A comparison of time domain boundary conditions for acoustic waves in wave guides

NASA Technical Reports Server (NTRS)

Banks, H. T.; Propst, G.; Silcox, R. J.

1991-01-01

Researchers consider several types of boundary conditions in the context of time domain models for acoustic waves. Experiments with four different duct terminations (hard wall, free radiation, foam, and wedge) were carried out in a wave duct from which reflection coefficients over a wide frequency range were measured. These reflection coefficients were used to estimate parameters in the time domain boundary conditions. A comparison of the relative merits of the models in describing the data is presented. Boundary conditions which yield a good fit of the model to the experimental data were found for all duct terminations except the wedge.

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

NASA Astrophysics Data System (ADS)

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

2017-11-01

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

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

PubMed

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

2017-11-01

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

Global and blowup solutions of a mixed problem with nonlinear boundary conditions for a one-dimensional semilinear wave equation

NASA Astrophysics Data System (ADS)

Kharibegashvili, S. S.; Jokhadze, O. M.

2014-04-01

A mixed problem for a one-dimensional semilinear wave equation with nonlinear boundary conditions is considered. Conditions of this type occur, for example, in the description of the longitudinal oscillations of a spring fastened elastically at one end, but not in accordance with Hooke's linear law. Uniqueness and existence questions are investigated for global and blowup solutions to this problem, in particular how they depend on the nature of the nonlinearities involved in the equation and the boundary conditions. Bibliography: 14 titles.

The surface-induced spatial-temporal structures in confined binary alloys

NASA Astrophysics Data System (ADS)

Krasnyuk, Igor B.; Taranets, Roman M.; Chugunova, Marina

2014-12-01

This paper examines surface-induced ordering in confined binary alloys. The hyperbolic initial boundary value problem (IBVP) is used to describe a scenario of spatiotemporal ordering in a disordered phase for concentration of one component of binary alloy and order parameter with non-linear dynamic boundary conditions. This hyperbolic model consists of two coupled second order differential equations for order parameter and concentration. It also takes into account effects of the “memory” on the ordering of atoms and their densities in the alloy. The boundary conditions characterize surface velocities of order parameter and concentration changing which is due to surface (super)cooling on walls confining the binary alloy. It is shown that for large times there are three classes of dynamic non-linear boundary conditions which lead to three different types of attractor’s elements for the IBVP. Namely, the elements of attractor are the limit periodic simple shock waves with fronts of “discontinuities” Γ. If Γ is finite, then the attractor contains spatiotemporal functions of relaxation type. If Γ is infinite and countable then we observe the functions of pre-turbulent type. If Γ is infinite and uncountable then we obtain the functions of turbulent type.

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

NASA Astrophysics Data System (ADS)

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.

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

NASA Astrophysics Data System (ADS)

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.

A characteristic based volume penalization method for general evolution problems applied to compressible viscous flows

NASA Astrophysics Data System (ADS)

Brown-Dymkoski, Eric; Kasimov, Nurlybek; Vasilyev, Oleg V.

2014-04-01

In order to introduce solid obstacles into flows, several different methods are used, including volume penalization methods which prescribe appropriate boundary conditions by applying local forcing to the constitutive equations. One well known method is Brinkman penalization, which models solid obstacles as porous media. While it has been adapted for compressible, incompressible, viscous and inviscid flows, it is limited in the types of boundary conditions that it imposes, as are most volume penalization methods. Typically, approaches are limited to Dirichlet boundary conditions. In this paper, Brinkman penalization is extended for generalized Neumann and Robin boundary conditions by introducing hyperbolic penalization terms with characteristics pointing inward on solid obstacles. This Characteristic-Based Volume Penalization (CBVP) method is a comprehensive approach to conditions on immersed boundaries, providing for homogeneous and inhomogeneous Dirichlet, Neumann, and Robin boundary conditions on hyperbolic and parabolic equations. This CBVP method can be used to impose boundary conditions for both integrated and non-integrated variables in a systematic manner that parallels the prescription of exact boundary conditions. Furthermore, the method does not depend upon a physical model, as with porous media approach for Brinkman penalization, and is therefore flexible for various physical regimes and general evolutionary equations. Here, the method is applied to scalar diffusion and to direct numerical simulation of compressible, viscous flows. With the Navier-Stokes equations, both homogeneous and inhomogeneous Neumann boundary conditions are demonstrated through external flow around an adiabatic and heated cylinder. Theoretical and numerical examination shows that the error from penalized Neumann and Robin boundary conditions can be rigorously controlled through an a priori penalization parameter η. The error on a transient boundary is found to converge as O(η), which is more favorable than the error convergence of the already established Dirichlet boundary condition.

Log corrections to entropy of three dimensional black holes with soft hair

NASA Astrophysics Data System (ADS)

Grumiller, Daniel; Perez, Alfredo; Tempo, David; Troncoso, Ricardo

2017-08-01

We calculate log corrections to the entropy of three-dimensional black holes with "soft hairy" boundary conditions. Their thermodynamics possesses some special features that preclude a naive direct evaluation of these corrections, so we follow two different approaches. The first one exploits that the BTZ black hole belongs to the spectrum of Brown-Henneaux as well as soft hairy boundary conditions, so that the respective log corrections are related through a suitable change of the thermodynamic ensemble. In the second approach the analogue of modular invariance is considered for dual theories with anisotropic scaling of Lifshitz type with dynamical exponent z at the boundary. On the gravity side such scalings arise for KdV-type boundary conditions, which provide a specific 1-parameter family of multi-trace deformations of the usual AdS3/CFT2 setup, with Brown-Henneaux corresponding to z = 1 and soft hairy boundary conditions to the limiting case z → 0+. Both approaches agree in the case of BTZ black holes for any non-negative z. Finally, for soft hairy boundary conditions we show that not only the leading term, but also the log corrections to the entropy of black flowers endowed with affine û (1) soft hair charges exclusively depend on the zero modes and hence coincide with the ones for BTZ black holes.

Well-posedness of the free boundary problem in compressible elastodynamics

NASA Astrophysics Data System (ADS)

Trakhinin, Yuri

2018-02-01

We study the free boundary problem for the flow of a compressible isentropic inviscid elastic fluid. At the free boundary moving with the velocity of the fluid particles the columns of the deformation gradient are tangent to the boundary and the pressure vanishes outside the flow domain. We prove the local-in-time existence of a unique smooth solution of the free boundary problem provided that among three columns of the deformation gradient there are two which are non-collinear vectors at each point of the initial free boundary. If this non-collinearity condition fails, the local-in-time existence is proved under the classical Rayleigh-Taylor sign condition satisfied at the first moment. By constructing an Hadamard-type ill-posedness example for the frozen coefficients linearized problem we show that the simultaneous failure of the non-collinearity condition and the Rayleigh-Taylor sign condition leads to Rayleigh-Taylor instability.

Evolution of the substructure of a novel 12% Cr steel under creep conditions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yadav, Surya Deo, E-mail: surya.yadav@tugraz.at; Kalácska, Szilvia, E-mail: kalacska@metal.elte.hu; Dománková, Mária, E-mail: maria.domankova@stuba.sk

2016-05-15

In this work we study the microstruture evolution of a newly developed 12% Cr martensitic/ferritic steel in as-received condition and after creep at 650 °C under 130 MPa and 80 MPa. The microstructure is described as consisting of mobile dislocations, dipole dislocations, boundary dislocations, precipitates, lath boundaries, block boundaries, packet boundaries and prior austenitic grain boundaries. The material is characterized employing light optical microscopy (LOM), scanning electron microscopy (SEM), transmission electron microscopy (TEM), X-ray diffraction (XRD) and electron backscatter diffraction (EBSD). TEM is used to characterize the dislocations (mobile + dipole) inside the subgrains and XRD measurements are used tomore » the characterize mobile dislocations. Based on the subgrain boundary misorientations obtained from EBSD measurements, the boundary dislocation density is estimated. The total dislocation density is estimated for the as-received and crept conditions adding the mobile, boundary and dipole dislocation densities. Additionally, the subgrain size is estimated from the EBSD measurements. In this publication we propose the use of three characterization techniques TEM, XRD and EBSD as necessary to characterize all type of dislocations and quantify the total dislocation densty in martensitic/ferritic steels. - Highlights: • Creep properties of a novel 12% Cr steel alloyed with Ta • Experimental characterization of different types of dislocations: mobile, dipole and boundary • Characterization and interpretation of the substructure evolution using unique combination of TEM, XRD and EBSD.« less

On Analysis of Stationary Viscous Incompressible Flow Through a Radial Blade Machine

NASA Astrophysics Data System (ADS)

Neustupa, Tomáš

2010-09-01

The paper is concerned with the analysis of the two dimensional model of incompressible, viscous, stationary flow through a radial blade machine. This type of turbine is sometimes called Kaplan's turbine. In the technical area the use is either to force some regular characteristic to the flow of the medium going through the turbine (flow of melted iron, air conditioning) or to gain some energy from the flowing medium (water). The inflow and outflow part of boundary are in general a concentric circles. The larger one represents an inflow part of boundary the smaller one the outflow part of boundary. Between them are regularly spaced the blades of the machine. We study the existence of the weak solution in the case of nonlinear boundary condition of the "do-nothing" type. The model is interesting for study the behavior of the flow when the boundary is formed by mutually disjoint and separated parts.

Iterative methods for plasma sheath calculations: Application to spherical probe

NASA Technical Reports Server (NTRS)

Parker, L. W.; Sullivan, E. C.

1973-01-01

The computer cost of a Poisson-Vlasov iteration procedure for the numerical solution of a steady-state collisionless plasma-sheath problem depends on: (1) the nature of the chosen iterative algorithm, (2) the position of the outer boundary of the grid, and (3) the nature of the boundary condition applied to simulate a condition at infinity (as in three-dimensional probe or satellite-wake problems). Two iterative algorithms, in conjunction with three types of boundary conditions, are analyzed theoretically and applied to the computation of current-voltage characteristics of a spherical electrostatic probe. The first algorithm was commonly used by physicists, and its computer costs depend primarily on the boundary conditions and are only slightly affected by the mesh interval. The second algorithm is not commonly used, and its costs depend primarily on the mesh interval and slightly on the boundary conditions.

Effect of boundary conditions on magnetocapacitance effect in a ring-type magnetoelectric structure

NASA Astrophysics Data System (ADS)

Zhang, Juanjuan

2017-12-01

By considering the nonlinear magneto-elastic coupling relationships of magnetostrictive materials, an analytical model is proposed. The resonance frequencies can be accurately predicted by this theoretical model, and they are in good agreement with experimental data. Subsequently, the magnetocapacitance effect in a ring-type magnetoelectric (ME) structure with different boundary conditions is investigated, and it is found that various mechanical boundaries, the frequency, the magnetic field, the geometric size, and the interface bonding significantly affect the capacitance of the ME structure. Further, additional resonance frequencies can be predicted by considering appropriate imperfect interface bonding. Finally, the influence of an external force on the capacitance is studied. The result shows that an external force on the boundary changes the capacitance, but has only a weak influence on the resonance frequency.

Nonreflective Conditions for Perfectly Matched Layer in Computational Aeroacoustics

NASA Astrophysics Data System (ADS)

Choung, Hanahchim; Jang, Seokjong; Lee, Soogab

2018-05-01

In computational aeroacoustics, boundary conditions such as radiation, outflow, or absorbing boundary conditions are critical issues in that they can affect the entire solution of the computation. Among these types of boundary conditions, the perfectly matched layer boundary condition, which has been widely used in computational fluid dynamics and computational aeroacoustics, is developed by augmenting the additional term in the original governing equations by an absorption function so as to stably absorb the outgoing waves. Even if the perfectly matched layer is analytically a perfectly nonreflective boundary condition, spurious waves occur at the interface, since the analysis is performed in discretized space. Hence, this study is focused on factors that affect numerical errors from perfectly matched layer to find the optimum conditions for nonreflective PML. Through a mathematical approach, a minimum width of perfectly matched layer and an optimum absorption coefficient are suggested. To validate the prediction of the analysis, numerical simulations are performed in a generalized coordinate system, as well as in a Cartesian coordinate system.

Integrable boundary value problems for elliptic type Toda lattice in a disk

DOE Office of Scientific and Technical Information (OSTI.GOV)

Guerses, Metin; Habibullin, Ismagil; Zheltukhin, Kostyantyn

The concept of integrable boundary value problems for soliton equations on R and R{sub +} is extended to regions enclosed by smooth curves. Classes of integrable boundary conditions in a disk for the Toda lattice and its reductions are found.

A finite-difference time-domain electromagnetic solver in a generalized coordinate system

NASA Astrophysics Data System (ADS)

Hochberg, Timothy Allen

A new, finite-difference, time-domain method for the simulation of full-wave electromagnetic wave propogation in complex structures is developed. This method is simple and flexible; it allows for the simulation of transient wave propogation in a large class of practical structures. Boundary conditions are implemented for perfect and imperfect electrically conducting boundaries, perfect magnetically conducting boundaries, and absorbing boundaries. The method is validated with the aid of several different types of test cases. Two types of coaxial cables with helical breaks are simulated and the results are discussed.

Assimilation of Wave and Current Data for Prediction of Inlet and River Mouth Dynamics

DTIC Science & Technology

2013-07-01

onto the Delft3D computational grid and the specification of Riemann -type boundary conditions for the boundary-normal velocity and surface elevation...conditions from time- history data from in situ tide gages. The corrections are applied to the surface-elevation contribution to the Riemann boundary...The algorithms described above are all of the strong-constraint variational variety, and make use of adjoint solvers corresponding to the various

Impact of interfacial imperfection on transverse wave in a functionally graded piezoelectric material structure with corrugated boundaries

NASA Astrophysics Data System (ADS)

Kumar Singh, Abhishek; Kumar, Santan; Kumari, Richa

2018-03-01

The propagation behavior of Love-type wave in a corrugated functionally graded piezoelectric material layered structure has been taken into account. Concretely, the layered structure incorporates a corrugated functionally graded piezoelectric material layer imperfectly bonded to a functionally graded piezoelectric material half-space. An analytical treatment has been employed to determine the dispersion relation for both cases of electrically open condition and electrically short condition. The phase velocity of the Love-type wave has been computed numerically and its dependence on the wave number has been depicted graphically for a specific type of corrugated boundary surfaces for both said conditions. The crux of the study lies in the fact that the imperfect bonding of the interface, the corrugated boundaries present in the layer, and the material properties of the layer and the half-space strongly influence the phase velocity of the Love-type wave. It can be remarkably noted that the imperfect bonding of the interface reduces the phase velocity of the Love-type wave significantly. As a special case of the problem, it is noticed that the procured dispersion relation for both cases of electrically open and electrically short conditions is in accordance with the classical Love wave equation.

On steady motion of viscoelastic fluid of Oldroyd type

DOE Office of Scientific and Technical Information (OSTI.GOV)

Baranovskii, E. S., E-mail: esbaranovskii@gmail.com

2014-06-01

We consider a mathematical model describing the steady motion of a viscoelastic medium of Oldroyd type under the Navier slip condition at the boundary. In the rheological relation, we use the objective regularized Jaumann derivative. We prove the solubility of the corresponding boundary-value problem in the weak setting. We obtain an estimate for the norm of a solution in terms of the data of the problem. We show that the solution set is sequentially weakly closed. Furthermore, we give an analytic solution of the boundary-value problem describing the flow of a viscoelastic fluid in a flat channel under a slipmore » condition at the walls. Bibliography: 13 titles. (paper)« less

Nonsteady Problem for an Elastic Half-Plane with Mixed Boundary Conditions

NASA Astrophysics Data System (ADS)

Kubenko, V. D.

2016-03-01

An approach to studying nonstationary wave processes in an elastic half-plane with mixed boundary conditions of the fourth boundary-value problem of elasticity is proposed. The Laplace and Fourier transforms are used. The sequential inversion of these transforms or the inversion of the joint transform by the Cagniard method allows obtaining the required solution (stresses, displacements) in a closed analytic form. With this approach, the problem can be solved for various types of loads

Ground state for a massive scalar field in the BTZ spacetime with Robin boundary conditions

NASA Astrophysics Data System (ADS)

Bussola, Francesco; Dappiaggi, Claudio; Ferreira, Hugo R. C.; Khavkine, Igor

2017-11-01

We consider a real, massive scalar field in Bañados-Teitelboim-Zanelli spacetime, a 2 +1 -dimensional black hole solution of Einstein's field equations with a negative cosmological constant. First, we analyze the space of classical solutions in a mode decomposition, and we characterize the collection of all admissible boundary conditions of Robin type which can be imposed at infinity. Second, we investigate whether, for a given boundary condition, there exists a ground state by constructing explicitly its two-point function. We demonstrate that for a subclass of the boundary conditions it is possible to construct a ground state that locally satisfies the Hadamard property. In all other cases, we show that bound state mode solutions exist and, therefore, such construction is not possible.

Experimental Investigation of Soil and Atmospheric Conditions on the Momentum, Mass, and Thermal Boundary Layers Above the Land Atmosphere Interface

NASA Astrophysics Data System (ADS)

Trautz, A.; Smits, K. M.; Illangasekare, T. H.; Schulte, P.

2014-12-01

The purpose of this study is to investigate the impacts of soil conditions (i.e. soil type, saturation) and atmospheric forcings (i.e. velocity, temperature, relative humidity) on the momentum, mass, and temperature boundary layers. The atmospheric conditions tested represent those typically found in semi-arid and arid climates and the soil conditions simulate the three stages of evaporation. The data generated will help identify the importance of different soil conditions and atmospheric forcings with respect to land-atmospheric interactions which will have direct implications on future numerical studies investigating the effects of turbulent air flow on evaporation. The experimental datasets generated for this study were performed using a unique climate controlled closed-circuit wind tunnel/porous media facility located at the Center for Experimental Study of Subsurface Environmental Processes (CESEP) at the Colorado School of Mines. The test apparatus consisting of a 7.3 m long porous media tank and wind tunnel, were outfitted with a sensor network to carefully measure wind velocity, air and soil temperature, relative humidity, soil moisture, and soil air pressure. Boundary layer measurements were made between the heights of 2 and 500 mm above the soil tank under constant conditions (i.e. wind velocity, temperature, relative humidity). The soil conditions (e.g. soil type, soil moisture) were varied between datasets to analyze their impact on the boundary layers. Experimental results show that the momentum boundary layer is very sensitive to the applied atmospheric conditions and soil conditions to a much less extent. Increases in velocity above porous media leads to momentum boundary layer thinning and closely reflect classical flat plate theory. The mass and thermal boundary layers are directly dependent on both atmospheric and soil conditions. Air pressure within the soil is independent of atmospheric temperature and relative humidity - wind velocity and soil moisture effects were observed. This data provides important insight into future work of accurately modeling the exchange processes associated with evaporation under various turbulent atmospheric conditions.

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

NASA Astrophysics Data System (ADS)

Kashiwabara, Takahito

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

Prediction of sound fields in acoustical cavities using the boundary element method. M.S. Thesis

NASA Technical Reports Server (NTRS)

Kipp, C. R.; Bernhard, R. J.

1985-01-01

A method was developed to predict sound fields in acoustical cavities. The method is based on the indirect boundary element method. An isoparametric quadratic boundary element is incorporated. Pressure, velocity and/or impedance boundary conditions may be applied to a cavity by using this method. The capability to include acoustic point sources within the cavity is implemented. The method is applied to the prediction of sound fields in spherical and rectangular cavities. All three boundary condition types are verified. Cases with a point source within the cavity domain are also studied. Numerically determined cavity pressure distributions and responses are presented. The numerical results correlate well with available analytical results.

Elastohydrodynamics of microfilament under distributed body actuation

NASA Astrophysics Data System (ADS)

Singh, T. Sonamani; Yadava, R. D. S.

2018-05-01

The dynamics of an active filament in low Reynolds (Re) number regime is analyzed under distributed body actuation represented by the sliding filament model. The governing elastohydrodynamic equations are formulated by assuming the resistive force theory (RFT). The effect of geometric nonlinearity in bending stiffness on the propulsive thrust has been analyzed where the former is introduced by cross-sectional tapering. Two types of boundary conditions (clamped-free and hinged-free) are analyzed. A comparison with the uniform filament dynamics reveals that the tapering enhances the thrust under both types of boundary conditions.

Nonequilibrium chemistry boundary layer integral matrix procedure

NASA Technical Reports Server (NTRS)

Tong, H.; Buckingham, A. C.; Morse, H. L.

1973-01-01

The development of an analytic procedure for the calculation of nonequilibrium boundary layer flows over surfaces of arbitrary catalycities is described. An existing equilibrium boundary layer integral matrix code was extended to include nonequilibrium chemistry while retaining all of the general boundary condition features built into the original code. For particular application to the pitch-plane of shuttle type vehicles, an approximate procedure was developed to estimate the nonequilibrium and nonisentropic state at the edge of the boundary layer.

Lichnerowicz-type equations with sign-changing nonlinearities on complete manifolds with boundary

NASA Astrophysics Data System (ADS)

Albanese, Guglielmo; Rigoli, Marco

2017-12-01

We prove an existence theorem for positive solutions to Lichnerowicz-type equations on complete manifolds with boundary (M , ∂ M , 〈 , 〉) and nonlinear Neumann conditions. This kind of nonlinear problems arise quite naturally in the study of solutions for the Einstein-scalar field equations of General Relativity in the framework of the so called Conformal Method.

Small scale effect on vibrational response of single-walled carbon nanotubes with different boundary conditions based on nonlocal beam models

NASA Astrophysics Data System (ADS)

Ansari, R.; Sahmani, S.

2012-04-01

The free vibration response of single-walled carbon nanotubes (SWCNTs) is investigated in this work using various nonlocal beam theories. To this end, the nonlocal elasticity equations of Eringen are incorporated into the various classical beam theories namely as Euler-Bernoulli beam theory (EBT), Timoshenko beam theory (TBT), and Reddy beam theory (RBT) to consider the size-effects on the vibration analysis of SWCNTs. The generalized differential quadrature (GDQ) method is employed to discretize the governing differential equations of each nonlocal beam theory corresponding to four commonly used boundary conditions. Then molecular dynamics (MD) simulation is implemented to obtain fundamental frequencies of nanotubes with different chiralities and values of aspect ratio to compare them with the results obtained by the nonlocal beam models. Through the fitting of the two series of numerical results, appropriate values of nonlocal parameter are derived relevant to each type of chirality, nonlocal beam model, and boundary conditions. It is found that in contrast to the chirality, the type of nonlocal beam model and boundary conditions make difference between the calibrated values of nonlocal parameter corresponding to each one.

Intrawellbore kinematic and frictional losses in a horizontal well in a bounded confined aquifer

NASA Astrophysics Data System (ADS)

Wang, Quanrong; Zhan, Hongbin

2017-01-01

Horizontal drilling has become an appealing technology for water resource exploration or aquifer remediation in recent decades, due to decreasing operational cost and many technical advantages over vertical wells. However, many previous studies on flow into horizontal wells were based on the Uniform Flux Boundary Condition (UFBC), which does not reflect the physical processes of flow inside the well accurately. In this study, we investigated transient flow into a horizontal well in an anisotropic confined aquifer laterally bounded by two constant-head boundaries. Three types of boundary conditions were employed to treat the horizontal well, including UFBC, Uniform-Head Boundary Condition (UHBC), and Mixed-Type Boundary Condition (MTBC). The MTBC model considered both kinematic and frictional effects inside the horizontal well, in which the kinematic effect referred to the accelerational and fluid-inflow effects. A new solution of UFBC was derived by superimposing the point sink/source solutions along the axis of a horizontal well with a uniform flux distribution. New solutions of UHBC and MTBC were obtained by a hybrid analytical-numerical method, and an iterative method was proposed to determine the well discretization required for achieving sufficiently accurate results. This study showed that the differences among the UFBC, UHBC, and MTBC solutions were obvious near the well screen, decreased with distance from the well, and became negligible near the constant-head boundary. The relationship between the flow rate and the drawdown was nonlinear for the MTBC solution, while it was linear for the UFBC and UHBC solutions.

Localized rotating convection with no-slip boundary conditions

NASA Astrophysics Data System (ADS)

Beaume, Cédric; Kao, Hsien-Ching; Knobloch, Edgar; Bergeon, Alain

2013-12-01

Localized patches of stationary convection embedded in a background conduction state are called convectons. Multiple states of this type have recently been found in two-dimensional Boussinesq convection in a horizontal fluid layer with stress-free boundary conditions at top and bottom, and rotating about the vertical. The convectons differ in their lengths and in the strength of the self-generated shear within which they are embedded, and exhibit slanted snaking. We use homotopic continuation of the boundary conditions to show that similar structures exist in the presence of no-slip boundary conditions at the top and bottom of the layer and show that such structures exhibit standard snaking. The homotopic continuation allows us to study the transformation from slanted snaking characteristic of systems with a conserved quantity, here the zonal momentum, to standard snaking characteristic of systems with no conserved quantity.

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

PubMed

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

2014-07-01

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

Techniques and computations for mapping plot clusters that straddle stand boundaries

Treesearch

Charles T. Scott; William A. Bechtold

1995-01-01

Many regional (extensive) forest surveys use clusters of subplots or prism points to reduce survey costs. Two common methods of handling clusters that straddle stand boundaries entail: (1) moving all subplots into a single forest cover type, or (2)"averaging" data across multiple conditions without regard to the boundaries. these methods result in biased...

Analysis of classical Fourier, SPL and DPL heat transfer model in biological tissues in presence of metabolic and external heat source

NASA Astrophysics Data System (ADS)

Kumar, Dinesh; Singh, Surjan; Rai, K. N.

2016-06-01

In this paper, the temperature distribution in a finite biological tissue in presence of metabolic and external heat source when the surface subjected to different type of boundary conditions is studied. Classical Fourier, single-phase-lag (SPL) and dual-phase-lag (DPL) models were developed for bio-heat transfer in biological tissues. The analytical solution obtained for all the three models using Laplace transform technique and results are compared. The effect of the variability of different parameters such as relaxation time, metabolic heat source, spatial heat source, different type boundary conditions on temperature distribution in different type of the tissues like muscle, tumor, fat, dermis and subcutaneous based on three models are analyzed and discussed in detail. The result obtained in three models is compared with experimental observation of Stolwijk and Hardy (Pflug Arch 291:129-162, 1966). It has been observe that the DPL bio-heat transfer model provides better result in comparison of other two models. The value of metabolic and spatial heat source in boundary condition of first, second and third kind for different type of thermal therapies are evaluated.

Cylindrical Couette flow of a rarefied gas: Effect of a boundary condition on the inverted velocity profile.

PubMed

Kosuge, Shingo

2015-07-01

The cylindrical Couette flow of a rarefied gas between a rotating inner cylinder and a stationary outer cylinder is investigated under the following two kinds of kinetic boundary conditions. One is the modified Maxwell-type boundary condition proposed by Dadzie and Méolans [J. Math. Phys. 45, 1804 (2004)] and the other is the Cercignani-Lampis condition, both of which have separate accommodation coefficients associated with the molecular velocity component normal to the boundary and with the tangential component. An asymptotic analysis of the Boltzmann equation for small Knudsen numbers and a numerical analysis of the Bhatnagar-Gross-Krook model equation for a wide range of the Knudsen number are performed to clarify the effect of each accommodation coefficient as well as of the boundary condition itself on the behavior of the gas, especially on the flow-velocity profile. As a result, the velocity-slip and temperature-jump conditions corresponding to the above kinetic boundary conditions are derived, which are necessary for the fluid-dynamic description of the problem for small Knudsen numbers. The parameter range for the onset of the velocity inversion phenomenon, which is related mainly to the decrease in the tangential momentum accommodation, is also obtained.

Numerical method for the solution of large systems of differential equations of the boundary layer type

NASA Technical Reports Server (NTRS)

Green, M. J.; Nachtsheim, P. R.

1972-01-01

A numerical method for the solution of large systems of nonlinear differential equations of the boundary-layer type is described. The method is a modification of the technique for satisfying asymptotic boundary conditions. The present method employs inverse interpolation instead of the Newton method to adjust the initial conditions of the related initial-value problem. This eliminates the so-called perturbation equations. The elimination of the perturbation equations not only reduces the user's preliminary work in the application of the method, but also reduces the number of time-consuming initial-value problems to be numerically solved at each iteration. For further ease of application, the solution of the overdetermined system for the unknown initial conditions is obtained automatically by applying Golub's linear least-squares algorithm. The relative ease of application of the proposed numerical method increases directly as the order of the differential-equation system increases. Hence, the method is especially attractive for the solution of large-order systems. After the method is described, it is applied to a fifth-order problem from boundary-layer theory.

Defect sink characteristics of specific grain boundary types in 304 stainless steels under high dose neutron environments

DOE PAGES

Field, Kevin G.; Yang, Ying; Busby, Jeremy T.; ...

2015-03-09

Radiation induced segregation (RIS) is a well-studied phenomena which occurs in many structurally relevant nuclear materials including austenitic stainless steels. RIS occurs due to solute atoms preferentially coupling to mobile point defect fluxes that migrate and interact with defect sinks. Here, a 304 stainless steel was neutron irradiated up to 47.1 dpa at 320 °C. Investigations into the RIS response at specific grain boundary types were utilized to determine the sink characteristics of different boundary types as a function of irradiation dose. A rate theory model built on the foundation of the modified inverse Kirkendall (MIK) model is proposed andmore » benchmarked to the experimental results. This model, termed the GiMIK model, includes alterations in the boundary conditions based on grain boundary structure and includes expressions for interstitial binding. This investigation, through experiment and modeling, found specific grain boundary structures exhibit unique defect sink characteristics depending on their local structure. Furthermore, such interactions were found to be consistent across all doses investigated and had larger global implications including precipitation of Ni-Si clusters near different grain boundary types.« less

The Forced Soft Spring Equation

ERIC Educational Resources Information Center

Fay, T. H.

2006-01-01

Through numerical investigations, this paper studies examples of the forced Duffing type spring equation with [epsilon] negative. By performing trial-and-error numerical experiments, the existence is demonstrated of stability boundaries in the phase plane indicating initial conditions yielding bounded solutions. Subharmonic boundaries are…

Comparability of ecological quality boundaries in the Mediterranean basin using freshwater benthic invertebrates. Statistical options and implications.

PubMed

Feio, M J; Ferreira, J; Buffagni, A; Erba, S; Dörflinger, G; Ferréol, M; Munné, A; Prat, N; Tziortzis, I; Urbanič, G

2014-04-01

Within the Mediterranean region each country has its own assessment method based on aquatic macroinvertebrates. However, independently of the classification system, quality assessments should be comparable across members of the European Commission, which means, among others, that the boundaries between classes should not deviate significantly. Here we check for comparability between High-Good and Good-Moderate classifications, through the use of a common metric. Additionally, we discuss the influence of the conceptual and statistical approaches used to calculate a common boundary within the Mediterranean countries participating in the Intercalibration Exercise (e.g., using individual national type-boundaries, one value for each common type or an average boundary by country; weighted average, median) in the overall outcome. All methods, except for the IBMWP (the Iberian BMWP) when applied to temporary rivers, were highly correlated (0.82

SPH for impact force and ricochet behavior of water-entry bodies

NASA Astrophysics Data System (ADS)

Omidvar, Pourya; Farghadani, Omid; Nikeghbali, Pooyan

The numerical modeling of fluid interaction with a bouncing body has many applications in scientific and engineering application. In this paper, the problem of water impact of a body on free-surface is investigated, where the fixed ghost boundary condition is added to the open source code SPHysics2D1 to rectify the oscillations in pressure distributions with the repulsive boundary condition. First, after introducing the methodology of SPH and the option of boundary conditions, the still water problem is simulated using two types of boundary conditions. It is shown that the fixed ghost boundary condition gives a better result for a hydrostatics pressure. Then, the dam-break problem, which is a bench mark test case in SPH, is simulated and compared with available data. In order to show the behavior of the hydrostatics forces on bodies, a fix/floating cylinder is placed on free surface looking carefully at the force and heaving profile. Finally, the impact of a body on free-surface is successfully simulated for different impact angles and velocities.

Stationary scalar clouds around a BTZ black hole

NASA Astrophysics Data System (ADS)

Ferreira, Hugo R. C.; Herdeiro, Carlos A. R.

2017-10-01

We establish the existence of stationary clouds of massive test scalar fields around BTZ black holes. These clouds are zero-modes of the superradiant instability and are possible when Robin boundary conditions (RBCs) are considered at the AdS boundary. These boundary conditions are the most general ones that ensure the AdS space is an isolated system, and include, as a particular case, the commonly considered Dirichlet or Neumann-type boundary conditions (DBCs or NBCs). We obtain an explicit, closed form, resonance condition, relating the RBCs that allow the existence of normalizable (and regular on and outside the horizon) clouds to the system's parameters. Such RBCs never include pure DBCs or NBCs. We illustrate the spatial distribution of these clouds, their energy and angular momentum density for some cases. Our results show that BTZ black holes with scalar hair can be constructed, as the non-linear realization of these clouds.

Chaotic Oscillations of Second Order Linear Hyperbolic Equations with Nonlinear Boundary Conditions: A Factorizable but Noncommutative Case

NASA Astrophysics Data System (ADS)

Li, Liangliang; Huang, Yu; Chen, Goong; Huang, Tingwen

If a second order linear hyperbolic partial differential equation in one-space dimension can be factorized as a product of two first order operators and if the two first order operators commute, with one boundary condition being the van der Pol type and the other being linear, one can establish the occurrence of chaos when the parameters enter a certain regime [Chen et al., 2014]. However, if the commutativity of the two first order operators fails to hold, then the treatment in [Chen et al., 2014] no longer works and significant new challenges arise in determining nonlinear boundary conditions that engenders chaos. In this paper, we show that by incorporating a linear memory effect, a nonlinear van der Pol boundary condition can cause chaotic oscillations when the parameter enters a certain regime. Numerical simulations illustrating chaotic oscillations are also presented.

Minimization of vibration in elastic beams with time-variant boundary conditions

NASA Technical Reports Server (NTRS)

Amirouche, F. M. L.; Xie, Mingjun

1992-01-01

This paper presents an innovative method for minimizing the vibration of structures with time-variant boundary conditions (supports). The elastic body is modeled in two ways: (1) the first model is a letter seven type beam with a movable mass not to exceed the lower tip; (2) the second model has an arm that is a hollow beam with an inside mass with adjustable position. The complete solutions to both problems are carried out where the body is undergoing large rotation. The quasi-static procedure is used for the time-variant boundary conditions. The method developed employs partial differential equations governing the motion of the beam, including the effects of rigid-body motion, time-variant boundary conditions, and calculus of variations. The analytical solution is developed using Laplace and Fourier transforms. Examples of elastic robotic arms are given to illustrate the effectiveness of the methods developed.

Inhomogeneous Radiation Boundary Conditions Simulating Incoming Acoustic Waves for Computational Aeroacoustics

NASA Technical Reports Server (NTRS)

Tam, Christopher K. W.; Fang, Jun; Kurbatskii, Konstantin A.

1996-01-01

A set of nonhomogeneous radiation and outflow conditions which automatically generate prescribed incoming acoustic or vorticity waves and, at the same time, are transparent to outgoing sound waves produced internally in a finite computation domain is proposed. This type of boundary condition is needed for the numerical solution of many exterior aeroacoustics problems. In computational aeroacoustics, the computation scheme must be as nondispersive ans nondissipative as possible. It must also support waves with wave speeds which are nearly the same as those of the original linearized Euler equations. To meet these requirements, a high-order/large-stencil scheme is necessary The proposed nonhomogeneous radiation and outflow boundary conditions are designed primarily for use in conjunction with such high-order/large-stencil finite difference schemes.

Initializing a Mesoscale Boundary-Layer Model with Radiosonde Observations

NASA Astrophysics Data System (ADS)

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.

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

NASA Astrophysics Data System (ADS)

Lau, Chun Sing

This thesis studies two types of problems in financial derivatives pricing. The first type is the free boundary problem, which can be formulated as a partial differential equation (PDE) subject to a set of free boundary condition. Although the functional form of the free boundary condition is given explicitly, the location of the free boundary is unknown and can only be determined implicitly by imposing continuity conditions on the solution. Two specific problems are studied in details, namely the valuation of fixed-rate mortgages and CEV American options. The second type is the multi-dimensional problem, which involves multiple correlated stochastic variables and their governing PDE. One typical problem we focus on is the valuation of basket-spread options, whose underlying asset prices are driven by correlated geometric Brownian motions (GBMs). Analytic approximate solutions are derived for each of these three problems. For each of the two free boundary problems, we propose a parametric moving boundary to approximate the unknown free boundary, so that the original problem transforms into a moving boundary problem which can be solved analytically. The governing parameter of the moving boundary is determined by imposing the first derivative continuity condition on the solution. The analytic form of the solution allows the price and the hedging parameters to be computed very efficiently. When compared against the benchmark finite-difference method, the computational time is significantly reduced without compromising the accuracy. The multi-stage scheme further allows the approximate results to systematically converge to the benchmark results as one recasts the moving boundary into a piecewise smooth continuous function. For the multi-dimensional problem, we generalize the Kirk (1995) approximate two-asset spread option formula to the case of multi-asset basket-spread option. Since the final formula is in closed form, all the hedging parameters can also be derived in closed form. Numerical examples demonstrate that the pricing and hedging errors are in general less than 1% relative to the benchmark prices obtained by numerical integration or Monte Carlo simulation. By exploiting an explicit relationship between the option price and the underlying probability distribution, we further derive an approximate distribution function for the general basket-spread variable. It can be used to approximate the transition probability distribution of any linear combination of correlated GBMs. Finally, an implicit perturbation is applied to reduce the pricing errors by factors of up to 100. When compared against the existing methods, the basket-spread option formula coupled with the implicit perturbation turns out to be one of the most robust and accurate approximation methods.

Generalized Forchheimer Flows of Isentropic Gases

NASA Astrophysics Data System (ADS)

Celik, Emine; Hoang, Luan; Kieu, Thinh

2018-03-01

We consider generalized Forchheimer flows of either isentropic gases or slightly compressible fluids in porous media. By using Muskat's and Ward's general form of the Forchheimer equations, we describe the fluid dynamics by a doubly nonlinear parabolic equation for the appropriately defined pseudo-pressure. The volumetric flux boundary condition is converted to a time-dependent Robin-type boundary condition for this pseudo-pressure. We study the corresponding initial boundary value problem, and estimate the L^∞ and W^{1,2-a} (with 0

On hydromagnetic oscillations in a rotating cavity.

NASA Technical Reports Server (NTRS)

Gans, R. F.

1971-01-01

Time-dependent hydromagnetic phenomena in a rotating spherical cavity are investigated in the framework of an interior boundary-layer expansion. The first type of wave is a modification of the hydrodynamic inertial wave, the second is a pseudo-geostrophic wave and is involved in spinup, and the third is related to the MAC waves of Braginskii (1967). It is shown that the MAC waves must satisfy more than the usual normal boundary conditions, and that reference must be made to the boundary-layer solution to resolve the ambiguity regarding which conditions are to be taken. The boundary-layer structure is investigated in detail to display the interactions between applied field, viscosity, electrical conductivity, frequency and latitu de.

Evolution of Grain Boundary Precipitates in an Al-Cu-Li Alloy During Aging

NASA Astrophysics Data System (ADS)

Ott, Noémie; Kairy, Shravan K.; Yan, Yuanming; Birbilis, Nick

2017-01-01

The grain boundary microstructure of Al-Cu-Li alloy AA2050 was investigated for different isothermal aging times to rationalize intergranular corrosion (IGC) characteristics. In the underaged condition, the dominant grain boundary precipitates are fine T1 (Al2CuLi). Extended aging revealed that grain boundaries were decorated by large T1 precipitates and S' phase (Al2CuMg), with S' growth not dimensionally constrained. Such a transition in the precipitate type at grain boundaries is a unique feature of the Al-Cu-Li system.

Homogenization of Winkler-Steklov spectral conditions in three-dimensional linear elasticity

NASA Astrophysics Data System (ADS)

Gómez, D.; Nazarov, S. A.; Pérez, M. E.

2018-04-01

We consider a homogenization Winkler-Steklov spectral problem that consists of the elasticity equations for a three-dimensional homogeneous anisotropic elastic body which has a plane part of the surface subject to alternating boundary conditions on small regions periodically placed along the plane. These conditions are of the Dirichlet type and of the Winkler-Steklov type, the latter containing the spectral parameter. The rest of the boundary of the body is fixed, and the period and size of the regions, where the spectral parameter arises, are of order ɛ . For fixed ɛ , the problem has a discrete spectrum, and we address the asymptotic behavior of the eigenvalues {β _k^ɛ }_{k=1}^{∞} as ɛ → 0. We show that β _k^ɛ =O(ɛ ^{-1}) for each fixed k, and we observe a common limit point for all the rescaled eigenvalues ɛ β _k^ɛ while we make it evident that, although the periodicity of the structure only affects the boundary conditions, a band-gap structure of the spectrum is inherited asymptotically. Also, we provide the asymptotic behavior for certain "groups" of eigenmodes.

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

PubMed

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

2014-09-01

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

Studies of the phase gradient at the boundary of the phase diffusion equation, motivated by peculiar wave patterns of rhythmic contraction in the amoeboid movement of Physarum polycephalum

NASA Astrophysics Data System (ADS)

Iima, Makoto; Kori, Hiroshi; Nakagaki, Toshiyuki

2017-04-01

The boundary of a cell is the interface with its surroundings and plays a key role in controlling the cell movement adaptations to different environments. We propose a study of the boundary effects on the patterns and waves of the rhythmic contractions in plasmodia of Physarum polycephalum, a tractable model organism of the amoeboid type. Boundary effects are defined as the effects of both the boundary conditions and the boundary shape. The rhythmicity of contraction can be modulated by local stimulation of temperature, light and chemicals, and by local deformation of cell shape via mechanosensitive ion channels as well. First, we examined the effects of boundary cell shapes in the case of a special shape resembling a tadpole, while requiring that the natural frequency in the proximity of the boundary is slightly higher and uniform. The simulation model reproduced the approximate propagated wave, from the tail to the head, while the inward waves were observed only near the periphery of the head section of the tadpole-shape. A key finding was that the frequency of the rhythmic contractions depended on the local shape of cell boundary. This implies that the boundary conditions of the phase were not always homogeneous. To understand the dependency, we reduced the two-dimensional model into a one-dimensional continuum model with Neumann boundary conditions. Here, the boundary conditions reflect the frequency distribution at the boundary. We described the analytic solutions and calculated the relationship between the boundary conditions and the wave propagation for a one-dimensional model of the continuous oscillatory field and a discrete coupled oscillator system. The results obtained may not be limited to cell movement of Physarum, but may be applicable to the other physical systems since the analysis used a generic phase diffusion equation.

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

NASA Astrophysics Data System (ADS)

Valchev, T. I.

2016-02-01

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

DOE Office of Scientific and Technical Information (OSTI.GOV)

Field, Kevin G.; Yang, Ying; Busby, Jeremy T.

Radiation induced segregation (RIS) is a well-studied phenomena which occurs in many structurally relevant nuclear materials including austenitic stainless steels. RIS occurs due to solute atoms preferentially coupling to mobile point defect fluxes that migrate and interact with defect sinks. Here, a 304 stainless steel was neutron irradiated up to 47.1 dpa at 320 °C. Investigations into the RIS response at specific grain boundary types were utilized to determine the sink characteristics of different boundary types as a function of irradiation dose. A rate theory model built on the foundation of the modified inverse Kirkendall (MIK) model is proposed andmore » benchmarked to the experimental results. This model, termed the GiMIK model, includes alterations in the boundary conditions based on grain boundary structure and includes expressions for interstitial binding. This investigation, through experiment and modeling, found specific grain boundary structures exhibit unique defect sink characteristics depending on their local structure. Furthermore, such interactions were found to be consistent across all doses investigated and had larger global implications including precipitation of Ni-Si clusters near different grain boundary types.« less

DOE Office of Scientific and Technical Information (OSTI.GOV)

Field, Kevin G.; Yang, Ying; Allen, Todd R.

Radiation induced segregation (RIS) is a well-studied phenomena which occurs in many structurally relevant nuclear materials including austenitic stainless steels. RIS occurs due to solute atoms preferentially coupling to mobile point defect fluxes that migrate and interact with defect sinks. Here, a 304 stainless steel was neutron irradiated up to 47.1 dpa at 320 °C. Investigations into the RIS response at specific grain boundary types were utilized to determine the sink characteristics of different boundary types as a function of irradiation dose. A rate theory model built on the foundation of the modified inverse Kirkendall (MIK) model is proposed andmore » benchmarked to the experimental results. This model, termed the GiMIK model, includes alterations in the boundary conditions based on grain boundary structure and includes expressions for interstitial binding. This investigation, through experiment and modeling, found specific grain boundary structures exhibit unique defect sink characteristics depending on their local structure. Such interactions were found to be consistent across all doses investigated and had larger global implications including precipitation of Ni-Si clusters near different grain boundary types.« less

Effect of geometric size on mechanical properties of dielectric elastomers based on an improved visco-hyperelastic film model

NASA Astrophysics Data System (ADS)

Chang, Mengzhou; Wang, Zhenqing; Tong, Liyong; Liang, Wenyan

2017-03-01

Dielectric polymers show complex mechanical behaviors with different boundary conditions, geometry size and pre-stress. A viscoelastic model suitable for inhomogeneous deformation is presented integrating the Kelvin-Voigt model in a new form in this work. For different types of uniaxial tensile test loading along the length direction of sample, single-step-relaxation tests, loading-unloading tests and tensile-creep-relaxation tests the improved model provides a quite favorable comparison with the experiment results. Moreover, The mechanical properties of test sample with several length-width ratios under different boundary conditions are also invested. The influences of the different boundary conditions are calculated with a stress applied on the boundary point and the result show that the fixed boundary will increase the stress compare with homogeneous deformation. In modeling the effect of pre-stress in the shear test, three pre-stressed mode are discussed. The model validation on the general mechanical behavior shows excellent predictive capability.

A shifted Jacobi collocation algorithm for wave type equations with non-local conservation conditions

NASA Astrophysics Data System (ADS)

Doha, Eid H.; Bhrawy, Ali H.; Abdelkawy, Mohammed A.

2014-09-01

In this paper, we propose an efficient spectral collocation algorithm to solve numerically wave type equations subject to initial, boundary and non-local conservation conditions. The shifted Jacobi pseudospectral approximation is investigated for the discretization of the spatial variable of such equations. It possesses spectral accuracy in the spatial variable. The shifted Jacobi-Gauss-Lobatto (SJ-GL) quadrature rule is established for treating the non-local conservation conditions, and then the problem with its initial and non-local boundary conditions are reduced to a system of second-order ordinary differential equations in temporal variable. This system is solved by two-stage forth-order A-stable implicit RK scheme. Five numerical examples with comparisons are given. The computational results demonstrate that the proposed algorithm is more accurate than finite difference method, method of lines and spline collocation approach

High order solution of Poisson problems with piecewise constant coefficients and interface jumps

NASA Astrophysics Data System (ADS)

Marques, Alexandre Noll; Nave, Jean-Christophe; Rosales, Rodolfo Ruben

2017-04-01

We present a fast and accurate algorithm to solve Poisson problems in complex geometries, using regular Cartesian grids. We consider a variety of configurations, including Poisson problems with interfaces across which the solution is discontinuous (of the type arising in multi-fluid flows). The algorithm is based on a combination of the Correction Function Method (CFM) and Boundary Integral Methods (BIM). Interface and boundary conditions can be treated in a fast and accurate manner using boundary integral equations, and the associated BIM. Unfortunately, BIM can be costly when the solution is needed everywhere in a grid, e.g. fluid flow problems. We use the CFM to circumvent this issue. The solution from the BIM is used to rewrite the problem as a series of Poisson problems in rectangular domains-which requires the BIM solution at interfaces/boundaries only. These Poisson problems involve discontinuities at interfaces, of the type that the CFM can handle. Hence we use the CFM to solve them (to high order of accuracy) with finite differences and a Fast Fourier Transform based fast Poisson solver. We present 2-D examples of the algorithm applied to Poisson problems involving complex geometries, including cases in which the solution is discontinuous. We show that the algorithm produces solutions that converge with either 3rd or 4th order of accuracy, depending on the type of boundary condition and solution discontinuity.

Non-CMC Solutions of the Einstein Constraint Equations on Compact Manifolds with Apparent Horizon Boundaries

NASA Astrophysics Data System (ADS)

Holst, Michael; Meier, Caleb; Tsogtgerel, G.

2018-01-01

In this article we continue our effort to do a systematic development of the solution theory for conformal formulations of the Einstein constraint equations on compact manifolds with boundary. By building in a natural way on our recent work in Holst and Tsogtgerel (Class Quantum Gravity 30:205011, 2013), and Holst et al. (Phys Rev Lett 100(16):161101, 2008, Commun Math Phys 288(2):547-613, 2009), and also on the work of Maxwell (J Hyperbolic Differ Eqs 2(2):521-546, 2005a, Commun Math Phys 253(3):561-583, 2005b, Math Res Lett 16(4):627-645, 2009) and Dain (Class Quantum Gravity 21(2):555-573, 2004), under reasonable assumptions on the data we prove existence of both near- and far-from-constant mean curvature (CMC) solutions for a class of Robin boundary conditions commonly used in the literature for modeling black holes, with a third existence result for CMC appearing as a special case. Dain and Maxwell addressed initial data engineering for space-times that evolve to contain black holes, determining solutions to the conformal formulation on an asymptotically Euclidean manifold in the CMC setting, with interior boundary conditions representing excised interior black hole regions. Holst and Tsogtgerel compiled the interior boundary results covered by Dain and Maxwell, and then developed general interior conditions to model the apparent horizon boundary conditions of Dainand Maxwell for compact manifolds with boundary, and subsequently proved existence of solutions to the Lichnerowicz equation on compact manifolds with such boundary conditions. This paper picks up where Holst and Tsogtgerel left off, addressing the general non-CMC case for compact manifolds with boundary. As in our previous articles, our focus here is again on low regularity data and on the interaction between different types of boundary conditions. While our work here serves primarily to extend the solution theory for the compact with boundary case, we also develop several technical tools that have potential for use for other cases.

Theoretical and Experimental Studies of the Transonic Flow Field and Associated Boundary Conditions near a Longitudinally-Slotted Wind-Tunnel Wall. Ph.D. Thesis

NASA Technical Reports Server (NTRS)

Everhart, Joel Lee

1988-01-01

A theoretical examination of the slotted-wall flow field is conducted to determine the appropriate wall pressure drop (or boundary condition) equation. This analysis improves the understanding of the fluid physics of these types of flow fields and helps in evaluating the uncertainties and limitations existing in previous mathematical developments. It is shown that the resulting slotted-wall boundary condition contains contributions from the airfoil-induced streamline curvature and the non-linear, quadratic, slot crossflow in addition to an often neglected linear term which results from viscous shearing in the slot. Existing and newly acquired experimental data are examined in the light of this formulation and theoretical developments.

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

NASA Astrophysics Data System (ADS)

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.

Numerical Treatment of Stokes Solvent Flow and Solute-Solvent Interfacial Dynamics for Nonpolar Molecules.

PubMed

Sun, Hui; Zhou, Shenggao; Moore, David K; Cheng, Li-Tien; Li, Bo

2016-05-01

We design and implement numerical methods for the incompressible Stokes solvent flow and solute-solvent interface motion for nonpolar molecules in aqueous solvent. The balance of viscous force, surface tension, and van der Waals type dispersive force leads to a traction boundary condition on the solute-solvent interface. To allow the change of solute volume, we design special numerical boundary conditions on the boundary of a computational domain through a consistency condition. We use a finite difference ghost fluid scheme to discretize the Stokes equation with such boundary conditions. The method is tested to have a second-order accuracy. We combine this ghost fluid method with the level-set method to simulate the motion of the solute-solvent interface that is governed by the solvent fluid velocity. Numerical examples show that our method can predict accurately the blow up time for a test example of curvature flow and reproduce the polymodal (e.g., dry and wet) states of hydration of some simple model molecular systems.

Numerical Treatment of Stokes Solvent Flow and Solute-Solvent Interfacial Dynamics for Nonpolar Molecules

PubMed Central

Sun, Hui; Zhou, Shenggao; Moore, David K.; Cheng, Li-Tien; Li, Bo

2015-01-01

We design and implement numerical methods for the incompressible Stokes solvent flow and solute-solvent interface motion for nonpolar molecules in aqueous solvent. The balance of viscous force, surface tension, and van der Waals type dispersive force leads to a traction boundary condition on the solute-solvent interface. To allow the change of solute volume, we design special numerical boundary conditions on the boundary of a computational domain through a consistency condition. We use a finite difference ghost fluid scheme to discretize the Stokes equation with such boundary conditions. The method is tested to have a second-order accuracy. We combine this ghost fluid method with the level-set method to simulate the motion of the solute-solvent interface that is governed by the solvent fluid velocity. Numerical examples show that our method can predict accurately the blow up time for a test example of curvature flow and reproduce the polymodal (e.g., dry and wet) states of hydration of some simple model molecular systems. PMID:27365866

Structural characteristics of the shock-induced boundary layer separation extended to the leading edge

NASA Astrophysics Data System (ADS)

Tao, Y.; Liu, W. D.; Fan, X. Q.; Zhao, Y. L.

2017-07-01

For a better understanding of the local unstart of supersonic/hypersonic inlet, a series of experiments has been conducted to investigate the shock-induced boundary layer separation extended to the leading edge. Using the nanoparticle-based planar laser scattering, we recorded the fine structures of these interactions under different conditions and paid more attention to their structural characteristics. According to their features, these interactions could be divided into four types. Specifically, Type A wave pattern is similar to the classic shock wave/turbulent boundary layer interaction, and Type B wave configuration consists of an overall Mach reflection above the large scale separation bubble. Due to the gradual decrease in the size of the separation bubble, the separation bubble was replaced by several vortices (Type C wave pattern). Besides, for Type D wave configuration which exists in the local unstart inlet, there appears to be some flow spillage around the leading edge.

Computer program for stress, stability, and vibration of complex branched shells of revolution: BOSOR 4

NASA Technical Reports Server (NTRS)

Bushnell, D.

1974-01-01

Code is easy to use yet is general with respect to: (a) type of analysis to be performed; (b) geometry of shell meridian; (c) type of wall construction; (d) type of boundary conditions, ring supports, and branching configuration; and (e) type of loading.

Boundary scattering in the Φ$$^{4}$$ model

DOE PAGES

Dorey, Patrick; Halavanau, Aliaksei; Mercer, James; ...

2017-05-19

Here, we study boundary scattering in themore » $$\\phi^4$$ model on a half-line with a one-parameter family of Neumann-type boundary conditions. A rich variety of phenomena is observed, which extends previously-studied behaviour on the full line to include regimes of near-elastic scattering, the restoration of a missing scattering window, and the creation of a kink or oscillon through the collision-induced decay of a metastable boundary state. We also study the decay of the vibrational boundary mode, and explore different scenarios for its relaxation and for the creation of kinks.« less

Generalized analytical solutions to sequentially coupled multi-species advective-dispersive transport equations in a finite domain subject to an arbitrary time-dependent source boundary condition

NASA Astrophysics Data System (ADS)

Chen, Jui-Sheng; Liu, Chen-Wuing; Liang, Ching-Ping; Lai, Keng-Hsin

2012-08-01

SummaryMulti-species advective-dispersive transport equations sequentially coupled with first-order decay reactions are widely used to describe the transport and fate of the decay chain contaminants such as radionuclide, chlorinated solvents, and nitrogen. Although researchers attempted to present various types of methods for analytically solving this transport equation system, the currently available solutions are mostly limited to an infinite or a semi-infinite domain. A generalized analytical solution for the coupled multi-species transport problem in a finite domain associated with an arbitrary time-dependent source boundary is not available in the published literature. In this study, we first derive generalized analytical solutions for this transport problem in a finite domain involving arbitrary number of species subject to an arbitrary time-dependent source boundary. Subsequently, we adopt these derived generalized analytical solutions to obtain explicit analytical solutions for a special-case transport scenario involving an exponentially decaying Bateman type time-dependent source boundary. We test the derived special-case solutions against the previously published coupled 4-species transport solution and the corresponding numerical solution with coupled 10-species transport to conduct the solution verification. Finally, we compare the new analytical solutions derived for a finite domain against the published analytical solutions derived for a semi-infinite domain to illustrate the effect of the exit boundary condition on coupled multi-species transport with an exponential decaying source boundary. The results show noticeable discrepancies between the breakthrough curves of all the species in the immediate vicinity of the exit boundary obtained from the analytical solutions for a finite domain and a semi-infinite domain for the dispersion-dominated condition.

Effects of thermal aging on the microstructure of Type-II boundaries in dissimilar metal weld joints

NASA Astrophysics Data System (ADS)

Yoo, Seung Chang; Choi, Kyoung Joon; Bahn, Chi Bum; Kim, Si Hoon; Kim, Ju Young; Kim, Ji Hyun

2015-04-01

In order to investigate the effects of long-term thermal aging on the microstructural evolution of Type-II boundary regions in the weld metal of Alloy 152, a representative dissimilar metal weld was fabricated from Alloy 690, Alloy 152, and A533 Gr.B. This mock-up was thermally aged at 450 °C to accelerate the effects of thermal aging in a nuclear power plant operation condition (320 °C). The microstructure of the Type-II boundary region of the weld root, which is parallel to and within 100 μm of the fusion boundary and known to be more susceptible to material degradation, was then characterized after different aging times using a scanning electron microscope equipped with an energy dispersive X-ray spectroscope for micro-compositional analysis, electron backscattered diffraction detector for grain and grain boundary orientation analysis, and a nanoindenter for measurement of mechanical properties. Through this, it was found that a steep compositional gradient and high grain average misorientation is created in the narrow zone between the Type-II and fusion boundaries, while the concentration of chromium and number of low-angle grain boundaries increases with aging time. A high average hardness was also observed in the same region of the dissimilar metal welds, with hardness peaking with thermal aging simulating an operational time of 15 years.

An iterative truncation method for unbounded electromagnetic problems using varying order finite elements

NASA Astrophysics Data System (ADS)

Paul, Prakash

2009-12-01

The finite element method (FEM) is used to solve three-dimensional electromagnetic scattering and radiation problems. Finite element (FE) solutions of this kind contain two main types of error: discretization error and boundary error. Discretization error depends on the number of free parameters used to model the problem, and on how effectively these parameters are distributed throughout the problem space. To reduce the discretization error, the polynomial order of the finite elements is increased, either uniformly over the problem domain or selectively in those areas with the poorest solution quality. Boundary error arises from the condition applied to the boundary that is used to truncate the computational domain. To reduce the boundary error, an iterative absorbing boundary condition (IABC) is implemented. The IABC starts with an inexpensive boundary condition and gradually improves the quality of the boundary condition as the iteration continues. An automatic error control (AEC) is implemented to balance the two types of error. With the AEC, the boundary condition is improved when the discretization error has fallen to a low enough level to make this worth doing. The AEC has these characteristics: (i) it uses a very inexpensive truncation method initially; (ii) it allows the truncation boundary to be very close to the scatterer/radiator; (iii) it puts more computational effort on the parts of the problem domain where it is most needed; and (iv) it can provide as accurate a solution as needed depending on the computational price one is willing to pay. To further reduce the computational cost, disjoint scatterers and radiators that are relatively far from each other are bounded separately and solved using a multi-region method (MRM), which leads to savings in computational cost. A simple analytical way to decide whether the MRM or the single region method will be computationally cheaper is also described. To validate the accuracy and savings in computation time, different shaped metallic and dielectric obstacles (spheres, ogives, cube, flat plate, multi-layer slab etc.) are used for the scattering problems. For the radiation problems, waveguide excited antennas (horn antenna, waveguide with flange, microstrip patch antenna) are used. Using the AEC the peak reduction in computation time during the iteration is typically a factor of 2, compared to the IABC using the same element orders throughout. In some cases, it can be as high as a factor of 4.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Dorey, Patrick; Halavanau, Aliaksei; Mercer, James

Here, we study boundary scattering in themore » $$\\phi^4$$ model on a half-line with a one-parameter family of Neumann-type boundary conditions. A rich variety of phenomena is observed, which extends previously-studied behaviour on the full line to include regimes of near-elastic scattering, the restoration of a missing scattering window, and the creation of a kink or oscillon through the collision-induced decay of a metastable boundary state. We also study the decay of the vibrational boundary mode, and explore different scenarios for its relaxation and for the creation of kinks.« less

Completed Beltrami-Michell formulation for analyzing mixed boundary value problems in elasticity

NASA Technical Reports Server (NTRS)

Patnaik, Surya N.; Kaljevic, Igor; Hopkins, Dale A.; Saigal, Sunil

1995-01-01

In elasticity, the method of forces, wherein stress parameters are considered as the primary unknowns, is known as the Beltrami-Michell formulation (BMF). The existing BMF can only solve stress boundary value problems; it cannot handle the more prevalent displacement of mixed boundary value problems of elasticity. Therefore, this formulation, which has restricted application, could not become a true alternative to the Navier's displacement method, which can solve all three types of boundary value problems. The restrictions in the BMF have been alleviated by augmenting the classical formulation with a novel set of conditions identified as the boundary compatibility conditions. This new method, which completes the classical force formulation, has been termed the completed Beltrami-Michell formulation (CBMF). The CBMF can solve general elasticity problems with stress, displacement, and mixed boundary conditions in terms of stresses as the primary unknowns. The CBMF is derived from the stationary condition of the variational functional of the integrated force method. In the CBMF, stresses for kinematically stable structures can be obtained without any reference to the displacements either in the field or on the boundary. This paper presents the CBMF and its derivation from the variational functional of the integrated force method. Several examples are presented to demonstrate the applicability of the completed formulation for analyzing mixed boundary value problems under thermomechanical loads. Selected example problems include a cylindrical shell wherein membrane and bending responses are coupled, and a composite circular plate.

Boundary States and Broken Bulk Symmetries in WAr Minimal Models

NASA Astrophysics Data System (ADS)

Caldeira, Alexandre F.; Wheater, J. F.

We review the free-field formalism for boundary states. The multi-component free-field formalism is then used to study the boundary states of (p',p) rational conformal field theories having a W symmetry of the type Ar. We show how the classification of primary fields for these models is obtained by demanding modular covariance of cylinder amplitudes and that the resulting modular S matrix satisfies all the necessary conditions. Basis states satisfying the boundary conditions are found in the form of coherent states and as expected we find that W violating states can be found for all these models. We construct consistent physical boundary states for all the rank 2 (p + 1,p) models (of which the already known case of the 3-state Potts model is the simplest example) and find that the W violating sector possesses a direct analogue of the Verlinde formula.

The effect of intra-wellbore head losses in a vertical well

NASA Astrophysics Data System (ADS)

Wang, Quanrong; Zhan, Hongbin

2017-05-01

Flow to a partially penetrating vertical well is made more complex by intra-wellbore losses. These are caused not only by the frictional effect, but also by the kinematic effect, which consists of the accelerational and fluid inflow effects inside a wellbore. Existing models of flow to a partially penetrating vertical well assume either a uniform-flux boundary condition (UFBC) or a uniform-head boundary condition (UHBC) for treating the flow into the wellbore. Neither approach considers intra-wellbore losses. In this study a new general solution, named the mixed-type boundary condition (MTBC) solution, is introduced to include intra-wellbore losses. It is developed from the existing solutions using a hybrid analytical-numerical method. The MTBC solution is capable of modeling various types of aquifer tests (constant-head tests, constant-rate tests, and slug tests) for partially or fully penetrating vertical wells in confined aquifers. Results show that intra-wellbore losses (both frictional and kinematic) can be significant in the early pumping stage. At later pumping times the UHBC solution is adequate because the difference between the MTBC and UHBC solutions becomes negligible.

Effect of Boundary Conditions on the Back Face Deformations of Flat UHMWPE Panels

DTIC Science & Technology

2014-12-01

Zhang [2] carried out a numerical study of the effects of clamping type and clamping pressure on the ballistic performance of woven Kevlar , and found...effects of composite size were also studied. Singletary [5] studied the effects of boundary conditions and panel sizes on V50 for Kevlar KM2 fabric. The...on the BFD in flat UHMWPE panels. UHMWPE possesses high tenacity and high strength compared to Kevlar , as a result of which it is the material of

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zanetti, F.M.; Vicentini, E.; Luz, M.G.E. da

It was proposed about a decade ago [M.G.E. da Luz, A.S. Lupu-Sax, E.J. Heller, Phys. Rev. E 56 (1997) 2496] a simple approach for obtaining scattering states for arbitrary disconnected open or closed boundaries C, with different boundary conditions. Since then, the so called boundary wall method has been successfully used to solve different open boundary problems. However, its applicability to closed shapes has not been fully explored. In this contribution we present a complete account of how to use the boundary wall to the case of billiard systems. We review the general ideas and particularize them to single connectedmore » closed shapes, assuming Dirichlet boundary conditions for the C's. We discuss the mathematical aspects that lead to both the inside and outside solutions. We also present a different way to calculate the exterior scattering S matrix. From it, we revisit the important inside-outside duality for billiards. Finally, we give some numerical examples, illustrating the efficiency and flexibility of the method to treat this type of problem.« less

46 CFR 171.080 - Damage stability standards for vessels with Type I or Type II subdivision.

Code of Federal Regulations, 2014 CFR

2014-10-01

... following conditions in the final stage of flooding: (1) On a vessel required to survive assumed damage with... in the final stage of flooding and to meet the conditions set forth in paragraphs (f) (8) and (9) of this section in each intermediate stage of flooding. For the purposes of establishing boundaries to...

46 CFR 171.080 - Damage stability standards for vessels with Type I or Type II subdivision.

Code of Federal Regulations, 2012 CFR

2012-10-01

... following conditions in the final stage of flooding: (1) On a vessel required to survive assumed damage with... in the final stage of flooding and to meet the conditions set forth in paragraphs (f) (8) and (9) of this section in each intermediate stage of flooding. For the purposes of establishing boundaries to...

46 CFR 171.080 - Damage stability standards for vessels with Type I or Type II subdivision.

Code of Federal Regulations, 2013 CFR

2013-10-01

... following conditions in the final stage of flooding: (1) On a vessel required to survive assumed damage with... in the final stage of flooding and to meet the conditions set forth in paragraphs (f) (8) and (9) of this section in each intermediate stage of flooding. For the purposes of establishing boundaries to...

Consistent lattice Boltzmann modeling of low-speed isothermal flows at finite Knudsen numbers in slip-flow regime: Application to plane boundaries

NASA Astrophysics Data System (ADS)

Silva, Goncalo; Semiao, Viriato

2017-07-01

The first nonequilibrium effect experienced by gaseous flows in contact with solid surfaces is the slip-flow regime. While the classical hydrodynamic description holds valid in bulk, at boundaries the fluid-wall interactions must consider slip. In comparison to the standard no-slip Dirichlet condition, the case of slip formulates as a Robin-type condition for the fluid tangential velocity. This makes its numerical modeling a challenging task, particularly in complex geometries. In this work, this issue is handled with the lattice Boltzmann method (LBM), motivated by the similarities between the closure relations of the reflection-type boundary schemes equipping the LBM equation and the slip velocity condition established by slip-flow theory. Based on this analogy, we derive, as central result, the structure of the LBM boundary closure relation that is consistent with the second-order slip velocity condition, applicable to planar walls. Subsequently, three tasks are performed. First, we clarify the limitations of existing slip velocity LBM schemes, based on discrete analogs of kinetic theory fluid-wall interaction models. Second, we present improved slip velocity LBM boundary schemes, constructed directly at discrete level, by extending the multireflection framework to the slip-flow regime. Here, two classes of slip velocity LBM boundary schemes are considered: (i) linear slip schemes, which are local but retain some calibration requirements and/or operation limitations, (ii) parabolic slip schemes, which use a two-point implementation but guarantee the consistent prescription of the intended slip velocity condition, at arbitrary plane wall discretizations, further dispensing any numerical calibration procedure. Third and final, we verify the improvements of our proposed slip velocity LBM boundary schemes against existing ones. The numerical tests evaluate the ability of the slip schemes to exactly accommodate the steady Poiseuille channel flow solution, over distinct wall slippage conditions, namely, no-slip, first-order slip, and second-order slip. The modeling of channel walls is discussed at both lattice-aligned and non-mesh-aligned configurations: the first case illustrates the numerical slip due to the incorrect modeling of slippage coefficients, whereas the second case adds the effect of spurious boundary layers created by the deficient accommodation of bulk solution. Finally, the slip-flow solutions predicted by LBM schemes are further evaluated for the Knudsen's paradox problem. As conclusion, this work establishes the parabolic accuracy of slip velocity schemes as the necessary condition for the consistent LBM modeling of the slip-flow regime.

Consistent lattice Boltzmann modeling of low-speed isothermal flows at finite Knudsen numbers in slip-flow regime: Application to plane boundaries.

PubMed

Silva, Goncalo; Semiao, Viriato

2017-07-01

The first nonequilibrium effect experienced by gaseous flows in contact with solid surfaces is the slip-flow regime. While the classical hydrodynamic description holds valid in bulk, at boundaries the fluid-wall interactions must consider slip. In comparison to the standard no-slip Dirichlet condition, the case of slip formulates as a Robin-type condition for the fluid tangential velocity. This makes its numerical modeling a challenging task, particularly in complex geometries. In this work, this issue is handled with the lattice Boltzmann method (LBM), motivated by the similarities between the closure relations of the reflection-type boundary schemes equipping the LBM equation and the slip velocity condition established by slip-flow theory. Based on this analogy, we derive, as central result, the structure of the LBM boundary closure relation that is consistent with the second-order slip velocity condition, applicable to planar walls. Subsequently, three tasks are performed. First, we clarify the limitations of existing slip velocity LBM schemes, based on discrete analogs of kinetic theory fluid-wall interaction models. Second, we present improved slip velocity LBM boundary schemes, constructed directly at discrete level, by extending the multireflection framework to the slip-flow regime. Here, two classes of slip velocity LBM boundary schemes are considered: (i) linear slip schemes, which are local but retain some calibration requirements and/or operation limitations, (ii) parabolic slip schemes, which use a two-point implementation but guarantee the consistent prescription of the intended slip velocity condition, at arbitrary plane wall discretizations, further dispensing any numerical calibration procedure. Third and final, we verify the improvements of our proposed slip velocity LBM boundary schemes against existing ones. The numerical tests evaluate the ability of the slip schemes to exactly accommodate the steady Poiseuille channel flow solution, over distinct wall slippage conditions, namely, no-slip, first-order slip, and second-order slip. The modeling of channel walls is discussed at both lattice-aligned and non-mesh-aligned configurations: the first case illustrates the numerical slip due to the incorrect modeling of slippage coefficients, whereas the second case adds the effect of spurious boundary layers created by the deficient accommodation of bulk solution. Finally, the slip-flow solutions predicted by LBM schemes are further evaluated for the Knudsen's paradox problem. As conclusion, this work establishes the parabolic accuracy of slip velocity schemes as the necessary condition for the consistent LBM modeling of the slip-flow regime.

The effect of thermal processing on microstructure and mechanical properties in a nickel-iron alloy

NASA Astrophysics Data System (ADS)

Yang, Ling

The correlation between processing conditions, resulted microstructure and mechanical properties is of interest in the field of metallurgy for centuries. In this work, we investigated the effect of thermal processing parameters on microstructure, and key mechanical properties to turbine rotor design: tensile yield strength and crack growth resistance, for a nickel-iron based superalloy Inconel 706. The first step of the designing of experiments is to find parameter ranges for thermal processing. Physical metallurgy on superalloys was combined with finite element analysis to estimate variations in thermal histories for a large Alloy 706 forging, and the results were adopted for designing of experiments. Through the systematic study, correlation was found between the processing parameters and the microstructure. Five different types of grain boundaries were identified by optical metallography, fractography, and transmission electron microscopy, and they were found to be associated with eta precipitation at the grain boundaries. Proportions of types of boundaries, eta size, spacing and angle respect to the grain boundary were found to be dependent on processing parameters. Differences in grain interior precipitates were also identified, and correlated with processing conditions. Further, a strong correlation between microstructure and mechanical properties was identified. The grain boundary precipitates affect the time dependent crack propagation resistance, and different types of boundaries have different levels of resistance. Grain interior precipitates were correlated with tensile yield strength. It was also found that there is a strong environmental effect on time dependent crack propagation resistance, and the sensitivity to environmental damage is microstructure dependent. The microstructure with eta decorated on grain boundaries by controlled processing parameters is more resistant to environmental damage through oxygen embrittlement than material without eta phase on grain boundaries. Effort was made to explore the mechanisms of improving the time dependent crack propagation resistance through thermal processing, several mechanisms were identified in both environment dependent and environment independent category, and they were ranked based on their contributions in affecting crack propagation.

SILICON CARBIDE GRAIN BOUNDARY DISTRIBUTIONS, IRRADIATION CONDITIONS, AND SILVER RETENTION IN IRRADIATED AGR-1 TRISO FUEL PARTICLES

DOE Office of Scientific and Technical Information (OSTI.GOV)

Lillo, T. M.; Rooyen, I. J.; Aguiar, J. A.

Precession electron diffraction in the transmission electron microscope was used to map grain orientation and ultimately determine grain boundary misorientation angle distributions, relative fractions of grain boundary types (random high angle, low angle or coincident site lattice (CSL)-related boundaries) and the distributions of CSL-related grain boundaries in the SiC layer of irradiated TRISO-coated fuel particles. Two particles from the AGR-1 experiment exhibiting high Ag-110m retention (>80%) were compared to a particle exhibiting low Ag-110m retention (<19%). Irradiated particles with high Ag-110m retention exhibited a lower fraction of random, high angle grain boundaries compared to the low Ag-110m retention particle. Anmore » inverse relationship between the random, high angle grain boundary fraction and Ag-110m retention is found and is consistent with grain boundary percolation theory. Also, comparison of the grain boundary distributions with previously reported unirradiated grain boundary distributions, based on SEM-based EBSD for similarly fabricated particles, showed only small differences, i.e. a greater low angle grain boundary fraction in unirradiated SiC. It was, thus, concluded that SiC layers with grain boundary distributions susceptible to Ag-110m release were present prior to irradiation. Finally, irradiation parameters were found to have little effect on the association of fission product precipitates with specific grain boundary types.« less

Vibration analysis of rotating functionally graded Timoshenko microbeam based on modified couple stress theory under different temperature distributions

NASA Astrophysics Data System (ADS)

Ghadiri, Majid; Shafiei, Navvab

2016-04-01

In this study, thermal vibration of rotary functionally graded Timoshenko microbeam has been analyzed based on modified couple stress theory considering temperature change in four types of temperature distribution on thermal environment. Material properties of FG microbeam are supposed to be temperature dependent and vary continuously along the thickness according to the power-law form. The axial forces are also included in the model as the thermal and true spatial variation due to the rotation. Governing equations and boundary conditions have been derived by employing Hamiltonian's principle. The differential quadrature method is employed to solve the governing equations for cantilever and propped cantilever boundary conditions. Validations are done by comparing available literatures and obtained results which indicate accuracy of applied method. Results represent effects of temperature changes, different boundary conditions, nondimensional angular velocity, length scale parameter, different boundary conditions, FG index and beam thickness on fundamental, second and third nondimensional frequencies. Results determine critical values of temperature changes and other essential parameters which can be applicable to design micromachines like micromotor and microturbine.

The numerical-analytical implementation of the cross-sections method to the open waveguide transition of the "horn" type

NASA Astrophysics Data System (ADS)

Divakov, Dmitriy; Malykh, Mikhail; Sevastianov, Leonid; Sevastianov, Anton; Tiutiunnik, Anastasiia

2017-04-01

In the paper we construct a method for approximate solution of the waveguide problem for guided modes of an open irregular waveguide transition. The method is based on straightening of the curved waveguide boundaries by introducing new variables and applying the Kantorovich method to the problem formulated in the new variables to get a system of ordinary second-order differential equations. In the method, the boundary conditions are formulated by analogy with the partial radiation conditions in the similar problem for closed waveguide transitions. The method is implemented in the symbolic-numeric form using the Maple computer algebra system. The coefficient matrices of the system of differential equations and boundary conditions are calculated symbolically, and then the obtained boundary-value problem is solved numerically using the finite difference method. The chosen coordinate functions of Kantorovich expansions provide good conditionality of the coefficient matrices. The numerical experiment simulating the propagation of guided modes in the open waveguide transition confirms the validity of the method proposed to solve the problem.

Weighted least-squares solver for determining pressure from particle image velocimetry data

NASA Astrophysics Data System (ADS)

de Kat, Roeland

2016-11-01

Currently, most approaches to determine pressure from particle image velocimetry data are Poisson approaches (e.g.) or multi-pass marching approaches (e.g.). However, these approaches deal with boundary conditions in their specific ways which cannot easily be changed-Poisson approaches enforce boundary conditions strongly, whereas multi-pass marching approaches enforce them weakly. Under certain conditions (depending on the certainty of the data or availability of reference data along the boundary) both types of boundary condition enforcement have to be used together to obtain the best result. In addition, neither of the approaches takes the certainty of the particle image velocimetry data (see e.g.) within the domain into account. Therefore, to address these shortcomings and improve upon current approaches, a new approach is proposed using weighted least-squares. The performance of this new approach is tested on synthetic and experimental particle image velocimetry data. Preliminary results show that a significant improvement can be made in determining pressure fields using the new approach. RdK is supported by a Leverhulme Trust Early Career Fellowship.

The asymptotic spectra of banded Toeplitz and quasi-Toeplitz matrices

NASA Technical Reports Server (NTRS)

Beam, Richard M.; Warming, Robert F.

1991-01-01

Toeplitz matrices occur in many mathematical, as well as, scientific and engineering investigations. This paper considers the spectra of banded Toeplitz and quasi-Toeplitz matrices with emphasis on non-normal matrices of arbitrarily large order and relatively small bandwidth. These are the type of matrices that appear in the investigation of stability and convergence of difference approximations to partial differential equations. Quasi-Toeplitz matrices are the result of non-Dirichlet boundary conditions for the difference approximations. The eigenvalue problem for a banded Toeplitz or quasi-Toeplitz matrix of large order is, in general, analytically intractable and (for non-normal matrices) numerically unreliable. An asymptotic (matrix order approaches infinity) approach partitions the eigenvalue analysis of a quasi-Toeplitz matrix into two parts, namely the analysis for the boundary condition independent spectrum and the analysis for the boundary condition dependent spectrum. The boundary condition independent spectrum is the same as the pure Toeplitz matrix spectrum. Algorithms for computing both parts of the spectrum are presented. Examples are used to demonstrate the utility of the algorithms, to present some interesting spectra, and to point out some of the numerical difficulties encountered when conventional matrix eigenvalue routines are employed for non-normal matrices of large order. The analysis for the Toeplitz spectrum also leads to a diagonal similarity transformation that improves conventional numerical eigenvalue computations. Finally, the algorithm for the asymptotic spectrum is extended to the Toeplitz generalized eigenvalue problem which occurs, for example, in the stability of Pade type difference approximations to differential equations.

A new analytical solution solved by triple series equations method for constant-head tests in confined aquifers

NASA Astrophysics Data System (ADS)

Chang, Ya-Chi; Yeh, Hund-Der

2010-06-01

The constant-head pumping tests are usually employed to determine the aquifer parameters and they can be performed in fully or partially penetrating wells. Generally, the Dirichlet condition is prescribed along the well screen and the Neumann type no-flow condition is specified over the unscreened part of the test well. The mathematical model describing the aquifer response to a constant-head test performed in a fully penetrating well can be easily solved by the conventional integral transform technique under the uniform Dirichlet-type condition along the rim of wellbore. However, the boundary condition for a test well with partial penetration should be considered as a mixed-type condition. This mixed boundary value problem in a confined aquifer system of infinite radial extent and finite vertical extent is solved by the Laplace and finite Fourier transforms in conjunction with the triple series equations method. This approach provides analytical results for the drawdown in a partially penetrating well for arbitrary location of the well screen in a finite thickness aquifer. The semi-analytical solutions are particularly useful for the practical applications from the computational point of view.

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

PubMed

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.

Early Warning Signals for Regime Transition in the Stable Boundary Layer: A Model Study

NASA Astrophysics Data System (ADS)

van Hooijdonk, I. G. S.; Moene, A. F.; Scheffer, M.; Clercx, H. J. H.; van de Wiel, B. J. H.

2017-02-01

The evening transition is investigated in an idealized model for the nocturnal boundary layer. From earlier studies it is known that the nocturnal boundary layer may manifest itself in two distinct regimes, depending on the ambient synoptic conditions: strong-wind or overcast conditions typically lead to weakly stable, turbulent nights; clear-sky and weak-wind conditions, on the other hand, lead to very stable, weakly turbulent conditions. Previously, the dynamical behaviour near the transition between these regimes was investigated in an idealized setting, relying on Monin-Obukhov (MO) similarity to describe turbulent transport. Here, we investigate a similar set-up, using direct numerical simulation; in contrast to MO-based models, this type of simulation does not need to rely on turbulence closure assumptions. We show that previous predictions are verified, but now independent of turbulence parametrizations. Also, it appears that a regime shift to the very stable state is signaled in advance by specific changes in the dynamics of the turbulent boundary layer. Here, we show how these changes may be used to infer a quantitative estimate of the transition point from the weakly stable boundary layer to the very stable boundary layer. In addition, it is shown that the idealized, nocturnal boundary-layer system shares important similarities with generic non-linear dynamical systems that exhibit critical transitions. Therefore, the presence of other, generic early warning signals is tested as well. Indeed, indications are found that such signals are present in stably stratified turbulent flows.

Matrix Sturm-Liouville equation with a Bessel-type singularity on a finite interval

NASA Astrophysics Data System (ADS)

Bondarenko, Natalia

2017-03-01

The matrix Sturm-Liouville equation on a finite interval with a Bessel-type singularity in the end of the interval is studied. Special fundamental systems of solutions for this equation are constructed: analytic Bessel-type solutions with the prescribed behavior at the singular point and Birkhoff-type solutions with the known asymptotics for large values of the spectral parameter. The asymptotic formulas for Stokes multipliers, connecting these two fundamental systems of solutions, are derived. We also set boundary conditions and obtain asymptotic formulas for the spectral data (the eigenvalues and the weight matrices) of the boundary value problem. Our results will be useful in the theory of direct and inverse spectral problems.

On flows of viscoelastic fluids under threshold-slip boundary conditions

NASA Astrophysics Data System (ADS)

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.

Spatially localized convection in a rotating layer

NASA Astrophysics Data System (ADS)

Knobloch, Edgar; Beaume, Cedric; Bergeon, Alain; Kao, Hsien-Ching

2014-11-01

We study two-dimensional stationary convection in a horizontal fluid layer heated from below and rotating about the vertical. With stress-free boundary conditions at top and bottom, spatially localized states can be found that are embedded in a self-generated background shear zone and lie on a pair of intertwined solution branches exhibiting ``slanted snaking.'' States of this type are present even in the absence of bistability between conduction and periodic convection - a consequence of the conservation of zonal momentum. With no-slip boundary conditions this quantity is no longer conserved but localized states continue to exist. These are no longer embedded in a background shear zone and exhibit standard snaking. Homotopic continuation from free-slip to no-slip boundary conditions is used to track the changes in the properties of the solutions and the associated bifurcation diagrams.

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

NASA Astrophysics Data System (ADS)

Gvozdic, Branko D.; Djurdjevic, Dusan Z.

2017-01-01

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

A nonreflecting upper boundary condition for anelastic nonhydrostatic mesoscale gravity-wave models

NASA Technical Reports Server (NTRS)

Kim, Young-Joon; Kar, Sajal K.; Arakawa, Akio

1993-01-01

A sponge layer is formulated to prevent spurious reflection of vertically propagating quasi-stationary gravity waves at the upper boundary of a two-dimensional numerical anelastic nonhydrostatic model. The sponge layer includes damping of both Newtonian-cooling type and Rayleigh-friction type, whose coefficients are determined in such a way that the reflectivity of wave energy at the bottom of the layer is zero. Unlike the formulations in earlier studies, our formulation includes the effects of vertical discretization, vertical mean density variation, and nonhydrostaticity. This sponge formulation is found effective in suppressing false downward reflection of waves for various types of quasi-stationary forcing.

General stability of memory-type thermoelastic Timoshenko beam acting on shear force

NASA Astrophysics Data System (ADS)

Apalara, Tijani A.

2018-03-01

In this paper, we consider a linear thermoelastic Timoshenko system with memory effects where the thermoelastic coupling is acting on shear force under Neumann-Dirichlet-Dirichlet boundary conditions. The same system with fully Dirichlet boundary conditions was considered by Messaoudi and Fareh (Nonlinear Anal TMA 74(18):6895-6906, 2011, Acta Math Sci 33(1):23-40, 2013), but they obtained a general stability result which depends on the speeds of wave propagation. In our case, we obtained a general stability result irrespective of the wave speeds of the system.

Effect of Boundary Conditions on the Axial Compression Buckling of Homogeneous Orthotropic Composite Cylinders in the Long Column Range

NASA Technical Reports Server (NTRS)

Mikulas, Martin M., Jr.; Nemeth, Michael P.; Oremont, Leonard; Jegley, Dawn C.

2011-01-01

Buckling loads for long isotropic and laminated cylinders are calculated based on Euler, Fluegge and Donnell's equations. Results from these methods are presented using simple parameters useful for fundamental design work. Buckling loads for two types of simply supported boundary conditions are calculated using finite element methods for comparison to select cases of the closed form solution. Results indicate that relying on Donnell theory can result in an over-prediction of buckling loads by as much as 40% in isotropic materials.

Numerical analysis of diffusion around a suspended expressway by a multi-scale CFD model

NASA Astrophysics Data System (ADS)

kondo, Hiroaki; Asahi, Kazutake; Tomizuka, Takayuki; Suzuki, Motoo

The diffusion of NO x around Ikegami-Shinmachi crossroads, which are among the most polluted roadside areas in Japan, was analyzed with a CFD model. This is a suspended four-lane express road with a six-lane ground-level road under the expressway and another four-lane ground-level road intersecting the two roads. Three types of boundary conditions for the CFD model were tested. In the first case, the boundary conditions were given with the results from the mesoscale meteorological model; in other words, the model was multi-scale. In the second case, the boundary conditions were given with the local one-point observation. In the third case, the conditions for the wind were given with the observation, and those for the turbulence were given with the mesoscale numerical model. All of the calculations indicated high concentrations in the morning and low ones in the afternoon, but they did not indicate high concentrations in the evening. The reasons for such time variations of NO x concentrations were investigated from the viewpoints of the wind direction, velocity, and boundary layer height. The results suggested that the extremely high concentration was generated by local sources and advection from the large source area of Tokyo. On the whole, the calculation with the boundary condition with the mesoscale model appears to be better than the other calculations.

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

NASA Astrophysics Data System (ADS)

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.

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

NASA Astrophysics Data System (ADS)

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.

Analytical and experimental investigation on transmission loss of clamped double panels: implication of boundary effects.

PubMed

Xin, F X; Lu, T J

2009-03-01

The air-borne sound insulation performance of a rectangular double-panel partition clamp mounted on an infinite acoustic rigid baffle is investigated both analytically and experimentally and compared with that of a simply supported one. With the clamped (or simply supported) boundary accounted for by using the method of modal function, a double series solution for the sound transmission loss (STL) of the structure is obtained by employing the weighted residual (Galerkin) method. Experimental measurements with Al double-panel partitions having air cavity are subsequently carried out to validate the theoretical model for both types of the boundary condition, and good overall agreement is achieved. A consistency check of the two different models (based separately on clamped modal function and simply supported modal function) is performed by extending the panel dimensions to infinite where no boundaries exist. The significant discrepancies between the two different boundary conditions are demonstrated in terms of the STL versus frequency plots as well as the panel deflection mode shapes.

Attractors for non-dissipative irrotational von Karman plates with boundary damping

NASA Astrophysics Data System (ADS)

Bociu, Lorena; Toundykov, Daniel

Long-time behavior of solutions to a von Karman plate equation is considered. The system has an unrestricted first-order perturbation and a nonlinear damping acting through free boundary conditions only. This model differs from those previously considered (e.g. in the extensive treatise (Chueshov and Lasiecka, 2010 [11])) because the semi-flow may be of a non-gradient type: the unique continuation property is not known to hold, and there is no strict Lyapunov function on the natural finite-energy space. Consequently, global bounds on the energy, let alone the existence of an absorbing ball, cannot be a priori inferred. Moreover, the free boundary conditions are not recognized by weak solutions and some helpful estimates available for clamped, hinged or simply-supported plates cannot be invoked. It is shown that this non-monotone flow can converge to a global compact attractor with the help of viscous boundary damping and appropriately structured restoring forces acting only on the boundary or its collar.

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

NASA Astrophysics Data System (ADS)

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

Displacement potential solution of a guided deep beam of composite materials under symmetric three-point bending

NASA Astrophysics Data System (ADS)

Rahman, M. Muzibur; Ahmad, S. Reaz

2017-12-01

An analytical investigation of elastic fields for a guided deep beam of orthotropic composite material having three point symmetric bending is carried out using displacement potential boundary modeling approach. Here, the formulation is developed as a single function of space variables defined in terms of displacement components, which has to satisfy the mixed type of boundary conditions. The relevant displacement and stress components are derived into infinite series using Fourier integral along with suitable polynomials coincided with boundary conditions. The results are presented mainly in the form of graphs and verified with finite element solutions using ANSYS. This study shows that the analytical and numerical solutions are in good agreement and thus enhances reliability of the displacement potential approach.

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

NASA Astrophysics Data System (ADS)

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.

What is Wrong with the Boundary Conditions in Column Tracer Tests

NASA Astrophysics Data System (ADS)

Zhan, H.

2007-12-01

Solute transport in a column is probably one of the most fundamental problems investigated in contaminant hydrology and soil physics because it serves as a benchmark for testing transport theories, for measuring dispersivities, etc. Despite its importance, there are still dispute and inconsistency on how to deal with the boundary conditions involved in such problems. The boundary condition could impose great influence upon transport in a column, particularly when the length of the column is relatively short, or the so-called Peclet number is not large. There are three types of boundary conditions to choose for transport in a column. Among these three types of boundary conditions, only the third-type boundary satisfies the mass balance requirement rigorously. The first type boundary, despite its frequent use in previous studies, could lead to serious mass balance problems. The most serious problem is on how to deal with the outlet boundary. Some studies have used a zero concentration gradient at the outlet (the so-called Danckwerts' boundary condition). This is named the model A. Another idea is to treat the finite length column as a part of an infinitely long column and to calculate the concentration at the outlet based on a formula developed for an infinitely long column. This is named the model B. The model A satisfies the mass balance requirement but was found to fit with the experimental data poorly. The model B does not satisfy the mass balance requirement, but usually agree well with the experimental data. So, the dilemma is: which model to choose? At present, most investigators prefer to choose the model B because of its close agreement with the experimental data, despite of its violation of the mass balance requirement. But the question is: why the model A, which satisfies the mass balance requirement, does not fit with the experimental data? It turns out that the advection-dispersion equation (ADE) that uses the Fick's first law to describe the hydrodynamic dispersion has some problems, particularly in the regions near the two boundaries. Taylor (1921) has pointed out that the dispersion coefficient varies linearly with time at the beginning and tends to its asymptotic, Fickian value after a travel time of a few correlation scales. Dagan and Bresler (1985) have further pointed out that the constant dispersivity is attained after the solute body has traveled tens of conductivity integral scales. For transport in a homogeneous column, the integral scale of the conductivity is probably around the pore scale or equivalent to the dispersivity value. Therefore, for a finite column whose length is not much greater than the dispersivity value, the transition zones in which solute transport is non-Fickian could consist of a significant portion of the column length. It is such non-Fickian transport in the column that is responsible for the failure of the model A. But still, why does the model B yield the right solution? There is no answer to this question based on a rigorous quantitative analysis yet. To resolve the dilemma, one must carry out a non-Fickian transport study to deal with the transition zones. It is my hypothesis that if the non-Fickian transport analysis succeeds, one will find that the mass balance requirement is indeed satisfied in the model B. Dagan and Bresler (1985) have pointed to the right direction, but a rigorous analysis has not followed. This is something interesting and worthwhile to investigate. REFERENCES CITED Dagan, G., and Bresler, E., 1985. Comment on ¡°Flux-averaged and volume-averaged concentration in continuum approaches to solute transport¡± by J.C. Parker and M.Th. van Genuchten. Water Resources Research, 21: 1299- 1300. Taylor, G.I., 1921. Diffusion by continuous movements. Proc. London Math Soc. 2: 196-212.

Mean Field Type Control with Congestion (II): An Augmented Lagrangian Method

DOE Office of Scientific and Technical Information (OSTI.GOV)

Achdou, Yves, E-mail: achdou@ljll.univ-paris-diderot.fr; Laurière, Mathieu

This work deals with a numerical method for solving a mean-field type control problem with congestion. It is the continuation of an article by the same authors, in which suitably defined weak solutions of the system of partial differential equations arising from the model were discussed and existence and uniqueness were proved. Here, the focus is put on numerical methods: a monotone finite difference scheme is proposed and shown to have a variational interpretation. Then an Alternating Direction Method of Multipliers for solving the variational problem is addressed. It is based on an augmented Lagrangian. Two kinds of boundary conditionsmore » are considered: periodic conditions and more realistic boundary conditions associated to state constrained problems. Various test cases and numerical results are presented.« less

A New Formulation of Time Domain Boundary Integral Equation for Acoustic Wave Scattering in the Presence of a Uniform Mean Flow

NASA Technical Reports Server (NTRS)

Hu, Fang; Pizzo, Michelle E.; Nark, Douglas M.

2017-01-01

It has been well-known that under the assumption of a constant uniform mean flow, the acoustic wave propagation equation can be formulated as a boundary integral equation, in both the time domain and the frequency domain. Compared with solving partial differential equations, numerical methods based on the boundary integral equation have the advantage of a reduced spatial dimension and, hence, requiring only a surface mesh. However, the constant uniform mean flow assumption, while convenient for formulating the integral equation, does not satisfy the solid wall boundary condition wherever the body surface is not aligned with the uniform mean flow. In this paper, we argue that the proper boundary condition for the acoustic wave should not have its normal velocity be zero everywhere on the solid surfaces, as has been applied in the literature. A careful study of the acoustic energy conservation equation is presented that shows such a boundary condition in fact leads to erroneous source or sink points on solid surfaces not aligned with the mean flow. A new solid wall boundary condition is proposed that conserves the acoustic energy and a new time domain boundary integral equation is derived. In addition to conserving the acoustic energy, another significant advantage of the new equation is that it is considerably simpler than previous formulations. In particular, tangential derivatives of the solution on the solid surfaces are no longer needed in the new formulation, which greatly simplifies numerical implementation. Furthermore, stabilization of the new integral equation by Burton-Miller type reformulation is presented. The stability of the new formulation is studied theoretically as well as numerically by an eigenvalue analysis. Numerical solutions are also presented that demonstrate the stability of the new formulation.

On the no-boundary proposal for ekpyrotic and cyclic cosmologies

DOE Office of Scientific and Technical Information (OSTI.GOV)

Battarra, Lorenzo; Lehners, Jean-Luc, E-mail: lorenzo.battarra@aei.mpg.de, E-mail: jlehners@aei.mpg.de

2014-12-01

The no-boundary proposal provides a compelling theory for the initial conditions of our universe. We study the implications of such initial conditions for ekpyrotic and cyclic cosmologies. These cosmologies allow for the existence of a new type of ''ekpyrotic instanton'', which describes the creation of a universe in the ekpyrotic contraction phase. Remarkably, we find that the ekpyrotic attractor can explain how the universe became classical. In a cyclic context, in addition to the ekpyrotic instantons there exist de Sitter-like instantons describing the emergence of the universe in the dark energy phase. Our results show that typically the ekpyrotic instantonsmore » yield a higher probability. In fact, in a potential energy landscape allowing both inflationary and cyclic cosmologies, the no-boundary proposal implies that the probability for ekpyrotic and cyclic initial conditions is vastly higher than that for inflationary ones.« less

Matched asymptotic expansion of the Hamilton-Jacobi-Bellman equation for aeroassisted plane-change maneuvers

NASA Technical Reports Server (NTRS)

Calise, Anthony J.; Melamed, Nahum

1993-01-01

In this paper we develop a general procedure for constructing a matched asymptotic expansion of the Hamilton-Jacobi-Bellman equation based on the method of characteristics. The development is for a class of perturbation problems whose solution exhibits two-time-scale behavior. A regular expansion for problems of this type is inappropriate since it is not uniformly valid over a narrow range of the independent variable. Of particular interest here is the manner in which matching and boundary conditions are enforced when the expansion is carried out to first order. Two cases are distinguished - one where the left boundary condition coincides with or lies to the right of the singular region and one where the left boundary condition lies to the left of the singular region. A simple example is used to illustrate the procedure, and its potential application to aeroassisted plane change is described.

A spectral approach for the stability analysis of turbulent open-channel flows over granular beds

NASA Astrophysics Data System (ADS)

Camporeale, C.; Canuto, C.; Ridolfi, L.

2012-01-01

A novel Orr-Sommerfeld-like equation for gravity-driven turbulent open-channel flows over a granular erodible bed is here derived, and the linear stability analysis is developed. The whole spectrum of eigenvalues and eigenvectors of the complete generalized eigenvalue problem is computed and analyzed. The fourth-order eigenvalue problem presents singular non-polynomial coefficients with non-homogenous Robin-type boundary conditions that involve first and second derivatives. Furthermore, the Exner condition is imposed at an internal point. We propose a numerical discretization of spectral type based on a single-domain Galerkin scheme. In order to manage the presence of singular coefficients, some properties of Jacobi polynomials have been carefully blended with numerical integration of Gauss-Legendre type. The results show a positive agreement with the classical experimental data and allow one to relate the different types of instability to such parameters as the Froude number, wavenumber, and the roughness scale. The eigenfunctions allow two types of boundary layers to be distinguished, scaling, respectively, with the roughness height and the saltation layer for the bedload sediment transport.

Evolution of passive movement in advective environments: General boundary condition

NASA Astrophysics Data System (ADS)

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.

Application of a low order panel method to complex three-dimensional internal flow problems

NASA Technical Reports Server (NTRS)

Ashby, D. L.; Sandlin, D. R.

1986-01-01

An evaluation of the ability of a low order panel method to predict complex three-dimensional internal flow fields was made. The computer code VSAERO was used as a basis for the evaluation. Guidelines for modeling internal flow geometries were determined and the effects of varying the boundary conditions and the use of numerical approximations on the solutions accuracy were studied. Several test cases were run and the results were compared with theoretical or experimental results. Modeling an internal flow geometry as a closed box with normal velocities specified on an inlet and exit face provided accurate results and gave the user control over the boundary conditions. The values of the boundary conditions greatly influenced the amount of leakage an internal flow geometry suffered and could be adjusted to eliminate leakage. The use of the far-field approximation to reduce computation time influenced the accuracy of a solution and was coupled with the values of the boundary conditions needed to eliminate leakage. The error induced in the influence coefficients by using the far-field approximation was found to be dependent on the type of influence coefficient, the far-field radius, and the aspect ratio of the panels.

The Boundary Function Method. Fundamentals

NASA Astrophysics Data System (ADS)

Kot, V. A.

2017-03-01

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

Emancipating traditional channel network types: quantification of topology and geometry, and relation to geologic boundary conditions

NASA Astrophysics Data System (ADS)

Temme, A.; Langston, A. L.

2017-12-01

Traditional classification of channel networks is helpful for qualitative geologic and geomorphic inference. For instance, a dendritic network indicates no strong lithological control on where channels flow. However, an approach where channel network structure is quantified, is required to be able to indicate for instance how increasing levels of lithological control lead, gradually or suddenly, to a trellis-type drainage network Our contribution aims to aid this transition to a quantitative analysis of channel networks. First, to establish the range of typically occurring channel network properties, we selected 30 examples of traditional drainage network types from around the world. For each of these, we calculated a set of topological and geometric properties, such as total drainage length, average length of a channel segment and the average angle of intersection of channel segments. A decision tree was used to formalize the relation between these newly quantified properties on the one hand, and traditional network types on the other hand. Then, to explore how variations in lithological and geomorphic boundary conditions affect channel network structure, we ran a set of experiments with landscape evolution model Landlab. For each simulated channel network, the same set of topological and geometric properties was calculated as for the 30 real-world channel networks. The latter were used for a first, visual evaluation to find out whether a simulated network that looked, for instance, rectangular, also had the same set of properties as real-world rectangular channel networks. Ultimately, the relation between these properties and the imposed lithological and geomorphic boundary conditions was explored using simple bivariate statistics.

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

NASA Astrophysics Data System (ADS)

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.

Analytical solutions with Generalized Impedance Boundary Conditions (GIBC)

NASA Technical Reports Server (NTRS)

Syed, H. H.; Volakis, John L.

1991-01-01

Rigorous uniform geometrical theory of diffraction (UTD) diffraction coefficients are presented for a coated convex cylinder simulated with generalized impedance boundary conditions. In particular, ray solutions are obtained which remain valid in the transition region and reduce uniformly to those in the deep lit and shadow regions. These involve new transition functions in place of the usual Fock-type integrals, characteristics to the impedance cylinder. A uniform asymptotic solution is also presented for observations in the close vicinity of the cylinder. The diffraction coefficients for the convex cylinder are obtained via a generalization of the corresponding ones for the circular cylinder.

On Raviart-Thomas and VMS formulations for flow in heterogeneous materials.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Turner, Daniel Zack

It is well known that the continuous Galerkin method (in its standard form) is not locally conservative, yet many stabilized methods are constructed by augmenting the standard Galerkin weak form. In particular, the Variational Multiscale (VMS) method has achieved popularity for combating numerical instabilities that arise for mixed formulations that do not otherwise satisfy the LBB condition. Among alternative methods that satisfy local and global conservation, many employ Raviart-Thomas function spaces. The lowest order Raviart-Thomas finite element formulation (RT0) consists of evaluating fluxes over the midpoint of element edges and constant pressures within the element. Although the RT0 element posesmore » many advantages, it has only been shown viable for triangular or tetrahedral elements (quadrilateral variants of this method do not pass the patch test). In the context of heterogenous materials, both of these methods have been used to model the mixed form of the Darcy equation. This work aims, in a comparative fashion, to evaluate the strengths and weaknesses of either approach for modeling Darcy flow for problems with highly varying material permeabilities and predominantly open flow boundary conditions. Such problems include carbon sequestration and enhanced oil recovery simulations for which the far-field boundary is typically described with some type of pressure boundary condition. We intend to show the degree to which the VMS formulation violates local mass conservation for these types of problems and compare the performance of the VMS and RT0 methods at boundaries between disparate permeabilities.« less

A discrete twin-boundary approach for simulating the magneto-mechanical response of Ni-Mn-Ga

NASA Astrophysics Data System (ADS)

Faran, Eilon; Shilo, Doron

2016-09-01

The design and optimization of ferromagnetic shape memory alloys (FSMA)-based devices require quantitative understanding of the dynamics of twin boundaries within these materials. Here, we present a discrete twin boundary modeling approach for simulating the behavior of an FSMA Ni-Mn-Ga crystal under combined magneto-mechanical loading conditions. The model is based on experimentally measured kinetic relations that describe the motion of individual twin boundaries over a wide range of velocities. The resulting calculations capture the dynamic response of Ni-Mn-Ga and reveal the relations between fundamental material parameters and actuation performance at different frequencies of the magnetic field. In particular, we show that at high field rates, the magnitude of the lattice barrier that resists twin boundary motion is the important property that determines the level of actuation strain, while the contribution of twinning stress property is minor. Consequently, type II twin boundaries, whose lattice barrier is smaller compared to type I, are expected to show better actuation performance at high rates, irrespective of the differences in the twinning stress property between the two boundary types. In addition, the simulation enables optimization of the actuation strain of a Ni-Mn-Ga crystal by adjusting the magnitude of the bias mechanical stress, thus providing direct guidelines for the design of actuating devices. Finally, we show that the use of a linear kinetic law for simulating the twinning-based response is inadequate and results in incorrect predictions.

The boundary integral theory for slow and rapid curved solid/liquid interfaces propagating into binary systems

NASA Astrophysics Data System (ADS)

Galenko, Peter K.; Alexandrov, Dmitri V.; Titova, Ekaterina A.

2018-01-01

The boundary integral method for propagating solid/liquid interfaces is detailed with allowance for the thermo-solutal Stefan-type models. Two types of mass transfer mechanisms corresponding to the local equilibrium (parabolic-type equation) and local non-equilibrium (hyperbolic-type equation) solidification conditions are considered. A unified integro-differential equation for the curved interface is derived. This equation contains the steady-state conditions of solidification as a special case. The boundary integral analysis demonstrates how to derive the quasi-stationary Ivantsov and Horvay-Cahn solutions that, respectively, define the paraboloidal and elliptical crystal shapes. In the limit of highest Péclet numbers, these quasi-stationary solutions describe the shape of the area around the dendritic tip in the form of a smooth sphere in the isotropic case and a deformed sphere along the directions of anisotropy strength in the anisotropic case. A thermo-solutal selection criterion of the quasi-stationary growth mode of dendrites which includes arbitrary Péclet numbers is obtained. To demonstrate the selection of patterns, computational modelling of the quasi-stationary growth of crystals in a binary mixture is carried out. The modelling makes it possible to obtain selected structures in the form of dendritic, fractal or planar crystals. This article is part of the theme issue `From atomistic interfaces to dendritic patterns'.

Sea Slugs, Subliminal Pictures, and Vegetative State Patients: Boundaries of Consciousness in Classical Conditioning

PubMed Central

Bekinschtein, Tristan A.; Peeters, Moos; Shalom, Diego; Sigman, Mariano

2011-01-01

Classical (trace) conditioning is a specific variant of associative learning in which a neutral stimulus leads to the subsequent prediction of an emotionally charged or noxious stimulus after a temporal gap. When conditioning is concurrent with a distraction task, only participants who can report the relationship (the contingency) between stimuli explicitly show associative learning. This suggests that consciousness is a prerequisite for trace conditioning. We review and question three main controversies concerning this view. Firstly, virtually all animals, even invertebrate sea slugs, show this type of learning; secondly, unconsciously perceived stimuli may elicit trace conditioning; and thirdly, some vegetative state patients show trace learning. We discuss and analyze these seemingly contradictory arguments to find the theoretical boundaries of consciousness in classical conditioning. We conclude that trace conditioning remains one of the best measures to test conscious processing in the absence of explicit reports. PMID:22164148

The charge conserving Poisson-Boltzmann equations: Existence, uniqueness, and maximum principle

DOE Office of Scientific and Technical Information (OSTI.GOV)

Lee, Chiun-Chang, E-mail: chlee@mail.nhcue.edu.tw

2014-05-15

The present article is concerned with the charge conserving Poisson-Boltzmann (CCPB) equation in high-dimensional bounded smooth domains. The CCPB equation is a Poisson-Boltzmann type of equation with nonlocal coefficients. First, under the Robin boundary condition, we get the existence of weak solutions to this equation. The main approach is variational, based on minimization of a logarithm-type energy functional. To deal with the regularity of weak solutions, we establish a maximum modulus estimate for the standard Poisson-Boltzmann (PB) equation to show that weak solutions of the CCPB equation are essentially bounded. Then the classical solutions follow from the elliptic regularity theorem.more » Second, a maximum principle for the CCPB equation is established. In particular, we show that in the case of global electroneutrality, the solution achieves both its maximum and minimum values at the boundary. However, in the case of global non-electroneutrality, the solution may attain its maximum value at an interior point. In addition, under certain conditions on the boundary, we show that the global non-electroneutrality implies pointwise non-electroneutrality.« less

Boundary lubrication of formulated C-ethers in air to 300 C. 2: Organic acid additives

NASA Technical Reports Server (NTRS)

Jones, W. R., Jr.

1973-01-01

Friction and wear measurements were made on CVM M-50 steel lubricated with three C-ether (modified polyphenyl ether) formulations in dry and moist air. Results were compared to those obtained with a formulated Type 2 ester and the C-ether base fluid. A ball-on-disk sliding friction apparatus was used. Experimental conditions were a 1-kilogram load, a 17-meter/minute surface speed, and a 25 to 300 C (77 to 572 F) disk temperature range. The three C-ether formulations yielded better boundary lubricating characteristics than the Type 2 ester under most test conditions. All C-ether formulations exhibited higher friction coefficients than the ester from 150 to 300 C (302 to 572 F) and similar or lower values from 25 to 150 C (77 to 302 F).

Intrinsic Compressive Stress in Polycrystalline Films is Localized at Edges of the Grain Boundaries.

PubMed

Vasco, Enrique; Polop, Celia

2017-12-22

The intrinsic compression that arises in polycrystalline thin films under high atomic mobility conditions has been attributed to the insertion or trapping of adatoms inside grain boundaries. This compression is a consequence of the stress field resulting from imperfections in the solid and causes the thermomechanical fatigue that is estimated to be responsible for 90% of mechanical failures in current devices. We directly measure the local distribution of residual intrinsic stress in polycrystalline thin films on nanometer scales, using a pioneering method based on atomic force microscopy. Our results demonstrate that, at odds with expectations, compression is not generated inside grain boundaries but at the edges of gaps where the boundaries intercept the surface. We describe a model wherein this compressive stress is caused by Mullins-type surface diffusion towards the boundaries, generating a kinetic surface profile different from the mechanical equilibrium profile by the Laplace-Young equation. Where the curvatures of both profiles differ, an intrinsic stress is generated in the form of Laplace pressure. The Srolovitz-type surface diffusion that results from the stress counters the Mullins-type diffusion and stabilizes the kinetic surface profile, giving rise to a steady compression regime. The proposed mechanism of competition between surface diffusions would explain the flux and time dependency of compressive stress in polycrystalline thin films.

Intrinsic Compressive Stress in Polycrystalline Films is Localized at Edges of the Grain Boundaries

NASA Astrophysics Data System (ADS)

Vasco, Enrique; Polop, Celia

2017-12-01

The intrinsic compression that arises in polycrystalline thin films under high atomic mobility conditions has been attributed to the insertion or trapping of adatoms inside grain boundaries. This compression is a consequence of the stress field resulting from imperfections in the solid and causes the thermomechanical fatigue that is estimated to be responsible for 90% of mechanical failures in current devices. We directly measure the local distribution of residual intrinsic stress in polycrystalline thin films on nanometer scales, using a pioneering method based on atomic force microscopy. Our results demonstrate that, at odds with expectations, compression is not generated inside grain boundaries but at the edges of gaps where the boundaries intercept the surface. We describe a model wherein this compressive stress is caused by Mullins-type surface diffusion towards the boundaries, generating a kinetic surface profile different from the mechanical equilibrium profile by the Laplace-Young equation. Where the curvatures of both profiles differ, an intrinsic stress is generated in the form of Laplace pressure. The Srolovitz-type surface diffusion that results from the stress counters the Mullins-type diffusion and stabilizes the kinetic surface profile, giving rise to a steady compression regime. The proposed mechanism of competition between surface diffusions would explain the flux and time dependency of compressive stress in polycrystalline thin films.

Finite volume solution of the compressible boundary-layer equations

NASA Technical Reports Server (NTRS)

Loyd, B.; Murman, E. M.

1986-01-01

A box-type finite volume discretization is applied to the integral form of the compressible boundary layer equations. Boundary layer scaling is introduced through the grid construction: streamwise grid lines follow eta = y/h = const., where y is the normal coordinate and h(x) is a scale factor proportional to the boundary layer thickness. With this grid, similarity can be applied explicity to calculate initial conditions. The finite volume method preserves the physical transparency of the integral equations in the discrete approximation. The resulting scheme is accurate, efficient, and conceptually simple. Computations for similar and non-similar flows show excellent agreement with tabulated results, solutions computed with Keller's Box scheme, and experimental data.

Correlation parameters for the study of leeside heating on a lifting body at hypersonic speeds

NASA Technical Reports Server (NTRS)

Vidal, R. J.

1974-01-01

Leeside heating was studied with the aim of gaining some insight into: (1) the magnitude of the leeside heating rates and (2) the methods to be used to extrapolate wind tunnel leeside heating rates to the full scale flight condition. This study was based on existing experimental data obtained in a hypersonic shock tunnel on lifting body configurations that are typical of shuttle orbiter vehicles. Heat transfer was first measured on the windward side to determine the boundary layer type. Then the leeside heating was investigated with the classified boundary layer. Correlation data are given on the windward turbulent boundary layer, the windward laminar boundary layer, and the leeside surfaces.

Automatic Control via Thermostats of a Hyperbolic Stefan Problem with Memory

DOE Office of Scientific and Technical Information (OSTI.GOV)

Colli, P.; Grasselli, M.; Sprekels, J.

1999-03-15

A hyperbolic Stefan problem based on the linearized Gurtin-Pipkin heat conduction law is considered. The temperature and free boundary are controlled by a thermostat acting on the boundary. This feedback control is based on temperature measurements performed by real thermal sensors located within the domain containing the two-phase system and/or at its boundary. Three different types of thermostats are analyzed: simple switch, relay switch, and a Preisach hysteresis operator. The resulting models lead to integrodifferential hyperbolic Stefan problems with nonlinear and nonlocal boundary conditions. Existence results are proved in all the cases. Uniqueness is also shown, except in the situationmore » corresponding to the ideal switch.« less

Heteroaggregation of oppositely charged particles in the presence of multivalent ions.

PubMed

Cao, Tianchi; Sugimoto, Takuya; Szilagyi, Istvan; Trefalt, Gregor; Borkovec, Michal

2017-06-14

Time-resolved dynamic light scattering is used to measure absolute heteroaggregation rate coefficients and the corresponding stability ratios for heteroaggregation between amidine and sulfate latex particles. These measurements are complemented by the respective quantities for the homoaggregation of the two systems and electrophoresis. Based on the latter measurements, the stability ratios are calculated using Derjaguin-Landau-Verwey-Overbeek (DLVO) theory. In monovalent salt solutions, the two types of particles investigated are oppositely charged. In the presence of multivalent ions, however, one particle type reverses its charge, while the charge of the other particle type is hardly affected. In this region, the heteroaggregation stability ratio goes through a pronounced maximum when plotted versus concentration. This region of slow aggregation is wider than the one observed in the corresponding homoaggregation process. One also finds that the onset of this region sensitively depends on the boundary conditions used to calculate the double layer force. The present results are more in line with constant potential boundary conditions.

Three-Dimensional Velocity Field De-Noising using Modal Projection

NASA Astrophysics Data System (ADS)

Frank, Sarah; Ameli, Siavash; Szeri, Andrew; Shadden, Shawn

2017-11-01

PCMRI and Doppler ultrasound are common modalities for imaging velocity fields inside the body (e.g. blood, air, etc) and PCMRI is increasingly being used for other fluid mechanics applications where optical imaging is difficult. This type of imaging is typically applied to internal flows, which are strongly influenced by domain geometry. While these technologies are evolving, it remains that measured data is noisy and boundary layers are poorly resolved. We have developed a boundary modal analysis method to de-noise 3D velocity fields such that the resulting field is divergence-free and satisfies no-slip/no-penetration boundary conditions. First, two sets of divergence-free modes are computed based on domain geometry. The first set accounts for flow through ``truncation boundaries'', and the second set of modes has no-slip/no-penetration conditions imposed on all boundaries. The modes are calculated by minimizing the velocity gradient throughout the domain while enforcing a divergence-free condition. The measured velocity field is then projected onto these modes using a least squares algorithm. This method is demonstrated on CFD simulations with artificial noise. Different degrees of noise and different numbers of modes are tested to reveal the capabilities of the approach. American Heart Association Award 17PRE33660202.

A computer program for the simulation of heat and moisture flow in soils

NASA Technical Reports Server (NTRS)

Camillo, P.; Schmugge, T. J.

1981-01-01

A computer program that simulates the flow of heat and moisture in soils is described. The space-time dependence of temperature and moisture content is described by a set of diffusion-type partial differential equations. The simulator uses a predictor/corrector to numerically integrate them, giving wetness and temperature profiles as a function of time. The simulator was used to generate solutions to diffusion-type partial differential equations for which analytical solutions are known. These equations include both constant and variable diffusivities, and both flux and constant concentration boundary conditions. In all cases, the simulated and analytic solutions agreed to within the error bounds which were imposed on the integrator. Simulations of heat and moisture flow under actual field conditions were also performed. Ground truth data were used for the boundary conditions and soil transport properties. The qualitative agreement between simulated and measured profiles is an indication that the model equations are reasonably accurate representations of the physical processes involved.

Homogenization models for thin rigid structured surfaces and films.

PubMed

Marigo, Jean-Jacques; Maurel, Agnès

2016-07-01

A homogenization method for thin microstructured surfaces and films is presented. In both cases, sound hard materials are considered, associated with Neumann boundary conditions and the wave equation in the time domain is examined. For a structured surface, a boundary condition is obtained on an equivalent flat wall, which links the acoustic velocity to its normal and tangential derivatives (of the Myers type). For a structured film, jump conditions are obtained for the acoustic pressure and the normal velocity across an equivalent interface (of the Ventcels type). This interface homogenization is based on a matched asymptotic expansion technique, and differs slightly from the classical homogenization, which is known to fail for small structuration thicknesses. In order to get insight into what causes this failure, a two-step homogenization is proposed, mixing classical homogenization and matched asymptotic expansion. Results of the two homogenizations are analyzed in light of the associated elementary problems, which correspond to problems of fluid mechanics, namely, potential flows around rigid obstacles.

Smoothed Particle Hydrodynamics Continuous Boundary Force method for Navier-Stokes equations subject to Robin boundary condition

DOE Office of Scientific and Technical Information (OSTI.GOV)

Pan, Wenxiao; Bao, Jie; Tartakovsky, Alexandre M.

2014-02-15

Robin boundary condition for the Navier-Stokes equations is used to model slip conditions at the fluid-solid boundaries. A novel Continuous Boundary Force (CBF) method is proposed for solving the Navier-Stokes equations subject to Robin boundary condition. In the CBF method, the Robin boundary condition at boundary is replaced by the homogeneous Neumann boundary condition at the boundary and a volumetric force term added to the momentum conservation equation. Smoothed Particle Hydrodynamics (SPH) method is used to solve the resulting Navier-Stokes equations. We present solutions for two-dimensional and three-dimensional flows in domains bounded by flat and curved boundaries subject to variousmore » forms of the Robin boundary condition. The numerical accuracy and convergence are examined through comparison of the SPH-CBF results with the solutions of finite difference or finite element method. Taken the no-slip boundary condition as a special case of slip boundary condition, we demonstrate that the SPH-CBF method describes accurately both no-slip and slip conditions.« less

Strain-Annealing Based Grain Boundary Engineering to Evaluate its Sole Implication on Intergranular Corrosion in Extra-Low Carbon Type 304L Austenitic Stainless Steel

NASA Astrophysics Data System (ADS)

Pradhan, S. K.; Bhuyan, P.; Kaithwas, C.; Mandal, Sumantra

2018-05-01

Strain-annealing based thermo-mechanical processing has been performed to promote grain boundary engineering (GBE) in an extra-low carbon type austenitic stainless steel without altering the grain size and residual strain to evaluate its sole influence on intergranular corrosion. Single-step processing comprising low pre-strain ( 5 and 10 pct) followed by annealing at 1273 K for 1 hour have resulted in a large fraction of Σ3 n boundaries and significant disruption in random high-angle grain boundaries (RHAGBs) connectivity. This is due to the occurrence of prolific multiple twinning in these specimens as confirmed by their large twin-related domain and twin-related grain size ratio. Among the iterative processing, the schedule comprising two cycles of 10 and 5 pct deformation followed by annealing at 1173 K for 1 hour has yielded the optimum GBE microstructure with the grain size and residual strain akin to the as-received condition. The specimens subjected to the higher number of iterations failed to realize GBE microstructures due to the occurrence of partial recrystallization. Owing to the optimum grain boundary character distribution, the GBE specimen has exhibited remarkable resistance against sensitization and intergranular corrosion as compared to the as-received condition. Furthermore, the lower depth of percolation in the GBE specimen is due to the significant disruption of RHAGBs connectivity as confirmed from its large twin-related domain and lower fractal dimension.

Strain-Annealing Based Grain Boundary Engineering to Evaluate its Sole Implication on Intergranular Corrosion in Extra-Low Carbon Type 304L Austenitic Stainless Steel

NASA Astrophysics Data System (ADS)

Pradhan, S. K.; Bhuyan, P.; Kaithwas, C.; Mandal, Sumantra

2018-07-01

Strain-annealing based thermo-mechanical processing has been performed to promote grain boundary engineering (GBE) in an extra-low carbon type austenitic stainless steel without altering the grain size and residual strain to evaluate its sole influence on intergranular corrosion. Single-step processing comprising low pre-strain ( 5 and 10 pct) followed by annealing at 1273 K for 1 hour have resulted in a large fraction of Σ3 n boundaries and significant disruption in random high-angle grain boundaries (RHAGBs) connectivity. This is due to the occurrence of prolific multiple twinning in these specimens as confirmed by their large twin-related domain and twin-related grain size ratio. Among the iterative processing, the schedule comprising two cycles of 10 and 5 pct deformation followed by annealing at 1173 K for 1 hour has yielded the optimum GBE microstructure with the grain size and residual strain akin to the as-received condition. The specimens subjected to the higher number of iterations failed to realize GBE microstructures due to the occurrence of partial recrystallization. Owing to the optimum grain boundary character distribution, the GBE specimen has exhibited remarkable resistance against sensitization and intergranular corrosion as compared to the as-received condition. Furthermore, the lower depth of percolation in the GBE specimen is due to the significant disruption of RHAGBs connectivity as confirmed from its large twin-related domain and lower fractal dimension.

Numerical modeling of coronal mass ejections based on various pre-event model atmospheres

NASA Technical Reports Server (NTRS)

Suess, S. T.; Wang, A. H.; Wu, S. T.; Poletto, G.

1994-01-01

We examine how the initial state (pre-event corona) affects the numerical MHD simulation for a coronal mass ejection (CME). Earlier simulations based on a pre-event corona with a homogeneous density and temperature distribution at lower boundary (i.e. solar surface) have been used to analyze the role of streamer properties in determining the characteristics of loop-like transients. The present paper extends these studies to show how a broader class of global coronal properties leads not only to different types of CME's, but also modifies the adjacent quiet corona and/or coronal holes. We consider four pre-event coronal cases: (1) Constant boundary conditions and a polytropic gas with gamma = 1.05; (2) Non-constant (latitude dependent) boundary conditions and a polytropic gas with gamma = 1.05; (3) Constant boundary conditions with a volumetric energy source and gamma = 1.67; (4) Non-constant (latitude dependent) boundary conditions with a volumetric energy source and gamma = 1.67. In all models, the pre-event magnetic fields separate the corona into closed field regions (streamers) and open field regions. The CME's initiation is simulated by introducing at the base of the corona, within the streamer region, a standard pressure pulse and velocity change. Boundary values are determined using MHD characteristic theory. The simulations show how different CME's, including loop-like transients, clouds, and bright rays, might occur. There are significant new features in comparison to published results. We conclude that the pre-event corona is a crucial factor in dictating CME's properties.

Numerical Modeling of Coronal Mass Ejections Based on Various Pre-event Model Atmospheres

NASA Technical Reports Server (NTRS)

Wang, A. H.; Wu, S. T.; Suess, S. T.; Poletto, G.

1995-01-01

We examine how the initial state (pre-event corona) affects the numerical MHD simulation for a coronal mass ejection (CME). Earlier simulations based on a pre-event corona with a homogeneous density and temperature distribution, at the lower boundary (i.e., solar surface) have been used to analyze the role of streamer properties in determining the characteristics of loop-like transients. The present paper extends these studies to show how a broader class of global coronal properties leads not only to different types of CME's, but also modifies the adjacent quiet corona and/or coronal holes. We consider four pre-event coronal cases: (1) constant boundary conditions and a polytropic gas with gamma = 1.05; (2) non-constant (latitude dependent) boundary conditions and a polytropic gas with gamma = 1.05; (3) constant boundary conditions with a volumetric energy source and gamma = 1.67; (4) non-constant (latitude dependent) boundary conditions with a volumetric energy source and gamma = 1.67. In all models, the pre-event magnetic fields separate the corona into closed field regions (streamers) and open field regions. The CME's initiation is simulated by introducing at the base of the corona, within the streamer region, a standard pressure pulse and velocity change. Boundary values are determined using magnetohydrodynamic (MHD) characteristic theory. The simulations show how different CME's, including loop-like transients, clouds and bright rays, might occur. There are significant new features in comparison to published results. We conclude that the pre-event corona is a crucial factor in dictating CME's properties.

Development of tf coil support concepts by design methodology in the case of a Bitter-type magnet. [Bitter-type magnets

DOE Office of Scientific and Technical Information (OSTI.GOV)

Brossmann, U.B.

1981-01-01

The application of the methodological design is demonstrated for the development of support concepts in the case of a Bitter-type magnet designed for a compact tokamak experimentat aiming at ignition of a DT plasma. With this methodology all boundary conditions and design criteria are more easily satisfied in a technical and economical way.

Transport phenomena in the micropores of plug-type phase separators

NASA Technical Reports Server (NTRS)

Fazah, M. M.

1995-01-01

This study numerically investigates the transport phenomena within and across a porous-plug phase separator. The effect of temperature differential across a single pore and of the sidewall boundary conditions, i.e., isothermal or linear thermal gradient, are presented and discussed. The effects are quantified in terms of the evaporation mass flux across the boundary and the mean surface temperature. A two-dimensional finite element model is used to solve the continuity, momentum, and energy equations for the liquid. Temperature differentials across the pore interface of 1.0, and 1.5 K are examined and their effect on evaporation flux and mean surface temperature is shown. For isothermal side boundary conditions, the evaporation flux across the pore is directly proportional and linear with Delta T. For the case of an imposed linear thermal gradient on the side boundaries, Biot numbers of 0.0, 0.15, and 0.5 are examined. The most significant effect of Biot number is to lower the overall surface temperature and evaporation flux.

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

PubMed Central

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

2016-01-01

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

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

PubMed

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

2015-08-15

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

Linear theory of boundary effects in open wind tunnels with finite jet lengths

NASA Technical Reports Server (NTRS)

Katzoff, S; Gardner, Clifford S; Diesendruck, Leo; Eisenstadt, Bertram J

1950-01-01

In the first part, the boundary conditions for an open wind tunnel (incompressible flow) are examined with special reference to the effects of the closed entrance and exit sections. Basic conditions are that the velocity must be continuous at the entrance lip and that the velocities in the upstream and downstream closed portions must be equal. In the second part, solutions are derived for four types of two-dimensional open tunnels, including one in which the pressures on the two free surfaces are not equal. Numerical results are given for every case. In general, if the lifting element is more than half the tunnel height from the inlet, the boundary effect at the lifting element is the same as for an infinitely long open tunnel. In the third part, a general method is given for calculating the boundary effect in an open circular wind tunnel of finite jet length. Numerical results are given for a lifting element concentrate at a point on the axis.

Uniqueness of boundary blow-up solutions on exterior domain of RN

NASA Astrophysics Data System (ADS)

Dong, Wei; Pang, Changci

2007-06-01

In this paper, we consider the existence and uniqueness of positive solutions of the degenerate logistic type elliptic equation where N[greater-or-equal, slanted]2, D[subset of]RN is a bounded domain with smooth boundary and a(x), b(x) are continuous functions on RN with b(x)[greater-or-equal, slanted]0, b(x)[not identical with]0. We show that under rather general conditions on a(x) and b(x) for large x, there exists a unique positive solution. Our results improve the corresponding ones in [W. Dong, Y. Du, Unbounded principal eigenfunctions and the logistic equation on RN, Bull. Austral. Math. Soc. 67 (2003) 413-427] and [Y. Du, L. Ma, Logistic type equations on RN by a squeezing method involving boundary blow-up solutions, J. London Math. Soc. (2) 64 (2001) 107-124].

Effect of initial densities in the lattice Boltzmann model for non-ideal fluid with curved interface

NASA Astrophysics Data System (ADS)

Gong, Jiaming; Oshima, Nobuyuki

2017-06-01

The effect of initial densities in a free energy based two-phase-flow lattice Boltzmann method for non-ideal fluids with a curved interface was investigated in the present work. To investigate this effect, the initial densities in the liquid and gas phases coming from the saturation points and the equilibrium state were adopted in the simulation of a static droplet in an open and a closed system. For the purpose of simplicity and easier comparison, the closed system is fabricated by the implementation of the periodic boundary condition at the inlet and outlet of a gas channel, and the open system is fabricated by the implementation of a constant flux boundary condition at the inlet and a free-out boundary condition at the outlet of the same gas channel. By comparing the simulation results from the two types of initial densities in the open and closed systems, it is proven that the commonly used saturation initial densities setting is the reason for droplet mass and volume variation which occurred in the simulation, particularly in the open system with a constant flux boundary condition. Such problems are believed to come from the curvature effect of the surface tension and can be greatly reduced by adopting the initial densities in the two phases from equilibrium state.

Inflow/Outflow Boundary Conditions with Application to FUN3D

NASA Technical Reports Server (NTRS)

Carlson, Jan-Renee

2011-01-01

Several boundary conditions that allow subsonic and supersonic flow into and out of the computational domain are discussed. These boundary conditions are demonstrated in the FUN3D computational fluid dynamics (CFD) code which solves the three-dimensional Navier-Stokes equations on unstructured computational meshes. The boundary conditions are enforced through determination of the flux contribution at the boundary to the solution residual. The boundary conditions are implemented in an implicit form where the Jacobian contribution of the boundary condition is included and is exact. All of the flows are governed by the calorically perfect gas thermodynamic equations. Three problems are used to assess these boundary conditions. Solution residual convergence to machine zero precision occurred for all cases. The converged solution boundary state is compared with the requested boundary state for several levels of mesh densities. The boundary values converged to the requested boundary condition with approximately second-order accuracy for all of the cases.

Finite size effects in ferroelectric-semiconductor thin films under open-circuit electric boundary conditions

DOE PAGES

Eliseev, Eugene A.; Kalinin, Sergei V.; Morozovska, Anna N.

2015-01-21

General features of finite size effects in the ferroelectric-semiconductor film under open-circuit electric boundary conditions are analyzed using Landau-Ginzburg-Devonshire theory and continuum media electrostatics. The temperature dependence of the film critical thickness, spontaneous polarization and depolarization field profiles of the open-circuited films are found to be significantly different from the characteristics of short-circuited ones. In particular, we predict the re-entrant type transition boundary between the mono-domain and poly-domain ferroelectric states due to reduced internal screening efficiency and analyzed possible experimental scenarios created by this mechanism. Performed analysis is relevant for the quantitative description of free-standing ferroelectric films phase diagrams andmore » polar properties. Also our results can be useful for the explanation of the scanning-probe microscopy experiments on free ferroelectric surfaces.« less

Hilbert complexes of nonlinear elasticity

NASA Astrophysics Data System (ADS)

Angoshtari, Arzhang; Yavari, Arash

2016-12-01

We introduce some Hilbert complexes involving second-order tensors on flat compact manifolds with boundary that describe the kinematics and the kinetics of motion in nonlinear elasticity. We then use the general framework of Hilbert complexes to write Hodge-type and Helmholtz-type orthogonal decompositions for second-order tensors. As some applications of these decompositions in nonlinear elasticity, we study the strain compatibility equations of linear and nonlinear elasticity in the presence of Dirichlet boundary conditions and the existence of stress functions on non-contractible bodies. As an application of these Hilbert complexes in computational mechanics, we briefly discuss the derivation of a new class of mixed finite element methods for nonlinear elasticity.

On the protein crystal formation as an interface-controlled process with prototype ion-channeling effect.

PubMed

Siódmiak, Jacek; Uher, Jan J; Santamaría-Holek, Ivan; Kruszewska, Natalia; Gadomski, Adam

2007-08-01

A superdiffusive random-walk action in the depletion zone around a growing protein crystal is considered. It stands for a dynamic boundary condition of the growth process and competes steadily with a quasistatic, curvature-involving (thermodynamic) free boundary condition, both of them contributing to interpret the (mainly late-stage) growth process in terms of a prototype ion-channeling effect. An overall diffusion function contains quantitative signatures of both boundary conditions mentioned and indicates whether the new phase grows as an orderly phase or a converse scenario occurs. This situation can be treated in a quite versatile way both numerically and analytically, within a generalized Smoluchowski framework. This study can help in (1) elucidating some dynamic puzzles of a complex crystal formation vs biomolecular aggregation, also those concerning ion-channel formation, and (2) seeing how ion-channel-type dynamics of non-Markovian nature may set properly the pace of model (dis)ordered protein aggregation.

Effect of Stress Relief Annealing on Microstructure & Mechanical Properties of Welded Joints Between Low Alloy Carbon Steel and Stainless Steel

NASA Astrophysics Data System (ADS)

Nivas, R.; Das, G.; Das, S. K.; Mahato, B.; Kumar, S.; Sivaprasad, K.; Singh, P. K.; Ghosh, M.

2017-01-01

Two types of welded joints were prepared using low alloy carbon steel and austenitic stainless steel as base materials. In one variety, buttering material and weld metal were Inconel 82. In another type, buttering material and weld metal were Inconel 182. In case of Inconel 82, method of welding was GTAW. For Inconel 182, welding was done by SMAW technique. For one set of each joints after buttering, stress relief annealing was done at 923 K (650 °C) for 90 minutes before further joining with weld metal. Microstructural investigation and sub-size in situ tensile testing in scanning electron microscope were carried out for buttered-welded and buttered-stress relieved-welded specimens. Adjacent to fusion boundary, heat-affected zone of low alloy steel consisted of ferrite-pearlite phase combination. Immediately after fusion boundary in low alloy steel side, there was increase in matrix grain size. Same trend was observed in the region of austenitic stainless steel that was close to fusion boundary between weld metal-stainless steel. Close to interface between low alloy steel-buttering material, the region contained martensite, Type-I boundary and Type-II boundary. Peak hardness was obtained close to fusion boundary between low alloy steel and buttering material. In this respect, a minimum hardness was observed within buttering material. The peak hardness was shifted toward buttering material after stress relief annealing. During tensile testing no deformation occurred within low alloy steel and failure was completely through buttering material. Crack initiated near fusion boundary between low alloy steel-buttering material for welded specimens and the same shifted away from fusion boundary for stress relieved annealed specimens. This observation was at par with the characteristics of microhardness profile. In as welded condition, joints fabricated with Inconel 82 exhibited superior bond strength than the weld produced with Inconel 182. Stress relief annealing reduced the strength of transition joints and the reduction was maximum for specimen welded with Inconel 82.

Recharge and groundwater models: An overview

USGS Publications Warehouse

Sanford, W.

2002-01-01

Recharge is a fundamental component of groundwater systems, and in groundwater-modeling exercises recharge is either measured and specified or estimated during model calibration. The most appropriate way to represent recharge in a groundwater model depends upon both physical factors and study objectives. Where the water table is close to the land surface, as in humid climates or regions with low topographic relief, a constant-head boundary condition is used. Conversely, where the water table is relatively deep, as in drier climates or regions with high relief, a specified-flux boundary condition is used. In most modeling applications, mixed-type conditions are more effective, or a combination of the different types can be used. The relative distribution of recharge can be estimated from water-level data only, but flux observations must be incorporated in order to estimate rates of recharge. Flux measurements are based on either Darcian velocities (e.g., stream base-flow) or seepage velocities (e.g., groundwater age). In order to estimate the effective porosity independently, both types of flux measurements must be available. Recharge is often estimated more efficiently when automated inverse techniques are used. Other important applications are the delineation of areas contributing recharge to wells and the estimation of paleorecharge rates using carbon-14.

Phase separation in the six-vertex model with a variety of boundary conditions

NASA Astrophysics Data System (ADS)

Lyberg, I.; Korepin, V.; Ribeiro, G. A. P.; Viti, J.

2018-05-01

We present numerical results for the six-vertex model with a variety of boundary conditions. Adapting an algorithm for domain wall boundary conditions, proposed in the work of Allison and Reshetikhin [Ann. Inst. Fourier 55(6), 1847-1869 (2005)], we examine some modifications of these boundary conditions. To be precise, we discuss partial domain wall boundary conditions, reflecting ends, and half turn boundary conditions (domain wall boundary conditions with half turn symmetry). Dedicated to the memory of Ludwig Faddeev

Matrix form of Legendre polynomials for solving linear integro-differential equations of high order

NASA Astrophysics Data System (ADS)

Kammuji, M.; Eshkuvatov, Z. K.; Yunus, Arif A. M.

2017-04-01

This paper presents an effective approximate solution of high order of Fredholm-Volterra integro-differential equations (FVIDEs) with boundary condition. Legendre truncated series is used as a basis functions to estimate the unknown function. Matrix operation of Legendre polynomials is used to transform FVIDEs with boundary conditions into matrix equation of Fredholm-Volterra type. Gauss Legendre quadrature formula and collocation method are applied to transfer the matrix equation into system of linear algebraic equations. The latter equation is solved by Gauss elimination method. The accuracy and validity of this method are discussed by solving two numerical examples and comparisons with wavelet and methods.

High frequency scattering by a smooth coated cylinder simulated with generalized impedance boundary conditions

NASA Technical Reports Server (NTRS)

Syed, Hasnain H.; Volakis, John L.

1991-01-01

Rigorous uniform geometrical theory of diffraction (UGTD) diffraction coefficients are presented for a coated convex cylinder simulated with generalized impedance boundary conditions. In particular, ray solutions are obtained which remain valid in the transition region and reduce uniformly to those in the deep lit and shadow regions. These involve new transition functions in place of the usual Fock-type integrals, characteristic to the impedance cylinder. A uniform asymptotic solution is also presented for observations in the close vicinity of the cylinder. As usual, the diffraction coefficients for the convex cylinder are obtained via a generalization of the corresponding ones for the circular cylinder.

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

NASA Astrophysics Data System (ADS)

Vaibhav, V.

2011-04-01

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

Electromagnetic MUSIC-type imaging of perfectly conducting, arc-like cracks at single frequency

NASA Astrophysics Data System (ADS)

Park, Won-Kwang; Lesselier, Dominique

2009-11-01

We propose a non-iterative MUSIC (MUltiple SIgnal Classification)-type algorithm for the time-harmonic electromagnetic imaging of one or more perfectly conducting, arc-like cracks found within a homogeneous space R2. The algorithm is based on a factorization of the Multi-Static Response (MSR) matrix collected in the far-field at a single, nonzero frequency in either Transverse Magnetic (TM) mode (Dirichlet boundary condition) or Transverse Electric (TE) mode (Neumann boundary condition), followed by the calculation of a MUSIC cost functional expected to exhibit peaks along the crack curves each half a wavelength. Numerical experimentation from exact, noiseless and noisy data shows that this is indeed the case and that the proposed algorithm behaves in robust manner, with better results in the TM mode than in the TE mode for which one would have to estimate the normal to the crack to get the most optimal results.

Subauroral polarization stream on the outer boundary of the ring current during an energetic ion injection event

NASA Astrophysics Data System (ADS)

Yuan, Zhigang; Qiao, Zheng; Li, Haimeng; Huang, Shiyong; Wang, Dedong; Yu, Xiongdong; Yu, Tao

2017-04-01

Subauroral polarization stream (SAPS) electric field can play an important role in the coupling between the inner magnetosphere and ionosphere; however, the production mechanism of SAPS has not been yet solved. During an energetic ion injection event on 26 March 2004, at latitudes lower than the equatorward boundaries of precipitating plasma sheet electrons and ions, the Defense Meteorological Satellite Program (DMSP) F13 satellite simultaneously observed a strong SAPS with the peak velocity of 1294 m/s and downward flowing field-aligned currents (FACs). Conjugate observations of DMSP F13 and NOAA 15 satellites have shown that FACs flowing into the ionosphere just lie in the outer boundary of the ring current (RC). The downward flowing FACs were observed in a region of positive latitudinal gradients of the ion energy density, implying that the downward flowing FACs are more likely linked to the azimuthal gradient than the radial gradient of the RC ion pressure. Our result demonstrates that RC ion pressure gradients on the outer boundary of the RC in the evening sector during energetic ion injection events can lead to downward flowing FACs so as to cause strong SAPS in condition of low ionospheric conductivities.