Passive Control of Plume Interference on Slender Axisymmetric Bodies

The physics of the plume-induced shock and separation, particularly at high plume to exit pressure ratios with and without shock-turbulent boundary-layer control methods, were studied using computational techniques. Mass-averaged Navier-Stokes equations with a two-equation turbulence model were solved by using a fully implicit finite volume scheme and time.marching algorithm. The control methodologies for shock interactions included a porous tail and a porous extension attached at the nozzle exit or trailing edge. The porous tail produced a weaker shock and fixed the shock position on the control surface. The effect of the porous extension on shock interactions was mainly to restrain the plume from strongly underexpanding during a change in flight conditions. These techniques could give an additional dimension to the design and control of supersonic missiles.

Validation study for prediction of iced aerofoil aerodynamics

Ice accretions can significantly change the aerodynamic performance of wings and rotor blades. Significant performance degradation can occur when ice accreations cause regions of separated flow, to predict this change implies, at a minimum, the solution of the Reynolds-Averaged Navier-Stokes equations. This paper presents validation for two generic cases involving the flow over aerofoil sections with added synthetic ice shapes. Results were obtained for two aerofoils, namely the NACA 23012 and a generic multi-element configuration. These results are compared with force and pressure coefficient measurements obtained in the NASA LTPT wind-tunnel for the NACA 23012, and force, PIV and boundary-layer measurements obtained at DNW for the multi-clement case. The level of agreement is assessed in the context of industrial requirements.

Multiphysics Simulation of Thermal Plasma Produced by Lightning Strike on Aircraft Structures

A novel numerical technique is proposed to model thermal plasma of microseconds/milliseconds time-scale effect. Modelling thermal plasma due to lightning strike will allow the estimation of electric current density, plasma pressure, and heat flux at the surface of the aircraft structure. These input data can then be used for better estimation of the mechanical/thermal induced damage on the aircraft structures for better protection systems design. Thermal plasma generated during laser cutting, electric (laser) welding and other plasma processing techniques have been the focus of many researchers. Thermal plasma is a gaseous state that consists from a mixture of electrons, ions, and natural particles. Thermal plasma can be assumed to be in local thermodynamic equilibrium, which means the electrons and the heavy species have equal temperature. Different numerical techniques have been developed using a coupled Navier Stokes – Heat transfer – Electromagnetic equations based on the assumption that the thermal plasma is a single laminar gas flow. These previous efforts focused on generating thermal plasma of time-scale in the range of seconds. Lighting strike on aircraft structures generates thermal plasma of time-scale of milliseconds/microseconds, which makes the previous physics used not applicable. The difficulty comes from the Navier-Stokes equations as the fluid is simulated under shock load, this introducing significant changes in the density and temperature of the fluid.

On the use of OpenFOAM to model Oscillating wave surge converters

The Computational Fluid Dynamic (CFD) toolbox OpenFOAM is used to assess the applicability of Reynolds-Averaged Navier-Stokes (RANS) solvers to the simulation of Oscillating Wave Surge Converters (OWSC) in significant waves. Simulation of these flap type devices requires the solution of the equations of motion and the representation of the OWSC’s motion in a moving mesh. A new way to simulate the sea floor inside a section of the moving mesh with a moving dissipation zone is presented. To assess the accuracy of the new solver, experiments are conducted in regular and irregular wave traces for a full three dimensional model. Results of acceleration and flow features are presented for numerical and experimental data. It is found that the new numerical model reproduces experimental results within the bounds of experimental accuracy.

Detection of improper installation from the sensor signal of vortex flowmeters

