Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “high order”

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 577 records · Page 32

Microgravity Processing and Photonic Applications of Organic and Polymeric Materials

Some of the primary purposes of this work are to study important technologies, particularly involving thin films, relevant to organic and polymeric materials for improving applicability to optical circuitry and devices and to assess the contribution of convection on film quality in unit and microgravity environments. Among the most important materials processing techniques of interest in this work are solution-based and by physical vapor transport, both having proven gravitational and acceleration dependence. In particular, PolyDiAcetylenes (PDA's) and PhthaloCyanines (Pc's) are excellent NonLinear Optical (NLO) materials with the promise of significantly improved NLO properties through order and film quality enhancements possible through microgravity processing. Our approach is to focus research on integrated optical circuits and optoelectronic devices relevant to solution-based and vapor processes of interest in the Space Sciences Laboratory at the Marshall Space Flight Center (MSFC). Modification of organic materials is an important aspect of achieving more highly ordered structures in conjunction with microgravity processing. Parallel activities include characterization of materials for particular NLO properties and determination of appropriation device designs consistent with selected applications. One result of this work is the determination, theoretically, that buoyancy-driven convection occurs at low pressures in an ideal gas in a thermalgradient from source to sink. Subsequent experiment supports the theory. We have also determined theoretically that buoyancy-driven convection occurs during photodeposition of PDA, an MSFC-patented process for fabricating complex circuits, which is also supported by experiment. Finally, the discovery of intrinsic optical bistability in metal-free Pc films enables the possibility of the development of logic gate technology on the basis of these materials.

Frazier, Donald O.↗

Order reduction for models of space structures using modal cost analysis

Modal cost analysis furnishes a promising methodology for developing dynamical models of space structures for use in control systems analysis. Economy and accuracy can be attained by only retaining vibration modes that contribute significantly to an appropriately defined cost function. Expressions for modal costs (especially simple for 'lightly damped' structures) are derived for attitude control, vibration suppression, and shape control. These techniques are illustrated through application to a high-order finite element model of a large platform-type structure.

Skelton, R. E.↗

Propagation of Computational Uncertainty Using the Modern Design of Experiments

This paper describes the use of formally designed experiments to aid in the error analysis of a computational experiment. A method is described by which the underlying code is approximated with relatively low-order polynomial graduating functions represented by truncated Taylor series approximations to the true underlying response function. A resource-minimal approach is outlined by which such graduating functions can be estimated from a minimum number of case runs of the underlying computational code. Certain practical considerations are discussed, including ways and means of coping with high-order response functions. The distributional properties of prediction residuals are presented and discussed. A practical method is presented for quantifying that component of the prediction uncertainty of a computational code that can be attributed to imperfect knowledge of independent variable levels. This method is illustrated with a recent assessment of uncertainty in computational estimates of Space Shuttle thermal and structural reentry loads attributable to ice and foam debris impact on ascent.

DeLoach, Richard↗

A Phase-Shifting Zernike Wavefront Sensor for the Palomar P3K Adaptive Optics System

A phase-shifting Zernike wavefront sensor has distinct advantages over other types of wavefront sensors. Chief among them are: 1) improved sensitivity to low-order aberrations and 2) efficient use of photons (hence reduced sensitivity to photon noise). We are in the process of deploying a phase-shifting Zernike wavefront sensor to be used with the realtime adaptive optics system for Palomar. Here we present the current state of the Zernike wavefront sensor to be integrated into the high-order adaptive optics system at Mount Palomar's Hale Telescope.

adaptive optics↗

Tunable Filter Made From Three Coupled WGM Resonators

