Results for “order” — Page 68

NASA NTRS2015

Atmospheric Turbulence Modeling for Aerospace Vehicles: Fractional Order Fit

An improved model for simulating atmospheric disturbances is disclosed. A scale Kolmogorov spectral may be scaled to convert the Kolmogorov spectral into a finite energy von Karman spectral and a fractional order pole-zero transfer function (TF) may be derived from the von Karman spectral. Fractional order atmospheric turbulence may be approximated with an integer order pole-zero TF fit, and the approximation may be stored in memory.

Kopasakis, George↗

NASA NTRS2016

Hyperbolic Method for Dispersive PDEs: Same High-Order of Accuracy for Solution, Gradient, and Hessian

In this paper, we introduce a new hyperbolic first-order system for general dispersive partial differential equations (PDEs). We then extend the proposed system to general advection-diffusion-dispersion PDEs. We apply the fourth-order RD scheme of Ref. 1 to the proposed hyperbolic system, and solve time-dependent dispersive equations, including the classical two-soliton KdV and a dispersive shock case. We demonstrate that the predicted results, including the gradient and Hessian (second derivative), are in a very good agreement with the exact solutions. We then show that the RD scheme applied to the proposed system accurately captures dispersive shocks without numerical oscillations. We also verify that the solution, gradient and Hessian are predicted with equal order of accuracy.

Mazaheri, Alireza↗

NASA NTRS2016

Tensor-Product Preconditioners for Higher-Order Space-Time Discontinuous Galerkin Methods

space-time discontinuous-Galerkin spectral-element discretization is presented for direct numerical simulation of the compressible Navier-Stokes equat ions. An efficient solution technique based on a matrix-free Newton-Krylov method is developed in order to overcome the stiffness associated with high solution order. The use of tensor-product basis functions is key to maintaining efficiency at high order. Efficient preconditioning methods are presented which can take advantage of the tensor-product formulation. A diagonalized Alternating-Direction-Implicit (ADI) scheme is extended to the space-time discontinuous Galerkin discretization. A new preconditioner for the compressible Euler/Navier-Stokes equations based on the fast-diagonalization method is also presented. Numerical results demonstrate the effectiveness of these preconditioners for the direct numerical simulation of subsonic turbulent flows.

Diosady, Laslo T.↗

NASA NTRS2018

Laboratory and On-Sky Validation of the Shaped Pupil Coronagraph's Sensitivity to Low-Order Aberrations with Active Wavefront Control

We present early laboratory simulations and extensive on-sky tests validating of the performance of a shaped pupil coronagraph (SPC) behind an extreme-AO (Adaptive Optics) corrected beam of the Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) system. In tests with the SCExAO internal source/wavefront error simulator, the normalized intensity profile for the SPC degrades more slowly than for the Lyot coronagraph as low-order aberrations reduce the Strehl ratio (S.R.) from extremely high values (S.R. approximately 0.93-0.99) to those characteristic of current ground-based extreme AO systems (S.R. approximately 0.74-0.93) and then slightly lower values down to S.R. approximately 0.57. On-sky SCExAO data taken with the SPC and other coronagraphs for brown dwarf/planet-hosting stars HD 1160 and HR 8799 provide further evidence for the SPC's robustness to low-order aberrations. From H-band Strehl ratios of 80 percent to 70 percent, the Lyot coronagraph's performance versus that of the SPC may degrade even faster on sky than is seen in our internal source simulations. The 5-sigma contrast also degrades faster (by a factor of two) for the Lyot than the SPC. The SPC we use was designed as a technology demonstrator only, with a contrast floor, throughput, and outer working angle poorly matched for SCExAO's current AO performance and poorly tuned for imaging the HR 8799 planets. Nevertheless, we detect HR 8799 cde with SCExAO/CHARIS using the SPC in broadband mode, where the signal-to-noise ratio for planet e is within 30 percent of that obtained using the vortex coronagraph. The shaped-pupil coronagraph is a promising design demonstrated to be robust in the presence of low-order aberrations and may be well-suited for future ground and space-based direct imaging observations, especially those focused on follow-up exoplanet characterization and technology demonstration of deep contrast within well-defined regions of the image plane.

