Bruno Savard

A computationally-efficient, semi-implicit, iterative method for the time-integration of reacting flows with stiff chemistry

A semi-implicit preconditioned iterative method is proposed for the time-integration of the stiff chemistry in simulations of unsteady reacting flows, such as turbulent flames, using detailed chemical kinetic mechanisms. Emphasis is placed on the simultaneous treatment of convection, diffusion, and chemistry, without using operator splitting techniques. The preconditioner corresponds to an approximation of the diagonal of the chemical Jacobian. Upon convergence of the sub-iterations, the fully-implicit, second-order time-accurate, Crank-Nicolson formulation is recovered. Performance of the proposed method is tested theoretically and numerically on one-dimensional laminar and three-dimensional high Karlovitz turbulent premixed n-heptane/air flames. The species lifetimes contained in the diagonal preconditioner are found to capture all critical small chemical timescales, such that the largest stable time step size for the simulation of the turbulent flame with the proposed method is limited by the convective CFL, rather than chemistry. The theoretical and numerical stability limits are in good agreement and are independent of the number of sub-iterations. The results indicate that the overall procedure is second-order accurate in time, free of lagging errors, and the cost per iteration is similar to that of an explicit time integration. The theoretical analysis is extended to a wide range of flames (premixed and non-premixed), unburnt conditions, fuels, and chemical mechanisms. In all cases, the proposed method is found (theoretically) to be stable and to provide good convergence rate for the sub-iterations up to a time step size larger than 1 µs. This makes the proposed method ideal for the simulation of turbulent flames.

Flight Control Using Wing-Tip Plasma Actuation

A concept for lift modification on a conventional aircraft wing for roll control at low angle of attack with dielectric barrier discharge plasma actuators is proposed and assessed through computational fluid dynamics simulations and preliminary wind-tunnel experiments. The concept consists of placing plasma actuators around the wing tip to add momentum in the direction opposite to that of the flow forming the tip vortex. Because of the limited strength of existing plasma actuators, the assessment is carried out for a relatively small two-dimensional wing (NACA 4418) with a rounded tip at zero angle of attack and 15 m/s for a Reynolds number in the range of 1.5 x 10 5 . Computational fluid dynamics simulations show a significant alteration of the vorticity field downstream of the trailing edge characterized by a more diffused vortex surrounded by zones of negative vorticity induced by the actuators and, but not necessarily, outboard displacement of the tip vortex. This leads to a reduced downwash that results in a change in lift of up to almost 20% for actuator strength levels that should be achievable in the short term with a new generation of dielectric barrier discharge actuators. The actuator placed on the suction side contributes the most to the lift increase, with its induced jet blocking the flow around the wind tip at the origin of the formation of the tip vortex. Wind-tunnel experimental results support the computational fluid dynamics predictions in both magnitude and trend. Furthermore, preliminary computational fluid dynamics simulations are carried out for a symmetric nonlifting wing (NACA 0018), representative of aircraft tail surfaces at zero angle of attack to generate lift for pitch and yaw control. Results indicate lift generation that increases and becomes larger than drag at higher actuator strengths. These promising results show a potential for the proposed concept to replace movable flight control surfaces on future aircraft wings and empennages.