A tunable third-order band-pass optical filter has been constructed as an assembly of three coupled, tunable, whispering-gallery-mode resonators similar to the one described in Whispering-Gallery-Mode Tunable Narrow-Band-Pass Filter (NPO-30896), NASA Tech Briefs, Vol. 28, No. 4 (April 2004), page 5a. This filter offers a combination of four characteristics that are desirable for potential applications in photonics: (1) wide real-time tunability accompanied by a high-order filter function, (2) narrowness of the passband, (3) relatively low loss between input and output coupling optical fibers, and (4) a sparse spectrum. In contrast, prior tunable band-pass optical filters have exhibited, at most, two of these four characteristics. As described in several prior NASA Tech Briefs articles, a whispering-gallery-mode (WGM) resonator is a spheroidal, disklike, or toroidal body made of a highly transparent material. It is so named because it is designed to exploit whispering-gallery electromagnetic modes, which are waveguide modes that propagate circumferentially and are concentrated in a narrow toroidal region centered on the equatorial plane and located near the outermost edge. Figure 1 depicts the optical layout of the present filter comprising an assembly of three coupled, tunable WGM resonators. Each WGM resonator is made from a disk of Z-cut LiNbO3 of 3.3-mm diameter and 50-m thickness. The perimeter of the disk is polished and rounded to a radius of curvature of 40 microns. The free spectral range of each WGM resonator is about 13.3 GHz. Gold coats on the flat faces of the disk serve as electrodes for exploiting the electro-optical effect in LiNbO3 for tuning. There is no metal coat on the rounded perimeter region, where the whispering-gallery modes propagate. Light is coupled from an input optical fiber into the whispering-gallery-modes of the first WGM resonator by means of a diamond prism. Another diamond prism is used to couple light from the whispering-gallery-modes of the third WGM resonator to an output optical fiber. The filter operates at a nominal wavelength of 1,550 nm and can be tuned over a frequency range of plus or minus 12 GHz by applying a potential in the range of plus or minus 150 V to the electrodes. The insertion loss (the loss between the input and output coupling optical fibers) was found to be repeatable at 6 dB. The resonance quality factor (Q) of the main sequence of resonator modes was found to be 5 x 10(exp 6), which corresponds to a bandwidth of 30 MHz. The filter can be shifted from one operating frequency to another within a tuning time less than or equal to 30 micro seconds. The transmission curve of the filter at frequencies near the middle of the passband closely approximates a theoretical third-order Butterworth filter profile, as shown in Figure 2.

Savchenkov, Anatoliy↗

A closer look at the induced drag of crescent-shaped wings

Recent interest in the induced drag characteristics of crescent-shaped wings has led to a closer look at the methods used for determination of induced drag from computational aerodynamic methods. Induced drag may be computed by integration of surface pressure, or by evaluation of a contour integral in the Trefftz plane. A high-order panel method was used to study the induced drag of crescent and elliptical wings using both techniques. Induced drag computations using surface-pressure integration were strongly affected by panel density and angle of attack. Accurate drag computations for the crescent wing were obtained only when the spanwise as well as chordwise panel density was extremely high. Trefftz-plane results for the two wing planforms do not show this sensitivity to panel density or angle of attack. Span efficiencies of 0.994 for the crescent wing and 0.987 for the elliptical wing were computed by the Trefftz-plane technique. Substitution of a force-free, rolled-up wake geometry on the crescent wing did not change the pressure-integrated drag significantly. The slightly higher span efficiency of the crescent wing is attributed to a more nearly elliptical spanwise lift distribution. The chord distribution of the elliptical wing was modified to produce an elliptical span-loading on a wing with an unswept quarter-chord line. This wing demonstrated a span efficiency equal to that of the crescent wing.

Smith, Stephen C.↗

Accurate Monotonicity - Preserving Schemes With Runge-Kutta Time Stepping

A new class of high-order monotonicity-preserving schemes for the numerical solution of conservation laws is presented. The interface value in these schemes is obtained by limiting a higher-order polynominal reconstruction. The limiting is designed to preserve accuracy near extrema and to work well with Runge-Kutta time stepping. Computational efficiency is enhanced by a simple test that determines whether the limiting procedure is needed. For linear advection in one dimension, these schemes are shown as well as the Euler equations also confirm their high accuracy, good shock resolution, and computational efficiency.

Suresh, A.↗

