Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “A posteriori analysis”

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 19 records

Slat Cove Noise Modeling: A Posteriori Analysis of Unsteady RANS Simulations

A companion paper by Khorrami et al demonstrates the feasibility of simulating the (nominally) self-sustained, large-scale unsteadiness within the leading-edge slat-cove region of multi-element airfoils using unsteady Reynolds-Averaged Navier-Stokes (URANS) equations, provided that the turbulence production term in the underlying two-equation turbulence model is switched off within the cove region. In conjunction with a FfowesWilliams-Hawkings solver, the URANS computations were shown to capture the dominant portion of the acoustic spectrum attributed to slat noise, as well as reproducing the increased intensity of slat cove motions (and, correspondingly, far-field noise as well) at the lower angles of attack. This paper examines that simulation database, augmented by additional simulations, with the objective of transitioning this apparent success to aeroacoustic predictions in an engineering context. As a first step towards this goal, the simulated flow and acoustic fields are compared with experiment and simplified analytical model. Rather intense near-field fluctuations in the simulated flow are found to be associated with unsteady separation along the slat bottom surface, relatively close to the slat cusp. Accuracy of the laminar-cove simulations in this near-wall region is raised to be an open issue. The adjoint Green's function approach is also explored in an attempt to identify the most efficient noise source locations.

Choudhari, Meelan↗

Shape Optimization by Bayesian-Validated Computer-Simulation Surrogates

A nonparametric-validated, surrogate approach to optimization has been applied to the computational optimization of eddy-promoter heat exchangers and to the experimental optimization of a multielement airfoil. In addition to the baseline surrogate framework, a surrogate-Pareto framework has been applied to the two-criteria, eddy-promoter design problem. The Pareto analysis improves the predictability of the surrogate results, preserves generality, and provides a means to rapidly determine design trade-offs. Significant contributions have been made in the geometric description used for the eddy-promoter inclusions as well as to the surrogate framework itself. A level-set based, geometric description has been developed to define the shape of the eddy-promoter inclusions. The level-set technique allows for topology changes (from single-body,eddy-promoter configurations to two-body configurations) without requiring any additional logic. The continuity of the output responses for input variations that cross the boundary between topologies has been demonstrated. Input-output continuity is required for the straightforward application of surrogate techniques in which simplified, interpolative models are fitted through a construction set of data. The surrogate framework developed previously has been extended in a number of ways. First, the formulation for a general, two-output, two-performance metric problem is presented. Surrogates are constructed and validated for the outputs. The performance metrics can be functions of both outputs, as well as explicitly of the inputs, and serve to characterize the design preferences. By segregating the outputs and the performance metrics, an additional level of flexibility is provided to the designer. The validated outputs can be used in future design studies and the error estimates provided by the output validation step still apply, and require no additional appeals to the expensive analysis. Second, a candidate-based a posteriori error analysis capability has been developed which provides probabilistic error estimates on the true performance for a design randomly selected near the surrogate-predicted optimal design.

Patera, Anthony T.↗

Local and global accuracy estimates for boundary element analysis

An a posteriori error estimation strategy for the boundary element solution of two-dimensional heat conduction problem is presented. The error estimate formulation is based on the boundary element equations and employs a least squares fit of the solutions obtained from these equations. A local error estimate is evaluated for the unknown quantity in each element. The global error estimates are found from local (element) quantities. The range and accuracy of the proposed strategy is very good as demonstrated by problems containing smooth and non-smooth solutions.

Rencis, Joseph J.↗

On the Hilbert-Huang Transform Theoretical Developments

One of the main heritage tools used in scientific and engineering data spectrum analysis is the Fourier Integral Transform and its high performance digital equivalent - the Fast Fourier Transform (FFT). Both carry strong a-priori assumptions about the source data, such as linearity, of being stationary, and of satisfying the Dirichlet conditions. A recent development at the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC), known as the Hilbert-Huang Transform (HHT), proposes a novel approach to the solution for the nonlinear class of spectrum analysis problems. Using a-posteriori data processing based on the Empirical Mode Decomposition (EMD) sifting process (algorithm), followed by the normalized Hilbert Transform of the decomposition data, the HHT allows spectrum analysis of nonlinear and nonstationary data. The EMD sifting process results in a non-constrained decomposition of a source real value data vector into a finite set of Intrinsic Mode Functions (IMF). These functions form a near orthogonal adaptive basis, a basis that is derived from the data. The IMFs can be further analyzed for spectrum interpretation by the classical Hilbert Transform. A new engineering spectrum analysis tool using HHT has been developed at NASA GSFC, the HHT Data Processing System (HHT-DPS). As the HHT-DPS has been successfully used and commercialized, new applications post additional questions about the theoretical basis behind the HHT and EMD algorithms. Why is the fastest changing component of a composite signal being sifted out first in the EMD sifting process? Why does the EMD sifting process seemingly converge and why does it converge rapidly? Does an IMF have a distinctive structure? Why are the IMFs near orthogonal? We address these questions and develop the initial theoretical background for the HHT. This will contribute to the developments of new HHT processing options, such as real-time and 2-D processing using Field Programmable Array (FPGA) computational resources, enhanced HHT synthesis, and broaden the scope of HHT applications for signal processing.

Kizhner, Semion↗

On Certain Theoretical Developments Underlying the Hilbert-Huang Transform

One of the main traditional tools used in scientific and engineering data spectral analysis is the Fourier Integral Transform and its high performance digital equivalent - the Fast Fourier Transform (FFT). Both carry strong a-priori assumptions about the source data, such as being linear and stationary, and of satisfying the Dirichlet conditions. A recent development at the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC), known as the Hilbert-Huang Transform (HHT), proposes a novel approach to the solution for the nonlinear class of spectral analysis problems. Using a-posteriori data processing based on the Empirical Mode Decomposition (EMD) sifting process (algorithm), followed by the normalized Hilbert Transform of the decomposed data, the HHT allows spectral analysis of nonlinear and nonstationary data. The EMD sifting process results in a non-constrained decomposition of a source real-value data vector into a finite set of Intrinsic Mode Functions (IMF). These functions form a nearly orthogonal derived from the data (adaptive) basis. The IMFs can be further analyzed for spectrum content by using the classical Hilbert Transform. A new engineering spectral analysis tool using HHT has been developed at NASA GSFC, the HHT Data Processing System (HHT-DPS). As the HHT-DPS has been successfully used and commercialized, new applications pose additional questions about the theoretical basis behind the HHT and EMD algorithms. Why is the fastest changing component of a composite signal being sifted out first in the EMD sifting process? Why does the EMD sifting process seemingly converge and why does it converge rapidly? Does an IMF have a distinctive structure? Why are the IMFs nearly orthogonal? We address these questions and develop the initial theoretical background for the HHT. This will contribute to the development of new HHT processing options, such as real-time and 2-D processing using Field Programmable Gate Array (FPGA) computational resources,

Kizhner, Semion↗

Adaptive p-version based finite element formulations for thermal modeling/analysis of structural configurations

Adaptive p-version based hierarchical finite element formulations in conjunction with a posteriori error estimation concepts are described with emphasis on applicability for thermal modeling/analysis of structural configurations. The basic concepts and formulations of hierarchical p-versiion finite element for thermal analysis are first described. A posteriori error estimation features are utililzed to steer the process of adaptive refinement. Several configurations comprised of one-dimensional structures are evaluated to validate the applicability of the proposed formulations and to demonstrate the potential of the p-version adaptive formulations for thermal modeling/analysis. The methodology offers potential and promises to be an attractive alternative to conventional finite element thermal modeling/analysis approaches.

Tamma, Kumar K.↗

A Survey of Autonomous Navigation Techniques Applicable to Lunar Surface Exploration

As humanity returns to the Moon, and more and more attention is being paid to lunar surface operations, there is a greater need than ever for methods of surface navigation. These could be methods of computer-assisted orienteering for astronauts exploring on foot during an Extra-Vehicular Activity (EVA), or methods of solving the Lost-on-the-Moon problem to initialize a crewed or autonomous rover’s state estimate. It may also be necessary to process navigation data associated with surface samples or other surface operations a posteriori to better understand where that analysis occurred. Autonomous rover operation will also require Hazard Detection and Avoidance (HDA) and terrain-aware pathfinding. While navigation on the surface of the Moon will likely rely on Earth-based assets such as the Deep Space Network (DSN) or communication with other spacecraft (e.g., LunaNet, LCRNS, pre-deployed moon beacons, a nearby lander) it may be necessary to navigate in a loss-of-communication scenario. This paper analyzes the methods of surface navigation used on other celestial bodies, such as those used during the Apollo missions and autonomous exploration of Mars, as well as novel methods which have been studied but not yet implemented which may prove useful. It is shown that the navigator has myriad options when processing data from an Inertial Measurement Unit (IMU), a star tracker, (rover) wheel encoders, optical cameras, and LIght Detection and Ranging (LIDAR) sensors. The intention of this paper is to provide a broad overview of what has been done and what could be done, to aid those designing vehicles and/or missions to the lunar surface.

Paul D Mckee↗

A Survey of Autonomous Navigation Techniques Applicable to Lunar Surface Exploration

As humanity returns to the Moon, and more and more attention is being paid to lunar surface operations, there is a greater need than ever for methods of surface navigation. These could be methods of computer-assisted orienteering for astronauts exploring on foot during an Extra-Vehicular Activity (EVA), or methods of solving the Lost-on-the-Moon problem to initialize a crewed or autonomous rover’s state estimate. It may also be necessary to process navigation data associated with surface samples or other surface operations a posteriori to better understand where that analysis occurred. Autonomous rover operation will also require Hazard Detection and Avoidance (HDA) and terrain-aware pathfinding. While navigation on the surface of the Moon will likely rely on Earth-based assets such as the Deep Space Network (DSN) or communication with other spacecraft (e.g., LunaNet, LCRNS, pre-deployed moon beacons, a nearby lander) it may be necessary to navigate in a loss-of-communication scenario. This paper analyzes the methods of surface navigation used on other celestial bodies, such as those used during the Apollo missions and autonomous exploration of Mars, as well as novel methods which have been studied but not yet implemented which may prove useful. It is shown that the navigator has myriad options when processing data from an Inertial Measurement Unit (IMU), a star tracker, (rover) wheel encoders, optical cameras, and LIght Detection and Ranging (LIDAR) sensors. The intention of this paper is to provide a broad overview of what has been done and what could be done, to aid those designing vehicles and/or missions to the lunar surface.

Paul McKee↗

Towards cheaper control centers

Today, any approach to the design of new space systems must take into consideration an important constraint, namely costs. This approach is our guideline for new missions and also applies to the ground segment, and particularly to the control center. CNES has carried out a study on a recent control center for application satellites in order to take advantage of the experience gained. This analysis, the purpose of which is to determine, a posteriori, the costs of architecture needs and choices, takes hardware and software costs into account and makes a number of recommendations.

Baize, Lionel↗

On the application of multifrequency polarimetric radar observations for sea-ice classification

The use of multifrequency polarimetric radar imagery to enhance the ability to separate different sea-ice types using single-frequency, single-polarization synthetic aperture radar (SAR) data is investigated. Backscatter characteristics of six radiometrically and polarimetrically distinct sea-ice types are selected in an unsupervised range-dependent analysis of multifrequency polarimetric SAR data using the maximum a posteriori (MAP) polarimetric classifier. Maximum ice discrimination is achieved with combined C- and L-band full polarimetry, and collocated passive microwave imagery suggests greater than 90 percent classification accuracy. C-band VV-pol alone achieves only 68 percent relative accuracy because it confuses multiyear and rough compressed first year ice. L-band, relative classification accuracy is 75 percent, 83 percent, and 85 percent, using HH-pol, HH- and VV-combined, or the full polarimetry, respectively. P-band is less accurate. Combinations of two frequencies at a single polarization show the greatest improvement over a single channel.

Rignot, Eric↗

On the Hilbert-Huang Transform Data Processing System Development

One of the main heritage tools used in scientific and engineering data spectrum analysis is the Fourier Integral Transform and its high performance digital equivalent - the Fast Fourier Transform (FFT). The Fourier view of nonlinear mechanics that had existed for a long time, and the associated FFT (fairly recent development), carry strong a-priori assumptions about the source data, such as linearity and of being stationary. Natural phenomena measurements are essentially nonlinear and nonstationary. A very recent development at the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC), known as the Hilbert-Huang Transform (HHT) proposes a novel approach to the solution for the nonlinear class of spectrum analysis problems. Using the Empirical Mode Decomposition (EMD) followed by the Hilbert Transform of the empirical decomposition data (HT), the HHT allows spectrum analysis of nonlinear and nonstationary data by using an engineering a-posteriori data processing, based on the EMD algorithm. This results in a non-constrained decomposition of a source real value data vector into a finite set of Intrinsic Mode Functions (IMF) that can be further analyzed for spectrum interpretation by the classical Hilbert Transform. This paper describes phase one of the development of a new engineering tool, the HHT Data Processing System (HHTDPS). The HHTDPS allows applying the "T to a data vector in a fashion similar to the heritage FFT. It is a generic, low cost, high performance personal computer (PC) based system that implements the HHT computational algorithms in a user friendly, file driven environment. This paper also presents a quantitative analysis for a complex waveform data sample, a summary of technology commercialization efforts and the lessons learned from this new technology development.

Kizhner, Semion↗

Adaptive finite element analysis based on p-convergence

The results of numerical experiments are presented in which a posteriori estimators of error in strain energy were examined on the basis of a typical problem in linear elastic fracture mechanics. Two estimators were found to give close upper and lower bounds for the strain energy error. The potential significance of this is that the same estimators may provide a suitable basis for adaptive redistribution of the degrees of freedom in finite element models.

Szabo, B. A.↗

Rotor systems research aircraft risk-reduction shake test

A shake test and an extensive analysis of results were performed to evaluate the possibility of and the method for dynamically calibrating the Rotor Systems Research Aircraft (RSRA). The RSRA airframe was subjected to known vibratory loads in several degrees of freedom and the responses of many aircraft transducers were recorded. Analysis of the transducer responses using the technique of dynamic force determination showed that the RSRA, when used as a dynamic measurement system, could predict, a posteriori, an excitation force in a single axis to an accuracy of about 5 percent and sometimes better. As the analysis was broadened to include multiple degrees of freedom for the excitation force, the predictive ability of the measurement system degraded to about 20 percent, with the error occasionally reaching 100 percent. The poor performance of the measurement system is explained by the nonlinear response of the RSRA to vibratory forces and the inadequacy of the particular method used in accounting for this nonlinearity.

Wellman, J. Brent↗

Eigenvalue error analysis of viscously damped structures using a Ritz reduction method

The efficient solution of the eigenvalue problem that results from inserting passive dampers with variable stiffness and damping coefficients into a structure is addressed. Eigenanalysis of reduced models obtained by retaining a number of normal modes augmented with Ritz vectors corresponding to the static solutions resulting from the load patterns introduced by the dampers has been empirically shown to yield excellent approximations to the full eigenvalue problem. An analysis of this technique in the case of a single damper is presented. A priori and a posteriori error estimates are generated and tested on numerical examples. Comparison theorems with modally truncated models and a Markov parameter matching reduced-order model are derived. These theorems corroborate the heuristic that residual flexibility methods improve low-frequency approximation of the system. The analysis leads to other techniques for eigenvalue approximation. Approximate closed-form solutions are derived that include a refinement to eigenvalue derivative methods for approximation. An efficient Newton scheme is also developed. A numerical example is presented demonstrating the effectiveness of each of these methods.

Chu, Cheng-Chih↗

Stress Recovery and Error Estimation for 3-D Shell Structures

The C1-continuous stress fields obtained from finite element analyses are in general lower- order accurate than are the corresponding displacement fields. Much effort has focussed on increasing their accuracy and/or their continuity, both for improved stress prediction and especially error estimation. A previous project developed a penalized, discrete least squares variational procedure that increases the accuracy and continuity of the stress field. The variational problem is solved by a post-processing, 'finite-element-type' analysis to recover a smooth, more accurate, C1-continuous stress field given the 'raw' finite element stresses. This analysis has been named the SEA/PDLS. The recovered stress field can be used in a posteriori error estimators, such as the Zienkiewicz-Zhu error estimator or equilibrium error estimators. The procedure was well-developed for the two-dimensional (plane) case involving low-order finite elements. It has been demonstrated that, if optimal finite element stresses are used for the post-processing, the recovered stress field is globally superconvergent. Extension of this work to three dimensional solids is straightforward. Attachment: Stress recovery and error estimation for shell structure (abstract only). A 4-node, shear-deformable flat shell element developed via explicit Kirchhoff constraints (abstract only). A novel four-node quadrilateral smoothing element for stress enhancement and error estimation (abstract only).

Riggs, H. R.↗

DES Y3 + KiDS-1000: Consistent Cosmology Combining Cosmic Shear Surveys

We present a joint cosmic shear analysis of the Dark Energy Survey (DES Y3) and the Kilo-Degree Survey (KiDS-1000) in a collaborative effort between the two survey teams. We find consistent cosmological parameter constraints between DES Y3 and KiDS-1000 which, when combined in a joint-survey analysis, constrain the parameter S 8 =σ 8 √(Ω m /0.3) with a mean value of 0.790 +0.018 −0.014. The mean marginal is lower than the maximum a posteriori estimate, S 8 =0.801, owing to skewness in the marginal distribution and projection effects in the multi-dimensional parameter space. Our results are consistent with S 8 constraints from observations of the cosmic microwave background by Planck, with agreement at the 1.7 σ level. We use a Hybrid analysis pipeline, defined from a mock survey study quantifying the impact of the different analysis choices originally adopted by each survey team. We review intrinsic alignment models, baryon feedback mitigation strategies, priors, samplers and models of the non-linear matter power spectrum.

cosmology↗

Approaches to the automatic generation and control of finite element meshes

The algorithmic approaches being taken to the development of finite element mesh generators capable of automatically discretizing general domains without the need for user intervention are discussed. It is demonstrated that because of the modeling demands placed on a automatic mesh generator, all the approaches taken to date produce unstructured meshes. Consideration is also given to both a priori and a posteriori mesh control devices for automatic mesh generators as well as their integration with geometric modeling and adaptive analysis procedures.

Shephard, Mark S.↗