Transient Dynamic System Behavior of Pressure Actuated Cellular Structures in a Morphing Wing

High aspect ratio aircraft have a significantly reduced induced drag, but have only limited installation space for control surfaces near the wingtip. This paper describes a multidisciplinary design methodology for a morphing aileron that is based on pressure-actuated cellular structures (PACS). The focus of this work is on the transient dynamic system behavior of the multi-functional aileron. Decisive design aspects are the actuation speed, the resistance against external loads, and constraints preparing for a future wind tunnel test. The structural stiffness under varying aerodynamic loads is examined while using a reduced-order truss model and a high-fidelity finite element analysis. The simulations of the internal flow investigate the transient pressurization process that limits the dynamic actuator response. The authors present a reduced-order model based on the Pseudo Bond Graph methodology enabling time-efficient flow simulation and compare the results to computational fluid dynamic simulations. The findings of this work demonstrate high structural resistance against external forces and the feasibility of high actuation speeds over the entire operating envelope. Future research will incorporate the fluid–structure interaction and the assessment of load alleviation capability.


Introduction
High aspect ratio (HAR) aircraft are in the focus of current research for their potential to significantly reduce induced drag [1]. The resulting long and slender wings provide little space for control surfaces, especially near the wingtip. The highest aspect ratios can be achieved with folding wingtips (FWT), which ensure the aircraft comply with airport gate size restrictions on ground. FWT are already installed on the Boeing 777X [2] and they are also under development in Europe [3,4]. The installation space in the FWT is particularly limited and the mass of the secondary systems must also be kept to a minimum. In HAR wings, conventional linear hydraulic and electro-mechanical actuators reach their limits regarding installation space and accessibility.
Shape-variable structures require a comparatively small installation space. Such structures integrate the actuator into the functional surface, which result in a small and lightweight component. Highly dynamic multi-functional movables on the trailing edge enable active load alleviation and, therefore, the design of lighter wings. Additional efficiency improvements arise by continuously adapting the wing geometry to changing various on-board pressure generation methods for the application of soft robotics [28]. Although, these pressure generation methods do not meet performance requirements, such as maximum pressure and continuous flow rate, which are essential for the intended aeronautical application. Vasista et al. evaluate different pneumatic and hydraulic system architectures for a fluid-driven morphing winglet trailing edge at an early conceptual stage [29]. However, none of the pressure-driven concepts has, so far, been investigated with regard to the flow inside the structure and transient dynamic behavior.
The overall aim of this paper is the multidisciplinary design of a shape-variable PACSdriven aileron, with particular focus on the transient flow within the cellular actuator. This leads to the following three main research objectives: (1) determination of the structural, dynamic, and systemic demands on the shape-variable aileron and associated actuator design requirements for a HAR wing with limited installation space; (2) investigation of the dynamic performance of the actuator limited by the internal flow during the fast filling process; and, (3) design of a wind tunnel experiment to study the structural resistance against transient external loads. First, the reduced-order PACS design methodology is described (Section 2.1) and then applied to the example of a morphing aileron (Section 2.2). The authors derive a two-dimensional geometry for a shape-variable aileron, which is structurally integrated into a NACA 0012 airfoil starting behind the wing rear spar. A reduced-order model that is based on the Pseudo Bond Graph methodology is developed to realize an efficient simulation of the internal flow (Section 2.3). Subsequenty, the results from the transient simulations on the pressurization process are presented and verified with state of the art computational fluid dynamics (CFD) simulations, and the effects of duct geometry and pressurization speed are investigated (Section 3.1). Finally, the influence of external pressure loads that are caused by the flow around the airfoil is evaluated for the entire operating envelope and the reduced-order results are compared to finite element analysis (FEA) results (Section 3.2). The article concludes with a discussion of the main findings and an outlook on realizing a wind tunnel test with a PACS-driven aileron as well as on further extending the reduced-order simulation models (Section 4).

