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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

The Parabolic Deformable Mirror Testbed at NASA Goddard Space Flight Center

One of the primary goals of future coronagraph space missions is to spectrally characterize Earth-like exoplanets. Observing efficiency is critical and directly related to how large of a field of view we can correct and the bandwidth of that correction. In addition to the coronagraph itself, the bandwidth and depth of contrast is related to the wavefront control system, which classically uses two deformable mirrors (DMs) in series. The minimum distance between these DMs is determined by the Talbot distance of the beam, meaning it scales as the square of the beam size. NASA Goddard is exploring an alternative design architecture to the two-DM approach in which the reimaging optics for the coronagraph are also deformable and we only retain a high-order pupil DM. With these off-axis parabolic elements being made deformable, the long propagation distance between the two DMs is eliminated. Goddard has procured an off-axis parabolic DM (PDM or OAP DM) from vendor ALPAO, France. A testbed is currently being built to test the PDM and demonstrate its utility for achieving high contrast in a coronagraph. The testbed is also being built with an integrated low-order wavefront sensing capability, allowing for integrated testing of a dynamic environment with multiple estimation and control loops.

Hari Subedi↗
NASA NTRS1994

Numerical Wave Propagation 211 Based on Wave Primitives

Compact higher order finite difference equations are applied to a sequence of problems in wave propagation and aeroacoustics. Systems of PDE's are reduced to a sequence of simple wave primitives using a local eigenvector decomposition. The wave primitives are first order PDE's in two independent variable and allow natural boundary conditions to be imposed for both single- and multidimensional problems. The method uses a "discrete dispersion relation" approach to obtain high order approximations to the wave primitives on a 3 spatial point / 2 time level computational molecule. The scheme is fourth order accurate for the class of system with constant coefficients, e.g., those that support exponential solutions. Weakly non-linear PDE's are solved in a similar manner using a variant of the "method of frozen coefficients." Experience with the new algorithm for linear, non-linear, and multi-dimensional test problems will be described.

Davis, Sanford S.↗
NASA NTRS1994

An HP Adaptive Discontinuous Galerkin Method for Hyperbolic Conservation Laws

This dissertation addresses various issues for model classes of hyperbolic conservation laws. The basic approach developed in this work employs a new family of adaptive, hp-version, finite element methods based on a special discontinuous Galerkin formulation for hyperbolic problems. The discontinuous Galerkin formulation admits high-order local approximations on domains of quite general geometry, while providing a natural framework for finite element approximations and for theoretical developments. The use of hp-versions of the finite element method makes possible exponentially convergent schemes with very high accuracies in certain cases; the use of adaptive hp-schemes allows h-refinement in regions of low regularity and p-enrichment to deliver high accuracy, while keeping problem sizes manageable and dramatically smaller than many conventional approaches. The use of discontinuous Galerkin methods is uncommon in applications, but the methods rest on a reasonable mathematical basis for low-order cases and has local approximation features that can be exploited to produce very efficient schemes, especially in a parallel, multiprocessor environment. The place of this work is to first and primarily focus on a model class of linear hyperbolic conservation laws for which concrete mathematical results, methodologies, error estimates, convergence criteria, and parallel adaptive strategies can be developed, and to then briefly explore some extensions to more general cases. Next, we provide preliminaries to the study and a review of some aspects of the theory of hyperbolic conservation laws. We also provide a review of relevant literature on this subject and on the numerical analysis of these types of problems.

Bey, Kim S.↗
NASA NTRS1974

Hypersonic flows generated by parabolic and paraboloidal shock waves

A computer algorithm has been developed to determine the blunt-body flowfields supporting symmetric parabolic and paraboloidal shock waves at infinite free-stream Mach number. Solutions are expressed in an analytic form as high-order power series, in the coordinate normal to the shock, whose coefficients can be determined exactly. Analytic continuation is provided by the use of Pade approximations. Test cases provide solutions of very high accuracy. In the axisymmetric case for gamma equals 715 the solution has been found far downstream, where it agrees with the modified blast-wave results. For plane flow, on the other hand, a limit line appears within the shock layer, a short distance past the sonic line, suggesting the presence of an imbedded shock. Local solutions in the downstream limit are discussed.

Schwartz, L. W.↗
NASA NTRS1981

The effects of missing data on global ozone estimates

