Aeroelastic simulations of a delta wing with a Chimera approach for deflected control surfaces

A numerical tool for the computation of aircraft control surface aerodynamics with flexibility effects is presented. The solution is based on coupled Computational Fluid Dynamics (CFD) and Computational Structural Mechanics (CSM) simulations embedded in the multidisciplinary simulation environment SimServer. In SimServer, the DLR-TAU Code is utilized to obtain the CFD solution by solving the Reynolds-Averaged Navier–Stokes (RANS) equations. Structural displacements are computed with a modal solver. The Chimera implementation of SimServer, suited for hybrid grids, is applied to model the control surfaces. Numerical simulations with the flexible Chimera method are performed for the Model53 wing configuration, which is a generic delta wing with a deployed slat as well as an inboard and outboard trailing edge flap. Aerodynamic and aeroelastic simulations at high dynamic pressure q=45\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q=45$$\end{document} kPa and transonic speed Ma=0.8\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text {Ma}} = 0.8$$\end{document} are performed for several angles of attack 10∘≤α≤25∘\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10^\circ \le \alpha \le 25^\circ$$\end{document} and flap deflection angles -30∘≤δ≤30∘\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-30^\circ \le \delta \le 30^\circ$$\end{document}. The effect of structural deformations on the flow field and control surface effectiveness are analyzed and compared to computations of components treated fully rigid. At the targeted freestream condition M=0.8\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M=0.8$$\end{document} and Re=15.1×107\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text {Re}}=15.1 \times 10^7$$\end{document}, the flow field around the Model53 configuration is characterized by the interaction of vortices and shock waves. The results of the lift and pitching moment coefficient for the rigid and flexible configuration revealed the importance of taking the structural flexibility into account in order to obtain more accurate results for the considered range of flap deflections. Furthermore, the computational effort of the aerodynamic and aeroelastic simulations are evaluated. The increase in computational effort is shown to be adequate for the given increase in accuracy.


Introduction
The influence of control surface deflections on the aerodynamic characteristics of an aircraft is of high interest in aerospace engineering. In particular, highly loaded flaps, e.g. for high agility aircraft [1][2][3], may cause large structural deformations, thus degrading the aerodynamic efficiency due to changes in spanwise flap gap adjustments, flap positions and wing twist. Aircraft maneuvers or trim state simulations with Computational Fluid Dynamics (CFD) solvers require the treatment of control surface deflections [4]. The movement of control surfaces in high-fidelity simulations is a challenging task. Re-meshing for each deflection angle is computationally expensive and not trivial for time-dependent motions [5]. Mesh deformation methods are problematic in combination with large deployments and do not allow for the opening of gaps between control surfaces and the wing. The Chimera technique, however, offers the possibility to model the deployment of control surfaces realistically, while enabling movements during a simulation without the necessity of re-meshing the computational grid. Thus, the efficiency of computing a series of deflection angles or an aerodynamic dataset for an aircraft configuration can be significantly improved with the Chimera method. To further increase the accuracy of such computations, structural stiffness needs to be considered due to their effect on control surface efficiency and the aerodynamic characteristics of the aircraft in general.
Over the last years, the modeling of control surface deflections with CFD solvers has been subject to several 1 3 projects, such as ComFlite [6], Digital-X [7], VicToria and Mephisto [8], as well as several publications, for instance by Vos et al. [9], Sørensen et al. [10] and Clifton et al. [11]. Similarly, aeroelastic simulations of aircraft using coupled iterative CFD-CSM simulations with high-fidelity CFD methods (RANS or higher precision) and detailed structural models have been investigated by several authors such as Voß et al. [12], Fornasier et al. [13] and Stickan et al. [14]. However, the combination of coupled high-fidelity CFD-CSM simulations with an accurate and efficient modeling of control surfaces is still an ongoing research topic. The focus of this work is the introduction and application of a numerical method based on the multidisciplinary simulation environment SimServer [15], which is capable of conducting coupled CFD-CSM simulation with control surface deflections. The modeling of control surfaces is performed with a Chimera approach for hybrid structured/unstructured grids, which ensures a robust, accurate and efficient handling of control surface deflections for aircraft configurations up to industrial complexity [10]. Due to the compatibility of the Chimera implementation with hybrid grids as well as intersecting surfaces, it offers the possibility to handle complex geometries as well as different types of control surfaces. Hence, the presented solution is considered to be advantageous compared to methods which are using sliding mesh approaches as described in [5], since the Chimera technique is more robust combined with mesh deformation, or to solutions which are limited to structured grids only, as presented in [9], where the mesh generation process can be time-consuming and difficult for complex geometries.
First, the numerical methods which are involved in the aeroelastic simulations are described. Subsequently, the Model53 wing geometry as well as the aerodynamic and structural grid is presented. Finally, the results obtained with the flexible Chimera method are analyzed for the Model53 at transonic speed and high dynamic pressure.