Effects of inappropriate installation can bias the measurements of flowmeters. For vortex flowmeters, a method is proposed to detect inappropriate installation of the flowmeter from the oscillatory signal of the vortex sensor. The method is based on assuming the process of vortex generation to be a generic, noisy, nonlinear oscillation, describable by a noisy Stuart-Landau equation, with a corresponding sensor signal that also contains higher harmonic excitations. By making use of the scaling properties of the Navier-Stokes Equation, the method was designed to be robust with respect to uncertainties in the fluid properties. The diagnostic functionality is demonstrated on measurement data. In the experiments, installation effects that lead to more than 0.5% error in the output of the flowmeter could clearly be detected. (C) 2003 Elsevier Ltd. All rights reserved.

Development of the kinetic model of platinum catalyzed ammonia oxidation in a microreactor

The ammonia oxidation reaction on supported polycrystalline platinum catalyst was investigated in an aluminum-based microreactor. An extensive set of reactions was included in the chemical reactor modeling to facilitate the construction of a kinetic model capable of satisfactory predictions for a wide range of conditions (NH3 partial pressure, 0.01-0.12 atm; O-2 partial pressure, 0.10-0.88 atm; temperature, 523-673 K; contact time, 0.3-0.7 ms). The elementary surface reactions used in developing the mechanism were chosen based on the literature data concerning ammonia oxidation on a Pt catalyst. Parameter estimates for the kinetic model were obtained using multi-response least squares regression analysis using the isothermal plug-flow reactor approximation. To evaluate the model, the behavior of a microstructured reactor was simulated by means of a complete Navier-Stokes model accounting for the reactions on the catalyst surface and the effect of temperature on the physico-chemical properties of the reacting mixture. In this way, the effect of the catalytic wall temperature non-uniformity and the effect of a boundary layer on the ammonia conversion and selectivity were examined. After further optimization of appropriate kinetic parameters, the calculated selectivities and product yields agree very well with the values actually measured in the microreactor. (C) 2002 Elsevier Science B.V. All rights reserved.

The opportunities and limitations of using CFD in the development of wave energy converters

Most Wave Energy Converters (WECs) being developed are fundamentally different from known marine structures. Limited experience is a fundamental challenge for the design, especially issues concerning load assumptions and power estimates. Reynolds-Averaged Navier-Stokes (RANS) CFD methods are being used successfully in many areas of marine engineering. They have been shown to accurately simulate many hydrodynamic effects and are a helpful tool for investigating complex flows. The major drawback is the significant computational power required and the associated overhead with pre and post-processing. This paper presents the challenges and advantages in the application of RANS CFD methods in the design process of a wave energy converter and compares the time, labour and ultimately financial requirements for obtaining practical results.

Modeling Multivariate Interest Rates Using Time-Varying Copulas and Reducible Nonlinear Stochastic Differential Equations

We propose a new approach for modeling nonlinear multivariate interest rate processes based on time-varying copulas and reducible stochastic differential equations (SDEs). In the modeling of the marginal processes, we consider a class of nonlinear SDEs that are reducible to Ornstein--Uhlenbeck (OU) process or Cox, Ingersoll, and Ross (1985) (CIR) process. The reducibility is achieved via a nonlinear transformation function. The main advantage of this approach is that these SDEs can account for nonlinear features, observed in short-term interest rate series, while at the same time leading to exact discretization and closed-form likelihood functions. Although a rich set of specifications may be entertained, our exposition focuses on a couple of nonlinear constant elasticity volatility (CEV) processes, denoted as OU-CEV and CIR-CEV, respectively. These two processes encompass a number of existing models that have closed-form likelihood functions. The transition density, the conditional distribution function, and the steady-state density function are derived in closed form as well as the conditional and unconditional moments for both processes. In order to obtain a more flexible functional form over time, we allow the transformation function to be time varying. Results from our study of U.S. and UK short-term interest rates suggest that the new models outperform existing parametric models with closed-form likelihood functions. We also find the time-varying effects in the transformation functions statistically significant. To examine the joint behavior of interest rate series, we propose flexible nonlinear multivariate models by joining univariate nonlinear processes via appropriate copulas. We study the conditional dependence structure of the two rates using Patton (2006a) time-varying symmetrized Joe--Clayton copula. We find evidence of asymmetric dependence between the two rates, and that the level of dependence is positively related to the level of the two rates. (JEL: C13, C32, G12) Copyright The Author 2010. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oxfordjournals.org, Oxford University Press.