PACS-Concept and Approach
This section gives a brief summary about the working principles of PACS and the underlying design methodology. The concept of PACS is based on a biomimetic approach that mimics the nastic movement of plants. The Venus flytrap (Dionaea muscipula) or the sensitive plant (Mimosa pudica), for example, are able to quickly and powerfully fold their leaves by varying the cells' turgor pressure. This change in turgor pressure in combination with flexible cell walls allows for the plant to move without the presence of muscles or rigid mechanisms. The transfer of this working principle into a technical system results in a cellular structure that is composed of polygonal cells with rigid cell walls and compliant hinge areas. An increase in cell pressure causes a deformation of each cell into a state of minimal internal energy, which occurs when all joints lie on a circular arc and the crosssection becomes maximum. Any further increase in pressure no longer leads to a change in geometry, but to a stiffening of the structure. The maximum angular deflection per cell is solely determined by the ratio of the cell wall lengths; the total deformation of the cellular structure is then described by the sum of the angular deflections of each individual cell.
Using an iterative form-finding algorithm allows for the design of arbitrary polygonal cell structures that can move between two predefined target shapes. A reduced-order PACS model is developed based on a truss structure that consists of rigid cell walls and discrete joints. The compensatory stiffness of each joint is derived from a bending beam with the flexural stiffness EI, the length l, and a constant thickness h. The approach of virtual work is then utilized to calculate the pressure-dependent state of equilibrium for the polygonal cell structure. For given external loads and maximum cell pressure, the length of each cell wall is iteratively adjusted until the deformed shape at equilibrium matches the predefined target shapes. Subsequently, the truss structure is translated into a cross-sectional design when considering flexure hinges and cell wall geometries. A comprehensive description of the holistic design approach of PACS, which is applied and extended in this work, is presented in [26].
When compared to previous publications, an enhanced geometric model is used in the work presented in this paper. For the first time, a tetragonal cell geometry is implemented for the outer cell row, whereas all previous realizations of PACS utilized a pentagonal geometry for these cells. The advantages of the tetragonal cell geometry are investigated in a preliminary study and they are pointed out by Gramüller et al. [30]: a significant increase in deformation capacity; a smooth and gapless outer surface for improved aerodynamic efficiency; lower deformation deviations between the reduced-order truss model and FEA results; and, local reduction of stress peaks due to the parallel alignment of the upper hinge and the resultant force vector.

Design of a Morphing PACS-Driven Aileron
In the following, the design of a shape-variable multi-functional aileron is described by applying the PACS design process that is presented in Section 2.1. The chosen application of a morphing aileron makes high structural and systemic demands on the actuator concept (high dynamic loads, high actuation speed, and fail-safe design). This application is intended to assess the conceptual limits of PACS with regard to the dynamic behavior under aeronautical operating conditions.

Operating Conditions and Design Constraints
The shape-variable PACS-driven aileron is designed in the context of a future wind tunnel test. The wind tunnel test will be conducted in a low speed closed-circuit wind tunnel (Göttingen-type) at the Aeronautics Research Centre Niedersachsen at the Technische Universität Braunschweig. The wind tunnel achieves a maximum speed of v ∞ = 60 m/s in the closed test section with the size of 1300 mm × 1300 mm × 3000 mm. Figure 1a shows the schematic integration of a two-dimensional wing model into the tunnel test section. A well-researched symmetrical NACA 0012 airfoil with 1000 mm chord length is chosen to investigate the principle effects of a PACS under the influence of external aerodynamic loads. The wing model with a total width of 1300 mm has two integrated trailing edge PACS actuators, each 500 mm wide. Figure 1b shows a cutout of one actuator section. The model consists of four basic structural elements: (1) the main wing, including leading edge, wing box, and the mounting interface to the wind tunnel walls; (2) the PACS actuator; (3) a rigid trailing edge; and, (4) an elastic skin at the lower surface. Performance-driven requirements from aileron operation are the maximum deflection angle ∆β max and the deflection rateβ. Vechtel et al. examined a modified Airbus A350 system architecture and identified a deflection rate of 16 • /s to fulfill requirements of the EASA certification specification CS-25 for roll control [31]. Multi-functional control capabilities demand faster actuation. Xu et al. investigated a gust load alleviation system for commercial aircraft, with rate limits of 30 • /s [32]. Gern et al. state that control surface deflection rates of 90 • /s are typically for high-performance fighter aircraft [33]. In general, few data on the deflection angles and rates for the latest generation of commercial aircraft are available in published literature. In this paper, the maximum aileron deflection angle is defined to ∆β max = ±21°. The deflection rate is assumed to beβ = 21 • /s for a baseline configuration andβ = 42 • /s for a faster actuation scenario. Higher actuation rates are beneficial for gust load alleviation and flutter suppression, and they are to be investigated in future research.

