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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

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↗
NASA NTRS2006

Adaptive Nulling for the Terrestrial Planet Finder Interferometer

A description of adaptive nulling for Terrestrial Planet Finder Interferometer (TPFI) is presented. The topics include: 1) Nulling in TPF-I; 2) Why Do Adaptive Nulling; 3) Parallel High-Order Compensator Design; 4) Phase and Amplitude Control; 5) Development Activates; 6) Requirements; 7) Simplified Experimental Setup; 8) Intensity Correction; and 9) Intensity Dispersion Stability. A short summary is also given on adaptive nulling for the TPFI.

adaptive nulling↗
NASA NTRS2004

An Integrated Optimal Estimation Approach to Spitzer Space Telescope Focal Plane Survey

This paper discusses an accurate and efficient method for focal plane survey that was used for the Spitzer Space Telescope. The approach is based on using a high-order 37-state Instrument Pointing Frame (IPF) Kalman filter that combines both engineering parameters and science parameters into a single filter formulation. In this approach, engineering parameters such as pointing alignments, thermomechanical drift and gyro drifts are estimated along with science parameters such as plate scales and optical distortions. This integrated approach has many advantages compared to estimating the engineering and science parameters separately. The resulting focal plane survey approach is applicable to a diverse range of science instruments such as imaging cameras, spectroscopy slits, and scanning-type arrays alike. The paper will summarize results from applying the IPF Kalman Filter to calibrating the Spitzer Space Telescope focal plane, containing the MIPS, IRAC, and the IRS science Instrument arrays.

optimal estimations↗
NASA NTRS2007

Revisiting and Extending Interface Penalties for Multi-Domain Summation-by-Parts Operators

General interface coupling conditions are presented for multi-domain collocation methods, which satisfy the summation-by-parts (SBP) spatial discretization convention. The combined interior/interface operators are proven to be L2 stable, pointwise stable, and conservative, while maintaining the underlying accuracy of the interior SBP operator. The new interface conditions resemble (and were motivated by) those used in the discontinuous Galerkin finite element community, and maintain many of the same properties. Extensive validation studies are presented using two classes of high-order SBP operators: 1) central finite difference, and 2) Legendre spectral collocation.

Carpenter, Mark H.↗
NASA NTRS2007

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↗
NASA NTRS2008

Compact Dielectric-Rod White-Light Delay Lines

Optical delay lines of a proposed type would be made from rods of such dielectric materials as calcium fluoride, fused silica, or sapphire. These would offer advantages over prior optical delay lines, as summarized below. Optical delay lines are key components of opto-electronic microwave oscillators, narrow-band opto-electronic microwave filters, evanescent-field optical biochemical detectors, and some Fourier-Transform spectrum analyzers. Heretofore, optical delay lines used in such applications have been of two types: resonators and coiled long optical fibers, both of which have disadvantages: Resonators are compact, but excitation must be provided by narrow-band lasers. Wide-band (including noisy) laser light cannot be coupled efficiently to narrow-band resonators. When light is coupled into a narrowband resonator from a source of reasonably high power, a significant amount of optical energy circulates within the resonator, causing nonlinear loss and significant noise. Typically, a coil-type optical delay line is made of fused-silica fiber, which exhibits fundamental loss. To overcome the limit imposed by the optical loss in fused silica, it would be necessary to use fibers having crystalline cores. Although space is saved by winding fibers into coils, fiber-coil delay lines are still inconveniently bulky. The proposed compact dielectric-rod delay lines would exploit the special class of non-diffracting light beams that are denoted Bessel beams because their amplitudes are proportional to Bessel functions of the radii from their central axes. High-order Bessel beams can have large values of angular momentum. They can be generated with the help of whispering-gallery-mode optical resonators, as described, for example, in "Simplified Generation of High-Angular-Momentum Light Beams" (NPO-42965) NASA Tech Briefs, Vol. 31, No. 3 (March 2007), page 8a. In a delay line according to the proposal, the dielectric rod would be dimensioned to function as a multimode waveguide. Suitably chosen high-angular-momentum modes in such a waveguide exhibit low group velocity (hence, long delay) and no resonance. Such a delay line could perform well at any wavelength or range of wavelengths within the transparency wavelength band of the dielectric material, and the maximum possible group delay achievable through suitable design would be limited only by the optical loss in the rod material.

Maleki, Lute↗
NASA NTRS2004

Spitzer Instrument Pointing Frame (IPF) Kalman Filter Algorithm

This paper discusses the Spitzer Instrument Pointing Frame (IPF) Kalman Filter algorithm. The IPF Kalman filter is a high-order square-root iterated linearized Kalman filter, which is parametrized for calibrating the Spitzer Space Telescope focal plane and aligning the science instrument arrays with respect to the telescope boresight. The most stringent calibration requirement specifies knowledge of certain instrument pointing frames to an accuracy of 0.1 arcseconds, per-axis, 1-sigma relative to the Telescope Pointing Frame. In order to achieve this level of accuracy, the filter carries 37 states to estimate desired parameters while also correcting for expected systematic errors due to: (1) optical distortions, (2) scanning mirror scale-factor and misalignment, (3) frame alignment variations due to thermomechanical distortion, and (4) gyro bias and bias-drift in all axes. The resulting estimated pointing frames and calibration parameters are essential for supporting on-board precision pointing capability, in addition to end-to-end 'pixels on the sky' ground pointing reconstruction efforts.

calibration↗
NASA NTRS2007

Simplified Generation of High-Angular-Momentum Light Beams

A simplified method of generating a beam of light having a relatively high value of angular momentum (see figure) involves the use of a compact apparatus consisting mainly of a laser, a whispering- gallery-mode (WGM) resonator, and optical fibers. The method also can be used to generate a Bessel beam. ( Bessel beam denotes a member of a class of non-diffracting beams, so named because their amplitudes are proportional to Bessel functions of the radii from their central axes. High-order Bessel beams can have high values of angular momentum.) High-angular-momentum light beams are used in some applications in biology and nanotechnology, wherein they are known for their ability to apply torque to make microscopic objects rotate. High-angular-momentum light beams could also be used to increase bandwidths of fiber-optic communication systems. The present simplified method of generating a high-angular-momentum light beam was conceived as an alternative to prior such methods, which are complicated and require optical setups that include, variously, holograms, modulating Fabry-Perot cavities, or special microstructures. The present simplified method exploits a combination of the complex structure of the electromagnetic field inside a WGM resonator, total internal reflection in the WGM resonator, and the electromagnetic modes supported by an optical fiber. The optical fiber used to extract light from the WGM resonator is made of fused quartz. The output end of this fiber is polished flat and perpendicular to the fiber axis. The input end of this fiber is cut on a slant and placed very close to the WGM resonator at an appropriate position and orientation. To excite the resonant whispering- gallery modes, light is introduced into the WGM resonator via another optical fiber that is part of a pigtailed fiber-optic coupler. Light extracted from the WGM resonator is transformed into a high-angular- momentum beam inside the extraction optical fiber and this beam is emitted from the polished flat output end. By adjusting the geometry of this apparatus, it is possible to generate a variety of optical beams characterized by a wide range of parameters. These beams generally have high angular momenta and can be of either Bessel or Bessel-related types.

Savchenkov, Anatoliy↗
NASA NTRS2007

Procedures for Tuning a Multiresonator Photonic Filter

Two procedures have been devised for tuning a photonic filter that comprises multiple whispering-gallery mode (WGM) disk resonators. As used here, tuning signifies locking the filter to a specific laser frequency and configuring the filter to obtain a desired high-order transfer function. The main problem in tuning such a filter is how to select the correct relative loading of the resonators to realize a prescribed filter function. The first of the two procedures solves this problem. As temperature gradients develop during operation, the spectra of individual resonators tend to drift, primarily because of the thermorefractive effect. Thus, there arises the additional problem of how to adjust the tuning during operation to maintain the desired transfer function. The second of the two procedures solves this problem. To implement the procedures, it is necessary to incorporate the resonators into an apparatus like that of Figure 1. In this apparatus, the spectrum of each resonator can be adjusted individually, via the electro-optical effect, by adjusting a bias voltage applied to that resonator. In addition, the positions of the coupling prisms and resonators can be adjusted to increase or reduce the gaps between them, thereby reducing or increasing, respectively, the optical coupling between them. The optical power (Pi) in resonator i is monitored by use of a tracking photodiode. Another tracking diode monitors the power reflected from the input terminal (Pr), and still others monitor the input power (Pin) and output power (Po). The readings of these photodiodes are used to guide the tuning adjustments described in this paper.

Matsko, Andrey↗
NASA NTRS2006

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↗
NASA NTRS2000

Cellular Responses to Mechanical Stress Selected Contribution: A Three-Dimensional Model for Assessment of in Vitro Toxicity in Balaena Mysticetus Renal Tissue

This study established two- and three-dimensional renal proximal tubular cell cultures of the endangered species bowhead whale (Balaena mysticetus), developed SV40-transfected cultures, and cloned the 61-amino acid open reading frame for the metallothionein protein, the primary binding site for heavy metal contamination in mammals. Microgravity research, modulations in mechanical culture conditions (modeled microgravity), and shear stress have spawned innovative approaches to understanding the dynamics of cellular interactions, gene expression, and differentiation in several cellular systems. These investigations have led to the creation of ex vivo tissue models capable of serving as physiological research analogs for three-dimensional cellular interactions. These models are enabling studies in immune function, tissue modeling for basic research, and neoplasia. Three-dimensional cellular models emulate aspects of in vivo cellular architecture and physiology and may facilitate environmental toxicological studies aimed at elucidating biological functions and responses at the cellular level. Marine mammals occupy a significant ecological niche (72% of the Earth's surface is water) in terms of the potential for information on bioaccumulation and transport of terrestrial and marine environmental toxins in high-order vertebrates. Few ex vivo models of marine mammal physiology exist in vitro to accomplish the aforementioned studies. Techniques developed in this investigation, based on previous tissue modeling successes, may serve to facilitate similar research in other marine mammals.

Goodwin, T. J.↗
NASA NTRS2004

Phase Contrast Wavefront Sensing for Adaptive Optics

Most ground-based adaptive optics systems use one of a small number of wavefront sensor technologies, notably (for relatively high-order systems) the Shack-Hartmann sensor, which provides local measurements of the phase slope (first-derivative) at a number of regularly-spaced points across the telescope pupil. The curvature sensor, with response proportional to the second derivative of the phase, is also sometimes used, but has undesirable noise propagation properties during wavefront reconstruction as the number of actuators becomes large. It is interesting to consider the use for astronomical adaptive optics of the "phase contrast" technique, originally developed for microscopy by Zemike to allow convenient viewing of phase objects. In this technique, the wavefront sensor provides a direct measurement of the local value of phase in each sub-aperture of the pupil. This approach has some obvious disadvantages compared to Shack-Hartmann wavefront sensing, but has some less obvious but substantial advantages as well. Here we evaluate the relative merits in a practical ground-based adaptive optics system.

adaptive optics↗
NASA NTRS2005

A Broad-Band Phase-Contrast Wave-Front Sensor

A broadband phase-contrast wave-front sensor has been proposed as a real-time wave-front sensor in an adaptive-optics system. The proposed sensor would offer an alternative to the Shack-Hartmann wave-front sensors now used in high-order adaptive-optics systems of some astronomical telescopes. Broadband sensing gives higher sensitivity than does narrow-band sensing, and it appears that for a given bandwidth, the sensitivity of the proposed phase-contrast sensor could exceed that of a Shack-Hartmann sensor. Relative to a Shack-Hartmann sensor, the proposed sensor may be optically and mechanically simpler. As described below, an important element of the principle of operation of a phase-contrast wave-front sensor is the imposition of a 90deg phase shift between diffracted and undiffracted parts of the same light beam. In the proposed sensor, this phase shift would be obtained by utilizing the intrinsic 90 phase shift between the transmitted and reflected beams in an ideal (thin, symmetric) beam splitter. This phase shift can be characterized as achromatic or broadband because it is 90deg at every wavelength over a broad wavelength range.

Bloemhof, Eric↗
NASA NTRS2004

Imaging of gamma-Irradiated Regions of a Crystal

A holographic technique has been devised for generating a visible display of the effect of exposure of a photorefractive crystal to gamma rays. The technique exploits the space charge that results from trapping of electrons in defects induced by gamma rays. The technique involves a three-stage process. In the first stage, one writes a holographic pattern in the crystal by use of the apparatus shown in Figure 1. A laser beam of 532-nm wavelength is collimated and split into signal and reference beams by use of a polarizing beam splitter. On its way to the crystal, the reference beam goes through a two-dimensional optical scanner that contains two pairs of lenses (L1y, L2y and L1x,L2x) and mirrors M1 and M2, which can be rotated by use of micrometer drives to make fine adjustments. The signal beam is sent through a spatial light modulator that imposes the holographic pattern, then through two imaging lenses L(sub img) on its way to the crystal. An aperture is placed at the common focus of lenses Limg to suppress high-order diffraction from the spatial light modulator. The hologram is formed by interference between the signal and reference beams. A camera lens focuses an image of the interior of the crystal onto a charge-coupled device (CCD). If the crystal is illuminated by only the reference beam once the hologram has been formed, then an image of the hologram is formed on the CCD: this phenomenon is exploited to make visible the pattern of gamma irradiation of the crystal, as described next. In the second stage of the process, the crystal is removed from the holographic apparatus and irradiated with rays at a dose of about 100 krad. In the third stage of the process, the crystal is remounted in the holographic apparatus in the same position as in the first stage and illuminated with only the reference beam to obtain the image of the hologram as modified by the effect of the rays. The orientations of M1 and M2 can be adjusted slightly, if necessary, to maximize the intensity of the image. Figure 2 shows such an image that was formed in a crystal of Fe:LiNbO3.

Dragoi, Danut↗
NASA NTRS2003

Parallel Subconvolution Filtering Architectures

These architectures are based on methods of vector processing and the discrete-Fourier-transform/inverse-discrete- Fourier-transform (DFT-IDFT) overlap-and-save method, combined with time-block separation of digital filters into frequency-domain subfilters implemented by use of sub-convolutions. The parallel-processing method implemented in these architectures enables the use of relatively small DFT-IDFT pairs, while filter tap lengths are theoretically unlimited. The size of a DFT-IDFT pair is determined by the desired reduction in processing rate, rather than on the order of the filter that one seeks to implement. The emphasis in this report is on those aspects of the underlying theory and design rules that promote computational efficiency, parallel processing at reduced data rates, and simplification of the designs of very-large-scale integrated (VLSI) circuits needed to implement high-order filters and correlators.

Gray, Andrew A.↗
DOE OSTI2025

Coupling Approaches with Non-matching Grids for Classical Linear Elasticity and Bond-based Peridynamic Models in 1D

Local-nonlocal coupling approaches provide a means to combine the computational efficiency of local models and the accuracy of nonlocal models. To facilitate the coupling of the two models, non-matching grids are often desirable as nonlocal grids usually require a finer resolution than local grids. In that case, it is often convenient to resort to interpolation operators so that models can exchange information in the overlap regions when nodes from the two grids do not coincide. This paper studies three existing coupling approaches, namely 1) a method that enforces matching displacements in an overlap region, 2) a variant that enforces a constraint on the stresses instead, and 3) a method that considers a variable horizon in the vicinity of the interfaces. Further, the effect of the interpolation order and of the grid ratio on the performance of the three coupling methods with non-matching grids is carefully studied on one-dimensional examples using polynomial manufactured solutions. The numerical results show that the degree of the interpolants should be chosen with care to avoid introducing additional modeling errors, or simply minimize these errors, in the coupling approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

An extension of the localized artificial diffusivity method for immiscible and high density ratio flows

The localized artificial diffusivity (LAD) method is widely regarded as the preferred multi-material regularization scheme for the compact finite difference method, because it is conservative, easy to implement, and generally robust for a wide range of multi-material problems. However, traditional LAD methods face significant challenges when applied to flows with large density ratios and when maintaining thermodynamic equilibrium across material interfaces. These limitations arise from the formulation of the artificial diffusivity flux and the reliance on enthalpy diffusion for interface regularization. Additionally, traditional LAD methods struggle to ensure stability under large density ratio conditions, fail to maintain a finite interface thickness, and are therefore unsuitable for modeling immiscible interfaces. Here, in this work, we discuss the origins of these issues in traditional LAD methods and propose modifications which enable the simulation of large density ratio and immiscible flows. The proposed method targets the artificial diffusion fluxes at gradients and ringing in the volume fraction, rather than the mass fraction in traditional methods, to consistently regularize large density ratio interfaces. Furthermore, the proposed method introduces an artificial bulk density diffusion term to enforce equilibrium conditions across interfaces. To address the challenge of modeling immiscible flows, a conservative diffuse interface term is incorporated into the formulation to ensure a finite interface thickness. Specific consideration is taken in the design of the method to ensure that these crucial properties are maintained for N -material flows. The effectiveness of the proposed method is demonstrated through a series of canonical test cases, and its accuracy is validated by comparison with experimental data on micro-bubble collapse in water. These results highlight the method’s robustness and its ability to overcome the limitations of traditional LAD approaches.

Artificial diffusivity↗