scholarly journals Development of an Aerostructural Analysis Tool for a Low-Sweep Parametric Wing

Author(s):  
Julian Bardin

An aerostructural analysis program was developed to predict the aerodynamic performance of a non-rigid, low-sweep wing. The wing planform was geometrically defined to have a rectangular section, and a trapezoidal section. The cross-section was further set to an airfoil shape which was consistent across the entire wingspan. Furthermore, to enable the inclusion of this multidisciplinary analysis module into an optimization scheme, the wing geometry was defined by a series of parameters: root chord, taper ratio, leading-edge sweep, semi-span length, and the kink location. Aerodynamic analysis was implemented through the quasi-three-dimensional approach, including a three-dimensional inviscid solution and a sectional two-dimensional viscous solution. The inviscid analysis was provided through the implementation of the vortex ring lifting surface method, which modelled the wing about its mean camber surface. The viscous aerodynamic solution was implemented through a sectional slicing of the wing. For each section, the effective angle of attack was determined and provided as an input to a two-dimensional airfoil solver. This airfoil solution was comprised of two subcomponents: a linear-strength vortex method inviscid solution, and a direct-method viscous boundary layer computation. The converged airfoil solution was developed by adjusting the effective airfoil geometry to account for the boundary layer displacement thickness, which in itself required the inviscid tangential speeds to compute. The structural solution was implemented through classical beam theory, with a torsion and bending calculator included. The torque and bending moment distribution along the wing were computed from the lift distribution, neglecting the effects of drag, and used to compute the twist and deflection of the wing. Interdisciplinary coupling was achieved through an iterative scheme. With the developed implementation, the inviscid lift loads were used to compute the deformation of the wing. This deformation was used to update the wing mesh, and the inviscid analysis was run again. This iteration was continued until the lift variation between computations was below 0.1%. Once the solution was converged upon by the inviscid and structural solutions, the viscous calculator was run to develop the parasitic drag forces. Once computation had completed, the aerodynamic lift and drag forces were output to mark the completion of execution.

2021 ◽  
Author(s):  
Julian Bardin

An aerostructural analysis program was developed to predict the aerodynamic performance of a non-rigid, low-sweep wing. The wing planform was geometrically defined to have a rectangular section, and a trapezoidal section. The cross-section was further set to an airfoil shape which was consistent across the entire wingspan. Furthermore, to enable the inclusion of this multidisciplinary analysis module into an optimization scheme, the wing geometry was defined by a series of parameters: root chord, taper ratio, leading-edge sweep, semi-span length, and the kink location. Aerodynamic analysis was implemented through the quasi-three-dimensional approach, including a three-dimensional inviscid solution and a sectional two-dimensional viscous solution. The inviscid analysis was provided through the implementation of the vortex ring lifting surface method, which modelled the wing about its mean camber surface. The viscous aerodynamic solution was implemented through a sectional slicing of the wing. For each section, the effective angle of attack was determined and provided as an input to a two-dimensional airfoil solver. This airfoil solution was comprised of two subcomponents: a linear-strength vortex method inviscid solution, and a direct-method viscous boundary layer computation. The converged airfoil solution was developed by adjusting the effective airfoil geometry to account for the boundary layer displacement thickness, which in itself required the inviscid tangential speeds to compute. The structural solution was implemented through classical beam theory, with a torsion and bending calculator included. The torque and bending moment distribution along the wing were computed from the lift distribution, neglecting the effects of drag, and used to compute the twist and deflection of the wing. Interdisciplinary coupling was achieved through an iterative scheme. With the developed implementation, the inviscid lift loads were used to compute the deformation of the wing. This deformation was used to update the wing mesh, and the inviscid analysis was run again. This iteration was continued until the lift variation between computations was below 0.1%. Once the solution was converged upon by the inviscid and structural solutions, the viscous calculator was run to develop the parasitic drag forces. Once computation had completed, the aerodynamic lift and drag forces were output to mark the completion of execution.


1976 ◽  
Vol 74 (1) ◽  
pp. 59-79 ◽  
Author(s):  
T. J. Pedley

The viscous boundary layer on a finite flat plate in a stream which reverses its direction once (at t = 0) is analysed using an improved version of the approximate method described earlier (Pedley 1975). Long before reversal (t < −t1), the flow at a point on the plate will be quasi-steady; long after reversal (t > t2), the flow will again be quasi-steady, but with the leading edge at the other end of the plate. In between (−t1 < t < t2) the flow is governed approximately by the diffusion equation, and we choose a simple solution of that equation which ensures that the displacement thickness of the boundary layer remains constant at t = −t1. The results of the theory, in the form of the wall shear rate at a point as a function of time, are given both for a uniformly decelerating stream, and for a sinusoidally oscillating stream which reverses its direction twice every cycle. The theory is further modified to cover streams which do not reverse, but for which the quasi-steady solution breaks down because the velocity becomes very small. The analysis is also applied to predict the wall shear rate at the entrance to a straight pipe when the core velocity varies with time as in a dog's aorta. The results show positive and negative peak values of shear very much larger than the mean. They suggest that, if wall shear is implicated in the generation of atherosclerosis because it alters the permeability of the wall to large molecules, then an appropriate index of wall shear at a point is more likely to be the r.m.s. value than the mean.


