Reinforced cables are usually installed on flexible airship structures to enhance their load-bearing capability. However, reinforced cables also increase the total weight of the airship. In order to find a balance between large loading-bear capability and light weight, a multi-objective optimisation scheme based on the genetic algorithm NSGA-II is put forward for the reinforced cable distribution on the airship. Firstly, different cable distribution schemes are presented according to engineering experience and the optimal one is determined by load analysis. Then, the CAE method and optimisation analysis are combined to achieve structure design optimisation. The parametric model of the airship structure with reinforced cables is established by ABAQUS secondary development and the load analysis is carried out. Parameter passing and optimisation algorithm are operated by Isight software and the optimisation analysis is conducted based on the NSGA-II algorithm. Finally, we draw some conclusions of the rules of optimised reinforcing cable distribution. The work of this paper has crucial engineering significance for improving performance of the airship structure design.

This study is focused on the development of longitudinal aerodynamic models for steady flight conditions. While several commercial solvers are available for this type of work, we seek to evaluate the accuracy of an open source software. This study aims to verify and demonstrate the accuracy of the OpenFoam solver when it is used on basic computers (32–64GB of RAM and eight cores). A new methodology was developed to show how an aerodynamic model of an aircraft could be designed using OpenFoam software. The mesh and the simulations were designed only using OpenFoam utilities, such as blockMesh, snappyHexMesh, simpleFoam and rhoSimpleFoam. For the methodology illustration, the process was applied to the Bombardier CRJ700 aircraft and simulations were performed for its flight envelope, up to M0.79. Forces and moments obtained with the OpenFoam model were compared with an accurate flight data source (level D flight simulator). Excellent results in data agreement were obtained with a maximum absolute error of 0.0026 for the drag coefficient, thus validating a high-fidelity aerodynamic model for the Bombardier CRJ-700 aircraft.

Orbiting objects in space are exposed to the risk of collision with space debris over their lifetime. Space debris orbiting in space experiences orbital decay due to various orbital perturbations. This work considers only orbital perturbations due to aerodynamic forces, which spacecraft experience due to the presence of a rarefied atmosphere, causing tumbling motion and orbital decay. Analysis of the orbital decay of a spacecraft is carried out by considering the variation of the drag coefficient as a function of its shape, motion and angle-of-attack. An in-house Direct Simulation Monte Carlo (DSMC) solver is modified for aerodynamic analysis of a spacecraft orbiting in the free molecular regime in low Earth orbit. In addition, an orbital dynamics model is developed to simulate the tumbling motion of a spacecraft and its orbital decay. The orbital decay trajectory is predicted for two sample spacecrafts using the aerodynamic coefficients obtained from the in-house DSMC solver as inputs to the orbital decay model. This study analyses and explores in detail the effects of the aerodynamic coefficients and shape of a spacecraft on its orbital decay.

The satellite attitude control system (SACS) is a complicated system. In order to reflect the relationship among different components in SACS and analyse the impact of component faults on system performance, a complete simulation platform of the SACS based on Simulink is built in this paper. With the embedding of the specific reaction flywheel, gyroscope and earth sensor model, and the design of the controller based on the quaternion feedback, the simulation platform can not only simulate the real SACS at the component level, but it can also realise the injection of component faults for analysing the system performance. Simulations are conducted to verify the performance of the simulation platform. Simulation results show that this simulation platform has the ability to accurately reflect the control performance of the SACS, and the output accuracy of the component model is high. The research results reveal that this simulation platform can provide model support for verifying the algorithm of fault diagnosis, prediction and tolerant control of the SACS. This simulation platform is easy to use and can be expanded and improved.

This paper introduces the Low Acoustic Noise and Turbulence (LANT) wind tunnel of the Sao Carlos School of Engineering, University of Sao Paulo (USP-EESC), Brazil. The closed-loop wind tunnel features several devices to improve flow uniformity, reduce swirl, and lower the background acoustic noise and turbulence, enabling stability and aeroacoustic experiments. The design criteria was based on the best practices reported, in particular for low turbulence wind tunnels. Yet, these criteria are conflicting and we discuss the decisions that had to be made and present flow quality results that were achieved. The 16-bladed axial fan with 13-blade stators is driven by a variable-speed electric motor. At the corners, 100 mm dense acoustic foam is installed on the vertical walls, floor and ceiling, and the turning vanes are filled with acoustic-absorbing material. The long settling chamber contains a 3.175 mm mesh hexagonal honeycomb and five fine mesh nylon screens, ending in a 7:1 area ratio short contraction. The 3-m long closed-working section has a $1\times 1\ {\rm m}^2$ cross-section area. At 15 m/s the working section wall boundary layer is less than 100 mm thick, providing an area of at least $800\times 800\ \mathrm{mm}^2$ where the streamwise flow uniformity was within 1%. In the 10–30 m/s flow speed range, the turbulence intensity ranged from 0.05% to 0.071% and the background acoustic noise level, obtained with an inflow microphone, ranged from 90 and 110 dB. A benchmark experiment on a flat plate boundary layer produced an almost perfect two-dimensional Blasius profile up to $Re_x \approx 2.5 \times 10^6$. A beamforming benchmark experiment on aeroacoustics accurately identified the sound emitted by a cylinder immersed in the flow.

