Prediction of Transonic LCOs using an Aeroelastic Harmonic Balance Method

This work proposes a novel approach to compute transonic Lim
it Cycle Oscillations using high fidelity analysis. CFD based Harmonic Balance methods have proven to be efficient tools to predict periodic phenomena. This paper’s contribution is to present a new methodology to determine the unknown frequency of oscillations, enabling HB methods to accurately capture Limit Cycle Oscillations (LCOs); this is achieved by defining a frequency updating procedure based on a coupled CFD/CSD Harmonic Balance formulation to find the LCO condition. A pitch/plunge aerofoil and delta wing aerodynamic and respective linear structural models are used to validate the new method against conventional time-domain simulations. Results show consistent agreement between the proposed and time-marching methods for both LCO amplitude and frequency, while producing at least one order of magnitude reduction in computational time.

Simulació mitjançant tècniques CFD d'una culata per al motor del vehicle de baix consum "Àliga-EPS-UdG"

El projecte pretén estudiar i quantificar les restriccions creades al fluid en circular pels conductes d’admissió i escapament de la culata del motor del vehicle Àliga. L’estudi consta de quatre etapes: estudi de les restriccions actuals dels sistemes d’admissió i escapament; anàlisi dels resultats de la culata de sèrie i proposta de millores aplicables al model real; càlcul de les restriccions creades pels models millorats, i finalment, estudi comparatiu dels resultats obtinguts, interpretant els resultats dels principals paràmetres a analitzar

Computational simulation of gas flow and heat transfer near an immersed object in fluidized beds

To improve the understanding of the heat transfer mechanism and to find a reliable and simple heat-transfer model, the gas flow and heat transfer between fluidized beds and the surfaces of an immersed object is numerically simulated based on a double particle-layer and porous medium model. The velocity field and temperature distribution of the gas and particles are analysed during the heat transfer process. The simulation shows that the change of gas velocity with the distance from immersed surface is consistent with the variation of bed voidage, and is used to validate approximately dimensional analysing result that the gas velocity between immersed surface and particles is 4.6Umf/εmf. The effects of particle size and particle residence time on the thermal penetration depth and the heat-transfer coefficients are also discussed.

Análise numérica de um coletor solar de tubo evacuado

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

CFD simulations of hydrocyclones with an air core - Comparison between large eddy simulations and a second moment closure

CFD simulations of the 75 mm, hydrocyclone of Hsieh (1988) have been conducted using Fluent TM. The simulations used 3-dimensional body fitted grids. The simulations were two phase simulations where the air core was resolved using the mixture (Manninen et al., 1996) and VOF (Hirt and Nichols, 1981) models. Velocity predictions from large eddy simulations (LES), using the Smagorinsky-Lilly sub grid scale model (Smagorinsky, 1963; Lilly, 1966) and RANS simulations using the differential Reynolds stress turbulence model (Launder et al., 1975) were compared with Hsieh's experimental velocity data. The LES simulations gave very good agreement with Hsieh's data but required very fine grids to predict the velocities correctly in the bottom of the apex. The DRSM/RANS simulations under predicted tangential velocities, and there was little difference between the velocity predictions using the linear (Launder, 1989) and quadratic (Speziale et al., 1991) pressure strain models. Velocity predictions using the DRSM turbulence model and the linear pressure strain model could be improved by adjusting the pressure strain model constants.

Global efficiency of innovative rotational mold directly heated by thermal fluid

The article is focused on analysis of global efficiency of new mold for rotational molding of plastic parts, being directly heated by thermal fluid. The overall efficiency is based on several items such as reduction of cycle time, better uniformity of heating-cooling and low energy consumption. The new tool takes advantage of additive fabrication and electroforming for making the optimal manifold and cavity shell of the mold. Experimental test of a prototype mold was carried out on an experimental rotational molding machine, developed for this purpose, measuring wall temperature, and internal air temperature, with and without plastic material inside. Results were compared with conventional mold heated into an oven and to theoretical simulations done by Computational Fluid Dynamic software (CFD). The analysis represents considerable improvement of cycle time related to conventional methods (heated by oven) and better thermal uniformity to conventional procedures by direct heating of oil with external channels. In addition to thermal analysis an energetic efficiency study was done. POLYM. ENG. SCI., 52:1998-2005, 2012. © 2012 Society of Plastics Engineers Copyright © 2012 Society of Plastics Engineers.

Coupled CFD Shape Optimization for aerodynamic profiles