1984 ◽  
Vol 106 (2) ◽  
pp. 260-267 ◽  
Author(s):  
R. J. Goldstein ◽  
J. Karni

A naphthalene sublimation technique is used to determine the circumferential and longitudinal variations of mass transfer from a smooth circular cylinder in a crossflow of air. The effect of the three-dimensional secondary flows near the wall-attached ends of a cylinder is discussed. For a cylinder Reynolds number of 19000, local enhancement of the mass transfer over values in the center of the tunnel are observed up to a distance of 3.5 cylinder diameters from the tunnel wall. In a narrow span extending from the tunnel wall to about 0.066 cylinder diameters above it (about 0.75 of the mainstream boundary layer displacement thickness), increases of 90 to 700 percent over the two-dimensional flow mass transfer are measured on the front portion of the cylinder. Farther from the wall, local increases of up to 38 percent over the two-dimensional values are measured. In this region, increases of mass transfer in the rear portion of the cylinder, downstream of separation, are, in general, larger and cover a greater span than the increases in the front portion of the cylinder.


1971 ◽  
Vol 22 (4) ◽  
pp. 346-362 ◽  
Author(s):  
J. F. Nash ◽  
R. R. Tseng

SummaryThis paper presents the results of some calculations of the incompressible turbulent boundary layer on an infinite yawed wing. A discussion is made of the effects of increasing lift coefficient, and increasing Reynolds number, on the displacement thickness, and on the magnitude and direction of the skin friction. The effects of the state of the boundary layer (laminar or turbulent) along the attachment line are also considered.A study is made to determine whether the behaviour of the boundary layer can adequately be predicted by a two-dimensional calculation. It is concluded that there is no simple way to do this (as is provided, in the laminar case, by the principle of independence). However, with some modification, a two-dimensional calculation can be made to give an acceptable numerical representation of the chordwise components of the flow.


2001 ◽  
Vol 426 ◽  
pp. 73-94 ◽  
Author(s):  
A. A. MASLOV ◽  
A. N. SHIPLYUK ◽  
A. A. SIDORENKO ◽  
D. ARNAL

Experimental investigations of the boundary layer receptivity, on the sharp leading edge of a at plate, to acoustic waves induced by two-dimensional and three- dimensional perturbers, have been performed for a free-stream Mach number M∞ = 5.92. The fields of controlled free-stream disturbances were studied. It was shown that two-dimensional and three-dimensional perturbers radiate acoustic waves and that these perturbers present a set of harmonic motionless sources and moving sources with constant amplitude. The disturbances excited in the boundary layer were measured. It was found that acoustic waves impinging on the leading edge generate Tollmien–Schlichting waves in the boundary layer. The receptivity coefficients were obtained for several radiation conditions and intensities. It was shown that there is a dependence of receptivity coefficients on the wave inclination angles.


2013 ◽  
Vol 724 ◽  
pp. 284-304 ◽  
Author(s):  
D. S. Pearson ◽  
P. J. Goulart ◽  
B. Ganapathisubramani

AbstractThe turbulent flow over a forward-facing step is studied using two-dimensional time-resolved particle image velocimetry. The structure and behaviour of the separation region in front of the step is investigated using conditional averages based on the area of reverse flow present. The relation between the position of the upstream separation and the two-dimensional shape of the separation region is presented. It is shown that when of ‘closed’ form, the separation region can become unstable resulting in the ejection of fluid over the corner of the step. The separation region is shown to grow simultaneously in both the wall-normal and streamwise directions, to a point where the maximum extent of the upstream position of separation is limited by the accompanying transfer of mass over the step corner. The conditional averages are traced backwards in time to identify the average behaviour of the boundary-layer displacement thickness leading up to such events. It is shown that these ejections are preceded by the convection of low-velocity regions from upstream, resulting in a three-dimensional interaction within the separation region. The size of the low-velocity regions, and the time scale at which the separation region fluctuates, is shown to be consistent with the large boundary layer structures observed in the literature. Instances of a highly suppressed separation region are accompanied by a steady increase in velocity in the upstream boundary layer.


2001 ◽  
Vol 439 ◽  
pp. 305-333 ◽  
Author(s):  
ZHIYIN YANG ◽  
PETER R. VOKE