PACS Topology and Initial Shape
The topology and dimensions of the PACS are determined by the available installation space that is prescribed by the cruise shape of the wing. The PACS is attached to the wing trailing edge. The rear spar is placed at 60% chord and no penetration of the wing box is allowed. Taking the mounting interface between actuator and rear spar into account, the actuated surface extents from 65% to 85% chord and the trailing edge behind 85% chord is rigid. The initial geometry is determined for the cruise configuration with zero aileron deflection, so that the outer surface follows the NACA 0012 contour. A suitable PACS topology with two cell rows results in n row 1 = 7 and n row 2 = 8 cells in upper and lower cell row respectively. A cell size of 32 mm is used for the first cell. Adjacent cells reduce in size with a tapering of 0.75, which describes the length ratio between first and last cell of the upper cell row. Figure 2 shows the resulting cross-sectional design and the integration into the NACA 0012 airfoil. The figure shows the final PACS geometry that is obtained from the optimization process described in the following sections. The demonstrated ability to scale the PACS into almost arbitrary design spaces enables the implementation of PACS as a control surface in HAR wings with limited installation space. Moreover, Figure 2 shows a smooth and gapless outer surface that results from the tetragonal cell design, which provides significant aerodynamic improvements when compared to any previously accomplished PACS design [30].

Target Shapes
The target shapes of the PACS actuator are derived based on the aileron operating conditions (∆β max = ±21°, see Section 2.2.1). The aileron deflection angle is assumed to be equal to the angular deflection at the tip of the PACS actuator, using the right mounting frame at 85% chord as a reference point for calculating the angular deflection. A pressurization of the upper cell row causes a downward deflection and a pressurization of the lower cell row an upward deflection. Figure 3a shows the deformed PACS cross section for both of the target shapes. The target shapes are obtained by superimposing a circular arc over the undeformed cruise shape. Under maximum cell pressure, each cell of the upper cell row deflects by the same angle ∆α. Combing the deformations of all n row 1 cells results in the total angular deflection of the PACS actuator. Subsequently, the angular deflection per cell used for the PACS design is ∆α = ∆β max /n row 1 = ±3°. ∆β ∆α

Loads and Boundary Conditions
The loads acting on the structure can be divided into internal and external loads. The primary source of internal loads results from the cell pressurization. The internal pressure acts on all inner surfaces and causes a deflection of the flexure hinges and, consequently, a global deformation of the structure. Simultaneously, a counteracting bending force occurs due to the inherent stiffness of the hinges. The maximum hinge deflection is mainly defined by the geometric dimensions of the cell walls and it will be further discussed in Section 2.2.6. The maximum cell pressure is set to p row 1,max = p row 2,max = 1 MPa, determined by the applied material (Section 2.2.5), and it is verified by FEA (Section 2.2.7).
External loads can act in the same or opposite direction as the deformation due to the internal pressurization of the PACS. In the presented application, the structure is exposed to aerodynamic loads that result from the external flow field, which counteract the deformation of the PACS actuator. Based on the wind tunnel test parameters, a free stream velocity of v ∞ = 60 m/s, an ambient pressure of p ∞ = 0.1013 MPa, and wind tunnel temperature of T ∞ = 40°C are defined. The pressure distribution around the airfoil is calculated with the subsonic analysis tool XFOIL that is based on a 2D panel method. The input parameters for XFOIL are a Mach number of Ma = 0.172 and a Reynolds number of Re = 3.754 × 10 6 that result for an ideal gas with the values of v ∞ , p ∞ and T ∞ shown above. The angle of attack is set to AOA = 0°to evaluate the performance of the PACS actuator while neglecting the influence of the specific airfoil characteristics. Geometric input data are the NACA 0012 airfoil coordinates that are shown in Figure 2 superimposed with the angular deflection ∆β acting between 65% to 85% chord length. Figure 3b shows the pressure distribution that is given by XFOIL for both target shapes. In the PACS design process, the aerodynamic pressure is applied as surface loads to the outer walls of the upper cell row. Additional aerodynamic loads act on the rigid trailing edge behind 85% chord and on the elastic skin. The bending moment that results from the accumulation of these pressure loads is applied as a discrete pair of forces at the right end of the PACS actuator.
A mounting interface is integrated into both edges of the PACS actuator. The connection to the rear spar of the wing box is made via a fixed bearing at the left end of the PACS. A rigid trailing edge is connected to the actuator on the right side.

Material Selection
Gramüller et al. show that the ratio of squared strength to stiffness R 2 /E provides a decisive criterion for the material selection and, thus, for the entire design [25]. Consequently, flexible but high-strength materials are demanded. Fiber reinforced plastics (FRP) offer the best ratio of R 2 /E, which makes them the most suitable materials for application in PACS. In this study, a glass-fiber reinforced plastic (GFRP)-HexPly913/28%/192/EC9756-is utilized, which has a Young's modulus of E = 42,000 MPa, a tensile strength of R = 1200 MPa, and a Poisson's ratio of ν = 0.29. This material is selected inside the group of FRP based on three reasons: (1) glass-fibers have a relatively high elongation at break and they are considerably less sensitive to compression than carbon or aramid fibers. Both is beneficial for the application as highly deflected flexure hinges; (2) the development of an automated manufacturing process for PACS made out of woven GFRP is part of a current research project [34]; and, (3) the usability of HexPly913 in the context of PACS has already been proven during the realization of the first functional demonstrator [25].