The present document deals with the optimization of shape of aerodynamic profiles -- The objective is to reduce the drag coefficient on a given profile without penalising the lift coefficient -- A set of control points defining the geometry are passed and parameterized as a B-Spline curve -- These points are modified automatically by means of CFD analysis -- A given shape is defined by an user and a valid volumetric CFD domain is constructed from this planar data and a set of user-defined parameters -- The construction process involves the usage of 2D and 3D meshing algorithms that were coupled into own- code -- The volume of air surrounding the airfoil and mesh quality are also parametrically defined -- Some standard NACA profiles were used by obtaining first its control points in order to test the algorithm -- Navier-Stokes equations were solved for turbulent, steady-state ow of compressible uids using the k-epsilon model and SIMPLE algorithm -- In order to obtain data for the optimization process an utility to extract drag and lift data from the CFD simulation was added -- After a simulation is run drag and lift data are passed to the optimization process -- A gradient-based method using the steepest descent was implemented in order to define the magnitude and direction of the displacement of each control point -- The control points and other parameters defined as the design variables are iteratively modified in order to achieve an optimum -- Preliminary results on conceptual examples show a decrease in drag and a change in geometry that obeys to aerodynamic behavior principles

Shear thinning and shear thickening non-Newtonian confined fluid flow over rotating cylinder

In present work, numerical solution is performed to study the confined flow of power-law non Newtonian fluids over a rotating cylinder. The main purpose is to evaluate drag and thermal coefficients as functions of the related governing dimensionless parameters, namely, power-law index (0.5 ≤ n ≤ 1.4), dimensionless rotational velocity (0 ≤ α ≤ 6) and the Reynolds number (100 ≤ Re ≤ 500). Over the range of Reynolds number, the flow is known to be steady. Results denoted that the increment of power law index and rotational velocity increases the drag coefficient due to momentum diffusivity improvement which is responsible for low rate of heat transfer, because the thicker the boundary layer, the lower the heat transfer is implemented.

Numerical study of lid-driven square cavity with heat generation using LBM

In this study, the mixed convection heat transfer and fluid flow behaviors in a lid-driven square cavity filled with high Prandtl number fluid (Pr = 5400, ν = 1.2×10-4 m2/s) at low Reynolds number is studied using thermal Lattice Boltzmann method (TLBM) where ν is the viscosity of the fluid. The LBM has built up on the D2Q9 model and the single relaxation time method called the Lattice-BGK (Bhatnagar-Gross-Krook) model. The effects of the variations of non dimensional mixed convection parameter called Richardson number(Ri) with and without heat generating source on the thermal and flow behavior of the fluid inside the cavity are investigated. The results are presented as velocity and temperature profiles as well as stream function and temperature contours for Ri ranging from 0.1 to 5.0 with other controlling parameters that present in this study. It is found that LBM has good potential to simulate mixed convection heat transfer and fluid flow problem. Finally the simulation results have been compared with the previous numerical and experimental results and it is found to be in good agreement.

Turbulence-induced melting of a nonequilibrium vortex crystal in a forced thin fluid film

To gain a better understanding of recent experiments on the turbulence-induced melting of a periodic array of vortices in a thin fluid film, we perform a direct numerical simulation of the two-dimensional Navier-Stokes equations forced such that, at low Reynolds numbers, the steady state of the film is a square lattice of vortices. We find that as we increase the Reynolds number, this lattice undergoes a series of nonequilibrium phase transitions, first to a crystal with a different reciprocal lattice and then to a sequence of crystals that oscillate in time. Initially, the temporal oscillations are periodic; this periodic behaviour becoming more and more complicated with increasing Reynolds number until the film enters a spatially disordered nonequilibrium statistical steady state that is turbulent. We study this sequence of transitions using fluid-dynamics measures, such as the Okubo-Weiss parameter that distinguishes between vortical and extensional regions in the flow, ideas from nonlinear dynamics, e.g. Poincare maps, and theoretical methods that have been developed to study the melting of an equilibrium crystal or the freezing of a liquid and that lead to a natural set of order parameters for the crystalline phases and spatial autocorrelation functions that characterize short- and long-range order in the turbulent and crystalline phases, respectively.

Study of Fluid Flow and Mixing Behaviour of a Vacuum Degasser

In recent times the demand of ultra-low carbon steel (ULCS) with improved mechanical properties such as good ductility and good workability has been increased as it is used to produce cold-rolled steel sheets for automobiles. For producing ULCS efficiently, it is necessary to improve the productivity of the vacuum degassers such as RH, DH and tank degasser. Recently, it has been claimed that using a new process, called REDA (revolutionary degassing activator), one can achieve the carbon content below 10 ppm in less time. As such, REDA process has not been studied thoroughly in terms of fluid flow and mass transfer which is a necessary precursor to understand and design this process. Therefore, momentum and mass transfer of the process has been studied by solving momentum and species balance equations along with k-epsilon turbulent model in two-dimension (2D) for REDA process. Similarly, computational fluid dynamic studies have been made in 2D for tank and RH degassers to compare them with REDA process. Computational results have been validated with published experimental and theoretical data. It is found that REDA process is the most efficient among all these processes in terms of mixing efficiency. Fluid flow phenomena have been studied in details for REDA process by varying gas flow rate, depth of immersed snorkel in the steel, diameter of the snorkel and change in vacuum pressure. It is found that design of snorkel affects the melt circulation in the bath significantly.

