Direct and adjoint global stability analysis of turbulent transonic flows over a NACA0012 profile

In this work, various turbulent solutions of the two-dimensional (2D) and three-dimensional compressible Reynolds averaged Navier?Stokes equations are analyzed using global stability theory. This analysis is motivated by the onset of flow unsteadiness (Hopf bifurcation) for transonic buffet conditions where moderately high Reynolds numbers and compressible effects must be considered. The buffet phenomenon involves a complex interaction between the separated flow and a shock wave. The efficient numerical methodology presented in this paper predicts the critical parameters, namely, the angle of attack and Mach and Reynolds numbers beyond which the onset of flow unsteadiness appears. The geometry, a NACA0012 profile, and flow parameters selected reproduce situations of practical interest for aeronautical applications. The numerical computation is performed in three steps. First, a steady baseflow solution is obtained; second, the Jacobian matrix for the RANS equations based on a finite volume discretization is computed; and finally, the generalized eigenvalue problem is derived when the baseflow is linearly perturbed. The methodology is validated predicting the 2D Hopf bifurcation for a circular cylinder under laminar flow condition. This benchmark shows good agreement with the previous published computations and experimental data. In the transonic buffet case, the baseflow is computed using the Spalart?Allmaras turbulence model and represents a mean flow where the high frequency content and length scales of the order of the shear-layer thickness have been averaged. The lower frequency content is assumed to be decoupled from the high frequencies, thus allowing a stability analysis to be performed on the low frequency range. In addition, results of the corresponding adjoint problem and the sensitivity map are provided for the first time for the buffet problem. Finally, an extruded three-dimensional geometry of the NACA0012 airfoil, where all velocity components are considered, was also analyzed as a Triglobal stability case, and the outcoming results were compared to the previous 2D limited model, confirming that the buffet onset is well detected.

Global Stability Analysis of Turbulent Transonic Flows on Airfoil geometries