Currie, Thayne↗

NASA NTRS2016

Low Order Wavefront Sensing and Control for WFIRST Coronagraph

To maintain the required WFIRST Coronagraph starlight suppression performance in a realistic space environment, a low order wavefront sensing and control (LOWFS/C) subsystem is necessary. The LOWFS/C uses 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 on 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 Coronagraph and describe the WFIRST Coronagraph LOWFS function, its design, and modeled performance. We will present experimental results on a dedicated LOWFS/C testbed that show that the LOWFS/C subsystem not only can sense pointing errors better than 0.2 mas but has also experimentally demonstrated closed loop pointing error suppression with residuals better than 0.4 mas rms per axis for the vast majority of observatory reaction wheel speeds.

wavefront sensing and control↗

NASA NTRS2018

Hybrid Reduced Order Model for N2-N2 Interactions for Application to Dissociation and Energy Transfer Processes

This work presents a general framework for model reduction of non-equilibrium energy transfer and dissociation processes. The multi-group maximum-entropy method is coupled with the quasi-classical trajectory method to directly construct a reduced order model for chemical non-equilibrium. Kinetic data is calculated by applying detailed balance at a microscopic level, overcoming the limitation of calculating recombination kinetic data. This approach enables the construction of a reduced order model for kinetics which bypasses the need to compute state-to-state kinetic data. This physics based reduced order model ensures that an equilibrium distribution is reached given infinite time, while allowing for non-equilibrium distributions during the relaxation and dissociation processes. A proof-of-concept test case demonstrates the applicability of this model by comparison with state-to-state kinetic data for the N2(Chi (sup 1) Sigma ((sub g), (sup plus sign)) minus N ((exp 4) S (sub u)) system with excellent agreement.

Maximum Entropy↗

NASA NTRS2018

Parametric Analysis of Aircraft WingWeight Using Low-Order Physics-Based Analysis

In the conceptual aircraft design phase, prediction of the empty weight typically relies on empirically-based regression equations which execute quickly and require little detailed information about the internal structural layout. Since they are based on existing aircraft, however, empirical methods can lose their validity for newer technologies and unconventional configurations. Designers can transition to higher-order, physics-based analysis methods to improve the accuracy of the weight prediction, but at the cost of complex model setup and increased computational time. This paper describes a methodology for low-order aero-structural analysis of conceptual aircraft configurations that increases the use of physics-based analysis in conceptual design, but is less complex and time-consuming than higher-order methods such as finite-element analysis. The methodology uses Vehicle Sketch Pad (OpenVSP) to model the aircraft geometry, and ASWING to perform the aero-structural analysis. The internal forces and moments from the ASWING analysis are post-processed to calculate the resulting direct and shear stresses in the structure, and the thickness distributions of the aircraft components are varied to match the maximum von Mises stress at each cross section to the material allowable. To offset the increased computational time relative to empirical weight equations, a process is studied which uses parametric variation to develop a regression equation relating the weight of the aircraft wing to major design variables. This new weight equation is similar to existing empirical equations, but is built using the more physics-based methodology; the new equation could be used to augment or replace portions of the empirical database to improve the validity of the wing weight prediction for unconventional configurations and advanced technologies.

Olson, Erik D.↗

NASA NTRS2019

Reduced Order Modeling for Transonic Aeroservoelastic Control Law Development

As aircraft become more flexible, aeroelastic considerations become increasingly important and complex, particularly for transonic flight where nonlinearities in the flow render linear analysis tools less effective. In order to analyze these aeroelastic interactions between the fluid and the structure efficiently, reduced order models (ROMs) are sometimes generated from and used in place of computational fluid dynamics solutions. In this paper, several aerodynamic ROMs are generated and coupled with structural models to form aeroelastic ROMs. The aerodynamic ROMs generated here include the effects of control surface motion. Hence, the aeroelastic ROMs presented here are appropriate for use in aeroservoelastic applications and are intended to be used for aeroservoelastic control law development. These ROMs are used to simulate a number of test cases with and without control surface involvement. Results show that several of the ROMs generated in the paper are able to predict results similar to solutions of higher-order computational methods.

Waite, Josiah M.↗

NASA NTRS2019

Aeroelastic Indicial Response Reduced-Order Modeling for Flexible Flight Vehicles

A reduced-order modeling method capable of providing computationally efficient predictions of the nonlinear, unsteady aerodynamics encountered by flexible flight vehicles under- going forced oscillations is presented. Models are developed using indicial response theory, which characterizes a vehicle’s dynamics through identification of time-accurate aerodynamic responses due to step changes in the vehicle-state parameters, e.g., angle-of-attack, pitch rate. A coupled computational fluid dynamics aeroelastic analysis is proposed for identifying step (indicial) responses of flexible vehicles. In this approach, aeroelastic indicial responses are simulated via prescribed rigid body motions, while fluid-structure interactions are captured at the subiterative level through coupling to a linear modal structural solver. A nonlinear extension of indicial response theory is applied through time-dependent linear interpolation of a database of locally linear aeroelastic step responses. Reduced-order models are then created using the mathematical principle of convolution applied to the interpolated aeroelastic indicial responses to predict the time-dependent aerodynamic response of a flexible vehicle to any arbitrary prescribed flight maneuver. The NASA FUN3D computational fluid dynamics solver is utilized for simulating full-order trajectories and indicial response functions. Aerodynamic predictions were generated for the X-56A aircraft undergoing a series of harmonic forced oscillations. The s are shown to provide a practical option for evaluating the unsteady aerodynamics of flexible vehicles using high-fidelity simulations.

Hiller, Brett↗

NASA NTRS2021

PatchworkWave: A Multipatch Infrastructure for Multiphysics/Multiscale/Multiframe/Multimethod Simulations at Arbitrary Order

We present an extension of the PatchworkMHD code [1], itself an MHD-capable extension of thePatch-workcode [2], for which several algorithms presented here were co-developed. Its purpose is to create a“multipatch” scheme compatible with numerical simulations of arbitrary equations of motion at any dis-cretization order in space and time. In thePatchworkframework, the global simulation is comprised of anarbitrary number of moving, local meshes, or “patches”, which are free to employ their own resolution, co-ordinate system/topology, physics equations, reference frame, and in our new approach, numerical method.Each local patch exchanges boundary data with a single global patch on which all other patches residethrough a client-router-server parallelization model. In generalizingPatchworkto be compatible witharbitrary order time integration,PatchworkMHDandPatchworkWavehave significantly improved theinterpatch interpolation accuracy by removing an interpolation of interpolated data feedback present in theoriginalPatchworkcode. Furthermore, we extendPatchworkto bemultimethodby allowing multiplestate vectors to be updated simultaneously, with each state vector providing its own interpatch interpolationand transformation procedures. As such, our scheme is compatible with nearly any set of hyperbolic partialdifferential equations. We demonstrate our changes through the implementation of a scalar wave toy-modelthat is evolved on arbitrary, time dependent patch configurations at 4th order accuracy.

Dennis B Bowen↗

NASA NTRS2018

Runtime Veriﬁcation: From Propositional to First-Order Temporal Logic

Runtime Veriﬁcation is a branch of formal methods concerned with analysis of execution traces for the purpose of determining the state or general quality of the executing system. The ﬁeld covers numerous approaches, one of which is speciﬁcation-based runtime veriﬁcation, where execution traces are checked against formal speciﬁcations. The paper presents syntax, semantics, and monitoring algorithms for respectively propositional and ﬁrst-order temporal logics. In propositional logics the observed events in the execution trace are represented using atomic propositions, while ﬁrst-order logic allows universal and existential quantiﬁcation over data occurring as arguments in events. Monitoring of the ﬁrst-order case is drastically more challenging than the propositional case, and we present a solution for this problem based on BDDs. We furthermore discuss monitorability of temporal properties by dividing them into different classes representing different degrees of monitorability.

Peled, Doron↗

NASA NTRS2016

Low order wavefront sensing and control for WFIRST coronagraph

To maintain the required WFIRST Coronagraph starlight suppression performance in a realistic space environment, a low order wavefront sensing and control (LOWFS/C) subsystem is necessary. The LOWFS/C uses 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 on 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 Coronagraph and describe the WFIRST Coronagraph LOWFS function, its design, and modeled performance. We will present experimental results on a dedicated LOWFS/C testbed that show that the LOWFS/C subsystem not only can sense pointing errors better than 0.2 mas but has also experimentally demonstrated closed loop pointing error suppression with residuals better than 0.4 mas rms per axis for the vast majority of observatory reaction wheel speeds.

Wilson, Dan↗

NASA NTRSDate not supplied

Efficient Preconditioning of a High-Order Solver for Multiple Physics

This work addresses preconditioning approaches for an implicit high-order solver frame-work applied to multiple physics. The solver is based on a space-time spectral element method and matrix-free Newton-Krylov solver developed at NASA over the recent years. Within this context, most preconditioning methods are impractical, as the computational time and memory requirements scale poorly with increasing polynomial orders. To improve computational efficiency, we first describe a novel entity-based Block Jacobi preconditioner for the continuous-Galerkin solution of the linear-elasticity and linear-shell equations. Second, we introduce a multigrid algorithm to further reduce time-to-solution on stiff cases arising from continuous-and discontinuous-Galerkin discretizations. Results obtained on relevant single-physics reference solutions, demonstrate the feasibility of the methods, paving the way for high-order solutions of fully coupled multi-physics problems.

STMD↗

NASA NTRS2021

Evaluation of Higher-order Quadrature Schemes in Improving Computational Efficiency for Orientation-averaged Single-Scattering Properties of Nonspherical Ice Particles

We evaluate several high-order quadrature schemes for accuracy and efficacy in obtaining orientation-averaged single-scattering properties (SSPs). We use the highly efficient MIDAS to perform electromagnetic scattering calculations to evaluate the gain in efficiency from these schemes. MIDAS is shown to be superior to DDSCAT, a popular discrete dipole approximation (DDA) method. This study is motivated by the fact that quality physical precipitation retrievals rely on using accurate orientation-averaged SSPs derived from realistic hydrometeors as input to radiative transfer Models (RTMs). The DDA has been a popular choice for single-scattering calculations, due to its versatility with respect to target geometry. However, being iterative-solver-based (ISB), the most used DDA codes, e.g. DDSCAT and ADDA, must solve the scattering problem for each orientation of the target separately. As the size parameter and geometric anisotropy of the hydrometeor increase, the number of orientations needed to obtain accurate orientation-averages can increase drastically and so does the computation cost incurred by the ISB-DDA methods. MIDAS is a Direct-Solver-Based (DSB) code, its decomposition of the original large matrix with a high rank into multiple more manageable smaller matrices of lower ranks makes it much more computationally efficient while maintaining excellent accuracy. In addition, direct solvers consider all requested orientations at once, giving MIDAS further advantage over popular ISB-DDA methods. MIDAS, when combined with high-order quadrature for orientation averaging, can be greater than three orders of magnitude more efficient in obtaining RTM-ready SSPs of complex-shaped hydrometeors than existing ISB-DDA methods, with the native quadrature schemes they offer.

Ines Fenni↗

NASA NTRSDate not supplied

Higher-Order Statistical Moments of Predicted Sonic Boom Waveforms Through Turbulence

Sonic booms generated by supersonic aircraft are affected by turbulence in the atmospheric boundary layer through which they propagate. Turbulence effects lead to random variability of the sonic boom waveforms measured on the ground, complicating the prediction of such waveforms. As an initial effort to predict the waveform variability, the solution of the coherent or mean sonic boom waveform has previously been formulated by the author. The current paper extends the formulation to derive an expression of the second-order statistical moment necessary to calculate the variance of the spectral amplitudes. Since the derivation uses a full wave equation, results using the derived expression will be compared with those obtained using a parabolic approximation. In order to fully quantify the uncertainty of our predictions, the higher order statistical moments are also formulated and are shown to be approximately zero. Consequently, the probability density function (pdf) of the spectral amplitudes is predicted to be Gaussian. The formulation is further extended to determine the pdf of the loudness of sonic booms and quantify the uncertainties associated with the loudness prediction. Results from the formulation are compared with available flight test data.

sonic booms↗

NASA NTRSDate not supplied

Shock Capturing via Limiting for High-Order Methods including Discontinuous Galerkin

High-order methods, such as discontinuous Galerkin (DG), spectral, and flux reconstruction (FR), are prone to generating unwanted oscillations near shocks and discontinuities. Conventional limiting techniques, while effective in suppressing oscillations near shocks, often compromise accuracy near extrema, where the solution is only first-order accurate. This paper introduces a novel limiting technique for these high-order schemes, aimed at effectively managing shocks while preserving accuracy. The key idea is to expand the standard monotonicity limits to provide “room” near smooth extrema, ensuring that limiting has no effect and thus preserving accuracy. Near a discontinuity, these expanded limits effectively reduce to the original monotonicity limits, suppressing oscillations. Additional motivation is drawn from a formula for the derivative of Radau polynomials, which depicts the behavior of oscillations resulting from discontinuities. This behavior leads to a simplification by applying the limits to the sum of magnitudes of all modes, linear and higher degree. Unlike typical approaches, which rely on successful detection to activate limiting, our limiter depends continuously on the data, there by avoiding potential issues if detection fails. To reduce computing time, efficient criteria for detecting smooth regions where limiting is unnecessary are presented. Combined with detection, the continuous dependence on the data is lost, but the method is more economical. A notable characteristic of the entire process is its simplicity in both concept and implementation. Numerical tests for advection and Euler equations are conducted to demonstrate the effectiveness of the proposed method.

numerical methods↗

NASA NTRSDate not supplied

Shock Capturing via Limiting for High-Order Methods Including Discontinuous Galerkin

As is well-known, popular high-order methods such as discontinuous Galerkin (DG) and flux reconstruction (FR) tend to generate oscillations near shocks and discontinuities, which can lead to negative pressure or density, ultimately causing code breakdown. For standard second-order methods, one approach to mitigating oscillations is to impose constraints on the calculated slopes, ensuring they do not become excessively steep. This limiting process is rooted in the idea of preserving monotonicity introduced by Van Leer (1974): when the data are monotone, limiting the slopes result in a monotone piecewise linear solution. However, a drawback of this approach is the loss of accuracy near extrema, where the non-monotone solution is somewhat flattened and achieves only first-order accuracy.

numerical methods↗

NASA NTRSDate not supplied

Shock Capturing via Limiting for High-Order Methods including Discontinuous Galerkin

High-order methods, such as discontinuous Galerkin (DG), spectral, and flux reconstruction (FR), are prone to generating unwanted oscillations near shocks and discontinuities. Conventional limiting techniques, while effective in suppressing oscillations near shocks, often compromise accuracy near extrema, where the solution is only first-order accurate. This paper introduces a novel limiting technique for these high-order schemes, aimed at effectively managing shocks while preserving accuracy. The key idea is to expand the standard monotonicity limits to provide “room” near smooth extrema, ensuring that limiting has no effect and thus preserving accuracy. Near a discontinuity, these expanded limits effectively reduce to the original monotonicity limits, suppressing oscillations. Additional motivation is drawn from a formula for the derivative of Radau polynomials, which depicts the behavior of oscillations resulting from discontinuities. This behavior leads to a simplification by applying the limits to the sum of magnitudes of all modes, linear and higher degree. Unlike typical approaches, which rely on successful detection to activate limiting, our limiter depends continuously on the data, there by avoiding potential issues if detection fails. To reduce computing time, efficient criteria for detecting smooth regions where limiting is unnecessary are presented. Combined with detection, the continuous dependence on the data is lost, but the method is more economical. A notable characteristic of the entire process is its simplicity in both concept and implementation. Numerical tests for advection and Euler equations are conducted to demonstrate the effectiveness of the proposed method.

numerical methods↗