Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

The a(3) Scheme--A Fourth-Order Space-Time Flux-Conserving and Neutrally Stable CESE Solver

The CESE development is driven by a belief that a solver should (i) enforce conservation laws in both space and time, and (ii) be built from a non-dissipative (i.e., neutrally stable) core scheme so that the numerical dissipation can be controlled effectively. To initiate a systematic CESE development of high order schemes, in this paper we provide a thorough discussion on the structure, consistency, stability, phase error, and accuracy of a new 4th-order space-time flux-conserving and neutrally stable CESE solver of an 1D scalar advection equation. The space-time stencil of this two-level explicit scheme is formed by one point at the upper time level and three points at the lower time level. Because it is associated with three independent mesh variables (the numerical analogues of the dependent variable and its 1st-order and 2ndorder spatial derivatives, respectively) and three equations per mesh point, the new scheme is referred to as the a(3) scheme. Through the von Neumann analysis, it is shown that the a(3) scheme is stable if and only if the Courant number is less than 0.5. Moreover, it is established numerically that the a(3) scheme is 4th-order accurate.

Chang, Sin-Chung↗
NASA NTRS2008

A Systematic Methodology for Constructing High-Order Energy-Stable WENO Schemes

A third-order Energy Stable Weighted Essentially Non-Oscillatory (ESWENO) finite difference scheme developed by Yamaleev and Carpenter (AIAA 2008-2876, 2008) was proven to be stable in the energy norm for both continuous and discontinuous solutions of systems of linear hyperbolic equations. Herein, a systematic approach is presented that enables \energy stable" modifications for existing WENO schemes of any order. The technique is demonstrated by developing a one-parameter family of fifth-order upwind-biased ESWENO schemes; ESWENO schemes up to eighth order are presented in the appendix. New weight functions are also developed that provide (1) formal consistency, (2) much faster convergence for smooth solutions with an arbitrary number of vanishing derivatives, and (3) improved resolution near strong discontinuities.

Yamaleev, Nail K.↗
NASA NTRS2009

High-Order Energy Stable WENO Schemes

A third-order Energy Stable Weighted Essentially Non-Oscillatory (ESWENO) finite difference scheme developed by Yamaleev and Carpenter was proven to be stable in the energy norm for both continuous and discontinuous solutions of systems of linear hyperbolic equations. Herein, a systematic approach is presented that enables 'energy stable' modifications for existing WENO schemes of any order. The technique is demonstrated by developing a one-parameter family of fifth-order upwind-biased ESWENO schemes; ESWENO schemes up to eighth order are presented in the appendix. New weight functions are also developed that provide (1) formal consistency, (2) much faster convergence for smooth solutions with an arbitrary number of vanishing derivatives, and (3) improved resolution near strong discontinuities.

Yamaleev, Nail K.↗
NASA NTRS2009

A Systematic Methodology for Constructing High-Order Energy Stable WENO Schemes

