Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary value problem”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Transformation of two and three-dimensional regions by elliptic systems

A reliable linear system is presented for grid generation in 2-D and 3-D. The method is robust in the sense that convergence is guaranteed but is not as reliable as other nonlinear elliptic methods in generating nonfolding grids. The construction of nonfolding grids depends on having reasonable approximations of cell aspect ratios and an appropriate distribution of grid points on the boundary of the region. Some guidelines are included on approximating the aspect ratios, but little help is offered on setting up the boundary grid other than to say that in 2-D the boundary correspondence should be close to that generated by a conformal mapping. It is assumed that the functions which control the grid distribution depend only on the computational variables and not on the physical variables. Whether this is actually the case depends on how the grid is constructed. In a dynamic adaptive procedure where the grid is constructed in the process of solving a fluid flow problem, the grid is usually updated at fixed iteration counts using the current value of the control function. Since the control function is not being updated during the iteration of the grid equations, the grid construction is a linear procedure. However, in the case of a static adaptive procedure where a trial solution is computed and used to construct an adaptive grid, the control functions may be recomputed at every step of the grid iteration.

Mastin, C. Wayne↗
NASA NTRS1994

Multilayer Thin Film Polarizer Design for Far Ultraviolet using Induced Transmission and Absorption Technique

Good theoretical designs of far ultraviolet polarizers have been reported using a MgF2/Al/MgF2 three layer structure on a thick Al layer as a substrate. The thicknesses were determined to induce transmission and absorption of p-polarized light. In these designs Al optical constants were used from films produced in ultrahigh vacuum (UHV: 10(exp -10) torr). Reflectance values for polarizers fabricated in a conventional high vacuum (p approx. 10(exp -6 torr)) using the UHV design parameters differed dramatically from the design predictions. Al is a highly reactive material and is oxidized even in a high vacuum chamber. In order to solve the problem other metals have been studied. It is found that a larger reflectance difference is closely related to higher amplitude and larger phase difference of Fresnel reflection coefficients between two polarizations at the boundary of MgF2/metal. It is also found that for one material a larger angle of incidence from the surface normal brings larger amplitude and phase difference. Be and Mo are found good materials to replace Al. Polarizers designed for 121.6 nm with Be at 60 deg and with Mo at 70 deg are shown as examples.

Kim, Jongmin↗
NASA NTRS2015

Analysis of Advanced Thermoelectric Materials and Their Functional Limits

The world's demand for energy is increasing dramatically, but the best energy conversion systems operate at approximately 30% efficiency. One way to decrease energy loss is in the recovery of waste heat using thermoelectric (TE) generators. A TE generator is device that generates electricity by exploiting heat flow across a thermal gradient. The efficiency of a TE material for power generation and cooling is determined by the dimensionless Figure of Merit (ZT): ZT = S(exp. 2)sigmaT/κ: where S is the Seebeck coefficient, sigma is the electrical conductivity, T is the absolute temperature, and κ is the thermal conductivity. The parameters are not physically independent, but intrinsically coupled since they are a function of the transport properties of electrons. Traditional research on TE materials has focused on synthesizing bulk semiconductor-type materials that have low thermal conductivity and high electrical conductivity affording ZT values of 1. The optimization of the σ/κ ratio is difficult to achieve using current material formats, as these material constants are complementary. Recent areas of research are focusing on using nanostructural artifacts that introduce specific dislocations and boundary conditions that scatter the phonons. This disrupts the physical link between thermal (phonon) and electrical (electron) transport. The result is that κ is decreased without decreasing σ. These material formats give ZT values of up to 2 which represent approximately 18% energy gain from waste heat recovery. The next challenge in developing the next generation of TE materials with superior performance is to tailor the interconnected thermoelectric physical parameters of the material system. In order to approach this problem, the fundamental physics of each parameter S, sigma, and κ need to be physically understood in their context of electron/phonon interaction for the construction of new high ZT thermoelectric devices. Is it possible to overcome the physical limit imposed by of the effect of phonon lattice oscillation and energetic electrons towards thermal conductivity? Is the Seebeck coefficient, based on the difference in voltage over temperature gradient ( deltaV/deltaT), an intrinsic parameter of each material? All these parameters were manipulated using nano-bridge and twin-lattice structural concepts at the NASA Langley Research Center. This talk will review the current trend of TE research to optimize the ZT and discuss about new approaches on increasing ZT within functional limits of each parameter.