Pilot-Induced Oscillation (PIO), although an old issue, still poses a significant threat to aviation safety. The introduction of new systems in modern aircraft modifies the human–machine interaction and makes it necessary for research to revisit the subject from time to time. Given the need of aircraft manufacturers to constantly perform PIO tests, this study analysed the feasibility of using three different computational pilot models (Tustin, Crossover and Precision) to simulate PIO conditions. Three aircraft models with different levels of propensity to PIO (original, low propensity and high propensity) were tested, as well as two pilot gain conditions (normal and high). Data were collected for a purely longitudinal synthetic task through simulations conducted in MATLAB®. PIO conditions were detect using a tuned PIO detection algorithm (ROVER). Data were analysed in terms of both whether the pilot models triggered a PIO condition and for how long the condition was sustained. The results indicated that the three pilot models only provoked PIO conditions when high gain inputs were applied. Additionally, Crossover was the only pilot model to trigger a PIO for the three aircraft models. There were also significant differences between the pilot models in the total PIO time, as the Tustin model typically sustained the oscillatory condition for longer.

Losses induced by tip clearance limit decisive improvements in the system efficiency and aerodynamic operational stability of aero-engine axial compressors. The tendency towards even lower blade heights to compensate for higher fluid densities aggravates their influence. Generally, it is emphasised that the tip clearance should be minimised but remain large enough to prevent collisions between the blade tip and the casing throughout the entire mission. The present work concentrates on the development of a preliminary aero-engine axial compressor casing design methodology involving meta-modelling techniques. Previous research work at the Institute for Turbomachinery and Flight Propulsion resulted in a Two-Dimensional (2D) axisymmetric finite element model for a generic multi-stage high-pressure axial compressor casing. Subsequent sensitivity studies led to the identification of significant parameters that are important for fine-tuning the tip clearance via specific flange design. This work is devoted to an exploration of the potential of surrogate modelling in preliminary compressor casing design with respect to rapid tip clearance assessments and its corresponding precision in comparison with finite element results. Reputed as data-driven mathematical approximation models and conceived for inexpensive numerical simulation result reproduction, surrogate models show even greater capacity when linked with extensive design space exploration and optimisation algorithms.

Compared with high-fidelity finite element simulations, the reductions obtained in computational time when using surrogate models amount to 99.9%. Validated via statistical methods and dependent on the size of the training database, the precision of surrogate models can reach down to the range of manufacturing tolerances. Subsequent inclusion of such surrogate models in a parametric optimisation process for tip clearance minimisation rapidly returned adaptions of the geometric design variables.

The commercial Computational Fluid Dynamics (CFD) software STAR-CCM+ was used to simulate the flow and breakup characteristics of a Liquid Jet Injected into the gaseous Crossflow (LJIC) under real engine operating conditions. The reasonable calculation domain geometry and flow boundary conditions were obtained based on a civil aviation engine performance model similar to the Leap-1B engine which was developed using the GasTurb software and the preliminary design results of its low-emission combustor. The Volume of Fluid (VOF) model was applied to simulate the breakup feature of the near field of LJIC. The numerical method was validated and calibrated through comparison with the public test data at atmospheric conditions. The results showed that the numerical method can capture most of the jet breakup structure and predict the jet trajectory with an error not exceeding ±5%. The verified numerical method was applied to simulate the breakup of LJIC at the real engine operating condition. The breakup mode of LJIC was shown to be surface shear breakup at elevated condition. The trajectory of the liquid jet showed good agreement with Ragucci’s empirical correlation.