Hinge and Cell Wall Sizing
The working principle of PACS is based on the combination of rigid cell walls with high thickness and flexible areas with reduced thickness. Since the entire cell deformation is concentrated on the latter, a key challenge is the design of these flexure hinges. The hinge design is driven by material properties and performance requirements. Each flexure hinge is described by its thickness h, length l, hinge angle u, and eccentricity ζ (see Figure 2). The hinge elements are modeled as bending beams with the deflection ∆u, whereby ∆u is the change in the hinge angle u due to the deformation of the pressurized cell. It is shown in previous publications that an optimum hinge thickness h = h opt exists that allows for applying the maximum pressure to a PACS cell [25]. By superimposing normal stresses that result from the internal pressure p and bending stresses caused by the hinge deflection ∆u, the following equation for the optimum hinge thickness can be derived: For the PACS actuator that is shown in Figure 2, the maximum hinge deflection appears in the bottom hinge in the lower cell row and it amounts to ∆u = 18.9°. Thus, using the material properties from Section 2.2.5 and a hinge length of l = 3 mm, the optimal hinge thickness for the critical hinge results in h opt = 0.26 mm. This constant thickness is applied to all hinges.
The dimensioning of the cell wall thickness is mainly done to fulfill the assumption of rigidity and allow for the axial sealing of the cells (the latter is not presented in this paper). A thickness ratio of about 1:12 between rigid cell walls and flexible hinge areas is chosen, which results in a wall thickness of h w = 3 mm.
The determination of the cell wall lengths is achieved while using an iterative optimization process. The deformability of each cell is geometrically defined by the ratio of its wall lengths. Under maximum internal pressure, the accumulated deflections of all cells must match the target shape. The cell wall lengths are determined applying the form-finding algorithm that is described in Section 2.1. During form-finding, the lengths are altered until a state of equilibrium is reached between internal and external loads. Table 1 presents a summary of all the design parameters used in the form-finding process. The resulting PACS geometry corresponds to the design that is shown in Figure 2. The pressure-dependent deformations derived from the reduced-order PACS design tool are numerically verified with FEA. The simulation software ANSYS is utilized to model the 2D solid structure ( Figure 2) with plane182-elements applying plain strain behavior. A MATLAB interface directly imports loads and boundary conditions for each load case from the PACS design tool into ANSYS (see Section 2.2.4). A mesh convergence study revealed that an element length that is equal to the hinge thickness h opt is sufficient to correctly predict the total angular deflection ∆β, which changes by less than 1% with increasing mesh density. However, previous analyses showed that a mesh refinement to at least six elements in the thickness direction of each flexure hinge is required to obtain correct stress values. Therefore, the element length is set to h opt /6 = 0.043 mm in the hinge regions, and it is increased in all other areas by a factor of ten. This results in a total number of 61,719 elements. Pressure loads are applied as a structural surface effect via surf153-elements, and the point loads are distributed to the right mounting frame by using MPC style contacts with targe169/conta172-elements.

Pressurization and Internal Flow Behavior of a PACS
The PACS actuator is a subsystem within the pneumatic system architecture of the wing. Further subsystems are the compressor as a compressed air generating component, tubes, or pipelines for distributing the compressed air and storage tanks for intermediate storage of the fluid during continuous pressurization and depressurization. The PACS actuator has one pressure port per cell row. Internal ducts in the cell walls realize the fluid transfer between the cells of each row (Figure 4a). For reasons of accessibility, the pressure ports are positioned at the mounting frame. Between adjacent cell, there are 12 internal ducts with a cross-sectional area of A Duct = 29.91 mm 2 each, which are equally distributed over the total depth of 500 mm. For a better understanding of the transient pressurization process, detailed studies of the flow within the PACS are performed using CFD. The results of the three-dimensional transient CFD simulations are used to verify a reduced-order model (ROM) based on Pseudo Bond Graph notation that is developed within this study. The derived ROM is used within the multidisciplinary PACS design approach to optimize the thermo-fluid dynamic components and gain a deeper understanding of the interaction between the transient internal flow and associated structural deformation.