The effects of missing data and model truncation on estimates of the global mean, zonal distribution, and global distribution of ozone are considered. It is shown that missing data can introduce biased estimates with errors that are not accounted for in the accuracy calculations of empirical modeling techniques. Data-fill techniques are introduced and used for evaluating error bounds and constraining the estimate in areas of sparse and missing data. It is found that the accuracy of the global mean estimate is more dependent on data distribution than model size. Zonal features can be accurately described by 7th order models over regions of adequate data distribution. Data variance accounted for by higher order models appears to represent climatological features of columnar ozone rather than pure error. Data-fill techniques can prevent artificial feature generation in regions of sparse or missing data without degrading high order estimates over dense data regions.

Drewry, J. W.↗
NASA NTRS1997

Runge-Kutta Methods for Linear Ordinary Differential Equations

Three new Runge-Kutta methods are presented for numerical integration of systems of linear inhomogeneous ordinary differential equations (ODES) with constant coefficients. Such ODEs arise in the numerical solution of the partial differential equations governing linear wave phenomena. The restriction to linear ODEs with constant coefficients reduces the number of conditions which the coefficients of the Runge-Kutta method must satisfy. This freedom is used to develop methods which are more efficient than conventional Runge-Kutta methods. A fourth-order method is presented which uses only two memory locations per dependent variable, while the classical fourth-order Runge-Kutta method uses three. This method is an excellent choice for simulations of linear wave phenomena if memory is a primary concern. In addition, fifth- and sixth-order methods are presented which require five and six stages, respectively, one fewer than their conventional counterparts, and are therefore more efficient. These methods are an excellent option for use with high-order spatial discretizations.

Zingg, David W.↗
NASA NTRS1973

Quasi-monochromatic measurements of homogeneous arc plasmas.

The refined diagnostic information obtainable by high-order spectrometry is illustrated by the results of quantitative measurements of a few rotational lines of OH in the ultraviolet spectrum of water-vapor plasmas generated in a wall-stabilized arc. Because of the high spectral and spatial resolution achieved in end-on measurements, the emission and also the absorption coefficients pertaining to homogeneous arc regions were obtained directly from measured line spectra - although the absorption was not measured explicitly - leading to the occupation of the upper and the lower state for the transition. The gas temperature was determined from the halfwidth of the Doppler-broadened rotational lines. The measured resolving power of the spectrometer was of the order of 400,000 in these measurements.

Klein, L.↗
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 NTRS1990

Numerical experiments on the accuracy of ENO and modified ENO schemes

Numerical experiments have been performed using different ENO schemes and different time discretizations, in order to assess accuracy-degeneracy phenomena of the type described by Rogerson and Meiburg (1990). A modified ENO scheme recovers the correct order of accuracy for all the test problems with smooth initial conditions and gives results comparable to the original ENO schemes for discontinuous problems. It is concluded that ENO schemes may lose the full high-order accuracy predicted by local truncation-error analysis for scalar linear conservation laws with smooth initial conditions. A modified ENO scheme can overcome this accuracy-degeneracy problem for the test problems without increasing the computational cost.

Shu, Chi-Wang↗
NASA NTRS1993

Protostellar hydrodynamics: Constructing and testing a spacially and temporally second-order accurate method. 2: Cartesian coordinates

In Boss & Myhill (1992) we described the derivation and testing of a spherical coordinate-based scheme for solving the hydrodynamic equations governing the gravitational collapse of nonisothermal, nonmagnetic, inviscid, radiative, three-dimensional protostellar clouds. Here we discuss a Cartesian coordinate-based scheme based on the same set of hydrodynamic equations. As with the spherical coorrdinate-based code, the Cartesian coordinate-based scheme employs explicit Eulerian methods which are both spatially and temporally second-order accurate. We begin by describing the hydrodynamic equations in Cartesian coordinates and the numerical methods used in this particular code. Following Finn & Hawley (1989), we pay special attention to the proper implementations of high-order accuracy, finite difference methods. We evaluate the ability of the Cartesian scheme to handle shock propagation problems, and through convergence testing, we show that the code is indeed second-order accurate. To compare the Cartesian scheme discussed here with the spherical coordinate-based scheme discussed in Boss & Myhill (1992), the two codes are used to calculate the standard isothermal collapse test case described by Bodenheimer & Boss (1981). We find that with the improved codes, the intermediate bar-configuration found previously disappears, and the cloud fragments directly into a binary protostellar system. Finally, we present the results from both codes of a new test for nonisothermal protostellar collapse.