Kim, Hyun Jung↗
NASA NTRS2021

Evaluation of Stereology for Snow Microstructure Measurement and Microwave Emission Modeling: A Case Study

Reliable microstructure measurement of snow is a requirement for microwave radiative transfer model validation. Snow specific surface area (SSA) can be measured using stereological methods, in which snow samples are cast in the field and photographed in the laboratory. Processing stereology photographs manually by counting intersections of test cycloids with air–ice boundaries reduces the problems in binary segmentation. This paper is a case study to evaluate the repeatability of the manually stereology interpretation by two independent research groups. We further assessed how uncertainty in snow SSA influences simulated brightness temperature(TB) driven by the Microwave Emission Model of Layered Snowpacks (MEMLS), and how stereology compares to Near Infrared (NIR) camera and hand lens. Data was obtained from two alpine snow profiles from Steamboat Springs, Colorado. Results showed that stereological SSA values measured by two groups are highly consistent, and the ground radiometer measured TBat 19 and 37 GHz was successfully predicted (RMSE<3.8 K);simulations using NIR SSA and hand-lens geometric grain size (Dg)measurements have larger errors. This conclusion was not sensitive to uncertainty in the free parameters of TB modeling.

Microwave radiometry↗
NASA NTRS1974

Finite element solution theory for three-dimensional boundary flows

A finite element algorithm is derived for the numerical solution of a three-dimensional flow field described by a system of initial-valued, elliptic boundary value partial differential equations. The familiar three-dimensional boundary layer equations belong to this description when diffusional processes in only one coordinate direction are important. The finite element algorithm transforms the original description into large order systems of ordinary differential equations written for the dependent variables discretized at node points of an arbitrarily irregular computational lattice. The generalized elliptic boundary conditions is piecewise valid for each dependent variable on boundaries that need not explicitly coincide with coordinate surfaces. Solutions for sample problems in laminar and turbulent boundary flows illustrate favorable solution accuracy, convergence, and versatility.

Baker, A. J.↗
NASA NTRS1944

On the general theory of thin airfoils for nonuniform motion

General thin-airfoil theory for a compressible fluid is formulated as boundary problem for the velocity potential, without recourse to the theory of vortex motion. On the basis of this formulation the integral equation of lifting-surface theory for an incompressible fluid is derived with the chordwise component of the fluid velocity at the airfoil as the function to be determined. It is shown how by integration by parts this integral equation can be transformed into the Biot-Savart theorem. A clarification is gained regarding the use of principal value definitions for the integral which occur. The integral equation of lifting-surface theory is used a s the starting point for the establishment of a theory for the nonstationary airfoil which is a generalization of lifting-line theory for the stationary airfoil and which might be called "lifting-strip" theory. Explicit expressions are given for section lift and section moment in terms of the circulation function, which for any given wing deflection is to be determined from an integral equation which is of the type of the equation of lifting-line theory. The results obtained are for airfoils of uniform chord. They can be extended to tapered airfoils. One of the main uses of the results should be that they furnish a practical means for the analysis of the aerodynamic span effect in the problem of wing flutter. The range of applicability of "lifting-strip" theory is the same as that of lifting-line theory so that its results may be applied to airfoils with aspect ratios as low as three.

BIBLIOGRAPHY-WINGS↗
NASA NTRS2000

Comparison of Analytical and CAA Solution(s): Unsteady Response of a Rectilinear Swept Cascade to an Incident Gust

This problem is aimed at testing the ability of computational aeroacoustic methods to calculate the acoustic pressure field generated as a result of the interaction of a convected harmonic gust with a rectilinear cascade of swept flat plates. All relevant geometric information is the same as the unsteady response of an isolated finite span swept airfoil to an incident gust problem with the exception of the sweep angle alpha, which is now taken to be a variable. Assume a cascade stagger angle of zero (with the x-axis aligned with the chord) and a gap-to-chord ratio of h/c = 1. The mean flow Mach number M(sub 0), gust frequency omega, gust amplitude A, and chordwise wavenumber k(sub x) are the same. But for the cascade problem take k(sub y) = pi and k(sub z) = 0. The appropriate physical scales are also the same, as is the governing equation and boundary conditions. Naturally, the impermeability condition now applies to the entire cascade. For this problem, determine the amplitude of upstream-radiated rms acoustic pressure as a function of the sweep angle at the specified frequency. Specifically, show the variations of rms acoustic pressure amplitude at the upstream location (-5c, 0, l/2) for sweep angle alpha in the range (0.0 deg, 30.0 deg). Use sweep angle increments no larger than 2.5 deg. Express the results in dB using the rms pressure value for alpha = 0.0 deg as the reference level. The appropriate FORTRAN output statement should read "WRITE(IUNIT,*) alpha, 20-log10 (p(sub rms)(alpha)/p(sub rms)(0))".