In this paper, the accuracy and practical capabilities of three different reduced-order models (ROMs) are explored: an enhanced implicit condensation and expansion (EnICE) model, a finite element beam model, and a finite volume beam model are compared for their capability to accurately predict the nonlinear structural response of geometrically nonlinear built-up wing structures. This work briefly outlines the different order reduction methods, highlighting the associated assumptions and computational effort. The ROMs are then used to calculate the wing deflection for different representative load cases and these results are compared with the global finite element model (GFEM) predictions when possible. Overall, the ROMs are found to be able to capture the nonlinear GFEM behaviour accurately, but differences are noticed at very large displacements and rotations due to local geometrical effects.

The design of a two-dimensional supersonic inlet with large bleed window at low Mach number was developed. Numerical simulation and wind tunnel experiments were carried out to investigate the aerodynamic performance and variable geometric rules of the inlet. The result indicates that the single-degree-of-freedom variable geometry scheme adopted in this paper guarantees the steady work of the inlet over a wide speed range. The large bleed window caused by rotation of the compression ramp appears near the throat at low Mach number. Low-pressure airflow near the bleed window neutralises the original high-pressure airflow behind the shock train, which decreases the overall pressure of the downstream region of the internal contraction section. To match the lower pressure, the structure of the shock train changes from strong $\lambda$-type to weak $\lambda$-type, and finally to a normal shock wave as backpressure increases at Mach number 2.5. Herein, the total pressure recovery coefficient of the inlet near the critical condition improves by 8.5% as the backpressure ratio (Pe/P0) adds from 13 to 14.6 at Mach number 2.5. It proves that the scheme is effective on terminal shock wave control and inlet performance improvement. In addition, due to the background wave and the bleed window, two kinds of shock wave oscillation occur when the backpressure ratio is 13.1.

The aerodynamic performance of conventional aircraft configurations are mainly affected by the wing and horizontal tail. Drag reduction by shape optimisation of the wing, while taking into account the aircraft trimmed constraint, has more benefit than focusing solely on the wing. So in order to evaluate this approach, the following study presents results of a single and multipoint aerodynamic shape optimisation of the wing-body-tail configuration, defined by the Aerodynamic Design Discussion Group (ADODG). Most of the aerodynamic shape optimisation problems published in the last years are focused mainly on the wing as the main driver for performance improvement, with no trim constraint and/or excess drag obtained from the fuselage, fins or other parts. This work partially fills this gap by an investigation of RANS-based aerodynamic optimisation for transonic trimmed flight. Mesh warping and geometry parametrisation is accomplished by fitting the multi-block structured grid to a B-spline volumes and performing the mesh movement by using surface control points embedded within the free-form deformation (FFD) volumes. A gradient-based optimisation algorithm is used with an adjoint method in order to compute the derivatives of the objective and constraint functions with respect to the design variables. In this work the aerodynamic shape optimisation of the CRM wing-body-tail configuration is investigated, including a trim constraint that is satisfied by rotating the horizontal tail. The shape optimisation is driven by 432 design variables that envelope the wing surface, and 120 shape variables for the tail, as well as the angle of attack and tail rotation angles. The constraints are the lift coefficient, wing’s thickness controlled by 1,000 control points, and the wing’s volume. For the untrimmed configuration the drag coefficient is reduced by 5.76%. Optimising the wing with a trim condition by tail rotation results in shock-free design with a considerably improved drag, even better than the untrimmed-optimised case. The second optimisation problem studied is a single and multi-point lift constraint drag minimisation of a gliding configuration wing in transonic viscous flow. The shock is eliminated, reducing the drag of the untrimmed configuration by more than 60%, using 192 design variables. Further robustness is achieved through a multi-point optimisation with more than 45% drag reduction.

Large birds have evolved an effective wing anatomy and mechanics, enabling airborne mastery of manoeuvres and endurance. For these very reasons, they are difficult to replicate and study. The aim of the present work is to achieve active wing articulations to mimic natural bird flapping towards efficient and agile Unmanned Aerial Vehicles (UAVs). The proposed design, bio-mimicking the black-headed gull, Larus ridibundus, has three active and independent servo-controlled wing joints at the shoulder, elbow and wrist to achieve complex controls. The construction of the wing is realised through a polymeric skin and carbon fibre–epoxy composite spars and ribs. The wing movements (flapping, span reduction and twisting) envelopes of the full-scale robotic gull (Robogull) are examined using the Digital Image Correlation (DIC) technique and laser displacement sensing. Its aerodynamic performance was evaluated in a wind tunnel at various flapping parameters, wind speeds and angles of attack. It is observed that a flapping amplitude of 45$^\circ$ is more favourable than 90$^\circ$ for generating higher lift and thrust, while also depending on the presence of span reduction, twist and wind speed. The model performs better at a flying velocity of 4m/s as compared with 8m/s. Both lift and thrust are high at a higher flapping frequency of 2.5Hz. Combined variation of the flapping frequency and stroke ratio should be considered for better aerodynamic performance. The combination of a lower stroke ratio of 0.5 with a flapping frequency of 2.5Hz generates higher lift and thrust than other combinations. Span reduction and wing twist notably and independently enhance lift and thrust, respectively. An increase in the angle-of-attack increases lift but decreases thrust. The model can also generate a significant rolling moment when set at a bank angle of 20$^\circ$ and operated with independently controlled flapping amplitudes for the wings (45$^\circ$ for the left wing and 90$^\circ$ for the right wing). Based on the optimal values for the flapping amplitude (45$^\circ$), flapping frequency (2.5Hz) and flying velocity (4m/s), the Strouhal number (St) of the Robogull model is 0.24, lying in the optimal range ($0.2 < \mathrm{St} < 0.4$) for natural flyers and swimmers.

