The 2D lumped wake vortex method applied to unsteady flapping

2020 ◽  
Vol 8 (3) ◽  
pp. 225-238
Author(s):  
David J. Talarico ◽  
Aaron Mazzeo ◽  
Mitsunori Denda

PurposeAdvancements in aerospace technologies, which rely on unsteady fluid dynamics, are being hindered by a lack of easy to use, computationally efficient unsteady computational fluid dynamics (CFD) software. Existing CFD platforms are capable of handling unsteady flapping, but the time, money and expertise required to run even a basic flapping simulation make design iteration and optimization prohibitively expensive for the average researcher.Design/methodology/approachIn the present paper, a remedy to model the effects of viscosity is introduced to the original vortex method, in which the pitching moment amplitude grew over time for simulations involving multiple flapping cycles. The new approach described herein lumps far-field wake vortices to mimic the vortex decay, which is shown to improve the accuracy of the solution while keeping the pitching moment amplitude under control, especially for simulations involving many flapping cycles.FindingsIn addition to improving the accuracy of the solution, the new method greatly reduces the computation time for simulations involving many flapping cycles. The solution of the original vortex method and the new method are compared to published Navier–Stokes solver data and show very good agreement.Originality/valueBy utilizing a novel unsteady vortex method, which has been designed specifically to handle the highly unsteady flapping wing problems, it has been shown that the time to compute a solution is reduced by several orders of magnitude (Denda et al., 2016). Despite the success of the vortex method, especially for a small number of flapping cycles, the solution deteriorates as the number of flapping cycles increases due to the inherent lack of viscosity in the vortex method.

2019 ◽  
Vol 92 (2) ◽  
pp. 145-155
Author(s):  
Kheireddine Choutri ◽  
Mohand Lagha ◽  
Laurent Dala

Purpose This paper aims to propose a new multi-layered optimal navigation system that jointly optimizes the energy consumption, improves the robustness and raises the performance of a quadrotor unmanned aerial vehicle (UAV). Design/methodology/approach The proposed system is designed as a multi-layered system. First, the control architecture layer links the input and the output spaces via quaternion-based differential flatness equations. Then, the trajectory generation layer determines the optimal reference path and avoids obstacles to secure the UAV from collisions. Finally, the control layer allows the quadrotor to track the generated path and guarantees the stability using a double loop non-linear optimal backstepping controller (OBS). Findings All the obtained results are confirmed using several scenarios in different situations to prove the accuracy, energy optimization and the robustness of the designed system. Practical implications The proposed controllers are easily implementable on-board and are computationally efficient. Originality/value The originality of this research is the design of a multi-layered optimal navigation system for quadrotor UAV. The proposed control architecture presents a direct relation between the states and their derivatives, which then simplifies the trajectory generation problem. Furthermore, the derived differentially flat equations allow optimization to occur within the output space as opposed to the control space. This is beneficial because constraints such as obstacle avoidance occur in the output space; hence, the computation time for constraint handling is reduced. For the OBS, the novelty is that all controller parameters are derived using the multi-objective genetic algorithm (MO-GA) that optimizes all the quadrotor state’s cost functions jointly.


2019 ◽  
Vol 91 (2) ◽  
pp. 309-316 ◽  
Author(s):  
Jernej Drofelnik ◽  
Andrea Da Ronch ◽  
Matteo Franciolini ◽  
Andrea Crivellini

Purpose This paper aims to present a numerical method based on computational fluid dynamics that allows investigating the buffet envelope of reference equivalent wings at the equivalent cost of several two-dimensional, unsteady, turbulent flow analyses. The method bridges the gap between semi-empirical relations, generally dominant in the early phases of aircraft design, and three-dimensional turbulent flow analyses, characterised by high costs in analysis setups and prohibitive computing times. Design/methodology/approach Accuracy in the predictions and efficiency in the solution are two key aspects. Accuracy is maintained by solving a specialised form of the Reynolds-averaged Navier–Stokes equations valid for infinite-swept wing flows. Efficiency of the solution is reached by a novel implementation of the flow solver, as well as by combining solutions of different fidelity spatially. Findings Discovering the buffet envelope of a set of reference equivalent wings is accompanied with an estimate of the uncertainties in the numerical predictions. Just over 2,000 processor hours are needed if it is admissible to deal with an uncertainty of ±1.0° in the angle of attack at which buffet onset/offset occurs. Halving the uncertainty requires significantly more computing resources, close to a factor 200 compared with the larger uncertainty case. Practical implications To permit the use of the proposed method as a practical design tool in the conceptual/preliminary aircraft design phases, the method offers the designer with the ability to gauge the sensitivity of buffet on primary design variables, such as wing sweep angle and chord to thickness ratio. Originality/value The infinite-swept wing, unsteady Reynolds-averaged Navier–Stokes equations have been successfully applied, for the first time, to identify buffeting conditions. This demonstrates the adequateness of the proposed method in the conceptual/preliminary aircraft design phases.