Envia, Edmane↗
NASA NTRS1979

Utility of a finite element solution algorithm for initial-value problems

The Galerkin criterion within a finite element Weighted Residuals formulation is employed to establish an implicit solution algorithm for an initial-value partial differential equation. Numerical solutions of a transient parabolic and a hyperbolic equation, obtained using linear, quadratic and two cubic finite element basis functions, are employed to quantize accuracy and confirm and refine theoretical convergence rate estimates. The linear basis algorithm for the hyperbolic equation displays excellent accuracy on a coarse computational grid and a high-order convergence rate with discretization refinement. Good accuracy and a strong convergence rate in surface flux are determined for a nonhomogeneous Neumann boundary constraint applied to a parabolic equation. The results amply demonstrate the impact of the nondiagonal finite element initial-value matrix structure on solution accuracy and/or convergence rate.

Baker, A. J.↗
NASA NTRS1988

Boundary element method for 3-D cracks in a plate

Fundamental solutions which automatically satisfy boundary conditions at the interfaces of an elastic plate perfectly bonded to two elastic halfspaces are implemented in a three-dimensional BEM for crack problems. The BEM features a new integration scheme for highly singular kernels. The capability is achieved through a part analytic and part numerical integration procedure, such that the analytic part of the integration is similar for all slip/opening variations. Part-through elliptic cracks in an elastic plate with traction-free surfaces are analyzed and the SIF values along the crack front are found to compare favorably with the numerical SIF results of Raju and Newman (1979).

Fares, N.↗
NASA NTRS2007

Convergent Close-Coupling Approach to Electron-Atom Collisions

It was with great pleasure and honour to accept the invitation to make a presentation at the symposium celebrating the life-long work of Aaron Temkin and Richard Drachman. The work of Aaron Temkin was particularly influential on our own during the development of the CCC method for electron-atom collisions. There are a number of key problems that need to be dealt with when developing a general computational approach to such collisions. Traditionally, the electron energy range was subdivided into the low, intermediate, and high energies. At the low energies only a finite number of channels are open and variational or close-coupling techniques could be used to obtain accurate results. At high energies an infinite number of discrete channels and the target continuum are open, but perturbative techniques are able to yield accurate results. However, at the intermediate energies perturbative techniques fail and computational approaches need to be found for treating the infinite number of open channels. In addition, there are also problems associated with the identical nature of electrons and the difficulty of implementing the boundary conditions for ionization processes. The beauty of the Temkin-Poet model of electron-hydrogen scattering is that it simplifies the full computational problem by neglecting any non-zero orbital angular momenta in the partial-wave expansion, without loosing the complexity associated with the above-mentioned problems. The unique nature of the problem allowed for accurate solution leading to benchmark results which could then be used to test the much more general approaches to electron-atom collision problems. The immense value of the Temkin-Poet model is readily summarised by the fact that the initial papers of Temkin and Poet have been collectively cited around 250 times to date and are still being cited in present times. Many of the citations came from our own work during the course of the development of the CCC method, which we now describe.

Bray, Igor↗
NASA NTRS1983

Modelling atmospheric aerosol backscatter at CO2 laser wavelengths. I - Aerosol properties, modeling techniques, and associated problems. II - Modeled values in the atmosphere. III - Effects of changes in wavelength and ambient conditions

The various methods of calculating the atmospheric aerosol backscattering function, beta(CO2), both from measured aerosol characteristics and from optical measurements made at other wavelengths, are discussed in detail, with limits placed on their accuracy. The most significant factor in determining beta(CO2) is found to be the aerosol size distribution and concentration; this should be known accurately for particle radii up to at least 1 micron for stratospheric particles and 5 microns for tropospheric particles. Results are then presented from the modeling of the aerosol backscattering function at a wavelength of 10.6 microns in the lowest 20 km of the atmosphere. It is found that beta(CO2) varies from 10 to the -6th per m per sr in the planetary boundary layer to less than 10 to the -11th per m per sr in the stratosphere. It is next shown that, with the exception of (NH4)2SO4-containing aerosols, whose size distributions have relatively large numbers of small particles, the variation of backscattering with CO2 wavelength is less than a factor of approximately 3. For such (NH4)2SO4 aerosol distributions, however, the variation of backscatter function with CO2 wavelengths between 9.1 and 11.1 microns may reach one order of magnitude.

