The interaction of riblets with wall-bounded turbulence

The purpose of this thesis is to give answer to the question: why do riblets stop working for a certain size? Riblets are small surface grooves aligned in the mean direction of an overlying turbulent flow, designed specifically to reduce the friction between the flow and the surface. They were inspired by biological surfaces, like the oriented denticles in the skin of fastswimming sharks, and were the focus of a significant amount of research in the late eighties and nineties. Although it was found that the drag reduction depends on the riblet size scaled in wall units, the physical mechanisms implicated have not been completely understood up to now. It has been explained how riblets of vanishing size interact with the turbulent flow, producing a change in the drag proportional to their size, but that is not the regime of practical interest. The optimum performance is achieved for larger sizes, once that linear behavior has broken down, but before riblets begin adopting the character of regular roughness and increasing drag. This regime, which is the most relevant from a technological perspective, was precisely the less understood, so we have focused on it. Our efforts have followed three basic directions. First, we have re-assessed the available experimental data, seeking to identify common characteristics in the optimum regime across the different existing riblet geometries. This study has led to the proposal of a new length scale, the square root of the groove crosssection, to substitute the traditional peak-to-peak spacing. Scaling the riblet dimension with this length, the size of breakdown of the linear behavior becomes roughly universal. This suggests that the onset of the breakdown is related to a certain, fixed value of the cross-section of the groove. Second, we have conducted a set of direct numerical simulations of the turbulent flow over riblets, for sizes spanning the full drag reduction range. We have thus been able to reproduce the gradual transition between the different regimes. The spectral analysis of the flows has proven particularly fruitful, since it has made possible to identify spanwise rollers immediately above the riblets, which begin to appear when the riblet size is close to the optimum. This is a quite surprising feature of the flow, not because of the uniqueness of the phenomenon, which had been reported before for other types of complex and porous surfaces, but because most previous studies had focused on the detail of the flow above each riblet as a unit. Our novel approach has provided the adequate tools to capture coherent structures with an extended spanwise support, which interact with the riblets not individually, but collectively. We have also proven that those spanwise structures are responsible for the increase in drag past the viscous breakdown. Finally, we have analyzed the stability of the flow with a simplified model that connects the appearance of rollers to a Kelvin–Helmholtz-like instability, as is the case also for the flow over plant canopies and porous surfaces. In spite of the model emulating the presence of riblets only in an averaged, general fashion, it succeeds to capture the essential attributes of the breakdown, and provides a theoretical justification for the scaling with the groove cross-section.

Vortical structures on three-dimensional shock control bumps

The vortical wake structure produced by a three-dimensional shock control bump (SCB) is thought to be useful for controlling transonic buffet on airfoils. However, at present the vorticity produced is relatively weak and the production mechanism is not well understood. Using a combined experimental and computational approach, a preliminary investigation on the wake vorticity for different bump geometries has been carried out. The structure of the wake for on and off-design conditions are considered, and the effects on the downstream boundary layer demonstrated. Three main vortical structures are observed: a primary vortex pair, weak inter-bump vortices and shear flow in the lambda-shock region. The effect of pressure gradients on vortex strength is examined and it is found that spanwise pressure gradients on the front section of the bump are the most significant parameter influencing vortex strength. © 2013 by S.P. Colliss et al.

A computational fluid-structure interaction study of inflation cycles of a tension cone decelerator

Inflatable aerodynamic decelerators present potential advantages for planetary entry in missions of robotic and human exploration. The design of these structures face many engineering challenges, including complex deformable geometries, anisotropic material response, and coupled shockturbulence interactions. In this paper, we describe a comprehensive computational fluid-structure interaction study of an inflation cycle of a tension cone decelerator in supersonic flow and compare the simulations with earlier published experimental results. The aeroshell design and flow conditions closely match recent experiments conducted at Mach 2.5. The structural model is a 16-sided polygonal tension cone with seams between each segment. The computational model utilizes adaptive mesh refinement, large-eddy simulation, and shell mechanics with self-contact modeling to represent the flow and structure interaction. This study focuses on the dynamics of the structure as the inflation pressure varies gradually, and the behavior of forces experienced by the flexible and rigid (the payload capsule) structures. © 2011 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Finite Element studies on indentation size effect using a higher order strain gradient theory