Numerical Simulation of the Flow inside a PACS Cell Row Using CFD
All the simulations of the present investigations are done with the commercially available three-dimensional (U)RANS solver ANSYS FLUENT. Computational grids of increasing number of cells are generated using the included grid generator FLUENT Meshing. Preliminary simulations are performed with different turbulence models, along with a grid sensitivity analysis. The grid sensitivity analysis shows that at least 9.48 million polyhedrons and prisms are required to achieve (i) a solid resolution of the boundary layers (y + ≈ 3 − 10) and (ii) the overall flow parameters accurately. The standard k-ω turbulence model turns out to be best suited for the present investigation. The computational domain that is regarded here covers the internal volume of the lower cell row of the PACS, including the pressure port and fluid channels that connect the single cells (Figure 4b). Transient simulations are calculated with a time step of 5 × 10 −5 s.
The pressurization is conducted as a fast filling process to simulate almost real operating conditions. Two pressurization scenarios are defined for increasing the cell pressure by p max = 1 MPa in order to deflect the aileron by ∆β max = ±21°: (i) full deflection in t = 1 s (corresponds toβ = 21 • /s), and (ii) full deflection in t = 0.5 s (corresponds tȯ β = 42 • /s). The filling process is realized using a linear pressure ramp of the total pressure from p t,start = 0.1 MPa to p t,end = 1.1 MPa at the inlet of the computational domain. Furthermore, a constant total temperature of T t,start = T t,end = 288.15 K and low turbulence intensity (Tu = 0.5%) are specified at the inlet, and no reverse flow is assumed.

PACS-Driven Aileron as a Reduced-Order Model by Pseudo Bond Graph Approach
State of the art high fidelity and high resolution tools are available for most of the physical domains that are involved in PACS. These include FEA for structural and thermal calculations as well as CFD methods for flow simulation inside the structure. While offering a high degree of detail, they require extensive sets of boundary conditions as well as computation time. Simulations of numerous PACS-setups during the design process are not currently feasible and there is a need for a reduced-order model. The Bond Graph theory features less level of detail but a higher degree of flexibility than the previously mentioned high-fidelity tools. The simulation tool ASTOR (Aircraft Engine Simulation for Transient Operation Research) that was developed at Technische Universität Braunschweig is based on the Pseudo Bond Graph theory, which is an extension of the Bond Graph theory. It can be used to develop reduced-order models for multiphysical systems, such as jet engines [35], electric air compressors [36], or other fluid dynamic systems [37]. An application towards PACS is also feasible, as the Pseudo Bond Graph theory on the one hand provides the ability to derive sketches of complex multiphysical systems with the advantage of giving clear instructions on the interconnections of the miscellaneous components. On the other hand, it enables the easy conversion of power flows in various physical domains and the derivation of a system of ordinary differential equations. ASTOR has been extensively described in the field of turbomachinery. The tool is used in this study to investigate transient phenomena in the field of PACS by deriving a ROM for the given structure.
In the following, the elements used in the Bond Graph theory are explained in more detail. Figure 5 shows a schematic drawing of a gas path element used in PACS. The product of effort e (e.g., a temperature T) and flow f (e.g., volume flowV) of a single bond have the physical dimension of a power, as per the Bond Graph formalism. However, this cannot be achieved easily in compressible flows. Therefore, the Pseudo Bond Graph notation introduces a second pair of effort and flow. With this extension, the Pseudo Bond Graph that is depicted at the bottom right of Figure 5 can fulfill the fundamental requirement of the Bond Graph theory. The upper bond provides mass flowṁ and pressure p, while the lower bond provides temperature T and energy flowĖ.
Subsequently, the simplified equations of motions, which are solved for the flow of process gas within the structure, are being described in more detail [37]. At the upper 1-junction, the conservation of momentum (Equation (2)) is enforced by using the inertia storage I.
At the following 0-junction, the balance of efforts is calculated, as shown in Equation (3) using the conservation of mass. The change in density can be obtained by balancing the incoming and outgoing mass flow. By incorporating the capacitive storage volume C, the volume dynamic effects can be modeled.
Finally, temperature changes are obtained at the lower 0-junction using the conservation of energy (Equation (4)) and the thermal capacity C T of the process gas that is currently stored within the components volume. By adding a heat flowQ from the inside of the component to the environment or vice versa, free or forced convective cooling can be simulated, but it is neglected in this study.
The temporal derivatives from Equations (2)-(4) are integrated and a temporal resolution of the quantities under investigation can be obtained. Currently, no fluid-structure interaction is modeled in ASTOR and the aileron's inertia is disregarded. Nevertheless, the coupling of internal pressure forces, structure deformation, and external aerodynamic loads will be included in future studies. This includes accurate aerodynamic loads (e.g., from XFOIL), together with the overall cellular structures represented as a truss model. The interfaces of the miscellaneous components and their interconnection are presented in, e.g., Karnopp [38].