Kent, G.S.↗
NASA NTRS2012

Anomaly Detection in Test Equipment via Sliding Mode Observers

Nonlinear observers were originally developed based on the ideas of variable structure control, and for the purpose of detecting disturbances in complex systems. In this anomaly detection application, these observers were designed for estimating the distributed state of fluid flow in a pipe described by a class of advection equations. The observer algorithm uses collected data in a piping system to estimate the distributed system state (pressure and velocity along a pipe containing liquid gas propellant flow) using only boundary measurements. These estimates are then used to further estimate and localize possible anomalies such as leaks or foreign objects, and instrumentation metering problems such as incorrect flow meter orifice plate size. The observer algorithm has the following parts: a mathematical model of the fluid flow, observer control algorithm, and an anomaly identification algorithm. The main functional operation of the algorithm is in creating the sliding mode in the observer system implemented as software. Once the sliding mode starts in the system, the equivalent value of the discontinuous function in sliding mode can be obtained by filtering out the high-frequency chattering component. In control theory, "observers" are dynamic algorithms for the online estimation of the current state of a dynamic system by measurements of an output of the system. Classical linear observers can provide optimal estimates of a system state in case of uncertainty modeled by white noise. For nonlinear cases, the theory of nonlinear observers has been developed and its success is mainly due to the sliding mode approach. Using the mathematical theory of variable structure systems with sliding modes, the observer algorithm is designed in such a way that it steers the output of the model to the output of the system obtained via a variety of sensors, in spite of possible mismatches between the assumed model and actual system. The unique properties of sliding mode control allow not only control of the model internal states to the states of the real-life system, but also identification of the disturbance or anomaly that may occur.

Solano, Wanda M.↗
NASA NTRSDate not supplied

Delta Method Application for the Correlation of IR Detector Thermal Parasitic Loads with Statistically Accurate Results

The Linear Etalon Imaging Spectral Array (LEISA) is a cryogenic IR detector that is passively cooled below 110K. LEISA is part of the L’Ralph instrument, which is going on the Lucy mission, with planned flyby’s of the Trojan asteroids around 5.5 AU. One of the primary thermal challenges is to successfully quantify each significant parasitic heat flow value with statistical certainty.Simulation of flight-like thermal environments and correlation of thermal models is an important part of this challenge, which can be critical for thermally sensitive systems such as cryogenic detectors that are passively cooled. It is common practice during TVAC to achieve several thermal balance points and then to correlate several design variables in order to achieve the closest match to all balance points tested. This approach lends itself to high levels of uncertainty in the correlated values due to the cumulative effect of absolute temperature uncertainty and differences between the sensors used for the measurement. Accuracy can be improved with the use of more accurate sensors, by calibrating all sensors with respect to each other, and by using Zero-Q methodologies, all of which add cost and complexity to the test.The Delta Method uses the difference in temperature caused by the changes in the boundary conditions to calculate heat flows and/or thermal resistances across temperature differentials. This eliminates the problem of high uncertainties with absolute temperature measurements, and yields statistically accurate results down to the sensor resolution of 1mK at < 100K, by taking the sensor’s calibration errors out of the equation. During testing, this method relies on parametric variations around an initial balance point, by varying the boundary conditions of the different parasitic heat sources. The Delta Method, as well as its application on the LEISA thermal test, will be discussed in detail.

Daniel G. Bae↗
NASA NTRS1999

Experimental Investigation of Pool Boiling Heat Transfer Enhancement in Microgravity in the Presence of Electric Fields