High-Fidelity Computational Aerodynamics of the Elytron 4S UAV

High-fidelity Computational Fluid Dynamics (CFD) have been carried out for the Elytron 4S Unmanned Aerial Vehicle (UAV), also known as the converticopter "proto12". It is the scaled wind tunnel model of the Elytron 4S, an Urban Air Mobility (UAM) concept, a tilt-wing, box-wing rotorcraft capable of Vertical Take-Off and Landing (VTOL). The three-dimensional unsteady Navier-Stokes equations are solved on overset grids employing high-order accurate schemes, dual-time stepping, and a hybrid turbulence model using NASA's CFD code OVERFLOW. The Elytron 4S UAV has been simulated in airplane mode and in helicopter mode.

High-Fidelity↗

Solving Nonlinear Euler Equations with Arbitrary Accuracy

A computer program that efficiently solves the time-dependent, nonlinear Euler equations in two dimensions to an arbitrarily high order of accuracy has been developed. The program implements a modified form of a prior arbitrary- accuracy simulation algorithm that is a member of the class of algorithms known in the art as modified expansion solution approximation (MESA) schemes. Whereas millions of lines of code were needed to implement the prior MESA algorithm, it is possible to implement the present MESA algorithm by use of one or a few pages of Fortran code, the exact amount depending on the specific application. The ability to solve the Euler equations to arbitrarily high accuracy is especially beneficial in simulations of aeroacoustic effects in settings in which fully nonlinear behavior is expected - for example, at stagnation points of fan blades, where linearizing assumptions break down. At these locations, it is necessary to solve the full nonlinear Euler equations, and inasmuch as the acoustical energy is of the order of 4 to 5 orders of magnitude below that of the mean flow, it is necessary to achieve an overall fractional error of less than 10-6 in order to faithfully simulate entropy, vortical, and acoustical waves.

Dyson, Rodger W.↗

Direct Numerical Simulations of Boundary Layer Transition on a Flat Plate

In recent years the techniques of computational fluid dynamics (CFD) have been used to compute flows associated with geometrically complex configurations. However, success in terms of accuracy and reliability has been limited to cases where the effects of turbulence and transition could be modeled in a straightforward manner. Even in simple flows, the accurate computation of skin friction and heat transfer using existing turbulence models has proved to be a difficult task, one that has required extensive fine-tuning of the turbulence models used. In more complex flows (for example, in turbomachinery flows in which vortices and wakes impinge on airfoil surfaces causing periodic transitions from laminar to turbulent flow) the development of a model that accounts for all scales of turbulence and predicts the onset of transition is an extremely difficult task. Fortunately, current trends in computing suggest that it may be possible to perform direct simulations of turbulence and transition at moderate Reynolds numbers in some complex cases in the near future. This presentation will focus on direct simulations of transition and turbulence using high-order accurate finite-difference methods. The advantage of the finite-difference approach over spectral methods is that complex geometries can be treated in a straightforward manner. Additionally, finite-difference techniques are the prevailing methods in existing application codes. An application of accurate finite-difference methods to direct simulations of transition and turbulence in a spatially evolving boundary layer subjected to high levels of freestream turbulence will be presented.

Rai, Man Mohan↗

3D Higher Order Modeling in the BEM/FEM Hybrid Formulation

