Comparison of DES predictions of turbulent flow and heat transfer in a 45° ribbed duct

Numerical predictions of the turbulent flow and heat transfer of a stationary duct with square ribs 45° angled to the main flow direction are presented. The rib height to channel hydraulic diameter is 0.1, the rib pitch to rib height is 10. The calculations have been carried out for a bulk Reynolds number of 50,000. The flows generated by ribs are dominated by separating and reattaching shear layers with vortex shedding and secondary flows in the cross-section. The hybrid RANS-LES approach is adopted to simulate such flows at a reasonable computation cost. The capability of the various versions of DES method, depending the RANS model, such as DES-SA, DES-RKE, DES-SST, have been compared and validated against the experiment. The significant effect of RANS model on the accuracy of the DES prediction has been shown. The DES-SST method, which was able to reproduce the correct physics of flow and heat transfer in a ribbed duct showed better performance than others.

Submillimeter wave frequency selective surface with polarisation independent spectral responses

This paper reports the design, construction and electromagnetic performance of a new freestanding frequency selective surface (FSS) structure which generates coincident spectral responses for dual polarisation excitation at oblique angles of incidence. The FSS is required to allow transmission of 316.5 - 325.5 GHz radiation with a loss = 0.6 dB and to achieve = 30 dB rejection from 349.5 - 358.5 GHz. It should also exhibit crosspolarisation levels below -25 dB, all criteria being satisfied simultaneously for TE and TM polarisations at 45° incidence. The filter consists of two identical, 30 mm diameter, 12.5 ?m thick, optically flat, perforated metal screens separated by 450 ?m. Each of the ˜5000 unit cells contains two nested, short circuited, rectangular loop slots and a rectangular dipole slot. The nested elements provide a passband spectral response centred at 320 GHz in the TE and TM planes; the dipole slot increases the filter roll-off above resonance. The FSS was fabricated from silicon-on-insulator wafers using precision micromachining and plating processes including the use of Deep Reactive Ion Etching (DRIE) to pattern the individual slots and remove the substrate under the periodic arrays. Quasi–optical transmission measurements in the 250 – 360 GHz range yielded virtually identical copolarised spectral responses, with the performance meeting or exceeding the above specifications. Experimental results are in excellent agreement with numerical predictions.

Viscoelastic Material Models of Polypropylene for Thermoforming Applications

Polypropylene (PP), a semi-crystalline material, is typically solid phase thermoformed at temperatures associated with crystalline melting, generally in the 150° to 160°Celsius range. In this very narrow thermoforming window the mechanical properties of the material rapidly decline with increasing temperature and these large changes in properties make Polypropylene one of the more difficult materials to process by thermoforming. Measurement of the deformation behaviour of a material under processing conditions is particularly important for accurate numerical modelling of thermoforming processes. This paper presents the findings of a study into the physical behaviour of industrial thermoforming grades of Polypropylene. Practical tests were performed using custom built materials testing machines and thermoforming equipment at Queen′s University Belfast. Numerical simulations of these processes were constructed to replicate thermoforming conditions using industry standard Finite Element Analysis software, namely ABAQUS and custom built user material model subroutines. Several variant constitutive models were used to represent the behaviour of the Polypropylene materials during processing. This included a range of phenomenological, rheological and blended constitutive models. The paper discusses approaches to modelling industrial plug-assisted thermoforming operations using Finite Element Analysis techniques and the range of material models constructed and investigated. It directly compares practical results to numerical predictions. The paper culminates discussing the learning points from using Finite Element Methods to simulate the plug-assisted thermoforming of Polypropylene, which presents complex contact, thermal, friction and material modelling challenges. The paper makes recommendations as to the relative importance of these inputs in general terms with regard to correlating to experimentally gathered data. The paper also presents recommendations as to the approaches to be taken to secure simulation predictions of improved accuracy.

Structural phase transitions in low-dimensional ion crystals

A chain of singly charged particles, confined by a harmonic potential, exhibits a sudden transition to a zigzag configuration when the radial potential reaches a critical value, depending on the particle number. This structural change is a phase transition of second order, whose order parameter is the crystal displacement from the chain axis. We study analytically the transition using Landau theory and find full agreement with numerical predictions by Schiffer [Phys. Rev. Lett. 70, 818 (1993)] and Piacente [Phys. Rev. B 69, 045324 (2004)]. Our theory allows us to determine analytically the system's behavior at the transition point.