Numerical methods
This chapter describes the applied numerical methods. The CFD-solver, the fluid-structure coupling and the CFD-grid deformation is introduced. Thereafter, the CSM-solver, the Chimera implementation and the CFD-CSM coupling process are described. All methods are implemented into the SimServer multidisciplinary simulation environment [15]. Emphasis is put on the integration of the Chimera method into the partitioned approach of SimServer for solving aeroelastic problems.

Computational fluid dynamics solver
The DLR-TAU flow code is employed in the SimServer environment to solve the compressible, three-dimensional Reynolds-Averaged Navier Stokes (RANS) equations [16]. It uses a vertex-centered dual cell method for hybrid, unstructured grids. Even though steady-state cases are targeted, the aerodynamic and aeroelastic simulations are performed as URANS (Unsteady-RANS) computations with a dual-time stepping and an implicit Backward-Euler scheme. A time step of Δt = 0.01 s is used, which corresponds to a characteristic time step of Δt * = (ΔtU ∞ ∕l ) = 0.5 , where U ∞ is the freestream velocity and l the mean aerodynamic chord. A comparison of CFD simulations with different time steps 0.001 s < Δt < 0.05 s and a steady-state RANS simulation revealed no difference with regard to the final results. The time step of Δt = 0.01 s is chosen, as it provided stable and fast convergence. The inner iterations are set to i = 50 but can be reduced to i = 25 by the solver if the convergence criterion is achieved for the lift and pitching moment coefficients, which is |C L | < 0.005 and |C My | < 0.001 within the last 500 iterations. After reaching the convergence criterion, the specified time steps for each deflection angle are completed but continued with the reduced inner iterations, so that a reduction in computational effort can be achieved. The inner iterations are reset to i = 50 after each flap rotation. In addition, convergence acceleration is achieved with the multigrid technique, for which a 3w cycle is utilized, using point-implicit relaxation. For the spatial discretization, a finite-volume formulation with second-order upwind scheme is applied. The RANS equations are closed with the Spalart-Allmaras one-equation turbulence model [17] in the Edward's modified version with rotational correction [18].

Fluid-structure coupling
For the computation of weakly coupled aeroelastic simulations with movable control surfaces, an important requirement to the fluid-structure coupling is the compatibility with the Chimera implementation. The existing coupling algorithms in SimServer did not fulfill this requirement, especially as problems arose in combination with intersecting geometries, that can occur in combination with Chimera. Thus, a new coupling algorithm has been developed for transferring forces and displacements between the surface nodes of the CFD-grid and the surface nodes of the structural finite-element (FE) model. The interpolation method is based on the point-element-relationship between the CFD-and the CSM-surface grid and has been validated with the BACT [19] and HIRENASD [20] test case. The point-element relationships are computed by a nearest neighbour search, which is based on finding the nearest CSM element for each CFD node.
The point-element-relationships are initialized with the CFD and CSM grid in their undeflected flap configuration during the first coupling exchange cycle. To reduce the computational effort for subsequent coupling cycles, the point-element-relationships are stored in registers. For the transfer of forces and displacements, the deflection of control surfaces is only considered for the CFD grid, while the CSM flaps are undeflected. The forces are evaluated on the deflected flaps of the CFD grid and then transferred to the undeflected CSM grid by maintaining the direction of the force vectors, see Fig. 1. Similarly, the displacements are computed on the undeflected FE model and coupled to the deflected CFD grid by maintaining the direction of the displacement vectors.
Compatibility with the Chimera implementation of Sim-Server is ensured by allowing the user to define the surface IDs of the structural and aerodynamic grid where data should be exchanged. Hence, the transfer of forces and displacements in Chimera regions, e.g. flaps, can be performed independently of other regions, such as the main wing or the fuselage.
Another important aspect in the case where the rotation of flaps are targeted in aeroelastic simulations is the adjustment of the hinge line axis according to the deformations of the wing. To take this into account, the coordinates of the two points defining the hinge axis of each flap are displaced by employing the same algorithm used for the transfer of displacements from the CSM to the CFD grid.