Author(s):  
M. R. Amiralaei ◽  
H. Alighanbari ◽  
S. M. Hashemi

The objective of the present study is to investigate the low Reynolds number (LRN) fluid dynamics of two airfoils in pitching oscillations. The airfoils are in a tandem configuration and perform in-phase oscillations. Navier-Stokes (NS) equations with Finite Volume Method (FVM) are used and the instantaneous aerodynamic force coefficients are analyzed. The effect of amplitude of pitching oscillations and Re are investigated on the fluid forces. It is found that the amplitude of pitching oscillations is of primary importance to the fluid forces, affecting them quantitatively and qualitatively. Re is found to be of secondary importance compared to the effects of the amplitude of pitching oscillations. It mainly affects the magnitude of the forces.


2016 ◽  
Vol 26 (8) ◽  
pp. 2462-2485 ◽  
Author(s):  
Juan Wen ◽  
Yinnian He ◽  
Xin Zhao

Purpose The purpose of this paper is to propose a new stabilized finite volume element method for the Navier-Stokes problem. Design/methodology/approach This new method is based on the multiscale enrichment and uses the lowest equal order finite element pairs P1/P1. Findings The stability and convergence of the optimal order in H1-norm for velocity and L2-norm for pressure are obtained. Originality/value Using a dual problem for the Navier-Stokes problem, the convergence of the optimal order in L2-norm for the velocity is obtained. Finally, numerical example confirms the theory analysis and validates the effectiveness of this new method.


2014 ◽  
Vol 24 (6) ◽  
pp. 1290-1300 ◽  
Author(s):  
Haiming Huang ◽  
Guo Huang ◽  
Xiaoliang Xu ◽  
Weijie Li

Purpose – Relevant analyses are presented on the base of the compressible vortex method for simulating the development of two or three co-rotating vortices with different characteristic Mach numbers. The paper aims to discuss this issue. Design/methodology/approach – In addition to having vorticity and dilatation properties, the vortex particles also carry density, enthalpy, and entropy. Taking co-rotating vortices in two-dimensional unsteady compressible flow for an example, truncation of unbounded domains via a nonreflecting boundary condition was considered in order to make the method computationally efficient. Findings – For two identical vortices, the effect of the vortex Mach number on merging process is not evident; if two vortices have the same circulation rather than different radiuses, the vorticity and dilatation fields of the vortex under a vortex Mach number will be absorbed by the vortex under a higher vortex Mach number. For three vortices, if the original arrangement of the vortices is changed, the evolvement of the vorticity and dilatation fields is different. Originality/value – The paper reveals new mechanism of the three co-rotating vortices by a feasible compressible vortex method.


2022 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Adam Targui ◽  
Wagdi George Habashi

Purpose Responsible for lift generation, the helicopter rotor is an essential component to protect against ice accretion. As rotorcraft present a smaller wing cross-section and a lower available onboard power compared to aircraft, electro-thermal heating pads are favored as they conform to the blades’ slender profile and limited volume. Their optimization is carried out here taking into account, for the first time, the highly three-dimensional (3D) nature of the flow and ice accretion, in contrast to the current state-of-the-art that is limited to two-dimensional (2D) airfoils. Design/methodology/approach Conjugate heat transfer simulation results are provided by the truly 3D finite element Navier–Stokes analysis package-ICE code, embedded in a proprietary rotorcraft simulation toolkit, with reduced-order modeling providing a time-efficient evaluation of the objective and constraint functions at every iteration. The proposed methodology optimizes heating pads extent and power usage and is versatile enough to address in a computationally efficient manner a wide variety of optimization formulations. Findings Low-error reduced-order modeling strategies are introduced to make the tackling of complex 3D geometries feasible in todays’ computers, with the developed framework applied to four problem formulations, demonstrating marked reductions to power consumption along with improved aerodynamics. Originality/value The present paper proposes a 3D framework for the optimization of electro-thermal rotorcraft ice protection systems, in hover and forward flight. The current state-of-the-art is limited to 2D airfoils.


2016 ◽  
Vol 73 (10) ◽  
pp. 4137-4149 ◽  
Author(s):  
Laurent Makké ◽  
Luc Musson-Genon ◽  
Bertrand Carissimo ◽  
Pierre Plion ◽  
Maya Milliez ◽  
...  

Abstract The atmospheric radiation field has seen the development of more accurate and faster methods to take into account absorption. Modeling fog formation, where infrared radiation is involved, requires accurate methods to compute cooling rates. Radiative fog appears under clear-sky conditions owing to a significant cooling during the night where absorption and emission are the dominant processes. Thanks to high-performance computing, high-resolution multispectral approaches to solving the radiative transfer equation are often used. Nevertheless, the coupling of three-dimensional radiative transfer with fluid dynamics is very computationally expensive. Radiation increases the computation time by around 50% over the pure computational fluid dynamics simulation. To reduce the time spent in radiation calculations, a new method using analytical absorption functions fitted by Sasamori on Yamamoto’s radiation chart has been developed to compute an equivalent absorption coefficient (spectrally integrated). Only one solution of the radiative transfer equation is needed against Nband × Ngauss for an Nband model with Ngauss quadrature points on each band. A comparison with simulation data has been done and the new parameterization of radiative properties proposed in this article shows the ability to handle variations of gas concentrations and liquid water.