Comparison of lagrangian and eulerian simulations of slurry flows in a sudden expansion

From a review of technical literature, it was not apparent if the Lagrangian or the Eulerian dispersed phase modeling approach was more valid to simulate dilute erosive slurry flow. In this study, both modeling approaches were employed and a comparative analysis of performances and accuracy between the two models was carried out. Due to an impossibility to define, for the Eulerian model already implemented in FLUENT, a set of boundary conditions consistent with the Lagrangian impulsive equations, an Eulerian dispersed phase model was integrated in the FLUENT code using subroutines and user-defined scalar equations. Numerical predictions obtained from the two different approaches for two-phase flow in a sudden expansion were compared with the measured data. Excellent agreement was attained between the predicted and observed fluid and particle velocity in the axial direction and for the kinetic energy. Erosion profiles in a sudden expansion computed using the Lagrangian scheme yielded good qualitative agreement with measured data and predicted a maximum impact angle of 29 deg at the fluid reattachment point. The Eulerian model was adversely affected by the reattachment of the fluid phase to the wall and the simulated erosion profiles were not in agreement with the Lagrangian or measured data. Furthermore, the Eulerian model under-predicted the Lagrangian impact angle at all locations except the reattachment point. © 2010 American Society of Mechanical Engineers.

Standing wave oscillations in an electrocatalytic reaction

Experimental standing wave oscillations of the interfacial potential across an electrode have been observed in the electrocatalytic oxidation of formic acid on a Pt ring working electrode. The instantaneous potential distribution was monitored by means of equispaced potential microprobes along the electrode. The oscillatory standing waves spontaneously arose from a homogeneous stationary state prior to a Hopf bifurcation if the reference electrode was placed close to the working electrode. Reduced electrolyte concentrations resulted in aperiodic potential patterns, while the presence of a sufficiently large ohmic resistance completely suppressed spatial inhomogeneities. The experimental findings confirm numerical predictions of a reaction-migration formalism: under the chosen geometry, a long-range negative potential coupling between distant points across the ring electrode can lead to oscillatory potential domains of distinct phase. It is further shown that the occurrence of oscillatory standing waves can be rationalized as the electrochemical equivalent of Turing's second bifurcation (wave bifurcation). In the presence of an external resistance, the coupling becomes positive throughout and leads to spatial synchronization.

Predicting the crushing behaviour of composite material using high-fidelity finite element modelling

The capability to numerically model the crushing behaviour of composite structures will enable the efficient design of structures with high specific energy absorption capacity. This is particularly relevant to the aerospace and automotive industries where cabin structures need to be shown to be crashworthy. In this paper, a three-dimensional damage model is presented, which accurately represents the behaviour of composite laminates under crush loading. Both intralaminar and interlaminar failure mechanisms are taken into account. The crush damage model was implemented in ABAQUS/Explicit as a VUMAT subroutine. Numerical predictions are shown to agree well with experimental results, accurately capturing the intralaminar and interlaminar damage for a range of stacking sequences, triggers and composite materials. The use of measured material parameters required by the numerical models, without the need to ‘calibrate’ this input data, demonstrates this computational tool's predictive capabilities

Finite Element Analysis of the Enhanced Strength of Laterally Restrained RC Slabs

This paper presents the numerical simulation of the ultimate behaviour of 85 one-way and two-way spanning laterally restrained concrete slabs of variable thickness, span, reinforcement ratio, strength and boundary conditions reported in literature by different authors. The developed numerical model was described and all the assumptions were illustrated. ABAQUS, a Finite Element Analysis suite of software, was employed. Non-linear implicit static general analysis method offered by ABAQUS was used. Other analysis methods were also discussed in general in terms of application such as Explicit Dynamic Analysis and Riks method. The aim is to demonstrate the ability and efficacy of FEA to simulate the ultimate load behaviour of slabs considering different material properties and boundary conditions. The authors intended to present a numerical model that provides consistent predictions of the ultimate behaviour of laterally restrained slabs that could be used as an alternative for expensive real life testing as well as for the design and assessment of new and existing structures respectively. The enhanced strength of laterally-restrained slabs compared with conventional design methods predictions is believed to be due to compressive membrane action (CMA). CMA is an inherent phenomenon of laterally restrained concrete beams/slabs. The numerical predictions obtained from the developed model were in good correlation with the experimental results and with those obtained from the CMA method developed at the Queen’s University Belfast, UK.