A third-order Energy Stable Weighted Essentially Non{Oscillatory (ESWENO) finite difference scheme developed by Yamaleev and Carpenter [1] was proven to be stable in the energy norm for both continuous and discontinuous solutions of systems of linear hyperbolic equations. Herein, a systematic approach is presented that enables "energy stable" modifications for existing WENO schemes of any order. The technique is demonstrated by developing a one-parameter family of fifth-order upwind-biased ESWENO schemes; ESWENO schemes up to eighth order are presented in the appendix. New weight functions are also developed that provide (1) formal consistency, (2) much faster convergence for smooth solutions with an arbitrary number of vanishing derivatives, and (3) improved resolution near strong discontinuities.

Yamaleev, Nail K.↗
NASA NTRS2009

The a(4) Scheme-A High Order Neutrally Stable CESE Solver

The CESE development is driven by a belief that a solver should (i) enforce conservation laws in both space and time, and (ii) be built from a nondissipative (i.e., neutrally stable) core scheme so that the numerical dissipation can be controlled effectively. To provide a solid foundation for a systematic CESE development of high order schemes, in this paper we describe a new high order (4-5th order) and neutrally stable CESE solver of a 1D advection equation with a constant advection speed a. The space-time stencil of this two-level explicit scheme is formed by one point at the upper time level and two points at the lower time level. Because it is associated with four independent mesh variables (the numerical analogues of the dependent variable and its first, second, and third-order spatial derivatives) and four equations per mesh point, the new scheme is referred to as the a(4) scheme. As in the case of other similar CESE neutrally stable solvers, the a(4) scheme enforces conservation laws in space-time locally and globally, and it has the basic, forward marching, and backward marching forms. Except for a singular case, these forms are equivalent and satisfy a space-time inversion (STI) invariant property which is shared by the advection equation. Based on the concept of STI invariance, a set of algebraic relations is developed and used to prove the a(4) scheme must be neutrally stable when it is stable. Numerically, it has been established that the scheme is stable if the value of the Courant number is less than 1/3

Chang, Sin-Chung↗
NASA NTRS2011

The Analysis of Rush Orders Risk in Supply Chain: A Simulation Approach

Satisfying customers by delivering demands at agreed time, with competitive prices, and in satisfactory quality level are crucial requirements for supply chain survival. Incidence of risks in supply chain often causes sudden disruptions in the processes and consequently leads to customers losing their trust in a company's competence. Rush orders are considered to be one of the main types of supply chain risks due to their negative impact on the overall performance, Using integrated definition modeling approaches (i.e. IDEF0 & IDEF3) and simulation modeling technique, a comprehensive integrated model has been developed to assess rush order risks and examine two risk mitigation strategies. Detailed functions sequence and objects flow were conceptually modeled to reflect on macro and micro levels of the studied supply chain. Discrete event simulation models were then developed to assess and investigate the mitigation strategies of rush order risks, the objective of this is to minimize order cycle time and cost.

Mahfouz, Amr↗
NASA NTRS2015

Projection-Based Reduced Order Modeling for Spacecraft Thermal Analysis

This paper presents a mathematically rigorous, subspace projection-based reduced order modeling (ROM) methodology and an integrated framework to automatically generate reduced order models for spacecraft thermal analysis. Two key steps in the reduced order modeling procedure are described: (1) the acquisition of a full-scale spacecraft model in the ordinary differential equation (ODE) and differential algebraic equation (DAE) form to resolve its dynamic thermal behavior; and (2) the ROM to markedly reduce the dimension of the full-scale model. Specifically, proper orthogonal decomposition (POD) in conjunction with discrete empirical interpolation method (DEIM) and trajectory piece-wise linear (TPWL) methods are developed to address the strong nonlinear thermal effects due to coupled conductive and radiative heat transfer in the spacecraft environment. Case studies using NASA-relevant satellite models are undertaken to verify the capability and to assess the computational performance of the ROM technique in terms of speed-up and error relative to the full-scale model. ROM exhibits excellent agreement in spatiotemporal thermal profiles (<0.5% relative error in pertinent time scales) along with salient computational acceleration (up to two orders of magnitude speed-up) over the full-scale analysis. These findings establish the feasibility of ROM to perform rational and computationally affordable thermal analysis, develop reliable thermal control strategies for spacecraft, and greatly reduce the development cycle times and costs.

0000↗
NASA NTRS2014

The Role of Higher-Order Modes on the Electromagnetic Whistler-Cyclotron Wave Fluctuations of Thermal and Non-Thermal Plasmas

Two fundamental challenging problems of laboratory and astrophysical plasmas are the understanding of the relaxation of a collisionless plasmas with nearly isotropic velocity distribution functions and the resultant state of nearly equipartition energy density with electromagnetic plasma turbulence. Here, we present the results of a study which shows the role that higher-order-modes play in limiting the electromagnetic whistler-like fluctuations in a thermal and non-thermal plasma. Our main results show that for a thermal plasma the magnetic fluctuations are confined by regions that are bounded by the least-damped higher order modes. We further show that the zone where the whistler-cyclotron normal modes merges the electromagnetic fluctuations shifts to longer wavelengths as the beta(sub e) increases. This merging zone has been interpreted as the beginning of the region where the whistler-cyclotron waves losses their identity and become heavily damped while merging with the fluctuations. Our results further indicate that in the case of nonthermal plasmas, the higher-order modes do not confine the fluctuations due to the effective higher-temperature effects and the excess of suprathermal plasma particles. The analysis presented here considers the second-order theory of fluctuations and the dispersion relation of weakly transverse fluctuations, with wave vectors parallel to the uniform background magnetic field, in a finite temperature isotropic bi-Maxwellian and Tsallis-kappa-like magnetized electron-proton plasma. Our results indicate that the spontaneously emitted electromagnetic fluctuations are in fact enhanced over these quasi modes suggesting that such modes play an important role in the emission and absorption of electromagnetic fluctuations in thermal or quasi-thermal plasmas.

Vinas, Adolfo F.↗
NASA NTRS2015

How Difficult is it to Reduce Low-Level Cloud Biases With the Higher-Order Turbulence Closure Approach in Climate Models?

Low-level clouds cover nearly half of the Earth and play a critical role in regulating the energy and hydrological cycle. Despite the fact that a great effort has been put to advance the modeling and observational capability in recent years, low-level clouds remains one of the largest uncertainties in the projection of future climate change. Low-level cloud feedbacks dominate the uncertainty in the total cloud feedback in climate sensitivity and projection studies. These clouds are notoriously difficult to simulate in climate models due to its complicated interactions with aerosols, cloud microphysics, boundary-layer turbulence and cloud dynamics. The biases in both low cloud coverage/water content and cloud radiative effects (CREs) remain large. A simultaneous reduction in both cloud and CRE biases remains elusive. This presentation first reviews the effort of implementing the higher-order turbulence closure (HOC) approach to representing subgrid-scale turbulence and low-level cloud processes in climate models. There are two HOCs that have been implemented in climate models. They differ in how many three-order moments are used. The CLUBB are implemented in both CAM5 and GDFL models, which are compared with IPHOC that is implemented in CAM5 by our group. IPHOC uses three third-order moments while CLUBB only uses one third-order moment while both use a joint double-Gaussian distribution to represent the subgrid-scale variability. Despite that HOC is more physically consistent and produces more realistic low-cloud geographic distributions and transitions between cumulus and stratocumulus regimes, GCMs with traditional cloud parameterizations outperform in CREs because tuning of this type of models is more extensively performed than those with HOCs. We perform several tuning experiments with CAM5 implemented with IPHOC in an attempt to produce the nearly balanced global radiative budgets without deteriorating the low-cloud simulation. One of the issues in CAM5-IPHOC is that cloud water content is much higher than in CAM5, which is combined with higher low-cloud coverage to produce larger shortwave CREs in some low-cloud prevailing regions. Thus, the cloud-radiative feedbacks are exaggerated there. The turning exercise is focused on microphysical parameters, which are also commonly used for tuning in climate models. The results will be discussed in this presentation.

Xu, Kuan-Man↗
NASA NTRS2017

A Fast Hybrid (3-D/1-D) Model for Thermal Radiative Transfer in Cirrus via Successive Orders of Scattering

We investigate the impact of cirrus cloud heterogeneity on the direct emission by cloud or surface and on the scattering by ice particles in the thermal infrared (TIR). Realistic 3-D cirri are modeled with the 3DCLOUD code, and top-of-atmosphere radiances are simulated by the 3-D Monte Carlo radiative transfer (RT) algorithm 3DMCPOL for two (8.65 micrometers and 12.05 micrometers) channels of the Imaging Infrared Radiometer on CALIPSO. At nadir, comparisons of 1-D and 3-D RT show that 3-D radiances are larger than their 1-D counterparts for direct emission but smaller for scattered radiation. For our cirrus cases, 99% of the 3-D total radiance is computed by the third scattering order, which corresponds to 90% of the total computational effort, but larger optical thicknesses need more scattering orders. To radically accelerate the 3-D RT computations (using only few percent of 3-D RT time with a Monte Carlo code), even in the presence of large optical depths, we develop a hybrid model based on exact 3-D direct emission, the first scattering order from 1-D in each homogenized column, and an empirical adjustment linearly dependent on the optical thickness to account for higher scattering orders. Good agreement is found between the hybrid model and the exact 3-D radiances for two very different cirrus models without changing the empirical parameters. We anticipate that a future deterministic implementation of the hybrid model will be fast enough to process multiangle thermal imagery in a practical tomographic reconstruction of 3-D cirrus fields.

radiative transfer↗
NASA NTRS2017

Impedance Eduction in Large Ducts Containing Higher-Order Modes and Grazing Flow

Impedance eduction test data are acquired in ducts with small and large cross-sectional areas at the NASA Langley Research Center. An improved data acquisition system in the large duct has resulted in increased control of the acoustic energy in source modes and more accurate resolution of higher-order duct modes compared to previous tests. Two impedance eduction methods that take advantage of the improved data acquisition to educe the liner impedance in grazing flow are presented. One method measures the axial propagation constant of a dominant mode in the liner test section (by implementing the Kumarsean and Tufts algorithm) and educes the impedance from an exact analytical expression. The second method solves numerically the convected Helmholtz equation and minimizes an objective function to obtain the liner impedance. The two methods are tested first on data synthesized from an exact mode solution and then on measured data. Results show that when the methods are applied to data acquired in the larger duct with a dominant higher-order mode, the same impedance spectra are educed as that obtained in the small duct where only the plane wave mode propagates. This result holds for each higher-order mode in the large duct provided that the higher-order mode is sufficiently attenuated by the liner.

Watson, Willie R.↗
NASA NTRS2015

Low Order Wavefront Sensing and Control for WFIRST-AFTA Coronagraph

To maintain the required WFIRST Coronagraph performance in a realistic space environment, a low order wavefront sensing and control (LOWFS/C) subsystem is necessary. The LOWFS/C use s the rejected stellar light from coronagraph to sense and suppress the telescope pointing drift and jitter as well as the low order wavefront errors due to changes in thermal loading of the telescope and the rest of the observatory. In this paper we will present an overview of the low order wavefront sensing and control subsystem for the WFIRST -AFTA Coronagraph. We will describe LOWFS/C's Zernike wavefront sensor concept and WFIRST LOWFS/C control design. We will present an overview of our analysis and modeling results on the Zernike wavefront sensor, the line -of-sight jitter suppression loop performance, as well as the low order wavefront error correction with the coronagraph's deformable mirror. In this paper we will also report the LOWFS/C testbed design and the preliminary in-air test results, which show a very promising performance of the Zernike wavefront sensor and FSM feedback loop.

exoplanets↗
NASA NTRS2019

Investigation of Reduced-Order Modeling for Aircraft Stability and Control Prediction

High-fidelity computational fluid dynamics tools offer the potential to approximate increments for ground-to-flight scaling effects, as well as to augment the dynamic damping derivative data for motion-based flight simulators. Unfortunately, the computational expense is currently prohibitive for populating a complete simulator database. This work investigates an existing surrogate-based, indicial response reduced-order model methodology as a means to efficiently augment a flight simulator database with high-fidelity nonlinear aerodynamic damping derivatives. Creation of the reduced-order model is based on the superposition integrals of the step response with the derivative of its corresponding input signal. Step responses are calculated using a computational grid motion approach that separates the effects of angle of attack and sideslip angle from angular rates, and rates from angle of attack and sideslip. It is demonstrated that the transients produced during the start of a forced-oscillation motion are captured by the reduced-order model to the level of fidelity of a comparable computational solution. Aerodynamic coefficients computed within minutes by the reduced-order model for an aircraft undergoing an 18-second half Lazy-8 maneuver and a 25-second Immelmann turn maneuver are compared with those from full computational flight solutions that required days to complete. Finally, a cost-benefit assessment is included that demonstrates a compelling advantage for this approach. d for maneuvering, flexible vehicles.

Frink, Neal T.↗
NASA NTRSDate not supplied

High-Order Shock Fitting with Finite Element Methods

A moving-grid, shock-tracking, finite element method has been implemented that can achieve high-order accuracy for flow simulations with shocks. In this approach, element edges in the computational mesh are fitted to the shock front and moved with the shock throughout the simulation. The Euler equations are solved on the moving mesh in an arbitrary Lagrangian-Eulerian framework. Three different methods for specifying the shock motion have been tested,and the order of accuracy of the resulting solutions was verified. It was found that one of the methods is applicable only for problems with a constant upstream flow state, whereas the other two methods converge to the expected solution for general shock problems. When using a finite element method with polynomial degree p, the two general shock motion methods converged with the expected order of accuracy of (p+1). Finally, accurate solutions are obtained for supersonic nozzle and blunt body bow shock problems. Based on the results, it is concluded that the shock tracking method can preserve high-order solution accuracy in the vicinity of a shock wave.

Luke D'Aquila↗
NASA NTRSDate not supplied

Efficient Basis Derivatives Evaluations for High-order Discontinuous Galerkin Schemes

Computational methods of evaluating a basis of Lagrange polynomials are developed for the purpose of implementing efficient discontinuous Galerkin conservation laws solvers. Special attention is payed to the computation of higher-order partial derivatives which may be required for certain applications. Two different approaches are considered; one involves hardcoding explicit one-line formulae into the source code in order to make it as simple as possible, whereas the other uses algorithms designed to minimize the asymptotic order of growth with respect to the order of the scheme. Timing experiments show that either approach can perform well if implemented effectively, and identify advantages and disadvantages of each.

Micaiah Smith-Pierce↗
NASA NTRS2022

Galaxy Correlation Function and Local Density from Photometric Redshifts Using the Stochastic Order Redshift Technique (SORT)

The stochastic order redshift technique (SORT) is a simple, efficient, and robust method to improve cosmological redshift measurements. The method relies upon having a small (∼10 per cent) reference sample of high-quality redshifts. Within pencil- beam-like sub-volumes surrounding each galaxy, we use the precise dN/dz distribution of the reference sample to recover new redshifts and assign them one-to-one to galaxies such that the original rank order of redshifts is preserved. Preserving the rank order is motivated by the fact that random variables drawn from Gaussian probability density functions with different means but equal standard deviations satisfy stochastic ordering. This process is repeated for sub-volumes surrounding each galaxy in the survey. This results in every galaxy being assigned multiple ‘recovered’ redshifts from which a new redshift estimate is determined. An earlier paper applied SORT to a mock Sloan Digital Sky Survey at z 0.2 and accurately recovered the two-point correlation function (2PCF) on scales > 4 h−1Mpc. In this paper, we test the performance of SORT in surveys spanning the redshift range 0.75 < z < 2.25. We used two mock surveys extracted from the Small MultiDark–Planck and Bolshoi–Planck N-body simulations with dark matter haloes that were populated by the Santa Cruz semi-analytic model. We find that SORT overall improves redshift estimates, accurately recovers the redshift-space 2PCF ξ (s) on scales > 2.5 h−1Mpc, and provides improved local density estimates in regions of average or higher density, which may allow for improved understanding of how galaxy properties relate to their environments.

James Kakos↗
NASA NTRS2024

Achieving High Efficiency in Reduced Order Modeling for Large Scale Polycrystal Plasticity Simulations

Reduced order models for the nonlinear response of heterogeneous microstructures typically require a construction (or training) stage to build the reduced order basis. In this manuscript, an efficient model construction strategy for the eigenstrain homogenization method (EHM) is presented. The proposed strategy relies on a parallel, element-by-element, conjugate gradient solver. Near linear scaling has been achieved with respect to the number of degrees of freedom used to resolve the microstructure. Linear scaling with respect to the number of pre-analyses required to construct the reduced order model (ROM) follows from the EHM formulation. Furthermore, a parallel implementation for fast evaluation of the constructed ROM has been developed using shared memory parallelization. It has been shown that for large microstructures with ≈ 10,000 grains, the total computational cost of evaluating the nonlinear response of a polycrystal could be reduced by approximately an order of magnitude using 32 cores with respect to serial ROM simulation. The present methodology has been verified using an additively manufactured polycrystalline microstructure of a nickel-based superalloy, Inconel 625. The capability of the developed framework to construct a ROM for such large microstructures, as well as the ability of the ROM to predict average and local quantities of interest has been demonstrated.

microscale↗
NASA NTRS1993

Spur-Reduced Digital Sinusoid Generation Using Higher-Order Phase Dithering

A higher-order phase dithering technique that reduces the complexity of digital sinusoid generation is presented and analyzed. M(sup th) order phase dithering denotes the addition of M independent, uniform variates to the sinusoid phase prior to word-length reduction. Spurious magnitudes due to quantization effects are accelerated from the usual -6 dBc per phase bit to -6(M + 1) dBc per phase bit, while the noise power increases linearly in M. For a given spurious specification, higher-order phase dithering permits the use of fewer phase bits than the non-dithered system, resulting in an exponential reduction in look-up table size, dramatically decreasing system complexity.

digital sinusoid generation higher-order phase dit↗