In this work, a Finite Element implementation of a higher order strain gradient theory (due to Fleck and Hutchinson, 2001) has been used within the framework of large deformation elasto-viscoplasticity to study the indentation of metals with indenters of various geometries. Of particular interest is the indentation size effect (ISE) commonly observed in experiments where the hardness of a range of materials is found to be significantly higher at small depths of indentation but reduce to a lower, constant value at larger depths. That the ISE can be explained by strain gradient plasticity is well known but this work aims to qualitatively compare a gamut of experimental observations on this effect with predictions from a higher order strain gradient theory. Results indicate that many of the experimental observations are qualitatively borne out by our simulations. However, areas exist where conflicting experimental results make assessment of numerical predictions difficult. © 2012 Elsevier Ltd. All rights reserved.

Manufacturing dish shaped rings on radial-axial ring rolling mills

Ring rolling is an established method to produce seamless rings of different cross-sectional geometries. For dish shaped rings, there are applications in different areas such as offshore, aeronautics or the energy sector. At the moment, dish shaped rings are produced by machining of rings with rectangular shaped cross section, by (open die) hollow forging on a conical mandrel or by using shaped ring rolling tools. These ways of manufacturing have the disadvantage of high material waste, additional costs for special tools, long process time and limited or inflexible geometries. Therefore, the manufacturing of dish shaped rings on conventional radial-axial ring rolling mills would expand the range of products for ring producers. The aim of this study is to investigate the feasibility of an alternative to the current manufacturing processes, without requiring additional tooling and material costs. Therefore, the intended formation of dish shaped rings-previously regarded as a form error-is investigated. Based on an analysis of geometrical requirements and metal flow mechanisms, a rolling strategy is presented, causing dishing and ring climbing by a large height reduction of the ring. Using this rolling strategy dish shaped rings with dishing angles up to 18° were achieved. In addition to the experiments finite element method (FEM)-simulations of the process have been successfully conducted, in order to analyze the local strain evolution. However, when the contact between ring and main roll is lost in the process the ring starts to oscillate around the mandrel and neither dishing nor ring climbing is observed. © 2013 German Academic Society for Production Engineering (WGP).

A fracture criterion for the notch strength of high strength steels in the presence of hydrogen

High strength steels can suffer from a loss of ductility when exposed to hydrogen, and this may lead to sudden failure. The hydrogen is either accommodated in the lattice or is trapped at defects, such as dislocations, grain boundaries and carbides. The challenge is to identify the effect of hydrogen located at different sites upon the drop in tensile strength of a high strength steel. For this purpose, literature data on the failure stress of notched and un-notched steel bars are re-analysed; the bars were tested over a wide range of strain rates and hydrogen concentrations. The local stress state at failure has been determined by the finite element (FE) method, and the concentration of both lattice and trapped hydrogen is predicted using Oriani's theory along with the stress-driven diffusion equation. The experimental data are rationalised in terms of a postulated failure locus of peak maximum principal stress versus lattice hydrogen concentration. This failure locus is treated as a unique material property for the given steel and heat treatment condition. We conclude that the presence of lattice hydrogen increases the susceptibility to hydrogen embrittlement whereas trapped hydrogen has only a negligible effect. It is also found that the observed failure strength of hydrogen charged un-notched bars is less than the peak local stress within the notched geometries. Weakest link statistics are used to account for this stressed volume effect. © 2013 Elsevier Ltd.

Joint experimental and numerical approach to three-dimensional shock control bump research