CFD volume grid deformation
The surface displacements are propagated to the CFD volume grid with a parallel mesh deformation method based on an advancing front algorithm. Since the repair of badly distorted elements was not available in the original code, a local re-meshing procedure according to Ref. [21] has been implemented in order to support larger local deformations. The repair procedure marks regions in the grid containing collapsed elements surrounded by valid cells only. Subsequently, the collapsed cells are rebuilt by employing a volume spline interpolation approach. More information on the mesh deformation method can be found in Refs. [15,21].

Chimera implementation
The Chimera technique is an efficient way to compute multiple body problems, such as control surface deflections, for hybrid grids [22]. In general, it assembles two or more grids into one computational domain through the application of interpolation boundary conditions. To apply the Chimera method in SimServer, one or more submeshes need to be contained within a background mesh. In case of a wing configuration with control surfaces, the wing mesh corresponds to the background mesh, whereas each control surface mesh corresponds to a submesh. Geometries belonging to the background mesh or submesh are allowed to intersect. A detailed explanation of the interpolation method and the hole-cutting procedure can be found in Ref. [10].

Structural solver
The structural deformations for the static aeroelastic simulations in this work are computed with a linear modal solver. The linear approach is justifiable for the investigated configuration, since the overall elastic deformations remain small with a wing tip deformation of approx. 4% relative to the semi-span. A mass-normalized modal matrix s , which contains the structural eigenmodes obtained by an eigenvalue analysis, is used to write the physical structural displacements in modal coordinates by

CFD-CSM equilibrium
The static aeroelastic equilibrium is found by an iterative process shown in Fig. 2, which is based on the numerical methods introduced above. The CFD-CSM loop is repeated until the convergence criterion or the specified number of CFD-CSM cycles is reached for the initial flap deflection angle. If a series of flap deflections are specified, the flap Chimera blocks are rotated accordingly and the CFD-CSM loop continues. The loop terminates, when the equilibrium of the final flap deflection angle is computed. This coupling process enables the efficient (1) s = s .
(2) s = s , computation of the aeroelastic equilibrium for a series a flap deflection angles, as less iterations are needed for subsequent flap angles compared to the initial angle due to the better initial guess of the flow field and deformations. For the aeroelastic simulations in this work, four to five CFD-CSM exchange cycles are specified for each flap deflection angle.

Model5delta wing
The numerical analysis in this work is performed for the Model53 test case, which is a cropped delta wing with a leading edge sweep angle of le = 53