A study on the usefulness of Support Vector Machines for the realtime computational simulation of soft biological organs

Realistic and realtime computational simulation of soft biological organs (e.g., liver, kidney) is necessary when one tries to build a quality surgical simulator that can simulate surgical procedures involving these organs. Since the realistic simulation of these soft biological organs should account for both nonlinear material behavior and large deformation, achieving realistic simulations in realtime using continuum mechanics based numerical techniques necessitates the use of a supercomputer or a high end computer cluster which are costly. Hence there is a need to employ soft computing techniques like Support Vector Machines (SVMs) which can do function approximation, and hence could achieve physically realistic simulations in realtime by making use of just a desktop computer. Present work tries to simulate a pig liver in realtime. Liver is assumed to be homogeneous, isotropic, and hyperelastic. Hyperelastic material constants are taken from the literature. An SVM is employed to achieve realistic simulations in realtime, using just a desktop computer. The code for the SVM is obtained from [1]. The SVM is trained using the dataset generated by performing hyperelastic analyses on the liver geometry, using the commercial finite element software package ANSYS. The methodology followed in the present work closely follows the one followed in [2] except that [2] uses Artificial Neural Networks (ANNs) while the present work uses SVMs to achieve realistic simulations in realtime. Results indicate the speed and accuracy that is obtained by employing the SVM for the targeted realistic and realtime simulation of the liver.

Micro-blast waves using detonation transmission tubing

Micro-blast waves emerging from the open end of a detonation transmission tube were experimentally visualized in this study. A commercially available detonation transmission tube was used (Nonel tube, M/s Dyno Nobel, Sweden), which is a small diameter tube coated with a thin layer of explosive mixture (HMX traces of Al) on its inner side. The typical explosive loading for this tube is of the order of 18 mg/m of tube length. The blast wave was visualized using a high speed digital camera (frame rate 1 MHz) to acquire time-resolved schlieren images of the resulting flow field. The visualization studies were complemented by computational fluid dynamic simulations. An analysis of the schlieren images showed that although the blast wave appears to be spherical, it propagates faster along the tube axis than along a direction perpendicular to the tube axis. Additionally, CFD analysis revealed the presence of a barrel shock and Mach disc, showing structures that are typical of an underexpanded jet. A theory in use for centered large-scale explosions of intermediate strength gave good agreement with the blast trajectory along the tube axis. The energy of these micro-blast waves was found to be J and the average TNT equivalent was found to be . The repeatability in generating these micro-blast waves using the Nonel tube was very good and this opens up the possibility of using this device for studying some of the phenomena associated with muzzle blasts in the near future.

NUMERICAL ANALYSIS OF THE FLOWFIELD IN A SUPERSONIC COIL WITH AN INTERLEAVED JET CONFIGURATION AND ITS EFFECT ON THE GAIN DISTRIBUTION

In a supersonic chemical oxygen-iodine laser (COIL) operating without primary buffer gas, the features of flowfield have significant effects on the Laser efficiency and beam quality. In this paper three-dimensional, multi-species, chemically reactive CFD technology was used to study the flowfield in mixing nozzle implemented with a supersonic interleaving jet configuration. The features of the flowfield as well as its effect on the spatial distribution of small signal gain were analyzed.

Vortex dynamics in the studies of looping in tropical cyclone tracks

This paper proposes criteria for predicting the tendency of looping in tropical cyclone tracks using the approach of vortex dynamics. We model the asymmetric structure of a cyclone by a system of vortex patches. The evolution of such system of vortices is simulated by the method of contour dynamics. A new set of exact analytic formulas for contour dynamics calculations is derived, which is shown to be more computationally effective. Based on point-vortex models, we derive analytic formulas for the criteria of looping in a cyclone track. From numerical experiments, the simulated trajectories obtained from the point-vortex system and vortex patch system agree quite well. Hence, the looping criteria obtained from the point-vortex system can be applied by forecasters to stay alert for tendency of looping in a cyclone track. To demonstrate the applicability of the proposed criteria, the trajectory of Typhoon Yancy (9012), whose field data are available from ''TCM-90'', is simulated. The case study shows that the asymmetric structure similar to the pattern of a beta gyre is responsible for its recurvature when Yancy landed Fujian Province, China on 20 August 1990.