STEADY-STATE SOLUTIONS IN NONLINEAR DIFFUSIVE SHOCK ACCELERATION

Stationary solutions to the equations of nonlinear diffusive shock acceleration play a fundamental role in the theory of cosmic-ray acceleration. Their existence usually requires that a fraction of the accelerated particles be allowed to escape from the system. Because the scattering mean free path is thought to be an increasing function of energy, this condition is conventionally implemented as an upper cutoff in energy space-particles are then permitted to escape from any part of the system, once their energy exceeds this limit. However, because accelerated particles are responsible for the substantial amplification of the ambient magnetic field in a region upstream of the shock front, we examine an alternative approach in which particles escape over a spatial boundary. We use a simple iterative scheme that constructs stationary numerical solutions to the coupled kinetic and hydrodynamic equations. For parameters appropriate for supernova remnants, we find stationary solutions with efficient acceleration when the escape boundary is placed at the point where growth and advection of strongly driven nonresonant waves are in balance. We also present the energy dependence of the distribution function close to the energy where it cuts off-a diagnostic that is in principle accessible to observation.

Field test of multiple 1/10 scale tidal turbines in steady flows

Queen's University Belfast and Wave Barrier Ltd have developed a tidal testing platform to test hydrokinetic turbines at medium scale. Multiple turbines can be pushed through still water conditions, in steady-state pushing tests. Experiments were conducted to evaluate the interactions between two identical, mono-strut, horizontal axis tidal turbines (HATTs) of 1.5 m diameter (D) rotor. Their relative performance when located individually, in-plane and in-line are investigated. The data shows a high consistency in the power curves at different flow speeds, which indicates high repeatability in this Reynolds range. For an individual turbine, there is no performance difference when the rotor is mounted either upstream or downstream of the supporting structure. When placed in-plane, the turbines have no adverse effect on one another. When spaced in-line with 2D separation, there is a 63% reduction in the performance of the downstream turbine. At 6D downstream this performance reduction is still 59%, indicating some wake recovery between 2D and 6D, though the influence from the upstream rotor persists to at least 6D downstream of the first device. In contrast the performance of the downstream turbine when placed at 1.5D offset of the upstream device at 6D downstream is approximately recovered to the individual turbine performance.

Nonequilibrium processes from generalized Langevin equations: Realistic nanoscale systems connected to two thermal baths

We extend the generalized Langevin equation (GLE) method [L. Stella, C. D. Lorenz, and L. Kantorovich, Phys. Rev. B 89, 134303 (2014)] to model a central classical region connected to two realistic thermal baths at two different temperatures. In such nonequilibrium conditions a heat flow is established, via the central system, in between the two baths. The GLE-2B (GLE two baths) scheme permits us to have a realistic description of both the dissipative central system and its surrounding baths. Following the original GLE approach, the extended Langevin dynamics scheme is modified to take into account two sets of auxiliary degrees of freedom corresponding to the mapping of the vibrational properties of each bath. These auxiliary variables are then used to solve the non-Markovian dissipative dynamics of the central region. The resulting algorithm is used to study a model of a short Al nanowire connected to two baths. The results of the simulations using the GLE-2B approach are compared to the results of other simulations that were carried out using standard thermostatting approaches (based on Markovian Langevin and Nosé-Hoover thermostats). We concentrate on the steady-state regime and study the establishment of a local temperature profile within the system. The conditions for obtaining a flat profile or a temperature gradient are examined in detail, in agreement with earlier studies. The results show that the GLE-2B approach is able to treat, within a single scheme, two widely different thermal transport regimes, i.e., ballistic systems, with no temperature gradient, and diffusive systems with a temperature gradient.