La aparición de inestabilidades en un flujo es un problema importante que puede afectar a algunas aplicaciones aerodinámicas. De hecho existen diferentes tipos de fenómenos no-estacionarios que actualmente son tema de investigación; casos como la separación a altos ángulos de ataque o el buffet transónico son dos ejemplos de cierta relevancia. El análisis de estabilidad global permite identificar la aparición de dichas condiciones inestables, proporcionando información importante sobre la región donde la inestabilidad es dominante y sobre la frecuencia del fenómeno inestable. La metodología empleada es capaz de calcular un flujo base promediado mediante una discretización con volúmenes finitos y posteriormente la solución de un problema de autovalores asociado a la linealización que aparece al perturbar el flujo base. El cálculo numérico se puede dividir en tres pasos: primero se calcula una solución estacionaria para las ecuaciones RANS, luego se extrae la matriz del Jacobiano que representa el problema linealizado y finalmente se deriva y se resuelve el problema de autovalores generalizado mediante el método iterativo de Arnoldi. Como primer caso de validación, la técnica descrita ha sido aplicada a un cilindro circular en condiciones laminares para detectar el principio de las oscilaciones de los vórtices de von Karman, y se han comparado los resultados con experimentos y cálculos anteriores. La parte más importante del estudio se centra en el análisis de flujos compresibles en régimen turbulento. La predicción de la aparición y la progresión de flujo separado a altos ángulos de ataque se han estudiado en el perfil NACA0012 en condiciones tanto subsónicas como supersónicas y en una sección del ala del A310 en condiciones de despegue. Para todas las geometrías analizadas, se ha podido observar que la separación gradual genera la aparición de un modo inestable específico para altos ángulos de ataque siempre mayores que el ángulo asociado al máximo coeficiente de sustentación. Además, se ha estudiado el problema adjunto para obtener información sobre la zona donde una fuerza externa provoca el máximo cambio en el campo fluido. El estudio se ha completado calculando el mapa de sensibilidad estructural y localizando el centro de la inestabilidad. En el presente trabajo de tesis se ha analizado otro importante fenómeno: el buffet transónico. En condiciones transónicas, la interacción entre la onda de choque y la capa límite genera una oscilación de la posición de la onda de choque y, por consiguiente, de las fuerzas aerodinámicas. El conocimiento de las condiciones críticas y su origen puede ayudar a evitar la oscilación causada por estas fuerzas. Las condiciones para las cuales comienza la inestabilidad han sido calculadas y comparadas con trabajos anteriores. Por otra parte, los resultados del correspondiente problema adjunto y el mapa de sensibilidad se han obtenido por primera vez para el buffet, indicando la región del dominio que sera necesario modificar para crear el mayor cambio en las propiedades del campo fluido. Dado el gran consumo de memoria requerido para los casos 3D, se ha realizado un estudio sobre la reducción del domino con la finalidad de reducirlo a la región donde está localizada la inestabilidad. La eficacia de dicha reducción de dominio ha sido evaluada investigando el cambio en la dimensión de la matriz del Jacobiano, no resultando muy eficiente en términos del consumo de memoria. Dado que el buffet es un problema en general tridimensional, el análisis TriGlobal de una geometría 3D podría considerarse el auténtico reto futuro. Como aproximación al problema, un primer estudio se ha realizado empleando una geometría tridimensional extruida del NACA00f2. El cálculo del flujo 3D y, por primera vez en casos tridimensionales compresibles y turbulentos, el análisis de estabilidad TriGlobal, se han llevado a cabo. La comparación de los resultados obtenidos con los resultados del anterior modelo 2D, ha permitido, primero, verificar la exactitud del cálculo 2D realizado anteriormente y también ha proporcionado una estimación del consumo de memoria requerido para el caso 3D. ABSTRACT Flow unsteadiness is an important problem in aerodynamic applications. In fact, there are several types of unsteady phenomena that are still at the cutting edge of research in the field; separation at high angles of attack and transonic buffet are two important examples. Global Stability Analysis can identify the unstable onset conditions, providing important information about the instability location in the domain and the frequency of the unstable phenomenon. The methodology computes a base flow averaged state based on a finite volume discretization and a solution for a generalized eigenvalue problem corresponding to the perturbed linearized equations. The numerical computation is then performed in three steps: first, a steady solution for the RANS equation is computed; second, the Jacobian matrix that represents the linearized problem is obtained; and finally, the generalized eigenvalue problem is derived and solved with an Arnoldi iterative method. As a first validation test, the technique has been applied on a laminar circular cylinder in order to detect the von Karman vortex shedding onset, comparing the results with experiments and with previous calculations. The main part of the study focusses on turbulent and compressible cases. The prediction of the origin and progression of separated flows at high angles of attack has been studied on the NACA0012 airfoil at subsonic and transonic conditions and for the A310 airfoil in take-off configuration. For all the analyzed geometries, it has been found that gradual separation generates the appearance of one specific unstable mode for angles of attack always greater than the ones related to the maximum lift coefficient. In addition, the adjoint problem has been studied to suggest the location of an external force that results in the largest change to the flow field. From the direct and the adjoint analysis the structural sensitivity map has been computed and the core of the instability has been located. The other important phenomenon analyzed in this work is the transonic buffet. In transonic conditions, the interaction between the shock wave and the boundary layer leads to an oscillation of the shock location and, consequently, of the aerodynamic forces. Knowing the critical operational conditions and its origin can be helpful in preventing such fluctuating forces. The instability onset has then been computed and compared with the literature. Moreover, results of the corresponding adjoint problem and a sensitivity map have been provided for the first time for the buffet problem, indicating the region that must be modified to create the biggest change in flow field properties. Because of the large memory consumption required when a 3D case is approached, a domain reduction study has been carried out with the aim of limiting the domain size to the region where the instability is located. The effectiveness of the domain reduction has been evaluated by investigating the change in the Jacobian matrix size, not being very efficient in terms of memory consumption. Since buffet is a three-dimensional problem, TriGlobal stability analysis can be seen as a future challenge. To approximate the problem, a first study has been carried out on an extruded three-dimensional geometry of the NACA0012 airfoil. The 3D flow computation and the TriGlobal stability analysis have been performed for the first time on a compressible and turbulent 3D case. The results have been compared with a 2D model, confirming that the buffet onset evaluated in the 2D case is well detected. Moreover, the computation has given an indication about the memory consumption for a 3D case.