Increasingly agile manoeuvre is an advantage in the flight of aircraft, missiles and aerial vehicles, but the principles of accelerating aerodynamics in the transonic regime are only now being fully investigated. This study contributes to the understanding of shock and separation effects on drag during axial acceleration, using a simple geometric configuration. Unsteady shock wave behavior was numerically investigated for an axisymmetric cone-cylinder using a commercial solver and the Moving Reference Frame acceleration technique. This acceleration technique was validated using unsteady numerical and experimental methods. The cone-cylinder was accelerated from Mach number 0.6 to Mach number 1.2 at 100g constant and deceleration was from Mach number 1.2 until Mach number 0.6 at –100g constant. Three cone angles were tested for the cone-cylinder with uniform cylinder diameter. Acceleration through the transonic Mach regime was characterised by a delayed and gradual shock wave development when compared to steady state, demonstrating a clear flow history effect. Deceleration through the transonic Mach regime was characterised by shock wave propagation from the base to the nose. New flow structures appeared during deceleration that do not have counterparts in the steady state, including shock interactions and propagating expansion-compression features. Gross changes in the unsteady drag coefficient curves for each cone-angle are explained with reference to unsteady shock wave behaviour for accelerating and decelerating motion.

Camber morphing is an effective way to control the lift generated by any aerofoil and potentially improve the range (as measured by the lift-to-drag ratio) and endurance (as measured by $C_l^{3/2}/C_d$). This can be especially useful for fixed-wing Unmanned Aerial Vehicles (UAVs) undergoing different flying manoeuvres and flight phases. This work investigates the aerodynamic characteristics of the NACA0012 aerofoil morphed using a Single Corrugated Variable-Camber (SCVC) morphing approach. Structural analysis and morphed shapes are obtained based on small-deformation beam theory using chain calculations and validated using finite-element software. The aerofoil is then reconstructed from the camber line using a Radial Basis Function (RBF)-based interpolation method (J.H.S. Fincham and M.I. Friswell, “Aerodynamic optimisation of a camber morphing aerofoil,” Aerosp. Sci. Technol., 2015). The aerodynamic analysis is done by employing two different finite-volume solvers (OpenFOAM and ANSYS-Fluent) and a panel method code (XFoil). Results reveal that the aerodynamic coefficients predicted by the two finite-volume solvers using a fully turbulent flow assumption are similar but differ from those predicted by XFoil. The aerodynamic efficiency and endurance factor of morphed aerofoils indicate that morphing is beneficial at moderate to high lift requirements. Further, the optimal morphing angle increases with an increase in the required lift. Finally, it is observed for a fixed angle-of-attack that an optimum morphing angle exists for which the aerodynamic efficiency becomes maximum.

Arc-jet facilities have been the norm for ablation experiments used to calibrate computational models to date. However, the arc jet has a few major limitations and challenges, including non-uniform enthalpy distribution, non-equilibrium state, change of surface quality during testing and the extent of oxidation, to name but a few. A novel plasma technique for preheating axisymmetric heatshield samples in hypersonic impulse facilities is presented herein. The major aim of this innovative work is to help reduce the large variations of ablation rate predictions, space vehicle materials and missile design/testing, obtain strongly coupled hypersonic boundary layers and achieve lower cost of aerothermodynamics experiments. This present work remains one of the most highly anticipated solutions to maximise payload success and replicate high surface temperatures identical to those experienced by real flight vehicles. This work makes a useful contribution to re-entry studies under conditions that replicate the characteristics of re-entry flights. Future applications for the technique are expected to be found in hypersonic impulse facilities that can simulate the true flow energy under re-entry conditions.