In boiling high heat fluxes are possible driven by relatively small temperature differences, which make its use increasingly attractive in aerospace applications. The objective of the research is to develop ways to overcome specific problems associated with boiling in the low gravity environment by substituting the buoyancy force with the electric force to enhance bubble removal from the heated surface. Previous studies indicate that in terrestrial applications nucleate boiling heat transfer can be increased by a factor of 50, as compared to values obtained for the same system without electric fields. The goal of our research is to experimentally explore the mechanisms responsible for EHD heat transfer enhancement in boiling in low gravity conditions, by visualizing the temperature distributions in the vicinity of the heated surface and around the bubble during boiling using real-time holographic interferometry (HI) combined with high-speed cinematography. In the first phase of the project the influence of the electric field on a single bubble is investigated. Pool boiling is simulated by injecting a single bubble through a nozzle into the subcooled liquid or into the thermal boundary layer developed along the flat heater surface. Since the exact location of bubble formation is known, the optical equipment can be aligned and focused accurately, which is an essential requirement for precision measurements of bubble shape, size and deformation, as well as the visualization of temperature fields by HI. The size of the bubble and the frequency of bubble departure can be controlled by suitable selection of nozzle diameter and mass flow rate of vapor. In this approach effects due to the presence of the electric field can be separated from effects caused by the temperature gradients in the thermal boundary layer. The influence of the thermal boundary layer can be investigated after activating the heater at a later stage of the research. For the visualization experiments a test cell was developed. All four vertical walls of the test cell are transparent, and they allow transillumination with laser light for visualization experiments by HI. The bottom electrode is a copper cylinder, which is electrically grounded. The copper block is heated with a resistive heater and it is equipped with 6 thermocouples that provide reference temperatures for the measurements with HI. The top electrode is a mesh electrode. Bubbles are injected with a syringe into the test cell through the bottom electrode. The working fluids presently used in the interferometric visualization experiments, water and PF 5052, satisfy requirements regarding thermophysical, optical and electrical properties. A 30kV power supply equipped with a voltmeter allows to apply the electric field to the electrodes during the experiments. The magnitude of the applied voltage can be adjusted either manually or through the LabVIEW data acquisition and control system connected to a PC. Temperatures of the heated block are recorded using type-T thermocouples, whose output is read by a data acquisition system. Images of the bubbles are recorded with 35mm photographic and 16mm high-speed cameras, scanned and analyzed using various software packages. Visualized temperature fields HI allows the visualization of temperature fields in the vicinity of bubbles during boiling in the form of fringes. Typical visualized temperature distributions around the air bubbles injected into the thermal boundary layer in PF5052 are shown. The temperature of the heated surface is 35 C. The temperature difference for a pair of fringes is approximately 0.05 C. The heat flux applied to the bottom surface is moderate, and the fringe patterns are regular. In the image a bubble penetrating the thermal boundary layer is visible. Because of the axial symmetry of the problem, simplified reconstruction techniques can be applied to recover the temperature field. The thermal plume developing above the heated surface for more intensive heating is shown. The temperature distribution in the liquid is clearly 3D, and tomographic techniques have to be applied to recover the temperature distribution in such a physical situation. A sequence of interferometric images showing the temperature distribution around the rising bubble, recorded with a high-speed camera is shown. Again, the temperature distribution is 3D, and a more complex approach to the evaluation, the tomographic reconstruction has to be taken. Measurement of the temperature distribution from the fringe pattern temperature distributions that yield important information regarding heat transfer are determined. Two algorithms that allow the quantitative evaluation of interferometric fringe patterns and the reconstruction of temperature fields during boiling have been developed at the Heat Transfer Laboratory of the Johns Hopkins University. In the first algorithm the bubble is assumed to be axially symmetrical, which significantly reduces the computational effort for quantifying the temperature distribution around the bubble. For this purpose the thermal boundary layer around the bubble is divided into equidistant concentric shells, and the refractive index is assumed to be constant in each of the shells. Since large temperature gradients are expected in the vicinity of the bubble during boiling, the deflection of the light beam cannot be neglected in boiling experiments. Since the exit angle of the light beam is known, this allows to account for the deflections and phase shifts outside the boundary layer (in the bulk fluid and in the windows of the test cell). Three dimensional temperature distributions in the vicinity of the bubble are reconstructed using tomographic techniques. In tomography, the measurement volume is sliced into 2D planes. In the present study these planes are parallel to the heated surface. The objective is to determine the values of the field parameter of interest in form of the field function in these 2D planes. The field parameter is the change of the refractive index of the liquid in the measurement volume caused by temperature changes. By superimposing data for many 2D planes recorded at the same time instant, the 3D temperature distribution in the measurement volume is recovered.

Herman, Cila↗
NASA NTRS1984

A model for the estimation of the surface fluxes of momentum, heat and moisture of the cloud topped marine atmospheric boundary layer from satellite measurable parameters