Ultra Thin Resistively Loaded FSS Absorber For Polarisation Independent Operation at Large Incident Angles

In this paper we report on a resistively loaded Frequency Selective Surface (FSS) absorber design which is insensitive to the polarization of microwave signals incident at angles of 45o ± 5o. The metal backed periodic structure is composed of an array of conductive rectangular loops, each loaded with a resistor at the center of the four sides. The geometry of the absorber and the resistance value of the vertical and horizontal resistor pairs are carefully chosen so that the structure presents a real impedance of 377 Ω at the center operating frequency for both TE and TM polarized waves incident at 45o. Numerical predictions of the electromagnetic scattering from three different absorbers, designed to work at X-band, are used to investigate the effect of thickness and resistance value on the reflectivity bandwidth and angular sensitivity.

A numerical investigation of the flow structures and losses for turbulent flow in 90 degrees elbow bends

Computational fluid dynamic modelling was carried out on a series of pipe bends having R/r values of 1.3, 5, and 20, with the purpose of determining the accuracy of numerical models in predicting pressure loss data from which to inform one-dimensional loss models. Four separate turbulence models were studied: the standard k-epsilon model, realizable k-epsilon model, k-omega model, and a Reynolds stress model (RSM). The results are presented for each bend in the form of upstream and downstream pressure profiles, pressure distributions along the inner and outer walls, detailed pressure and velocity fields as well as overall loss values. In each case, measured data were presented to evaluate the predictive ability of each model. The RSM was found to perform the best, producing accurate pressure loss data for bends with R/r values of 5 and 20. For the tightest bend with an R/r value of 1.3, however, predictions were significantly worse due to the presence of flow separation, stronger pressure gradients, and high streamline curvature.

Numerical and Experimental Investigation of Aircraft Panel Deformations During Riveting Process

In collaboration with Airbus-UK, the dimensional growth of aircraft panels while being riveted with stiffeners is investigated. Small panels are used in this investigation. The stiffeners have been fastened to the panels with rivets and it has been observed that during this operation the panels expand in the longitudinal and transverse directions. It has been observed that the growth is variable and the challenge is to control the riveting process to minimize this variability. In this investigation, the assembly of the small panels and longitudinal stiffeners has been simulated using static stress and nonlinear explicit finite element models. The models have been validated against a limited set of experimental measurements; it was found that more accurate predictions of the riveting process are achieved using explicit finite element models. Yet, the static stress finite element model is more time efficient, and more practical to simulate hundreds of rivets and the stochastic nature of the process. Furthermore, through a series of numerical simulations and probabilistic analyses, the manufacturing process control parameters that influence panel growth have been identified. Alternative fastening approaches were examined and it was found that dimensional growth can be controlled by changing the design of the dies used for forming the rivets.

Nonlinear numerical modelling of lightning strike effect on composite panels with temperature dependent material properties

This paper presents a physics based modelling procedure to predict the thermal damage of composite material when struck by lightning. The procedure uses the Finite Element Method with non-linear material models to represent the extreme thermal material behaviour of the composite material (carbon/epoxy) and an embedded copper mesh protection system. Simulation predictions are compared against published experimental data, illustrating the potential accuracy and computational cost of virtual lightning strike tests and the requirement for temperature dependent material modelling. The modelling procedure is then used to examine and explain a number of practical solutions to minimize thermal material damage. © 2013 Elsevier Ltd.

Thermodynamic Properties of Dichloromethane, Bromochloromethane, and Dibromomethane under Elevated Pressure: Experimental Results and SAFT-VR Mie Predictions

The speeds of sound in dibromomethane, bromochloromethane, and dichloromethane have been measured in the temperature range from 293.15 to 313.15 K and at pressures up to 100 MPa. Densities and isobaric heat capacities at atmospheric pressure have been also determined. Experimental results were used to calculate the densities and isobaric heat capacities as the function of temperature and pressure by means of a numerical integration technique. Moreover, experimental data at atmospheric pressure were then used to determine the SAFT-VR Mie molecular parameters for these liquids. The accuracy of the model has been then evaluated using a comparison of derived experimental high-pressure data with those predicted using SAFT. It was found that the model provide the possibility to predict also the isobaric heat capacity of all selected haloalkanes within an error up to 6%.