Extension of equivalent angle-of-attack method for nonlinear flowfields
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
Numerical accuracy for the boundary and domain methods of the material derivative approach to shape design sensitivities is investigated through the use of mesh refinement. The results show that the domain method is generally more accurate than the boundary method, using the finite element technique. It is also shown that the domain method is equivalent, under certain assumptions, to the implicit differentiation approach not only theoretically but also numerically.
Nonlinear equations of motion of symmetrically laminated anisotropic plates are derived accounting for von Karman strains. The effect of transverse shear is included in the formulation. Using a single-mode Galerkin procedure, the nonlinear modal equation is obtained. The direct equivalent linearization method is employed for solution of this equation. The response to acoustic excitation of moderately thick composite panels is studied. Further, the effects of transverse shear on large deflection vibration of laminates under random excitation are studied. Mean-square deflections and mean-square inplane stresses are obtained for symmetric graphite-epoxy laminates. Using equilibrium equations and the continuity requirements, the mean-square transverse shear stresses are calculated. The results obtained will be useful in the sonic fatigue design of composite aircraft panels. The analysis is presented in detail for simply supported plates. The analogous equations for a clamped case are also obtained.
The backscattering from a circular disk is analyzed using the geometrical theory of diffraction. First-, second-, and third-order diffractions are included in the hard polarization analysis, while first-, second-, and third-order slope diffractions are included for soft polarization. Improvements in the prediction of the monostatic radar cross section over previous works are noted. For hard polarization, an excellent agreement is exhibited between experimental and theoretical results, while a very good agreement is noted for soft polarization. To further improve the soft polarization results for wide angles, a model for the creeping wave or circulating current on the edge of the disk is obtained and used to find an additional component of the backscattered field. The addition of this component significantly improves the results for wide angles, leading to excellent agreement for soft polarization also. An axial-caustic correction method using equivalent currents is also included in the analysis.
Nonlinear equations of motion of symmetrically laminated anisotropic plates are derived considering von Karman strains and transverse shear effects. Using a single-mode Galerkin procedure, the nonlinear modal equation is obtained. The direct equivalent linearization method is employed for solution of this equation. The effects of transverse shear on large deflection vibration of laminated plates with mixed boundary support conditions under random excitation are studied. Mean-square deflections and mean-square inplane strains are obtained for symmetric graphite-epoxy laminates. Using equilibrium equations and the continuity requirements, the mean-square transverse shear stresses are calculated. The results obtained will be useful in the sonic fatigue design of composite aircraft panels. The analysis is presented in detail for plates with two opposing edges simply supported and the remaining edges clamped.
By using the linear mode shapes of a thermally buckled structure, a governing system of finite element equations for large-deflection random vibrations is reduced to a system of nonlinear modal equations. This system is then solved using the method of equivalent linearization. The methodology developed in this study is applicable to the design of aerospace vehicles subjected to combined thermal/acoustic loads.
An efficient high-frequency ray analysis, which combines the geometrical theory of diffraction (GTD) with the equivalent current method, is used for evaluating modal reflection at the open ends of three-dimensional waveguides. After reviewing a method used for the two-dimensional case and introducing some extensions necessary for evaluating modal reflection of open-ended three-dimensional waveguides, the modal reflection coefficients are then given for open-ended circular and rectangular waveguides. Numerical results calculated using these expressions agree favorably with those obtained by a different analytical technique and by experiment.
The problem of electromagnetic coupling between two horns is of interest for the Microwave Reflectometer Ionization Sensor (MRIS) that will be used in the Aeroassist Flight Experiment (AFE). Laboratory measurements of mutual coupling between conical horns (using a flat metallic reflector to simulate a critically dense plasma outside) have shown a strong dependence on the finite dimensions of the shuttle tile over the apertures. Since both, the dielectric tile and the plasma outside the tile reflect microwaves, a study should be done to isolate the two mechanisms so that the MRIS reentry flight data can be interpreted correctly. Once the coupling due to the tile itself is determined then the location of the critial electron number density layers can be determined. As a first attempt to tackle this problem the Geometrical Theory of Diffraction was used to modify the existing solution to mutual coupling between apertures with infinite dielectric sheets. By using the equivalent current method, aperture theory to determine the radiated fields inside the dielectric tiles, and ray tracing the contributions to mutual coupling were determined. Results from two cases with different tile thicknesses have indicated that the main contribution to mutual coupling is due to diffraction from the bottom and top (back and front) wedges.
Several different high-frequency methods for modeling the radar cross sections (RCSs) of plate geometries are examined. The Method of Equivalent Currents and a numerically derived corner diffraction coefficient are used to model the RCS of a rectangular, perfectly conducting plate in nonprincipal planes. The Uniform Theory of Diffraction is used to model the RCS of a rectangular, perfectly conducting plate in principal planes. For the soft polarization case, first-order and slope-diffraction terms are included. For the hard polarization case, up to four orders of diffraction are included. Finally, the Uniform Theory of Diffraction for impedance wedges and the Impedance Boundary Condition are used to model the RCS of a coated, rectangular plate in principal planes. In most of the cases considered, comparisons are made between theoretical and experimental results.
Radar cross section (RCS) prediction of several rectangular plate geometries is discussed using high-frequency techniques such as the Uniform Theory of Diffraction (UTD) for perfectly conducting and impedance wedges and the Method of Equivalent Currents (MEC). Previous reports have presented detailed solutions to the principal-plane scattering by a perfectly conducting and a coated rectangular plate and nonprincipal-plane scattering by a perfectly conducting plate. These solutions are briefly reviewed and a modified model is presented for the coated plate. Theoretical and experimental data are presented for the perfectly conducting geometries. Agreement between theory and experiment is very good near and at normal incidence. In regions near and at grazing incidence, the disagreement between the data vary according to diffraction distances and angles involved. It is these areas of disagreement which are of extreme interest as an explanation for the disagreement will yield invaluable insight into scattering mechanisms which are not yet identified as major contributors near and at grazing incidence. Areas of disagreement between theory and experiment are identified and examined in an attempt to better understand and predict near-grazing incidence, grazing incidence, and nonprincipal-plane diffractions.
A high-accuracy numerical method for the solution of a 1D Schroedinger equation that is suitable for a diatomic molecule, obtained by combining a finite-difference method with iterative extrapolation to the limit, is presently shown to have several advantages over more conventional methods. Initial guesses for the term values are obviated, and implementation of the algorithm is straightforward. The method is both less sensitive to round-off error, and faster than conventional methods for equivalent accuracy. These advantages are illustrated through the solution of Schroedinger's equation for a Morse potential function suited for HCl and a numerically derived Rydberg-Klein-Rees potential function for the X 1Sigma(+) state of CO.
Cross-stream models are defined as solutions to the Orr-Sommerfeld equation which are propagating normal to the flow direction. These models are utilized as a basis for a Hilbert space to approximate the spectrum of the Orr-Sommerfeld equation with plane Poiseuille flow. The cross-stream basis leads to a standard eigenvalue problem for the frequencies of Poiseuille flow instability waves. The coefficient matrix in the eigenvalue problem is shown to be the sum of a real matrix and a negative-imaginary diagonal matrix which represents the frequencies of the cross-stream modes. The real coefficient matrix is shown to approach a Toeplitz matrix when the row and column indices are large. The Toeplitz matrix is diagonally dominant, and the diagonal elements vary inversely in magnitude with diagonal position. The Poiseuille flow eigenvalues are shown to lie within Gersgorin disks with radii bounded by the product of the average flow speed and the axial wavenumber. It is shown that the eigenvalues approach the Gersgorin disk centers when the mode index is large, so that the method may be used to compute spectra with an essentially unlimited number of elements. When the mode index is large, the real part of the eigenvalue is the product of the axial wavenumber and the average flow speed, and the imaginary part of the eigen value is identical to the corresponding cross-stream mode frequency. The cross-stream method is numerically well-conditioned in comparison to Chebyshev based methods, providing equivalent accuracy for small mode indices and superior accuracy for large indices.
The governing nonlinear equations of motion are derived for a rectangular isotropic plate subjected to a combined thermal-acoustic loading. The thermal loading is taken to be a steady state temperature distribution, and the acoustic loading is considered to be a stationary Gaussian pressure with zero mean and uniform magnitude and phase over the plate surface. With these types of loads, the thermal-acoustic analysis can be separated into two problems: (1) thermal post-buckling; and (2) random vibration of the thermally buckled plate. Since both analyses are nonlinear, iterative techniques are used to solve each. For the random vibration analysis, the method of equivalent linearization is used; both auto- and cross-correlation terms are included. Results obtained using the present formulation are compared with previous classical equivalent linearization results.
We have used recent observational data for the Cepheids in the rich, young LMC clusters NGC 1866 and NGC 2031 to constrain the cluster distances with the mass-equivalency (ME) method. The basis of this approach is to fix the cluster distances by requiring the Cepheid evolutionary and pulsational masses to be equal. Using evolutionary models incorporating a mild amount of core and envelope overshooting along with recent pulsational models, we derive distance moduli of 18.51 +/- 0.21 and 18.32 +/- 0.20 for NGC 1866 and NGC 2031, respectively. The quoted errors are dominated by the uncertainties in the heavy element abundances of the clusters (assumed to be 0.3 dex for both clusters), with a smaller contribution due to the apparently intrinsic spread in the masses of the Cepheids in each cluster. For the ME method, we find that Delta(m-M)(0)/Delta sub Z(1) = 0.69, where Z(1) = log (Z/0.016). This result implies that the cluster distances can be determined to better than +/- 5 percent if the cluster abundances can be measured to better than about +/- 0.15 dex.
Most computational fluid dynamics (CFD) schemes are not adequately accurate for solving aeroacoustics problems, which have wave amplitudes several orders of magnitude smaller yet with frequencies larger than the flow field variations generating the sound. Hence, a computational aeroacoustics (CAA) algorithm should have minimal dispersion and dissipation features. A dispersion relation preserving (DRP) scheme is, therefore, applied to solve the linearized Euler equations in order to simulate the propagation of three types of waves, namely: acoustic, vorticity, and entropy waves. The scheme is derived using an optimization procedure to ensure that the numerical derivatives preserve the wave number and angular frequency of the partial differential equations being discretized. Consequently, simulated waves propagate with the correct wave speeds and exhibit their appropriate properties. A set of radiation and outflow boundary conditions, compatible with the DRP scheme and derived from the asymptotic solutions of the governing equations, are also implemented. Numerical simulations are performed to test the effectiveness of the DRP scheme and its boundary conditions. The computed solutions are shown to agree favorably with the exact solutions. The major restriction appears to be that the dispersion relations can be preserved only for waves with wave lengths longer than four or five spacings. The boundary conditions are found to be transparent to the outgoing disturbances. However, when the disturbance source is placed closer to a boundary, small acoustic reflections start appearing. CAA generates enormous amounts of temporal data which needs to be reduced to understand the physical problem being simulated. Spectral analysis is one approach that helps us in extracting information which often can not be easily interpreted in the time domain. Thus, three different methods for the spectral analysis of numerically generated aeroacoustic data are studied. First, the capabilities of two traditional methods for spectral analysis, namely, the Blackman-Tukey method and periodogram method, are compared in estimating the spectra of a simple-periodic process. The periodogram is then applied to analyze transitory-deterministic processes. Finally, these two methods are compared with a more recent method, referred as the Weighted-Overlapped-Segment-Averaging (WOSA) method, in estimating the spectra of a chaotic (random-like) process. From the demonstrative case for the spectral analyses of data generated by simple-periodic process, the periodogram method is found to give a better estimate of the steep-sloped spectra than the Blackman-Tukey method. Also, for this problem, the Hanning window is found to perform better with the periodogram method than with the Blackman-Tukey method. Finally, for the spectral analysis of data generated by the chaotic process, the periodogram method does not perform well, whereas, the WOSA and Blackman-Tukey methods give equivalently good results.
A process for extracting long-equivalent wavelength interferometric information from a two-wavelength polychromatic or achromatic interferometer. The process comprises the steps of simultaneously recording a non-linear sum of two different frequency visible light interferograms on a high resolution film and then placing the developed film in an optical train for Fourier transformation, low pass spatial filtering and inverse transformation of the film image to produce low spatial frequency fringes corresponding to a long-equivalent wavelength interferogram. The recorded non-linear sum irradiance derived from the two-wavelength interferometer is obtained by controlling the exposure so that the average interferogram irradiance is set at either the noise level threshold or the saturation level threshold of the film.
The Fast Scattering Code (FSC) is a frequency domain noise prediction program developed at the NASA Langley Research Center (LaRC) to simulate the acoustic field produced by the interaction of known, time harmonic incident sound with bodies of arbitrary shape and surface impedance immersed in a potential flow. The code uses the equivalent source method (ESM) to solve an exterior 3-D Helmholtz boundary value problem (BVP) by expanding the scattered acoustic pressure field into a series of point sources distributed on a fictitious surface placed inside the actual scatterer. This work provides additional code validation studies and illustrates the range of code parameters that produce accurate results with minimal computational costs. Systematic noise prediction studies are presented in which monopole generated incident sound is scattered by simple geometric shapes - spheres (acoustically hard and soft surfaces), oblate spheroids, flat disk, and flat plates with various edge topologies. Comparisons between FSC simulations and analytical results and experimental data are presented.
The computational fluid dynamics (CFD) prediction workshops sponsored by the AIAA have created invaluable opportunities in which to discuss the predictive capabilities of CFD in areas in which it has struggled, e.g., cruise drag, high-lift, and sonic boom pre diction. While there are many factors that contribute to disagreement between simulated and experimental results, such as modeling or discretization error, quantifying the errors contained in a simulation is important for those who make decisions based on the computational results. The linearized error transport equations (ETE) combined with a truncation error estimate is a method to quantify one source of errors. The ETE are implemented with a complex-step method to provide an exact linearization with minimal source code modifications to CFD and multidisciplinary analysis methods. The equivalency of adjoint and linearized ETE functional error correction is demonstrated. Uniformly refined grids from a series of AIAA prediction workshops demonstrate the utility of ETE for multidisciplinary analysis with a connection between estimated discretization error and (resolved or under-resolved) flow features.