Left ventricular mass in patients with type 2 diabetes is independently associated with central but not peripheral pulse pressure

An increase in left ventricular mass (LVM) occurs in the presence of type 2 diabetes, apparently independent of hypertension (1), but the determinants of this process are unknown. Brachial blood pressure is not representative of that at the ascending aorta (2) because the pressure wave is amplified from central to peripheral arteries. Central blood pressure is probably more clinically important since local pulsatile pressure determines adverse arterial and myocardial remodeling (3,4). Thus, an inaccurate assessment of the contribution of arterial blood pressure to LVM may occur if only brachial blood pressure is taken into consideration. In this study we sought the contribution of central blood pressure (and other interactive factors known to affect wave reflection, e.g., glycemic control and total arterial compliance) to LVM in patients with type 2 diabetes.

A systematic review and meta-analysis of clinical trials on physical interventions for lateral epicondylalgia

A systematic review of the literature on the effectiveness of physical interventions for lateral epicondylalgia ( tennis elbow) was carried out. Seventy six randomised controlled trials were identified, 28 of which satisfied the minimum criteria for meta-analysis. The evidence suggests that extracorporeal shock wave therapy is not beneficial in the treatment of tennis elbow. There is a lack of evidence for the long term benefit of physical interventions in general. However, further research with long term follow up into manipulation and exercise as treatments is indicated.

The effect of exercise on large artery haemodynamics in healthy young men

Background Brachial blood pressure predicts cardiovascular outcome at rest and during exercise. However, because of pulse pressure amplification, there is a marked difference between brachial pressure and central (aortic) pressure. Although central pressure is likely to have greater clinical importance, very little data exist regarding the central haemodynamic response to exercise. The aim of the present study was to determine the central and peripheral haemodynamic response to incremental aerobic exercise. Materials and methods Twelve healthy men aged 31 +/- 1 years (mean +/- SEM) exercised at 50%, 60%, 70% and 80% of their maximal heart rate (HRmax) on a bicycle ergometer. Central blood pressure and estimated aortic pulse wave velocity, assessed by timing of the reflected wave (T-R), were obtained noninvasively using pulse wave analysis. Pulse pressure amplification was defined as the ratio of peripheral to central pulse pressure and, to assess the influence of wave reflection on amplification, the ratio of peripheral pulse pressure to nonaugmented central pulse pressure (PPP : CDBP-P-1) was also calculated. Results During exercise, there was a significant, intensity-related, increase in mean arterial pressure and heart rate (P < 0.001). There was also a significant increase in pulse pressure amplification and in PPP : CDBP-P-1 (P < 0.001), but both were independent of exercise intensity. Estimated aortic pulse wave velocity increased during exercise (P < 0.001), indicating increased aortic stiffness. There was also a positive association between aortic pulse wave velocity and mean arterial pressure (r = 0.54; P < 0.001). Conclusions Exercise significantly increases pulse pressure amplification and estimated aortic stiffness.

Validation of a generalized transfer function to noninvasively derive central blood pressure during exercise

Exercise brachial blood pressure ( BP) predicts mortality, but because of wave reflection, central ( ascending aortic) pressure differs from brachial pressure. Exercise central BP may be clinically important, and a noninvasive means to derive it would be useful. The purpose of this study was to test the validity of a noninvasive technique to derive exercise central BP. Ascending aortic pressure waveforms were recorded using a micromanometer-tipped 6F Millar catheter in 30 patients (56 +/- 9 years; 21 men) undergoing diagnostic coronary angiography. Simultaneous recordings of the derived central pressure waveform were acquired using servocontrolled radial tonometry at rest and during supine cycling. Pulse wave analysis of the direct and derived pressure signals was performed offline (SphygmoCor 7.01). From rest to exercise, mean arterial pressure and heart rate were increased by 20 +/- 10 mm Hg and 15 +/- 7 bpm, respectively, and central systolic BP ranged from 77 to 229 mm Hg. There was good agreement and high correlation between invasive and noninvasive techniques with a mean difference (+/- SD) for central systolic BP of -1.3 +/- 3.2 mm Hg at rest and -4.7 +/- 3.3 mm Hg at peak exercise ( for both r=0.995; P < 0.001). Conversely, systolic BP was significantly higher peripherally than centrally at rest (155 +/- 33 versus 138 +/- 32mm Hg; mean difference, -16.3 +/- 9.4mm Hg) and during exercise (180 +/- 34 versus 164 +/- 33 mm Hg; mean difference, -15.5 +/- 10.4 mm Hg; for both P < 0.001). True myocardial afterload is not reliably estimated by peripheral systolic BP. Radial tonometry and pulse wave analysis is an accurate technique for the noninvasive determination of central BP at rest and during exercise.