Previous studies of transonic shock control bumps have often been either numerical or experimental. Comparisons between the two have been hampered by the limitations of either approach. The present work aims to bridge the gap between computational fluid dynamics and experiment by planning a joint approach from the outset. This enables high-quality validation data to be produced and ensures that the conclusions of either aspect of the study are directly relevant to the application. Experiments conducted with bumps mounted on the floor of a blowdown tunnel were modified to include an additional postshock adverse pressure gradient through the use of a diffuser as well as introducing boundary-layer suction ahead of the test section to enable the in-flow boundary layer to be manipulated. This has the advantage of being an inexpensive and highly repeatable method. Computations were performed on a standard airfoil model, with the flight conditions as free parameters. The experimental and computational setups were then tuned to produce baseline conditions that agree well, enabling confidence that the experimental conclusions are relevant. The methods are then applied to two different shock control bumps: a smoothly contoured bump, representative of previous studies, and a novel extended geometry featuring a continuously widening tail, which spans the wind-tunnel width at the rear of the bump. Comparison between the computational and experimental results for the contour bump showed good agreement both with respect to the flow structures and quantitative analysis of the boundary-layer parameters. It was seen that combining the experimental and numerical data could provide valuable insight into the flow physics, which would not generally be possible for a one-sided approach. The experiments and computational fluid dynamics were also seen to agree well for the extended bump geometry, providing evidence that, even though thebumpinteracts directly with the wind-tunnel walls, it was still possible to observe the key flow physics. The joint approach is thus suitable even for wider bump geometries. Copyright © 2013 by S. P. Colliss, H. Babinsky, K. Nubler, and T. Lutz. Published by the American Institute of Aeronautics and Astronautics, Inc.

High performance bilayer-graphene Terahertz detectors

© 2014 AIP Publishing LLC. We report bilayer-graphene field effect transistors operating as Terahertz (THz) broadband photodetectors based on plasma-waves excitation. By employing wide-gate geometries or buried gate configurations, we achieve a responsivity ∼1.2 V/W (1.3 mA/W) and a noise equivalent power ∼2 × 10-9 W/√Hz in the 0.29-0.38 THz range, in photovoltage and photocurrent mode. The potential of this technology for scalability to higher frequencies and the development of flexible devices makes our approach competitive for a future generation of THz detection systems.

The effect of carbon nanotube orientation on erosive wear resistance of CNT-epoxy based composites

Aligned carbon nanotube (CNT) polymer composites are envisioned as the next-generation composite materials for a wide range of applications. In this work, we investigate the erosive wear behavior of epoxy matrix composites reinforced with both randomly dispersed and aligned carbon nanotube (CNT) arrays. The aligned CNT composites are prepared in two different configurations, where the sidewalls and ends of nanotubes are exposed to the composite surface. Results have shown that the composite with vertically aligned CNT-arrays exhibits superior erosive wear resistance compared to any of the other types of composites, and the erosion rate reaches a similar performance level to that of carbon steel at 20° impingement angle. The erosive wear mechanism of this type of composite, at various impingement angles, is studied by Scanning Electron Microscopy (SEM). We report that the erosive wear performance shows strong dependence on the alignment geometries of CNTs within the epoxy matrix under identical nanotube loading fractions. Correlations between the eroded surface roughness and the erosion rates of the CNT composites are studied by surface profilometry. This work demonstrates methods to fabricate CNT based polymer composites with high loading fractions of the filler, alignment control of nanotubes and optimized erosive wear properties. © 2014 Elsevier Ltd. All rights reserved.

Efficient strategies for process design and optimization in ring rolling

Ring rolling is an incremental bulk forming process for the near-net-shape production of seamless rings. This paper shows how nowadays the process design and optimization can be efficiently supported by simulation methods. For reliable predictions of the material flow and the microstructure evolution it's necessary to include a real ring rolling mill's control algorithm into the model. Furthermore an approach for the online measurement of the profile evolution during the process is presented by means of axial profiling in ring rolling. Hence the definition of new ring rolling strategies is possible even for advanced geometries.

Characterizing the class of local metal sheet thickening processes

