Fluid-structure interaction simulations of a tension-cone inflatable aerodynamic decelerator

A series of fluid-structure interaction simulations of an aerodynamic tension-cone supersonic decelerator prototype intended for large mass payload deployment in planetary explorations are discussed. The fluid-structure interaction computations combine large deformation analysis of thin shells with large-eddy simulation of compressible turbulent flows using a loosely coupled approach to enable quantification of the dynamics of the vehicle. The simulation results are compared with experiments carried out at the NASA Glenn Research Center. Reasonably good agreement between the simulations and the experiment is observed throughout a deflation cycle. The simulations help to illuminate the details of the dynamic progressive buckling of the tension-cone decelerator that ultimately results in the collapse of the structure as the inflation pressure is decreased. Furthermore, the tension-cone decelerator exhibits a transient oscillatory behavior under impulsive loading that ultimately dies out. The frequency of these oscillations was determined to be related to the acoustic time scale in the compressed subsonic region between the bow shock and the structure. As shown, when the natural frequency of the structure and the frequency of the compressed subsonic region approximately match, the decelerator exhibits relatively large nonaxisymetric oscillations. The observed response appears to be a fluid-structure interaction resonance resulting from an acoustic chamber (pistonlike) mode exciting the structure. Copyright © 2013 by Christopher Porter, R. Mark Rennie, Eric J. Jumper.

Large eddy simulations in low-pressure turbines: Effect of wakes at elevated free-stream turbulence

The transition of a separated shear layer over a flat plate, in the presence of periodic wakes and elevated free-stream turbulence (FST), is numerically investigated using Large Eddy Simulation (LES). The upper wall of the test section is inviscid and specifically contoured to impose a streamwise pressure distribution over the flat plate to simulate the suction surface of a low-pressure turbine (LPT) blade. Two different distributions representative of a 'high-lift' and an 'ultra high-lift' turbine blade are examined. Results obtained from the current LES compare favourably with the extensive experimental data previously obtained for these configurations. The LES results are then used to further investigate the flow physics involved in the transition process.In line with experimental experience, the benefit of wakes and FST obtained by suppressing the separation bubble, is more pronounced in 'ultra high-lift' design when compared to the 'high-lift' design. Stronger 'Klebanoff streaks' are formed in the presence of wakes when compared to the streaks due to FST alone. These streaks promoted much early transition. The weak Klebanoff streaks due to FST continued to trigger transition in between the wake passing cycles.The experimental inference regarding the origin of Klebanoff streaks at the leading edge has been confirmed by the current simulations. While the wake convects at local free-stream velocity, its impression in the boundary layer in the form of streaks convects much slowly. The 'part-span' Kelvin-Helmholtz structures, which were observed in the experiments when the wake passes over the separation bubble, are also captured. The non-phase averaged space-time plots manifest that reattachment is a localized process across the span unlike the impression of global reattachment portrayed by phase averaging. © 2013 Elsevier Inc.

Numerical simulations supporting the process design of ring rolling processes

In conventional Finite Element Analysis (FEA) of radial-axial ring rolling (RAR) the motions of all tools are usually defined prior to simulation in the preprocessing step. However, the real process holds up to 8 degrees of freedom (DOF) that are controlled by industrial control systems according to actual sensor values and preselected control strategies. Since the histories of the motions are unknown before the experiment and are dependent on sensor data, the conventional FEA cannot represent the process before experiment. In order to enable the usage of FEA in the process design stage, this approach integrates the industrially applied control algorithms of the real process including all relevant sensors and actuators into the FE model of ring rolling. Additionally, the process design of a novel process 'the axial profiling', in which a profiled roll is used for rolling axially profiled rings, is supported by FEA. Using this approach suitable control strategies can be tested in virtual environment before processing. © 2013 AIP Publishing LLC.

Thermo-hydro-geomechanical simulations of methane gas production from deep sea methane hydrate formations

A new constitutive model called Methane Hydrate Critical State (MHCS) model was conducted to investigate the geomechanical response of the gas-hydrate-bearing sediments at the Nankai Trough during the wellbore construction process. The strength and dilatancy of gas-hydrate-bearing soil would gradually disappear when the bonds are destroyed because of excessively shearing, which are often observed in dense soils and also in bonded soils such as cemented soil and unsaturated soil. In this study, the MHCS model, which presents such softening features, would be incorporated into a staged-finite-element model in ABAQUS, which mainly considered the loading history of soils and the interaction between cement-casing-formation. This model shows the influence of gas-hydrate-bearing soil to the deformation and stability of a wellbore and the surrounding sediments during wellbore construction. At the same time, the conventional Mohr-Coulomb model was used in the model to show the advantages of MHCS model by comparing the results of the two models.

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.

Simulations of autoignition and laminar premixed flames in methane/air mixtures diluted with hot products

This article considers constant-pressure autoignition and freely propagating premixed flames of cold methane/air mixtures mixed with equilibrium hot products at high enough dilution levels to burn within the moderate to intense low oxygen dilution (MILD) combustion regime. The analysis is meant to provide further insight on MILD regime boundaries and to identify the effect of hot products speciation. As the mass fraction of hot products in the reactants mixture increases, autoignition occurs earlier. Species profiles show that the products/reactants mixture approximately equilibrates to a new state over a quick transient well before the main autoignition event, but as dilution becomes very high, this equilibration transient becomes more prominent and eventually merges with the primary ignition event. The dilution level at which these two reactive zones merge corresponds well with that marking the transition into the MILD regime, as defined according to conventional criteria. Similarly, premixed flame simulations at high dilutions show evidence of significant reactions involving intermediate species prior to the flame front. Since the premixed flame governing equations system demands that the species and temperature gradients be zero at the "cold" boundary, flame speed cannot be calculated above a certain dilution level. Up to this point, which again agrees reasonably well with the transition into the MILD regime according to convention, the laminar burning velocity was found to increase with hot product dilution while flame thickness remained largely unchanged. Some comments on the MILD combustion regime boundary definition for gas turbine applications are included. Copyright © Taylor & Francis Group, LLC.

Monitoring and preventing numerical oscillations in 3D simulations with coupled Monte Carlo codes

Previous studies have reported that different schemes for coupling Monte Carlo (MC) neutron transport with burnup and thermal hydraulic feedbacks may potentially be numerically unstable. This issue can be resolved by application of implicit methods, such as the stochastic implicit mid-point (SIMP) methods. In order to assure numerical stability, the new methods do require additional computational effort. The instability issue however, is problem-dependent and does not necessarily occur in all cases. Therefore, blind application of the unconditionally stable coupling schemes, and thus incurring extra computational costs, may not always be necessary. In this paper, we attempt to develop an intelligent diagnostic mechanism, which will monitor numerical stability of the calculations and, if necessary, switch from simple and fast coupling scheme to more computationally expensive but unconditionally stable one. To illustrate this diagnostic mechanism, we performed a coupled burnup and TH analysis of a single BWR fuel assembly. The results indicate that the developed algorithm can be easily implemented in any MC based code for monitoring of numerical instabilities. The proposed monitoring method has negligible impact on the calculation time even for realistic 3D multi-region full core calculations. © 2014 Elsevier Ltd. All rights reserved.