Hybrid approaches to estimating freeway travel times using point traffic detector data

The accurate and reliable estimation of travel time based on point detector data is needed to support Intelligent Transportation System (ITS) applications. It has been found that the quality of travel time estimation is a function of the method used in the estimation and varies for different traffic conditions. In this study, two hybrid on-line travel time estimation models, and their corresponding off-line methods, were developed to achieve better estimation performance under various traffic conditions, including recurrent congestion and incidents. The first model combines the Mid-Point method, which is a speed-based method, with a traffic flow-based method. The second model integrates two speed-based methods: the Mid-Point method and the Minimum Speed method. In both models, the switch between travel time estimation methods is based on the congestion level and queue status automatically identified by clustering analysis. During incident conditions with rapidly changing queue lengths, shock wave analysis-based refinements are applied for on-line estimation to capture the fast queue propagation and recovery. Travel time estimates obtained from existing speed-based methods, traffic flow-based methods, and the models developed were tested using both simulation and real-world data. The results indicate that all tested methods performed at an acceptable level during periods of low congestion. However, their performances vary with an increase in congestion. Comparisons with other estimation methods also show that the developed hybrid models perform well in all cases. Further comparisons between the on-line and off-line travel time estimation methods reveal that off-line methods perform significantly better only during fast-changing congested conditions, such as during incidents. The impacts of major influential factors on the performance of travel time estimation, including data preprocessing procedures, detector errors, detector spacing, frequency of travel time updates to traveler information devices, travel time link length, and posted travel time range, were investigated in this study. The results show that these factors have more significant impacts on the estimation accuracy and reliability under congested conditions than during uncongested conditions. For the incident conditions, the estimation quality improves with the use of a short rolling period for data smoothing, more accurate detector data, and frequent travel time updates.

Theoretical Performance Model and Initial Experimentation of a Baffled-Tube Ram Accelerator

Thesis (Master's)--University of Washington, 2016-08

P-wave velocities of samples from the Penninic of the western Alps along NFP20-west seismic-reflection profiles

Using the transit pulse method, we have determined compressional wave velocities of rocks from various geological units belonging to the Penninic zone along the NFP20-West profiles of the Swiss western Alps. The velocities have been measured at confining pressures up to 400 MPa, along three orthogonal axes defined by the macrostructure of the rocks. The samples analysed show a degree of metamorphism ranging from greenschist to eclogite facies. This collection includes schists, dolomites, gneisses and ophiolitic rocks. The mean velocities range from 5.9 km/s for a quartzitic calcschist to 7.9 km/s for an eclogitic metagabbro. The velocity anisotropy is as high as 20 %. The range of acoustic impedance is wide, from 15 to 27 10(6) kg/m2s. From these measurements, normal incident reflection coefficients for likely rock assemblages within and between geological units were estimated in order to interpret zone of the strong reflections recorded along the seismic profiles. Reflection coefficients as high as 0.17 could be determined.

A similarity solution for a shock reflection problem in isentropic gas dynamics

A one-dimensional shock-reflection test problem in the case of slab, cylindrical, or spherical symmetry is discussed. The differential equations for a similarity solution are derived and solved numerically in conjunction with the Rankie-Hugoniot shock relations.

Reflection characteristics of screen wave absorbers : hydraulic model investigation /

Bibliography: p. 42.

Water wave transmission through and reflection by pervious coastal structures : hydraulic laboratory investigation /

"October 1969."

Two-dimensional tests of wave transmission and reflection characteristics of laboratory breakwaters /

"June 1980."