Studies on the Transient Internal Flow
The following results describe the transient system behavior of the case of maximum aileron deflection for one PACS cell row. Only the lower cell row is pressurized and the upper cell row is kept at an ambient pressure. Fluid-structure interaction is neglected here. A goal of the work presented in this paper is to create a ROM for representing the threedimensional time-dependent flow within the cell rows of the PACS that can be integrated into the holistic design process. First, the three-dimensional time-dependent behavior of a fast filling process of the lower cell row is investigated while using CFD. Subsequently, the ability of a ROM that is based on the Pseudo Bond Graph approach for representing the transient system behavior of the PACS is described. ROM results are shown for a 2D section and compared to the three-dimensional time-dependent CFD results.

Flow Behavior Through Ducts and Cells
CFD simulations of the PACS are performed for a pressurization process corresponding to a deflection rate ofβ = 21 • /s, as described in Section 2.3. The fast filling process is prescribed here as a linear pressure ramp from p t,start to p t,end within t = 1 s (cf. Figure 8a). The time-dependent three-dimensional flow behavior for the x-y and x-z plane can be seen in Figures 6 and 7, respectively. The normalized static pressure, including streamlines, is shown in Figure 6 as contour plots at mid-span of the lower cell row for three different time steps. It can be seen that the fast filling process within the first time steps t = 0.001 s to 0.005 s results in an uneven distribution of the normalized static pressure throughout the single cells. The pressurization needs some time to proceed through the PACS, resulting in a delayed deformation of the cells and, hence, deflection of the cell row. The delay in pressure within the last cells is overcome after a certain time step of t = 0.010 s. Here, the velocity direction of the main flow indicated by the streamlines also shows the end of a homogenization process resulting in recirculating zones in each cell, where the dynamic pressure is transferred in static pressure. Figure 7 shows the normalized static pressure contour plot in a x-z view through the lower cell row for three different time steps. It can be seen that the fast filling process, initializing mass flow at the inlet duct, as seen in the top left of the figures, results in a second uneven pressurization of the cell rows. Within the first time steps t = 0.001 s to 0.005 s, not much fluid is transported through the whole PACS system. The transient pressurization results in a diagonal filling, which can result in an uneven deflection of the HAR wing aileron.

Variation of Pressurization Time
Two scenarios with different pressurization times are analyzed. The baseline scenario (Section 3.1.1) is compared to a scenario with doubled actuation speed (β = 42 • /s), i.e., half pressurization time (t = 0.5 s). The transient flow behavior turned out to be comparable for both of the pressurization scenarios. For faster filling speeds, locally increased flow velocities up to critical flow conditions inside the internal ducts lead to increased total pressure loss. The faster filling leads to higher temperatures, hence having an impact on material selection. Overall, the results show a potential for improving the internal flow by finding an optimal duct geometry (see Section 3.1.4).

Pseudo Bond Graph Methodology vs. CFD
The three-dimensional time-dependent flow behavior is represented by a ROM using the Pseudo Bond Graph methodology. Figure 8 compares the filling process that was calculated with the spatially and time resolved CFD and the time resolved ROM simulation. The results have been normalized with the steady-state value at the end of the simulation duration. The overall performance of the filling process presented in Figure 8a is similar throughout the majority of the simulation time. However, the predicted behavior during the beginning (Figure 8b) and end (Figure 8c) of the filling process deviates significantly. On the one hand, the responsiveness of the ROM is much better when compared to CFD. On the other hand, a series of over-and undershoots occurs at the end of the simulation. These phenomena are both due to the reduced damping of the ROM. The highly threedimensional flow characteristics with large vortices, as visible in Figure 7, are not accurately captured in the ROM. This accounts for large deviations at the start of the filling process as well as at the end. Furthermore, the ROM approach, which only considers 1/12th of the 3D-system as a simplification, cannot capture the highly three-dimensional character of the flow propagation visible in Figure 7a. Therefore, further efforts are required in the future to assess how to correct the single-passage ROM for 3D filling phenomena.