Higher order divergence- and curl-conforming bases have been shown to provide significant benefits, in both convergence rate and accuracy, in the 2D hybrid finite element/boundary element formulation (P. Fink and D. Wilton, National Radio Science Meeting, Boulder, CO, Jan. 2000). A critical issue in achieving the potential for accuracy of the approach is the accurate evaluation of all matrix elements. These involve products of high order polynomials and, in some instances, singular Green's functions. In the 2D formulation, the use of a generalized Gaussian quadrature method was found to greatly facilitate the computation and to improve the accuracy of the boundary integral equation self-terms. In this paper, a 3D, hybrid electric field formulation employing higher order bases and higher order elements is presented. The improvements in convergence rate and accuracy, compared to those resulting from lower order modeling, are established. Techniques developed to facilitate the computation of the boundary integral self-terms are also shown to improve the accuracy of these terms. Finally, simple preconditioning techniques are used in conjunction with iterative solution procedures to solve the resulting linear system efficiently. In order to handle the boundary integral singularities in the 3D formulation, the parent element- either a triangle or rectangle-is subdivided into a set of sub-triangles with a common vertex at the singularity. The contribution to the integral from each of the sub-triangles is computed using the Duffy transformation to remove the singularity. This method is shown to greatly facilitate t'pe self-term computation when the bases are of higher order. In addition, the sub-triangles can be further divided to achieve near arbitrary accuracy in the self-term computation. An efficient method for subdividing the parent element is presented. The accuracy obtained using higher order bases is compared to that obtained using lower order bases when the number of unknowns is approximately equal. Also, convergence rates obtained using higher order bases are compared to those obtained with lower order bases for selected sample

Fink, P. W.↗

Sigma Filter

The GLA Fourth-Order model is needed to smooth the topography. This is to remove the Gibbs phenomenon. The Gibbs phenomenon occurs whenever we truncate a Fourier Series. The Sigma factors were introduced to reduce the Gibbs phenomenon. It is found that the smooth Fourier series is nothing but the original Fourier series with its coefficients multiplied by corresponding sigma factors. This operator can be applied many times to obtain high order sigma filtered field and is easily applicable using FFT. It is found that this filter is beneficial in deriving the topography.

Balgovind, R. C.↗

Direct simulations of turbulent flow using finite-difference schemes

A high-order accurate finite-difference approach is presented for calculating incompressible turbulent flow. The methods used include a kinetic energy conserving central difference scheme and an upwind difference scheme. The methods are evaluated in test cases for the evolution of small-amplitude disturbances and fully developed turbulent channel flow. It is suggested that the finite-difference approach can be applied to complex geometries more easilty than highly accurate spectral methods. It is concluded that the upwind scheme is a good candidate for direct simulations of turbulent flows over complex geometries.

Rai, Man Mohan↗

Improving Pixel Level Cloud Optical Property Retrieval using Monte Carlo Simulations

The accurate pixel-by-pixel retrieval of cloud optical properties from space is influenced by radiative smoothing due to high order photon scattering and radiative roughening due to low order scattering events. Both are caused by cloud heterogeneity and the three-dimensional nature of radiative transfer and can be studied with the aid of computer simulations. We use Monte Carlo simulations on variable 1-D and 2-D model cloud fields to seek for dependencies of smoothing and roughening phenomena on single scattering albedo, solar zenith angle, and cloud characteristics. The results are discussed in the context of high resolution satellite (such as Landsat) retrieval applications. The current work extends the investigation on the inverse NIPA (Non-local Independent Pixel Approximation) as a tool for removing smoothing and improving retrievals of cloud optical depth. This is accomplished by: (1) Delineating the limits of NIPA applicability; (2) Exploring NIPA parameter dependences on cloud macrostructural features, such as mean cloud optical depth and geometrical thickness, degree of extinction and cloud top height variability. We also compare parameter values from empirical and theoretical considerations; (3) Examining the differences between applying NIPA on radiation quantities vs direct application on optical properties; (4) Studying the radiation budget importance of the NIPA corrections as a function of scale. Finally, we discuss fundamental adjustments that need to be considered for successful radiance inversion at non-conservative wavelengths and oblique Sun angles. These adjustments are necessary to remove roughening signatures which become more prominent with increasing absorption and solar zenith angle.

Oreopoulos, Lazaros↗

Constrained Field Construction Using Bernstein Polynomial Derived Operators

This method recasts the mass-conservation equation as a Sylvester equation employing high-order derivative operators derived from modified Bernstein polynomial expansions. Given prescribed velocity fields, the algorithm yields a constrained, least-squared solution for the associated density field. To demonstrate the practical utility of this methodology, it is applied to a computationally-derived two-dimensional isolator dataset. This reconstruction method is envisioned to be used in conjunction with diagnostic techniques to aid in the quantification of isolator flow fields, since obtaining highly characterized datasets within this engine component is exceedingly difficult.