Finite element model
The FE model of the Model53, Fig. 5a, consists of 7928 grid points, 17,147 triangle elements and 24 bar elements. For the modeling of the ribs, spars, stiffening elements (shown in turquoise in Fig. 5a) and skin, shell elements are used. Stiffening elements are explicitly modeled and are employed to keep buckling fields sufficiently small as well as to reduce local eigenmodes. The inner and outer flap are structurally modeled with skin, ribs and spars and are attached to the wing with a hinge concept introduced in Ref. [24], which is shown in Fig. 5b for the Model53. Each attachment consists of two triangle elements, four bar elements and a double node connecting the bar and triangle elements of the wing and flap side with a rigid body element. The inner and outer flap are each connected to the main wing with three attachments, of which one attachment can transfer translational and rotational degrees of freedom and two attachments only translational degrees of freedom. Each triangle element of the flap attachments is either connected to a spar of the main wing or the flap. The stiffness of the attachment can be controlled with the bar's material properties. This hinge concept enables a realistic behavior and physically meaningful analysis of control surface loads [8].
For the coupling, the 25 lowest mode shapes are used for the modal basis, since higher modes are dominated by buckling effects. The first four modes are shown in Fig. 6. The first bending mode occurs at f 1 = 16.9 Hz and the first torsional mode is split in two frequencies f 2 = 33.4 Hz and f 3 = 34.3 Hz, where the lower frequency contains a larger contribution of the inner flap movement and the higher frequency a larger contribution of the outer flap movement around their hinge axis. The second bending is observable at the fourth natural frequency f 4 = 51.7 Hz. Clearly, the individual movements of the inner and outer flap with respect to the main wing are contained in the set of 25 modes. This is important for a meaningful analysis of the influence of flexibility effects from the control surfaces. The wing is clamped at the root and inertial forces are not considered. For all structural components, aluminum has been chosen as material. A mesh independence study for the structural grid has been performed and revealed deviations smaller than 1% compared to the finer grid (18,499 nodes, 38,823 triangle elements) in terms of the natural frequencies and corresponding eigenvectors.

Flow field description
Before investigating the influence of flexibility on the aerodynamic characteristics of the Model53, the flow field is computed assuming the wing to be fully rigid. To verify the robustness of the presented flexible Chimera  At = 20 • , the complexity of the flow field around the Model53 wing increases. Compared to = 10 • , a stronger apex vortex located further inboard is observable from the more pronounced suction footprint in the pressure distribution as well as higher vorticity values in the slice contours, Fig. 7b. In addition, large regions of supersonic flow are present above the upper wing surface in comparison to Fig. 7c. The highest Mach number is reached in Fig. 6 First four eigenmodes of the Model53 wing configuration the apex vortex core with Ma max ≈ 2 . Further downstream, a weaker suction footprint and lower vorticity values are present, which can be related to the bursting of the apex vortex at x∕c r ≈ 0.2 . The breakdown of the apex vortex is influenced by complex shock/vortex interactions and occurs downstream of an initial shock, visible in the pressure distribution, Fig. 7c. Downstream of the apex vortex breakdown, the vortical flow is subsonic and additional shocks occur inboard of the vortex, which are clearly visible through the pressure distribution on the symmetry plane, Fig. 7c.
A second vortex separates at the deployed slat. Due to the smaller local angle of attack resulting from the slat deployment, the slat vortex has a different, rather wake-type, topology compared to the apex vortex. The breakdown position of the slat vortex is located in the rear half of the wing and is indicated by low vorticity magnitudes as well as reverse axial flow u∕U ∞ < 0 , Fig. 7b. Similar to the apex vortex, the formation and breakdown of the slat vortex is significantly influenced by shock/vortex interactions. Additionally, an interaction between both vortices is observable from the vorticity contours in Fig. 7b. Figure 8 shows the lift coefficient C L and the lift coefficient deviation ΔC L = C L,flexible − C L,rigid over several flap deflection angles. Accordingly, in  Fig. 8a. The gain in lift coefficient gets smaller for each Δ = 5 • increase in angle of attack, which means that C L has a sub-linear behaviour with respect to for the given configuration and flow conditions. The pitching moment coefficient C My decreases from = −30 • to = 30 • for all angles of attack, except for = 25 • from = −30 • to = −20 • , Fig. 9a, b. Starting from the baseline of = 0 • , the highest lift and pitching moment coefficient increments, or highest aerodynamic control surface efficiency, can be achieved within −10 • < < 10 • for all angles of attack. Both curves flatten towards the minimum and maximum deflection angle. The wing's flexibility has a varying influence on the lift and pitching moment coefficient characteristics over the flap deflection angles dependent on the angle of attack. In general, the trend is that the lift coefficient decreases for ≥ 0 • and increases for < 0 • compared to the rigid configuration, whereas the pitching moment increases for ≥ 0 • and decreases for < 0 • , Figs. 8 and 9. Furthermore, the offset in lift and pitching moment coefficient decreases with higher angle of attack for upward deflections < 0 • . For downward deflections and no deflection ≥ 0 • , the offsets increase for = 25 • compared to = 10 • , however, in between no clear trend can be identified. To understand, which changes in the flow field lead to the differences in the lift and pitching moment coefficient characteristics between the rigid and flexible case, the surface pressure distribution is analyzed for both. Figure 10 depicts the rigid wing on the left and the flexible wing in a mirrored view on the right for selected angles of attack and flap deflection angles. The z-displacement component of the same selection of angles of attack and flap deflection angles is shown in Fig. 11.