The need to create high-value products for specialist applications, and the search for efficient forming routes that obviate the need for some machining steps, is driving Interest In a novel class of forming processes aiming to create locally thickened features within sheet work- pieces. A number of novel forming processes have been proposed to meet this need, but it is as yet unclear which processes will be most effective in creating local thickening of various geometries, and many process configurations have yet to be tried. This paper aims to provide some basic principles for designing and characterising process behaviour. A simplified generic description of sheet thickening processes is provided, with two tools of variable operating on a sheet workpiece in plane strain, with different tool separations and motions parameterised. A comprehensive numerical study of the behaviour of this class of processes is conducted in Abaqus to predict the main characteristics of the material flow in each configuration. The results are used to classify the different basic behaviours that can be achieved by the sheet-bulk thickening processes and to give guidance on future process development, capability and applicability. © 2011 Wiley-VCH Verlag GmbH & Co. KGaA. Weinheim.

The forming of variable cross-section I-beams by hot rolling

A process is presented for the forming of variable cross-section I-beams by hot rolling. Optimized I-beams with variable cross-section offer a significant weight advantage over prismatic beams. By tailoring the cross-section to the bending moment experienced within the beam, around 30% of the material can be saved compared to a standard section. Production of such beams by hot rolling would be advantageous, as It combines high volume capacity with high material yields. Through controlled variation of the roll gap during multiple passes, beams with a variable cross-section have been created using shaped rolls similar to those used for conventional I-beam rolling. The process was tested experimentally on a small scale rolling mill, using plasticine as the modelling material. These results were then compared to finite element simulations of individual stages of the process conducted using Abaqus/Standard. Results here show that the process can successfully form a beam with a variable depth web. The main failure modes of the process, and the limitations on the achievable variations In geometry are also presented. Finally, the question of whether or not optimal beam geometries can be created by this process Is discussed. © 2011 Wiley-VCH Verlag GmbH & Co. KGaA. Weinheim.

Differential equation based length scales to improve DES and RANS simulations

Surprisingly expensive to compute wall distances are still used in a range of key turbulence and peripheral physics models. Potentially economical, accuracy improving differential equation based distance algorithms are considered. These involve elliptic Poisson and hyperbolic natured Eikonal equation approaches. Numerical issues relating to non-orthogonal curvilinear grid solution of the latter are addressed. Eikonal extension to a Hamilton-Jacobi (HJ) equation is discussed. Use of this extension to improve turbulence model accuracy and, along with the Eikonal, enhance Detached Eddy Simulation (DES) techniques is considered. Application of the distance approaches is studied for various geometries. These include a plane channel flow with a wire at the centre, a wing-flap system, a jet with co-flow and a supersonic double-delta configuration. Although less accurate than the Eikonal, Poisson method based flow solutions are extremely close to those using a search procedure. For a moving grid case the Poisson method is found especially efficient. Results show the Eikonal equation can be solved on highly stretched, non-orthogonal, curvilinear grids. A key accuracy aspect is that metrics must be upwinded in the propagating front direction. The HJ equation is found to have qualitative turbulence model improving properties. © 2003 by P. G. Tucker.

Fixed-rank matrix factorizations and Riemannian low-rank optimization

Motivated by the problem of learning a linear regression model whose parameter is a large fixed-rank non-symmetric matrix, we consider the optimization of a smooth cost function defined on the set of fixed-rank matrices. We adopt the geometric framework of optimization on Riemannian quotient manifolds. We study the underlying geometries of several well-known fixed-rank matrix factorizations and then exploit the Riemannian quotient geometry of the search space in the design of a class of gradient descent and trust-region algorithms. The proposed algorithms generalize our previous results on fixed-rank symmetric positive semidefinite matrices, apply to a broad range of applications, scale to high-dimensional problems, and confer a geometric basis to recent contributions on the learning of fixed-rank non-symmetric matrices. We make connections with existing algorithms in the context of low-rank matrix completion and discuss the usefulness of the proposed framework. Numerical experiments suggest that the proposed algorithms compete with state-of-the-art algorithms and that manifold optimization offers an effective and versatile framework for the design of machine learning algorithms that learn a fixed-rank matrix. © 2013 Springer-Verlag Berlin Heidelberg.