A model is developed for the estimation of the surface fluxes of momentum, heat, and moisture of the cloud topped marine atmospheric boundary layer by use of satellite remotely sensed parameters. The parameters chosen for the problem are the integrated liquid water content, q sub li, the integrated water vapor content, q sub vi, the cloud top temperature, and either a measure of the 10 meter neutral wind speed or the friction velocity at the surface. Under the assumption of a horizontally homogeneous, well-mixed boundary layer, the model calculates the equivalent potential temperature and total water profiles of the boundary layer along with the boundary layer height from inputs of q sub li, q sub vi, and cloud top temperature. These values, along with the 10m neutral wind speed or friction velocity and the sea surface temperature are then used to estimate the surface fluxes. The development of a scheme to parameterize the integrated water vapor outside of the boundary layer for the cases of cold air outbreak and California coastal stratus is presented.

Allison, D. E.↗
NASA NTRS1993

Noniterative implicit method for tracking particles in mixed Lagrangian-Eulerian formulations

The existing implicit methods for the current initial value problems (IVPs) concerning particle-laden flows are complicated and iterative in nature. This paper presents a noniterative implicit method which can be used with pressure-based as well as with density-based algorithms. The method is illustrated by analyzing a dilute dispersion of noninteracting solid particles in an isothermal flow in a passage bounded by one straight wall and one wavy wall, in which all particles are spherical and have a finite velociy relative to the continuum phase at the inflow boundary.

Shih, T. I.-P.↗
NASA NTRS2005

Updated Chemical Kinetics and Sensitivity Analysis Code

An updated version of the General Chemical Kinetics and Sensitivity Analysis (LSENS) computer code has become available. A prior version of LSENS was described in "Program Helps to Determine Chemical-Reaction Mechanisms" (LEW-15758), NASA Tech Briefs, Vol. 19, No. 5 (May 1995), page 66. To recapitulate: LSENS solves complex, homogeneous, gas-phase, chemical-kinetics problems (e.g., combustion of fuels) that are represented by sets of many coupled, nonlinear, first-order ordinary differential equations. LSENS has been designed for flexibility, convenience, and computational efficiency. The present version of LSENS incorporates mathematical models for (1) a static system; (2) steady, one-dimensional inviscid flow; (3) reaction behind an incident shock wave, including boundary layer correction; (4) a perfectly stirred reactor; and (5) a perfectly stirred reactor followed by a plug-flow reactor. In addition, LSENS can compute equilibrium properties for the following assigned states: enthalpy and pressure, temperature and pressure, internal energy and volume, and temperature and volume. For static and one-dimensional-flow problems, including those behind an incident shock wave and following a perfectly stirred reactor calculation, LSENS can compute sensitivity coefficients of dependent variables and their derivatives, with respect to the initial values of dependent variables and/or the rate-coefficient parameters of the chemical reactions.

Radhakrishnan, Krishnan↗
NASA NTRS1992

Stability investigations of airfoil flow by global analysis

As the result of global, non-parallel flow stability analysis the single value of the disturbance growth-rate and respective frequency is obtained. This complex value characterizes the stability of the whole flow configuration and is not referred to any particular flow pattern. The global analysis assures that all the flow elements (wake, boundary and shear layer) are taken into account. The physical phenomena connected with the wake instability are properly reproduced by the global analysis. This enhances the investigations of instability of any 2-D flows, including ones in which the boundary layer instability effects are known to be of dominating importance. Assuming fully 2-D disturbance form, the global linear stability problem is formulated. The system of partial differential equations is solved for the eigenvalues and eigenvectors. The equations, written in the pure stream function formulation, are discretized via FDM using a curvilinear coordinate system. The complex eigenvalues and corresponding eigenvectors are evaluated by an iterative method. The investigations performed for various Reynolds numbers emphasize that the wake instability develops into the Karman vortex street. This phenomenon is shown to be connected with the first mode obtained from the non-parallel flow stability analysis. The higher modes are reflecting different physical phenomena as for example Tollmien-Schlichting waves, originating in the boundary layer and having the tendency to emerge as instabilities for the growing Reynolds number. The investigations are carried out for a circular cylinder, oblong ellipsis and airfoil. It is shown that the onset of the wake instability, the waves in the boundary layer, the shear layer instability are different solutions of the same eigenvalue problem, formulated using the non-parallel theory. The analysis offers large potential possibilities as the generalization of methods used till now for the stability analysis.

Morzynski, Marek↗