General trends
Angle of attack = • : At = 10 • and a flap deflection of = 0 • , the flexible wing provides a lower lift coefficient value compared to the rigid wing, Fig. 8. This results from the reduced local angle of attack due to increased wing twist    Fig. 10a. From the displacement contours, it is apparent that the wing is deformed downwards, Fig. 11a, as the overall configuration generates negative lift due to the high downforce from the flaps, Fig. 8a. The maximum displacement for = −30 • is d z ≈ −0.05 m at the outer flap's wing tip.
For a downward deflection of = 30 • , the most apparent difference in pressure distribution is the increased suction footprint of the wing-tip vortex present on the flexible configuration, which is a result of the deformed wing geometry that enables the flow to maintain supersonic speeds further aft, Fig. 10b. The effective flexible flap deflection angle is eff = 25.8 • for a rigid downward deflection of = 30 • . In contrast to the upward flap deflection = −30 • , a larger contribution of wing twist can be seen for = 30 • , Fig. 11b. A reduction in local angle of attack due to wing twist as well as lower effective flap deflection angles contribute to lower lift and higher pitching moment coefficient values for downward deflected flaps. Even though the largest displacements and biggest difference in effective deflection angle Δ occurs at = 30 • , the highest offset ΔC L and ΔC My is reached at = 20 • , because the flow around the main wing is more sensitive to effective flap deflection angle changes at = 20 • . A similar course for the absolute pitching moment deviations over the flap deflection angles can be seen in Fig. 9c.
Angle of attack = • With a higher angle of attack = 15 • , the offset in lift and pitching moment coefficient between the rigid and flexible configuration rises at the zero flap deflection angle = 0 • compared to = 10 • , Figs. 8 and 9. This is caused by a higher wind loading which leads to larger displacements and hence to a greater reduction in lift and increase in pitching moment coefficient.
Compared to the other angles of attack, an interesting tendency can be seen at = 15 • for downward deflections, where the deviation in lift and pitching moment coefficient changes its sign at = 15 • from negative ΔC L and positive ΔC My at = 20 • , to positive ΔC L and negative ΔC My at = 30 • , Figs. 8b and 9c. From the surface pressure distribution in Fig. 10c, d, it is evident that the decrease of lift for the rigid wing is caused by an upstream movement of the breakdown position of the slat vortex from = 20 • to = 30 • . Concerning the influence of flexibility, a further downstream breakdown position of the slat vortex can be seen for = 20 • and = 30 • , even though the displacement contours for both flap deflection angles do not differ significantly, Fig. 11c, d. Angle of attack = • In accordance with the decreasing offset in lift coefficient between the rigid and flexible configuration from = 10 • to = 20 • for upward flap deflection angles < 0 • , Figs. 8 and 9, the pressure distribution contours differ only marginally for = 20 • , Fig. 10e, f. The reason is that for < 0 • and increasing angles of attack, the downward force or efficiency of the flaps decreases, since the flap's upper sides encounter wake type flow conditions with decreased pressure levels. Less down-force acting on the flaps leads to a smaller difference between the effective flexible and rigid deflection angle and hence a smaller offset in lift and pitching moment coefficient. The effective angle of attack is eff = −29 • for a rigid upward deflection of = −30 • . From the displacement contours, Fig. 11e, it is apparent that the contour lines are more aligned with the x-axis, which relates to a small amount of twist.