Isolator↗

Higher Order Modulation Intersymbol Interference Caused by Traveling-wave Tube Amplifiers

For the first time, a time-dependent, physics-based computational model has been used to provide a direct description of the effects of the traveling wave tube amplifier (TWTA) on modulated digital signals. The TWT model comprehensively takes into account the effects of frequency dependent AM/AM and AM/PM conversion; gain and phase ripple; drive-induced oscillations; harmonic generation; intermodulation products; and backward waves, Thus, signal integrity can be investigated in the presence of these sources of potential distortion as a function of the physical geometry and operating characteristics of the high power amplifier and the operational digital signal. This method promises superior predictive fidelity compared to methods using TWT models based on swept-amplitude and/or swept-frequency data. First, the TWT model using the three dimensional (3D) electromagnetic code MAFIA is presented. Then, this comprehensive model is used to investigate approximations made in conventional TWT black-box models used in communication system level simulations, To quantitatively demonstrate the effects these approximations have on digital signal performance predictions, including intersymbol interference (ISI), the MAFIA results are compared to the system level analysis tool, Signal Processing, Workstation (SPW), using high order modulation schemes including 16 and 64-QAM.

Kory, Carol L.↗

Convergence Analysis of Turbulent Flow Solutions

Data from the "Turbulence Modeling Resource" website for turbulent flow over an NACA-0012 airfoil is analyzed to determine the convergence behavior of three second-order CFD (Computational Fluid Dynamics) codes: CFL3D (Computational Fluids Lab 3 Dimensional flow solver), FUN3D (Fully Unstructured Navier-stokes flow solver), and TAU (German Aerospace Center (DLR) 2 dimensional code for unstructured hybrid grids solving the Reynolds-Averaged Navier-Stokes equations or the Euler equations). The convergence of both integrated properties and pointwise data are examined. Several different methods for estimating errors and computing convergence rates are compared. A high-order extension to the Richardson extrapolation is developed that improves the accuracy of the mesh limit values and provides a quantitative estimate of the threshold of the asymptotic regime. The coefficient of total drag exhibits second-order convergence for all three codes, and convergence is monotone over a sequence of 7 grids. Other force coefficients are not so well behaved. The convergence rates of the viscous component of drag on the three nest grids ranges from 3:0 for CFL3D to 1:0 for FUN3D. The three codes are converging to similar but not identical solutions. The largest differences between the codes are in the coefficient of lift for which the difference between CFL3D and FUN3D is greater than 10 (sup minus 4). The best agreement occurs in the viscous component of drag, which is the only force component for which all three codes are converging towards each other at a rate of second-order. The agreement between the two unstructured grid codes is good with all properties except lift converging towards common values at a rate of second-order. No one code was universally better than the other. The TAU code has the lowest error in total drag, FUN3D has the lowest error in lift, and CFL3D has the lowest error in the viscous component of drag.

Atkins, Harold L.↗

Non-Gaussian approach for parametric random vibration of non-linear structures

The dynamic response of a nonlinear, single degree of freedom structural system subjected to a physically white noise parametric excitation is investigated. The Ito stochastic calculus is employed to derive a general differential equation for the moments of the response coordinates. The differential equations of moments of any order are found to be coupled with higher order moments. A non-Gaussian closure scheme is developed to truncate the moment equations up to fourth order. The statistical of the stationary response are computed numerically and compared with analytical solutions predicted by a Gaussian closure scheme and the stochastic averaging method. It is found that the computed results exhibit the jump phenomenon which is typical of the characteristics of deterministic nonlinear systems. In addition, the numerical algorithm leads to multiple solutions all of which give positive mean squares. However, two of these solutions are found to violate the properties of high order moments. One solution preserves the moments properties and demonstrates that the system achieves a stationary response.

Ibrahim, R. A.↗