Vibration of a high-altitude tethered balloon with application to climate engineering

A balloon tethered at an altitude of 20 km could deliver a particulate cloud leading to global cooling. Tethering a balloon at this altitude poses significant problems with respect to vibration and stability, especially in regions of high wind. No-one has ever proposed, yet alone launched, a balloon at an altitude of 20 km tethered to the ground. Owing to wind, the tether needs to be 23 km in length and is to be fixed to a ship at sea or on land in equatorial regions. Whilst the balloon at 20 km is subject to relatively modest wind conditions, at jet stream altitudes (10km) the tether will experience much higher wind loadings, not only because of the high wind speeds of up to 300 km / hr but also because of the high air density. A tether of circular cross section in these high winds will be subject to horizontal and downward drag forces that would bring the aerostat down. For this reason it is advantageous to consider a self-aligning tether of an aerodynamic cross section whereby it is possible to reduce the drag substantially. One disadvantage of a non-circular tether is the possibility of flutter and galloping instabilities. It is reasonably straightforward to model these phenomena for short lengths of aerofoil, but the situation becomes more complex for a 20 km tensioned tether with large deflection and curvature, variable wind speed, variable air density and variable tension. Analysis using models of infinite length are used to establish the stability at a local scale where the tension, aerodynamic and geometric properties are considered constant. Dispersion curve analysis is useful here. But for dynamics on a long-wavelength scale (several km) then a full non-linear analysis is required. This non-linear model can be used to establish the local values of tension appropriate for the dispersion analysis. This keynote presentation will give some insight into these issues.

HTGR fuel element depletion benchmark: Stage three results

Recently, a new numerical benchmark exercise for High Temperature Gas Cooled Reactor (HTGR) fuel depletion was defined. The purpose of this benchmark is to provide a comparison basis for different codes and methods applied to the burnup analysis of HTGRs. The benchmark specifications include three different models: (1) an infinite lattice of tristructural isotropic (TRISO) fuel particles, (2) an infinite lattice of fuel pebbles, and (3) a prismatic fuel including fuel and coolant channels. In this paper, we present the results of the third stage of the benchmark obtained with MCNP based depletion code BGCore and deterministic lattice code HELIOS 1.9. The depletion calculations were performed for three-dimensional model of prismatic fuel with explicitly described TRISO particles as well as for two-dimensional model, in which double heterogeneity of the TRISO particles was eliminated using reactivity equivalent physical transformation (RPT). Generally, good agreement in the results of the calculations obtained using different methods and codes was observed.

Comparison among MCNP-based depletion codes applied to burnup calculations of pebble-bed HTR lattices

The double-heterogeneity characterising pebble-bed high temperature reactors (HTRs) makes Monte Carlo based calculation tools the most suitable for detailed core analyses. These codes can be successfully used to predict the isotopic evolution during irradiation of the fuel of this kind of cores. At the moment, there are many computational systems based on MCNP that are available for performing depletion calculation. All these systems use MCNP to supply problem dependent fluxes and/or microscopic cross sections to the depletion module. This latter then calculates the isotopic evolution of the fuel resolving Bateman's equations. In this paper, a comparative analysis of three different MCNP-based depletion codes is performed: Montburns2.0, MCNPX2.6.0 and BGCore. Monteburns code can be considered as the reference code for HTR calculations, since it has been already verified during HTR-N and HTR-N1 EU project. All calculations have been performed on a reference model representing an infinite lattice of thorium-plutonium fuelled pebbles. The evolution of k-inf as a function of burnup has been compared, as well as the inventory of the important actinides. The k-inf comparison among the codes shows a good agreement during the entire burnup history with the maximum difference lower than 1%. The actinide inventory prediction agrees well. However significant discrepancy in Am and Cm concentrations calculated by MCNPX as compared to those of Monteburns and BGCore has been observed. This is mainly due to different Am-241 (n,γ) branching ratio utilized by the codes. The important advantage of BGCore is its significantly lower execution time required to perform considered depletion calculations. While providing reasonably accurate results BGCore runs depletion problem about two times faster than Monteburns and two to five times faster than MCNPX. © 2009 Elsevier B.V. All rights reserved.

Effects of cooling time on a closed LWR fuel cycle

In this study, the effects of cooling time prior to reprocessing spent LWR fuel has on the reactor physics characteristics of a PWR fully loaded with homogeneously mixed U-Pu or U-TRU oxide (MOX) fuel is examined. A reactor physics analysis was completed using the CASM04e code. A void reactivity feedback coefficient analysis was also completed for an infinite lattice of fresh fuel assemblies. Some useful conclusions can be made regarding the effect that cooling time prior to reprocessing spent LWR fuel has on a closed homogeneous MOX fuel cycle. The computational analysis shows that it is more neutronically efficient to reprocess cooled spent fuel into homogeneous MOX fuel rods earlier rather than later as the fissile fuel content decreases with time. Also, the number of spent fuel rods needed to fabricate one MOX fuel rod increases as cooling time increases. In the case of TRU MOX fuel, with time, there is an economic tradeoff between fuel handling difficulty and higher throughput of fuel to be reprocessed. The void coefficient analysis shows that the void coefficient becomes progressively more restrictive on fuel Pu content with increasing spent fuel cooling time before reprocessing.

Aerodynamic noise from a poroelastic trailing edge with implications for the silent flight of owls

The interaction of a turbulent eddy with a semi-infinite, poroelastic edge is examined with respect to the effects of both elasticity and porosity on the efficiency of scattered aerodynamic noise. The scattering problem is solved using the Wiener-Hopf technique for constant plate properties to identify their scaling dependence on the resulting aerodynamic noise, including the dependence on flight velocity, where special attention is paid to the limiting cases of rigid, porous and elastic, impermeable plate conditions. Results from these analyses attempt to address how trailing edge noise may be mitigated by porosity and seek to deepen the understanding of how owls hunt in acoustic stealth. © 2012 by Justin W. Jaworski and Nigel Peake. Published by the American Institute of Aeronautics and Astronautics, Inc.

An analytical study of the low frequency structural responses of a submerged hull

This paper studies the low frequency vibrational behaviour of a submerged hull. The submerged hull is modelled as a finite fluid-loaded cylindrical shell closed at each end by circular plates. The external pressure acting on the hull due to the fluid loading is analytically calculated using an infinite model. Three excitation cases of the hull are considered. In the first model, an axial point force is applied at the centre of one end plate, giving rise to an axisymmetric case in which only the zeroth circumferential shell modes are excited. In the second model, an axial point force is applied at the edge of the end plate. In the third model, a radial point force is applied also at the edge of the end plate. In the second and third load cases, all cylindrical shell circumferential modes are excited. The effects of fluid loading and different excitation locations are studied. A more complex hull model including stiffeners and bulkheads is then examined. A smeared approach is used to analytically model the ring stiffeners. All load cases are again considered and the effects of the various influencing factors on the low frequency responses are described.

Lock-in and quasiperiodicity in a forced hydrodynamically self-excited jet

The ability of hydrodynamically self-excited jets to lock into strong external forcing is well known. Their dynamics before lock-in and the specific bifurcations through which they lock in, however, are less well known. In this experimental study, we acoustically force a low-density jet around its natural global frequency. We examine its response leading up to lock-in and compare this to that of a forced van der Pol oscillator. We find that, when forced at increasing amplitudes, the jet undergoes a sequence of two nonlinear transitions: (i) from periodicity to T{double-struck}2 quasiperiodicity via a torus-birth bifurcation; and then (ii) from T{double-struck}2 quasiperiodicity to 1:1 lock-in via either a saddle-node bifurcation with frequency pulling, if the forcing and natural frequencies are close together, or a torus-death bifurcation without frequency pulling, but with a gradual suppression of the natural mode, if the two frequencies are far apart. We also find that the jet locks in most readily when forced close to its natural frequency, but that the details contain two asymmetries: the jet (i) locks in more readily and (ii) oscillates more strongly when it is forced below its natural frequency than when it is forced above it. Except for the second asymmetry, all of these transitions, bifurcations and dynamics are accurately reproduced by the forced van der Pol oscillator. This shows that this complex (infinite-dimensional) forced self-excited jet can be modelled reasonably well as a simple (three-dimensional) forced self-excited oscillator. This result adds to the growing evidence that open self-excited flows behave essentially like low-dimensional nonlinear dynamical systems. It also strengthens the universality of such flows, raising the possibility that more of them, including some industrially relevant flames, can be similarly modelled. © 2013 Cambridge University Press.

Parametric guidance for turbulent noise reduction from poroelastic trailing edges and owls

The interaction of a turbulent eddy with a semi-infinite, poroelastic edge is examined with respect to the effects of both elasticity and porosity on the efficiency of aerodynamic noise generation. The edge is modelled as a thin plate poroelastic plate, which is known to admit fifth-, sixth-, and seventh-power noise dependences on a characteristic velocity U of the turbulent eddy. The associated acoustic scattering problem is solved using the Wiener-Hopf technique for the case of constant plate properties. For the special cases of porous-rigid and impermeable-elastic plate conditions, asymptotic analysis of the Wiener- Hopf kernel function furnishes the parameter groups and their ranges where U5, U6, and U7 behaviours are expected to occur. Results from this analysis attempt to help guide the search for passive edge treatments to reduce trailing-edge noise that are inspired by the wing features of silently flying owls. Furthermore, the appropriateness of the present model to the owl noise problem is discussed with respect to the acoustic frequencies of interest, wing chord-lengths, and foraging behaviour across a representative set of owl species.

Models for acoustic propagation through turbofan exhaust flows

Turbomachinery noise radiating into the rearward arc is an important problem. This noise is scattered by the trailing edges of the nacelle and the jet exhaust, and interacts with the shear layers between the external flow, bypass stream and jet, en route to the far field. In the past a range of relevant model problems involving semi-infinite cylinders have been solved. However, one limitation of previous solutions is that they do not allow for the jet nozzle to protrude a finite distance beyond the end of the nacelle (or in certain configurations being buried a finite distance upstream). In this paper we use the matrix Wiener-Hopf technique, which will allow precisely the finite nacelle-jet nozzle separation to be included. The crucial step in our work is to factorise a certain matrix as a product of terms analytic and invertible in the upper/lower halves of the complex plane. The way we do this matrix factorisation is quite different in the buried and protruding nozzle cases. In the buried case our solution method is the so-called pole-removal technique. In the technically more demanding protruding case, however, we must first use Pade approximants to generate a uniformly-valid, meromorphic representation of a certain function, before the same pole-removal method can be applied. Sample results are presented, investigating in particular the effects of exit plane stagger. © 2007 by B Veitch and N Peake.

Predicting the response of structures to transient shock loading.

This work concerns the prediction of the response of an uncertain structure to a load of short duration. Assuming an ensemble of structures with small random variations about a nominal form, a mean impulse response can be found using only the modal density of the structure. The mean impulse response turns out to be the same as the response of an infinite structure: the response is calculated by taking into account the direct field only, without reflections. Considering the short duration of an impulsive loading, the approach is reasonable before the effect of the reverberant field becomes important. The convolution between the mean impulse response and the shock loading is solved in discrete time to calculate the response at the driving point and at remote points. Experimental and numerical examples are presented to validate the theory presented for simple structures such as beams, plates, and cylinders.

Climb-enabled discrete dislocation plasticity

A small strain two-dimensional discrete dislocation plasticity framework coupled to vacancy diffusion is developed wherein the motion of edge dislocations is by a combination of glide and climb. The dislocations are modelled as line defects in a linear elastic medium and the mechanical boundary value problem is solved by the superposition of the infinite medium elastic fields of the dislocations and a complimentary non-singular solution that enforces the boundary conditions. Similarly, the climbing dislocations are modelled as line sources/sinks of vacancies and the vacancy diffusion boundary value problem is also solved by a superposition of the fields of the line sources/sinks in an infinite medium and a complementary non-singular solution that enforces the boundary conditions. The vacancy concentration field along with the stress field provides the climb rate of the dislocations. Other short-range interactions of the dislocations are incorporated via a set of constitutive rules. We first employ this formulation to investigate the climb of a single edge dislocation in an infinite medium and illustrate the existence of diffusion-limited and sink-limited climb regimes. Next, results are presented for the pure bending and uniaxial tension of single crystals oriented for single slip. These calculations show that plasticity size effects are reduced when dislocation climb is permitted. Finally, we contrast predictions of this coupled framework with an ad hoc model in which dislocation climb is modelled by a drag-type relation based on a quasi steady-state solution. © 2013 Elsevier Ltd. All rights reserved.

Grain boundary interface mechanics in strain gradient crystal plasticity

Interactions between dislocations and grain boundaries play an important role in the plastic deformation of polycrystalline metals. Capturing accurately the behaviour of these internal interfaces is particularly important for applications where the relative grain boundary fraction is significant, such as ultra fine-grained metals, thin films and microdevices. Incorporating these micro-scale interactions (which are sensitive to a number of dislocation, interface and crystallographic parameters) within a macro-scale crystal plasticity model poses a challenge. The innovative features in the present paper include (i) the formulation of a thermodynamically consistent grain boundary interface model within a microstructurally motivated strain gradient crystal plasticity framework, (ii) the presence of intra-grain slip system coupling through a microstructurally derived internal stress, (iii) the incorporation of inter-grain slip system coupling via an interface energy accounting for both the magnitude and direction of contributions to the residual defect from all slip systems in the two neighbouring grains, and (iv) the numerical implementation of the grain boundary model to directly investigate the influence of the interface constitutive parameters on plastic deformation. The model problem of a bicrystal deforming in plane strain is analysed. The influence of dissipative and energetic interface hardening, grain misorientation, asymmetry in the grain orientations and the grain size are systematically investigated. In each case, the crystal response is compared with reference calculations with grain boundaries that are either 'microhard' (impenetrable to dislocations) or 'microfree' (an infinite dislocation sink). © 2013 Elsevier Ltd. All rights reserved.

Efficient statistical electromagnetic analysis in reverberant environments

This paper presents a statistical approach to the electromagnetic analysis of a system that lies within a reverberant cavity that has random or uncertain properties. The need to solve Maxwell's equations within the cavity is avoided by employing a relation known as the diffuse field reciprocity principle, which leads directly to the ensemble mean squared response of the system; all that is required is the impedance matrix of the system associated with radiation into infinite space. The general theoretical approach is presented, and the analysis is then applied to a five-cable bundle in a reverberation room © 2013 EMC Europe Foundation.

Low-energy exciton states in a ZnO cylindrical nanodisk

We consider the electron-hole pair confined in a simplified infinite potential. The low-lying excition states in a ZnO cylindrical nanodisk are calculated based on effective-mass theory. To further understand the optical properties, we calculate the linear optical susceptibilities chi(w) and the radiative recombination lifetime tau of excitons in a ZnO nanodisk. The exciton radiative lifetime in a cylindrical nanodisk is of the order of tens of picoseconds, which is small compared with the lifetime of bulk ZnO material. (C) 2008 American Institute of Physics. [DOI: 10.1063/1.3006134]

Resonant tunneling through double-bended graphene nanoribbons

We investigate theoretically resonant tunneling through double-bended graphene nanoribbon (GNR) structures, i.e., armchair-edged GNRs (AGNRs) in between two semi-infinite zigzag GNR leads. Our numerical results demonstrate that the resonant tunneling can be tuned dramatically by the Fermi energy and the length and/or widths of the AGNR for both the metallic and semiconductorlike AGNRs. The structure can also be use to control the valley polarization of the tunneling currents and could be useful for potential application in valleytronics devices. (C) 2008 American Institute of Physics.