Computational effort
The completion of the aerodynamic simulation took 1500 CPU hours on Intel Xeon processors with a core frequency of 3.1 GHz. For the aeroelastic simulation, 2250 CPU hours have been used, which is 1.5 times more compared to the rigid case. Concerning the aeroelastic case, 92% of the computational time has been used for CFD solver iterations, 3 % for the volume grid deformation, 3 % for the fluid-structure coupling and the remaining 2 % for pre-processing and the creation of log files. The increase in computational time for the aeroelastic simulation is considered acceptable given the improvement in the accuracy of the results. Figure 12 depicts the convergence history of the lift and pitching moment coefficient for the aerodynamic and aeroelastic simulation. In addition, the flap deflection angles for the corresponding iteration range is shown.
The iteration at which the rotation of the flap deflection angle is performed, is also clearly visible by identifying strong changes in the course of the lift or pitching moment coefficient. Compared to the initial flap deflection angle, a decrease in the amount of iterations can be seen for subsequent deflection angles. This highlights the advantage of using the Chimera technique to reduce computational time for the simulation of several flap deflections.

Conclusion
In this work, the developed numerical method for performing aeroelastic simulations including deformable control surfaces modeled through the Chimera technique was introduced and described. This method is based on the SimServer multi-disciplinary simulation environment. A fluid-structure coupling was developed which is suited to work in combination with the Chimera technique and the deflection of control surfaces. Furthermore, the aeroelastic simulation process has been adapted to include the rotation of flaps.
Aerodynamic and aeroelastic computations involving the Model53 delta wing configuration were performed at transonic speed and high dynamic pressure. The wing features a deployed slat and an inboard and outboard trailing edge flap. The flaps were deflected from = −30 • to = 30 • using a Chimera mesh in the aerodynamic model. In the structural model, the flaps are attached to the main wing with a hinge concept that ensures a realistic modeling of the control surfaces. The complex flow field around the Model53 configuration at M = 0.8 and Re = 15.1 × 10 7 is characterized by the interaction of vortices and shock waves.
Aerodynamic computations under-predicted lift and over-predicted pitching moment values for upward flap deflections compared to aeroelastic computations. For each angle of attack, the offset between the rigid and flexible wing increased with higher upward flap deflections. With increasing angle of attack, the differences between the aerodynamic and aeroelastic results diminished as a result of decreasing flap down-force. The highest difference for upward flap deflections was observed at = 10 • , = −30 • with ΔC L = +0.047 and ΔC My = −0.014 compared to the rigid wing ( C L,rigid = −0.07 and C My,rigid = 0.204 ). The differences in lift and pitching moment are mainly driven by a change in the effective flap deflection angle, which for upward deflected flaps reaches the highest difference of For downward flap deflections, lift was generally overpredicted and pitching moment under-predicted by the aerodynamic simulations. The deviation in lift and pitching moment coefficient generally increased with higher downward deflection angles, but decreased again after a certain deflection angle. Dependent on the angle of attack, the maximum offset in lift and pitching moment coefficient occurred at differing flap deflection angles, but within the range of 10 • ≤ ≤ 20 • . Starting from = 10 • the offset in lift coefficient values between the rigid and flexible configuration initially decreased until = 20 • , but increased considerably at = 25 • even though the overall deformations and effective  The analysis of the lift and pitching moment characteristics of the rigid and flexible configuration revealed the importance of taking the structural flexibility into account to obtain more accurate results for the considered range of flap deflections. Compared to the aerodynamic simulations, the computational time increased by a factor of 1.5 for the aeroelastic simulations. Thus, in case of high dynamic pressure, coupled CFD-CSM simulations are recommended for delta wings instead of CFD simulations to gain more accurate results.