3-D inelastic analysis methods for hot section components. Volume 2: Advanced special functions models

This Annual Status Report presents the results of work performed during the third year of the 3-D Inelastic Analysis Methods for Hot Sections Components program (NASA Contract NAS3-23697). The objective of the program is to produce a series of computer codes that permit more accurate and efficient three-dimensional analyses of selected hot section components, i.e., combustor liners, turbine blades, and turbine vanes. The computer codes embody a progression of mathematical models and are streamlined to take advantage of geometrical features, loading conditions, and forms of material response that distinguish each group of selected components

Gravity-driven draining of a thin rivulet with constant width down a slowly varying substrate

The locally unidirectional gravity-driven draining of a thin rivulet with constant width but slowly varying contact angle down a slowly varying substrate is considered. Specifically, the flow of a rivulet in the azimuthal direction from the top to the bottom of a large horizontal cylinder is investigated. In particular, it is shown that, despite behaving the same locally, this flow has qualitatively different global behaviour from that of a rivulet with constant contact angle but slowly varying width. For example, whereas in the case of constant contact angle there is always a rivulet that runs all the way from the top to the bottom of the cylinder, in the case of constant width this is possible only for sufficiently narrow rivulets. Wider rivulets with constant width are possible only between the top of the cylinder and a critical azimuthal angle on the lower half of the cylinder. Assuming that the contact lines de-pin at this critical angle (where the contact angle is zero) the rivulet runs from the critical angle to the bottom of the cylinder with zero contact angle, monotonically decreasing width and monotonically increasing maximum thickness. The total mass of fluid on the cylinder is found to be a monotonically increasing function of the value of the constant width

The Stokes boundary layer for a thixotropic or antithixotropic fluid

We present a mathematical investigation of the oscillatory boundary layer (‘Stokes layer’) in a semi-infinite fluid bounded by an oscillating wall (the socalled ‘Stokes problem’), when the fluid has a thixotropic or antithixotropic rheology. We obtain asymptotic solutions in the limit of small-amplitude oscillations, and we use numerical integration to validate the asymptotic solutions and to explore the behaviour of the system for larger-amplitude oscillations. The solutions that we obtain differ significantly from the classical solution for a Newtonian fluid. In particular, for antithixotropic fluids the velocity reaches zero at a finite distance from the wall, in contrast to the exponential decay for a thixotropic or a Newtonian fluid. For small amplitudes of oscillation, three regimes of behaviour are possible: the structure parameter may take values defined instantaneously by the shear rate, or by a long-term average; or it may behave hysteretically. The regime boundaries depend on the precise specification of structure build-up and breakdown rates in the rheological model, illustrating the subtleties of complex fluid models in non-rheometric settings. For larger amplitudes of oscillation the dominant behaviour is hysteretic. We discuss in particular the relationship between the shear stress and the shear rate at the oscillating wall

Quantum Random Walks do not need a Coin Toss

Classical randomized algorithms use a coin toss instruction to explore different evolutionary branches of a problem. Quantum algorithms, on the other hand, can explore multiple evolutionary branches by mere superposition of states. Discrete quantum random walks, studied in the literature, have nonetheless used both superposition and a quantum coin toss instruction. This is not necessary, and a discrete quantum random walk without a quantum coin toss instruction is defined and analyzed here. Our construction eliminates quantum entanglement from the algorithm, and the results match those obtained with a quantum coin toss instruction.Comment: 6 pages, 4 figures, RevTeX (v2) Expanded to include relation to quantum walk with a coin. Connection with Dirac equation pointed out. Version to be published in Phys. Rev.

Three-dimensional coating and rimming flow: a ring of fluid on a rotating horizontal cylinder

The steady three-dimensional flow of a thin, slowly varying ring of Newtonian fluid on either the outside or the inside of a uniformly rotating large horizontal cylinder is investigated. Specifically, we study “full-ring” solutions, corresponding to a ring of continuous, finite and non-zero thickness that extends all the way around the cylinder. In particular, it is found that there is a critical solution corresponding to either a critical load above which no full-ring solution exists (if the rotation speed is prescribed) or a critical rotation speed below which no full-ring solution exists (if the load is prescribed). We describe the behaviour of the critical solution and, in particular, show that the critical flux, the critical load, the critical semi-width and the critical ring profile are all increasing functions of the rotation speed. In the limit of small rotation speed, the critical flux is small and the critical ring is narrow and thin, leading to a small critical load. In the limit of large rotation speed, the critical flux is large and the critical ring is wide on the upper half of the cylinder and thick on the lower half of the cylinder, leading to a large critical load.\ud \ud We also describe the behaviour of the non-critical full-ring solution, and, in particular, show that the semi-width and the ring profile are increasing functions of the load but, in general, non-monotonic functions of the rotation speed. In the limit of large rotation speed, the ring approaches a limiting non-uniform shape, whereas in the limit of small load, the ring is narrow and thin with a uniform parabolic profile. Finally, we show that, while for most values of the rotation speed and the load the azimuthal velocity is in the same direction as the rotation of the cylinder, there is a region of parameter space close to the critical solution for sufficiently small rotation speed in which backflow occurs in a small region on the right-hand side of the cylinder

Thermoviscous Coating and Rimming Flow

A comprehensive description is obtained of steady thermoviscous (i.e. with temperature-dependent viscosity) coating and rimming flow on a uniformly rotating horizontal cylinder that is uniformly hotter or colder than the surrounding atmosphere. It is found that, as in the corresponding isothermal problem, there is a critical solution with a corresponding critical load (which depends, in general, on both the Biot number and the thermoviscosity number) above which no ``full-film'' solutions corresponding to a continuous film of fluid covering the entire outside or inside of the cylinder exist. The effect of thermoviscosity on both the critical solution and the full-film solution with a prescribed load is described. In particular, there are no full-film solutions with a prescribed load M for any value of the Biot number when M is greater than or equal to M_{c0} divided by the square root of f for positive thermoviscosity number and when M is greater than M_{c0} for negative thermoviscosity number, where f is a monotonically decreasing function of the thermoviscosity number and M_{c0} = 4.44272 is the critical load in the constant-viscosity case. It is also found that when the prescribed load M is less than 1.50315 there is a narrow region of the Biot number - thermoviscosity number parameter plane in which backflow occurs

Interfacing peripheral nerve with macro-sieve electrodes following spinal cord injury

Macro-sieve electrodes were implanted in the sciatic nerve of five adult male Lewis rats following spinal cord injury to assess the ability of the macro-sieve electrode to interface regenerated peripheral nerve fibers post-spinal cord injury. Each spinal cord injury was performed via right lateral hemisection of the cord at the T9–10 site. Five months post-implantation, the ability of the macro-sieve electrode to interface the regenerated nerve was assessed by stimulating through the macro-sieve electrode and recording both electromyography signals and evoked muscle force from distal musculature. Electromyography measurements were recorded from the tibialis anterior and gastrocnemius muscles, while evoked muscle force measurements were recorded from the tibialis anterior, extensor digitorum longus, and gastrocnemius muscles. The macro-sieve electrode and regenerated sciatic nerve were then explanted for histological evaluation. Successful sciatic nerve regeneration across the macro-sieve electrode interface following spinal cord injury was seen in all five animals. Recorded electromyography signals and muscle force recordings obtained through macro-sieve electrode stimulation confirm the ability of the macro-sieve electrode to successfully recruit distal musculature in this injury model. Taken together, these results demonstrate the macro-sieve electrode as a viable interface for peripheral nerve stimulation in the context of spinal cord injury