2019 ◽  
Vol 13 (2) ◽  
pp. 174-180
Author(s):  
Poonam Sharma ◽  
Ashwani Kumar Dubey ◽  
Ayush Goyal

Background: With the growing demand of image processing and the use of Digital Signal Processors (DSP), the efficiency of the Multipliers and Accumulators has become a bottleneck to get through. We revised a few patents on an Application Specific Instruction Set Processor (ASIP), where the design considerations are proposed for application-specific computing in an efficient way to enhance the throughput. Objective: The study aims to develop and analyze a computationally efficient method to optimize the speed performance of MAC. Methods: The work presented here proposes the design of an Application Specific Instruction Set Processor, exploiting a Multiplier Accumulator integrated as the dedicated hardware. This MAC is optimized for high-speed performance and is the application-specific part of the processor; here it can be the DSP block of an image processor while a 16-bit Reduced Instruction Set Computer (RISC) processor core gives the flexibility to the design for any computing. The design was emulated on a Xilinx Field Programmable Gate Array (FPGA) and tested for various real-time computing. Results: The synthesis of the hardware logic on FPGA tools gave the operating frequencies of the legacy methods and the proposed method, the simulation of the logic verified the functionality. Conclusion: With the proposed method, a significant improvement of 16% increase in throughput has been observed for 256 steps iterations of multiplier and accumulators on an 8-bit sample data. Such an improvement can help in reducing the computation time in many digital signal processing applications where multiplication and addition are done iteratively.


Author(s):  
Claudio Miccoli ◽  
Alessandro Turchi ◽  
Pierre Schrooyen ◽  
Domenic D’Ambrosio ◽  
Thierry Magin

AbstractThis work deals with the analysis of the cork P50, an ablative thermal protection material (TPM) used for the heat shield of the qarman Re-entry CubeSat. Developed for the European Space Agency (ESA) at the von Karman Institute (VKI) for Fluid Dynamics, qarman is a scientific demonstrator for Aerothermodynamic Research. The ability to model and predict the atypical behavior of the new cork-based materials is considered a critical research topic. Therefore, this work is motivated by the need to develop a numerical model able to respond to this demand, in preparation to the post-flight analysis of qarman. This study is focused on the main thermal response phenomena of the cork P50: pyrolysis and swelling. Pyrolysis was analyzed by means of the multi-physics Computational Fluid Dynamics (CFD) code argo, developed at Cenaero. Based on a unified flow-material solver, the Volume Averaged Navier–Stokes (VANS) equations were numerically solved to describe the interaction between a multi-species high enthalpy flow and a reactive porous medium, by means of a high-order Discontinuous Galerkin Method (DGM). Specifically, an accurate method to compute the pyrolysis production rate was implemented. The modeling of swelling was the most ambitious task, requiring the development of a physical model accounting for this phenomenon, for the purpose of a future implementation within argo. A 1D model was proposed, mainly based on an a priori assumption on the swelling velocity and the resolution of a nonlinear advection equation, by means of a Finite Difference Method (FDM). Once developed, the model was successfully tested through a matlab code, showing that the approach is promising and thus opening the way to further developments.


2019 ◽  
Vol 72 (5) ◽  
pp. 557-565
Author(s):  
Dilek Bulut ◽  
Tatjana Krups ◽  
Gerhard Poll ◽  
Ulrich Giese

Purpose Elastomer seals are used in many applications. They are exposed to lubricants and additives at elevated temperatures, as well as mechanical stresses. They can only provide good sealing function when they have resistance to those factors. There are many elastomer-lubricant compatibility tests based on DIN ISO 1817 in industry. However, they are insufficient and costly. Correlations between the tests and the applications are inadequate. The purpose of this study is investigating lubricant compatibility of fluoroelastomers (FKM) seals in polyethylene-glycol (PG)- and polyalphaolefin (PAO)- based synthetic oils and developing a methodology to predict seal service life. Design/methodology/approach A new compatibility test which is more sufficient in terms of time and cost was developed and compared with a standard test, currently used in industry. Compatibility of FKM radial lip seals with PG- and PAO-based synthetic oils with different additives was investigated chemically and dynamically. Failure mechanisms were examined. Findings The new method and the Freudenberg Flender Test FB 73 11 008 showed similar results concerning damages and similar tendencies regarding wear. The additive imidazole derivative was the most critical. Static tests give indications of possible chemically active additives, but alone they are insufficient to simulate the dynamic applications. Originality/value The paper describes a new method to investigate elastomer-lubricant compatibility and gives first results with a variety of lubricants.


Sign in / Sign up

Export Citation Format

Share Document