Comparison of Duct Cross-Section Geometries
A study on the influence of cross-sectional area of the connecting ducts is performed using the ROM. With short simulation run-times, this type of model is of good use for parameter studies. Therefore, the area of the connecting ducts has been doubled and halved in order to study the impact on the pressurization performance ( Figure 9). The duct geometry significantly influences the performance of the pressure build up, both during the beginning and the end of the filling process. At the beginning, an increased area leads to almost immediate response of the pressure inside the cells with very small deviations between the first and the last cell. However, with reduced area, the initial delay is approximately doubled and the overall pressure gradient is smaller. At the end of the filling process, an increased duct area leads to less oscillation in the last cell, while a reduced area increases them significantly. Because this study has been conducted without the feedback loop from a structural perspective, evaluations of the feasibility regarding structural requirements will have to be conducted in the future.  Figure 10a shows the deformation of the PACS for a pressurization of only one cell row at a time. The figure demonstrates that the reduced-order truss model exactly predicts both of the design points. A comparison between the reduced-order truss model and the high-fidelity finite element model shows a difference in deflection of 12.8% for the first design point and 14.7% for the second design point. A smaller deflection is shown than in the FEA, as the reduced-order truss model neglects the deformation of the cell walls and does not consider the hinge transition zones.
Furthermore, the resistance of the structure against external aerodynamic loads is evaluated. Three load scenarios are considered. The baseline scenario with design loads from XFOIL is compared to a scenario with loads increased by a factor of 1.5 and to a scenario with zero aerodynamic loads. For each deformation state, the subsonic analysis tool XFOIL iteratively recalculates the aerodynamic pressure distribution. Neglecting the external loads leads to an increased deformation (dashed lines), since the aerodynamic loads counteract the deflection. Increasing the external loads by a factor of 1.5 reduces the deflection (dotted lines). The difference in deflection due to the 50% load increase amounts to 3.5%, respectively, 10.6% for both of the design points. The influence of the external loads is less pronounced at higher cell pressures, because higher pressures lead to a stiffening of the structure. Figure 10b presents the angular deflection ∆β covering the entire operating envelope of the PACS. The surface plot is obtained for different pressure combinations in both cell rows from a total of 625 data points per calculation model. The figure indicates a deformation neutral line formed by the intersection with a plane through ∆β = 0°. The cruise configuration is specified by an aileron deflection of 0°and, thus, can be achieved by any pressure combination on this line. A shift of this deformation neutral line towards small cell pressures in the lower cell row can be justified by a larger lever arm as compared to the bending neutral axis. In order to achieve a high resistance to external loads and, at the same time, allow for a fast deflection into both target shapes, the cruise design point is set to p row 1 = 0.8 MPa and p row 2 = 0.35 MPa.  Figure 10c shows an evaluation of the structural stresses for a pressurization of only the lower cell row, since this is the critical load condition. All other pressure combinations cause lower stresses, because pressurizing both cell rows counteracts the deformation, and less bending stresses are experienced by the flexural hinges. For the same reasons, an increase in aerodynamic loads results in lower structural stress. The maximum equivalent stress σ max that is caused by stress concentrations in the hinge transition regions is evaluated by means of FEA and normalized to the material strength R. A safety factor of 1.14 is achieved at the design point at p row 1 = 0 MPa and p row 2 = 1 MPa. Figure 10d presents the angular deflection ∆β over time for the pressurization of the lower cell row with the upper row remaining at ambient conditions. In addition, the linear pressure rise at the inlet is plotted in black. A highly nonlinear deflection characteristic of the structure can be observed, even though the pressure rise is almost linear according to Figure 8a. Deviations between the reduced-order model and CFD simulation are limited to the initial deformation, where the ROM predicts a higher deflection due to the increased responsiveness. Howeber, at the end of the maneuver, discrepancies diminish due to the almost identical pressure rise as well as the stiffening characteristics of the PACS cells at near maximum pressure. The overshoots that are predicted by the ROM do not appear to have great influence on the deflection, as they are in the stiffening zone of the PACS, where angular deflections due to pressure changes diminish.