Myhill, Elizabeth A.↗
NASA NTRS2003

Functional nanocomposites prepared by self-assembly and polymerization of diacetylene surfactants and silicic acid

Conjugated polymer/silica nanocomposites with hexagonal, cubic, or lamellar mesoscopic order were synthesized by self-assembly using polymerizable amphiphilic diacetylene molecules as both structure-directing agents and monomers. The self-assembly procedure is rapid and incorporates the organic monomers uniformly within a highly ordered, inorganic environment. By tailoring the size of the oligo(ethylene glycol) headgroup of the diacetylene-containing surfactant, we varied the resulting self-assembled mesophases of the composite material. The nanostructured inorganic host altered the diacetylene polymerization behavior, and the resulting nanocomposites show unique thermo-, mechano-, and solvatochromic properties. Polymerization of the incorporated surfactants resulted in polydiacetylene (PDA)/silica nanocomposites that were optically transparent and mechanically robust. Molecular modeling and quantum calculations and (13)C spin-lattice relaxation times (T(1)) of the PDA/silica nanocomposites indicated that the surfactant monomers can be uniformly organized into precise spatial arrangements prior to polymerization. Nanoindentation and gas transport experiments showed that these nanocomposite films have increased hardness and reduced permeability as compared to pure PDA. Our work demonstrates polymerizable surfactant/silica self-assembly to be an efficient, general approach to the formation of nanostructured conjugated polymers. The nanostructured inorganic framework serves to protect, stabilize, and orient the polymer, mediate its performance, and provide sufficient mechanical and chemical stability to enable integration of conjugated polymers into devices and microsystems.

NASA Discipline Life Support Systems↗
NASA NTRS2011

CFD Study of the Hydrocarbon Boost Low-Pressure Inducer and Kicker in the Presence of a Circumferential Groove

Results are presented of a computational fluid dynamics (CFD) study done in support of Marshall Space Flight Center's (MSFC) sub-scale water flow experiments of the Hydrocarbon Boost (HCB) Oxidizer Turbopump (OTP) being developed by the Air Force Research Laboratory (AFRL) and Aerojet. A circumferential groove may be added to the pump to reduce synchronous cavitation and subsequent bearing loads at a minimal performance cost. However, the energy may reappear as high order cavitation (HOC) that spans a relatively large frequency range. Thus, HOC may have implications for the full-scale OTP inducer in terms of reduced structural margin at higher mode frequencies. Simulations using the LOCI/Stream CFD program were conducted in order to explore the fluid dynamical impact of the groove on the low-pressure inducer and kicker. It was found that the circumferential groove has minimal head performance impact, but causes back-flowing high-swirl fluid to interact with the nearly-axial incoming fluid just above the inducer blades. The high-shear interface between the fluids is Kelvin-Helmholtz unstable, resulting in trains of low pressure regions or 'pearls' forming near the upstream edge of the groove. When the static pressure in these regions becomes low enough and they get cut by the blade leading edge, HOC is thought to occur. Although further work is required, the numerical models indicate that HOC will occur in the runbox of the AFRL/Aerojet HCB OTP. Comparisons to the ongoing water flow experiments will be discussed, as well as possible designs that may mitigate HOC while continuing to reduce synchronous cavitation. December 2011 MSS/LPS/SPS Joint Subcommittee Meeting ABSTRACT SUBMITTAL FORM

Coker, Robert↗
NASA NTRS2017

Comparison of High-Fidelity Computational Tools for Wing Design of a Distributed Electric Propulsion Aircraft

A variety of tools, from fundamental to high order, have been used to better understand applications of distributed electric propulsion to aid the wing and propulsion system design of the Leading Edge Asynchronous Propulsion Technology (LEAPTech) project and the X-57 Maxwell airplane. Three high-fidelity, Navier-Stokes computational fluid dynamics codes used during the project with results presented here are FUN3D, STAR-CCM+, and OVERFLOW. These codes employ various turbulence models to predict fully turbulent and transitional flow. Results from these codes are compared for two distributed electric propulsion configurations: the wing tested at NASA Armstrong on the Hybrid-Electric Integrated Systems Testbed truck, and the wing designed for the X-57 Maxwell airplane. Results from these computational tools for the high-lift wing tested on the Hybrid-Electric Integrated Systems Testbed truck and the X-57 high-lift wing presented compare reasonably well. The goal of the X-57 wing and distributed electric propulsion system design achieving or exceeding the required 𝐶 (sub L) = 3.95 for stall speed was confirmed with all of the computational codes.

Deere, Karen A.↗
NASA NTRS2011

LES of Temporally Evolving Mixing Layers by an Eighth-Order Filter Scheme

An eighth-order filter method for a wide range of compressible flow speeds (H.C. Yee and B. Sjogreen, Proceedings of ICOSAHOM09, June 22-26, 2009, Trondheim, Norway) are employed for large eddy simulations (LES) of temporally evolving mixing layers (TML) for different convective Mach numbers (Mc) and Reynolds numbers. The high order filter method is designed for accurate and efficient simulations of shock-free compressible turbulence, turbulence with shocklets and turbulence with strong shocks with minimum tuning of scheme parameters. The value of Mc considered is for the TML range from the quasi-incompressible regime to the highly compressible supersonic regime. The three main characteristics of compressible TML (the self similarity property, compressibility effects and the presence of large-scale structure with shocklets for high Mc) are considered for the LES study. The LES results using the same scheme parameters for all studied cases agree well with experimental results of Barone et al. (2006), and published direct numerical simulations (DNS) work of Rogers & Moser (1994) and Pantano & Sarkar (2002).

Hadjadj, A↗
NASA NTRS2017

Navier-Stokes Simulation of UH-60A Rotor/Wake Interaction Using Adaptive Mesh Refinement

Time-dependent Navier-Stokes simulations have been carried out for a flexible UH-60A rotor in forward flight, where the rotor wake interacts with the rotor blades. These flow conditions involved blade vortex interaction and dynamic stall, two common conditions that occur as modern helicopter designs strive to achieve greater flight speeds and payload capacity. These numerical simulations utilized high-order spatial accuracy and delayed detached eddy simulation. Emphasis was placed on understanding how improved rotor wake resolution affects the prediction of the normal force, pitching moment, and chord force of the rotor. Adaptive mesh refinement was used to highly resolve the turbulent rotor wake in a computationally efficient manner. Moreover, blade vortex interaction was found to trigger dynamic stall. Time-dependent flow visualization was utilized to provide an improved understanding of the numerical and physical mechanisms involved with three-dimensional dynamic stall.

UH-60A RotorWake↗
NASA NTRS1999

Computation of Steady and Unsteady Laminar Flames: Theory

In this paper we describe the numerical analysis underlying our efforts to develop an accurate and reliable code for simulating flame propagation using complex physical and chemical models. We discuss our spatial and temporal discretization schemes, which in our current implementations range in order from two to six. In space we use staggered meshes to define discrete divergence and gradient operators, allowing us to approximate complex diffusion operators while maintaining ellipticity. Our temporal discretization is based on the use of preconditioning to produce a highly efficient linearly implicit method with good stability properties. High order for time accurate simulations is obtained through the use of extrapolation or deferred correction procedures. We also discuss our techniques for computing stationary flames. The primary issue here is the automatic generation of initial approximations for the application of Newton's method. We use a novel time-stepping procedure, which allows the dynamic updating of the flame speed and forces the flame front towards a specified location. Numerical experiments are presented, primarily for the stationary flame problem. These illustrate the reliability of our techniques, and the dependence of the results on various code parameters.

Hagstrom, Thomas↗
NASA NTRS1990

Time-optimal Aircraft Pursuit-evasion with a Weapon Envelope Constraint

The optimal pursuit-evasion problem between two aircraft including a realistic weapon envelope is analyzed using differential game theory. Six order nonlinear point mass vehicle models are employed and the inclusion of an arbitrary weapon envelope geometry is allowed. The performance index is a linear combination of flight time and the square of the vehicle acceleration. Closed form solution to this high-order differential game is then obtained using feedback linearization. The solution is in the form of a feedback guidance law together with a quartic polynomial for time-to-go. Due to its modest computational requirements, this nonlinear guidance law is useful for on-board real-time implementation.

Menon, P. K. A.↗