Transition arising from a separated region of flow is quite common and plays an important role in engineering. It is difficult to predict using conventional models and the transition mechanism is still not fully understood. We report the results of a numerical simulation to study the physics of separated boundary-layer transition induced by a change of curvature of the surface. The geometry is a flat plate with a semicircular leading edge. The Reynolds number based on the uniform inlet velocity and the leading-edge diameter is 3450. The simulated mean and turbulence quantities compare well with the available experimental data.The numerical data have been comprehensively analysed to elucidate the entire transition process leading to breakdown to turbulence. It is evident from the simulation that the primary two-dimensional instability originates from the free shear in the bubble as the free shear layer is inviscidly unstable via the Kelvin–Helmholtz mechanism. These initial two-dimensional instability waves grow downstream with a amplification rate usually larger than that of Tollmien–Schlichting waves. Three-dimensional motions start to develop slowly under any small spanwise disturbance via a secondary instability mechanism associated with distortion of two-dimensional spanwise vortices and the formation of a spanwise peak–valley wave structure. Further downstream the distorted spanwise two-dimensional vortices roll up, leading to streamwise vorticity formation. Significant growth of three-dimensional motions occurs at about half the mean bubble length with hairpin vortices appearing at this stage, leading eventually to full breakdown to turbulence around the mean reattachment point. Vortex shedding from the separated shear layer is also observed and the ‘instantaneous reattachment’ position moves over a distance up to 50% of the mean reattachment length. Following reattachment, a turbulent boundary layer is established very quickly, but it is different from an equilibrium boundary layer.


1990 ◽  
Vol 217 ◽  
pp. 367-390 ◽  
Author(s):  
Philip Hall ◽  
Sharon O. Seddougui

The three-dimensional boundary layer on a swept wing can support different types of hydrodynamic instability. Here attention is focused on the so-called ‘spanwise instability’ problem which occurs when the attachment-line boundary layer on the leading edge becomes unstable to Tollmien–Schlichting waves. In order to gain insight into the interactions that are important in that problem a simplified basic state is considered. This simplified flow corresponds to the swept attachment-line boundary layer on an infinite flat plate. The basic flow here is an exact solution of the Navier–Stokes equations and its stability to two-dimensional waves propagating along the attachment line can be considered exactly at finite Reynolds number. This has been done in the linear and weakly nonlinear regimes by Hall, Malik & Poll (1984) and Hall & Malik (1986). Here the corresponding problem is studied for oblique waves and their interaction with two-dimensional waves is investigated. In fact oblique modes cannot be described exactly at finite Reynolds number so it is necessary to make a high-Reynolds-number approximation and use triple-deck theory. It is shown that there are two types of oblique wave which, if excited, cause the destabilization of the two-dimensional mode and the breakdown of the disturbed flow at a finite distance from the leading edge. First a low-frequency mode closely related to the viscous stationary crossflow mode discussed by Hall (1986) and MacKerrell (1987) is a possible cause of breakdown. Secondly a class of oblique wave with frequency comparable with that of the two-dimensional mode is another cause of breakdown. It is shown that the relative importance of the modes depends on the distance from the attachment line.


2014 ◽  
Vol 741 ◽  
pp. 228-251 ◽  
Author(s):  
Mochamad Dady Ma’mun ◽  
Masahito Asai ◽  
Ayumu Inasawa

AbstractThe effects of surface corrugation with small amplitude on the growth of Tollmien–Schlichting (T–S) waves were examined experimentally in a zero-pressure-gradient boundary layer. Two- and three-dimensional corrugations of sinusoidal geometry with wavelengths of the same order as that of the two-dimensional T–S wave were considered. The corrugation amplitudes were one order of magnitude smaller than the boundary-layer displacement thickness. Streamwise growth of T–S waves on the corrugated walls was compared with that in the boundary layer on the smooth surface. A distinct difference was found in the destabilizing effect between the two- and three-dimensional corrugations. The two-dimensional corrugation significantly enhanced the growth of two-dimensional T–S waves even when the corrugation amplitude was only ∼10% of the displacement thickness. On decreasing the corrugation amplitude, the growth rate of two-dimensional T–S waves asymptotically approached that in the smooth-wall case. On the other hand, the three-dimensional corrugation had only a small influence on the growth of two-dimensional T–S waves even when the corrugation amplitude was as large as 20% of the displacement thickness. For three-dimensional corrugations, however, a pair of oblique waves was generated and developed by an interaction between the two-dimensional T–S wave and the corrugation-induced mean-flow distortion for the corrugation wavelength considered. On increasing the corrugation amplitude, the oblique waves generated were increased in amplitude and thus significantly influenced the secondary instability process.


Sign in / Sign up

Export Citation Format

Share Document