Discussion and Conclusions
The objective of this work is developing a multidisciplinary design methodology for PACS-driven morphing ailerons, with a focus on the transient dynamic behavior of the cellular structure. A challenge with HAR wings is the small installation space, especially near the wingtip or in the FWT. In the scope of this paper, the almost arbitrary scalability of PACS to the limited installation space at a wing trailing edge is demonstrated. Decisive design aspects for an overall PACS system are pointed out in the context of a future wind tunnel test. The resistance of the structure to changing external loads on the one hand and the dynamic performance of the flight control surface on the other are the main research questions. Two reduced-order models are developed to assess both questions: one for the structural design and one for the simulation of the flow inside the structure.
A characterization of the PACS deflection under varying external loads shows a high structural stiffness, which is more pronounced with increasing internal pressure. Using the PACS design tool that is based on a reduced-order truss model, feasible deformations are calculated for the entire operating envelope. For each deformation state, the subsonic analysis tool XFOIL iteratively recalculates the aerodynamic loads. A load increase of 50% results in a change of deflection angle of 3.5% to 10.6% for the two target shapes. For the cruise configuration, a pressure state of p row 1 = 0.8 MPa and p row 2 = 0.35 MPa is identified, as it provides high structural stiffness and simultaneously allows for rapid actuation into both target shapes.
Numerical verification utilizing FEA shows that the reduced-order truss model predicts 12.8% to 14.7% lower deflection. The deviations result from the model assumptions of rigid cell walls and the neglect of hinge transition zones. Comparing the results presented here with the variable-camber PACS wing of Gramüller et al. shows an increased model accuracy resulting from the improved hinge alignment in the tetragonal cell design. The relative deviation in the publication of Gramüller et al. ranges from 15.1% to 79.9% between the reduced-order truss model and FEA for both target shapes [26]. A similar improvement due to the tetragonal cell design over the pentagonal one is also shown in [30]. An extended calculation approach that is presented by Pagitz et al. further increases the prediction accuracy by taking the cell wall compliance into account [39]. The subject of the current work is the design of flexure hinges considering manufacturing constraints that are specific to automated textile fabrication [34]. Based on these findings, an improved modeling approach can be implemented to account for the hinge transition geometry.
The presented results from the flow simulation inside the pressure-driven structure allow for the assessment of its dynamic performance. A multi-functional aileron with the capability of active load alleviation and flutter suppression places high demands on actuating speeds. The response time of the PACS actuator is mainly limited by the transient pressurization process, which is investigated with state of the art CFD simulations.
Calculations of the most unfavorable operating points show that a deflection rate oḟ β = 42 • /s is achievable in the entire operating range. In all cases, inflation from ambient pressure to the maximum cell pressure is considered and, thus, an actuator deflection from the neutral position to the maximum deflection angle of ∆β = 21°. However, during flight operations, there is always a certain internal pressure to maintain structural stiffness (cf. cruise configuration). In addition, full actuator deflection is not necessary for the purpose of gust load alleviation and deflections angles are therefore limited to ±10°in this case [32]. With these adjustments and while taking improvements in the flow channel geometry into account, the authors consider the operation of PACS as a highly dynamic gust load alleviation system as feasible. Further investigations on higher deflection rates for load alleviation and flutter suppression are part of the current work.
The results from CFD calculations point out the transient behavior of the pressurization, leading to a transient, non-uniform aileron deflection. Local critical flow conditions can be observed in the inlet duct between the piping and the first cell for the fastest filling process investigated. Because the study presented here has so far neglected the fluid-structure interaction, the influence of dynamic pressure from the airflow around the model, and the connection to the upper cell row, the results should be verified using an enhanced CFD model. Dynamic effects from structural damping and inertia are also to be investigated in the future.
The authors introduce a reduced-order model that is based on the Pseudo Bond Graph methodology enabling efficient flow simulation. Thus, the application of the Pseudo Bond Graph method allows the calculation of a large number of conceptual variants for the optimization of the internal flow using minimal computational resources. In future research, the flow simulation will be extended to include fluid-structure interaction. The Bond Graph methodology will provide a good basis, as it clearly specifies the interconnections of the fluid and structural domains via corresponding 1-and 0-junctions. Thus, the consideration of variable volumes in the Pseudo Bond Graph model with corresponding models of structural stiffness, like beams representing the PACS structure, are going to be incorporated. A more detailed discretization in the third dimension allows for representing the diagonal flow component during the filling process. In addition, losses due to three-dimensional flow effects encountered in the CFD simulations have to be modeled. Finally, the enhanced model accounting for fluid-structure interaction has to be extended by external aerodynamic forces. These forces can either be generated by XFOIL for steady-state loads or by unsteady CFD for time-varying loads, like gusts.
The realization of the wind tunnel model requires the design of further components. These include the elastic skin at the lower surface and the main wing, together with sensor technology, like static pressure taps and force measurement. A review on approaches for elastic skins can be found in [6]. Moreover, the axial cell closure is to be designed for the given PACS structure. A method for determining the shape of a deformation supportive end cap for PACS is described in [27]. Future research includes the development of a holistic system architecture that considers the pressure supply in the overall aircraft system. Next to performance-related studies, the backlash between PACS and the pressure supply system is investigated by extending the Pseudo Bond Graph method to include these components.

Acknowledgments:
The authors thank Jan Göing for sharing his insight into the Pseudo Bond Graph theory and Hendrik Traub for helpful discussions on the topic of PACS, which greatly improved this manuscript.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations
The following abbreviations are used in this manuscript: