Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Methods”

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 1,639 records · Page 91

Invariant Imbedded T-Matrix Method for Axial Symmetric Hydrometeors with Extreme Aspect Ratios

The single-scattering properties (SSPs) of hydrometeors are the fundamental quantities for physics-based precipitation retrievals. Thus, efficient computation of their electromagnetic scattering is of great value. Whereas the semi-analytical T-Matrix methods are likely the most efficient for nonspherical hydrometeors with axial symmetry, they are not suitable for arbitrarily shaped hydrometeors absent of any significant symmetry, for which volume integral methods such as those based on Discrete Dipole Approximation (DDA) are required. Currently the two leading T-matrix methods are the Extended Boundary Condition Method (EBCM) and the Invariant Imbedding T-matrix Method incorporating Lorentz-Mie Separation of Variables (IITM+SOV). EBCM is known to outperform IITM+SOV for hydrometeors with modest aspect ratios. However, in cases when aspect ratios become extreme, such as needle-like particles with large height to diameter values, EBCM fails to converge. Such hydrometeors with extreme aspect ratios are known to be present in solid precipitation and their SSPs are required to model the radiative responses accurately. In these cases, IITM+SOV is shown to converge. An efficient, parallelized C++ implementation for both EBCM and IITM+SOV has been developed to conduct a performance comparison between EBCM, IITM+SOV, and DDSCAT (a popular implementation of DDA). We present the comparison results and discuss details. Our intent is to release the combined ECBM IITM+SOV software to the community under an open source license.

Pelissier, Craig↗

A Fast Monte Carlo Method for Model-Based Prognostics Based on Stochastic Calculus

This work proposes a fast Monte Carlo method to solve differential equations utilized in model-based prognostics. The methodology is derived from the theory of stochastic calculus, and the goal of such a method is to speed up the estimation of the probability density functions describing the independent variable evolution over time. In the prognostic scenarios presented in this paper, the stochastic differential equations describe variables directly or indirectly related to the degradation of a monitored system. The method allows the estimation of the probability density functions by solving the deterministic equation and approximating the stochastic integrals using samples of the model noise. By so doing, the prognostic problem is solved without the Monte Carlo simulation based on Euler's forward method, which is typically the most time consuming task of the prediction stage. Three different prognostic scenarios are presented as proof of concept: (i) life prediction of electrolytic capacitors, (ii) remaining time to discharge of Lithium-ion batteries, and (iii) prognostic of cracked structures under fatigue loading. The paper shows how the method produces probability density functions that are statistically indistinguishable from the distributions estimated with Euler's forward Monte Carlo simulation. However, the proposed solution is orders of magnitude faster when computing the time-to-failure distribution of the monitored system. The approach may enable complex real-time prognostics and health management solutions with limited computing power.

Corbetta, M.↗

A Fast Monte Carlo Method for Model-Based Prognostics Based on Stochastic Calculus

This work proposes a fast Monte Carlo method to solve differential equations utilized in model-based prognostics. The methodology is derived from the theory of stochastic calculus, and the goal of such a method is to speed up the estimation of the probability density functions describing the independent variable evolution over time. In the prognostic scenarios presented in this paper, the stochastic differential equations describe variables directly or indirectly related to the degradation of a monitored system. The method allows the estimation of the probability density functions by solving the deterministic equation and approximating the stochastic integrals using samples of the model noise. By so doing, the prognostic problem is solved without the Monte Carlo simulation based on Euler's forward method, which is typically the most time consuming task of the prediction stage. Three different prognostic scenarios are presented as proof of concept: (i) life prediction of electrolytic capacitors, (ii) remaining time to discharge of Lithium-ion batteries, and (iii) prognostic of cracked structures under fatigue loading. The paper shows how the method produces probability density functions that are statistically indistinguishable from the distributions estimated with Euler's forward Monte Carlo simulation. However, the proposed solution is orders of magnitude faster when computing the time-to-failure distribution of the monitored system. The approach may enable complex real-time prognostics and health management solutions with limited computing power.

stochastic calculus↗

Comparison of Cross Correlation and Optical Flow Methods for Processing Retroreflective and Natural Background BOS Data

Background oriented schlieren images have historically been generated by calculating the observed pixel displacement between an image pair using normalized cross-correlation methods. This work uses optical flow data reduction methods to solve the displacement fields. A well established method used in the computer vision community, optical flow is the apparent brightness motion in an image sequence. The regularization method of Horn and Schunck is used to create schlieren images using two data sets: a supersonic shockjet plume interaction at the NASA Ames Unitary Plan Wind Tunnel, and a transonic flight test of a T-38 using a naturally occurring background performed in conjunction with NASA Ames and Armstrong research centers. Results are presented and contrasted with those using normalized cross-correlation methods. The optical flow images are found to provided significantly more detail at a decreased computational time.

schlieren↗

Cleanroom Contamination Identification Method Development

During fabrication, assembly, and testing of spacecraft and flight hardware it is vital to avoid contaminants that can cause degradation and could result in significant failure. Yet, there is no existing contamination monitoring method that provides the identity of airborne particles in a cleanroom facility. Knowing the particle identities, would allow scientists and engineers to determine the source of the contaminants and prevent setbacks before they occur or cause damage. Current cleanliness monitoring methods include airborne particle counters (APCs), fallout filters, and visual inspections. Particle counts from APCs are the primary metric used to define a cleanroom class and hence its level of cleanliness, but do not provide identification nor can they differentiate between large and small sizes of particles. In addition, using fallout filters is not a proactive, timely, or representative approach to cleanroom contamination monitoring because these samples are only retrieved after 30 days and are placed away from spacecraft processing to avoid interference with operations. In contrast, the forced air sampling method can collect a sample within an hour at any location required and provide results in less than a day. This system uses a cassette and filter sample medium to capture airborne particles which are then taken to a scanning electron microscope with energy dispersive spectroscopy (SEM/EDS) to identify and size the captured particles. Development of forced air sampling into an established laboratory capability will allow for fast sampling and routine identification of unknown contamination sources within the cleanroom. The test method development required market research for an air sampling cassette that increases sample collection efficiency and a filter with low enough background contamination to allow differentiation between a blank (control) and the collected sample. It was determined that a conductive black cassette and a polycarbonate filter were the best options. Conductive black cassettes, in comparison to the standard styrene, are manufactured using polypropylene filled with carbon. This makes the cassette conductive and minimizes the tendency of particles to stick to the wall of the cassette due to electrostatic force. In previous trials a mixed cellulose ester (MCE) filter was used to capture the contaminants, however the rougher surface of the filter contributed to entrapment of the particles within the filter structure and made it harder to identify the particles. In comparison, track etched polycarbonate filters have random cylindrical pores and a smooth surface which contributes to uniform sample distribution on the surface of the filter. Future work includes: testing the system using control samples to determine the efficiency and suitability of the medium, performing sample collection in various environments to establish ideal operating parameters and analyzing contaminant particles using SEM/EDS and assistant characterization techniques. Once fully developed, employing the forced air sampling method will help to prevent damage to spacecraft, avoid schedule delays, and allow for mission success.

Hernandez Melendez, Jailyn M.↗

Entropy Stable Method for the Euler Equations Revisited: Central Differencing via Entropy Splitting and SBP

The two decades old high order central differencing via entropy splitting and summation-by-parts (SBP) difference boundary closure of Ols- son & Oliger (1994), Gerritsen & Olsson (1996), and Yee et al. (2000) is revisited. The entropy splitting is a form of skew-symmetric splitting of the nonlinear Euler flux derivatives. Central differencing applied to the entropy splitting form of the Euler flux derivatives together with SBP difference operators will, hereafter, be referred to as entropy split schemes. This study is prompted by the recent growing interest in numerical methods for which a discrete entropy conservation law holds, a discrete global entropy conservation can be proved and/or the numerical method possesses a stable entropy in the framework of SBP difference operators and L2-energy norm estimate. The objective of this paper is to recast the entropy split scheme as the re- cent definition of an entropy stable method for central differencing with SBP operators for both periodic and non-periodic boundary conditions for non- linear Euler equations. Standard high order spatial central differencing as well as high order central spatial DRP (dispersion relation preserving) spatial differencing is part of the entropy stable methodology framework. Long time integration of 2D and 3D test cases is included to show the comparison of this efficient entropy stable method with the Tadmor-type of entropy conservative methods. Studies also include the comparison among the three skew-symmetric splittings on their nonlinear stability and accuracy performance without added numerical dissipations for smooth flows. These are, namely, entropy splitting, Ducros et al. splitting and the Kennedy & Grub- ber splitting.

Sjogreen, Bjorn↗

Weighted Least-Squares Cell-Average Gradient Construction Methods for the VULCAN-CFD Second-Order Accurate Unstructured-Grid Cell-Centered Finite-Volume Solver

The ability to solve the equations governing the hypersonic turbulent flow of a real gas on unstructured grids using a spatially-elliptic, 2nd-order accurate, cell-centered, finite-volume method has been recently implemented in the VULCAN-CFD code. The construction of cell-average gradients using a weighted linear least-squares method and the use of these gradients in the construction of the inviscid fluxes is the focus of this paper. A comparison of least-squares stencil construction methodologies is presented and approaches to augment the number of cells participating in the stencil while preserving accuracy are explored. Due to our interest in hypersonic flow, a robust multidimensional cell-average gradient limiter procedure that is consistent with the stencil used to construct the cell-average gradients is described and investigated. Canonical problems are computed to illustrate the challenges and investigate the accuracy, robustness and convergence behavior of the cell-average gradient methods on unstructured cell-centered finite-volume grids. Finally, thermally perfect, chemically frozen, Mach 8 turbulent flow of air around a blunt wedge is computed to demonstrate the robustness and convergence behavior of the new method for constructing stencils of use in a weighted linear least-squares gradient method for a hypersonic flow.

White, Jeffery A.↗

Developments on Low Cost Manufacturing Methods for Cylindrical Launch Vehicle Structures

Due to the increasing international competition and growing pressure to lower the cost of launch vehicles, enhanced technologies and competences for efficient, low cost manufacturing have to be developed. In this paper, an innovative manufacturing technology for launch vehicle applications is presented that was developed over the last few years at MT Aerospace AG (MTA) in collaboration with NASA Langley Research Center, the European Space Agency (ESA), Lockheed Martin Corporation (LMCO), Leifeld Metal Spinning and Scot Forge. The Integrally Stiffened Cylinder (ISC) manufacturing method was first investigated in 2010. Early results indicated a high mass saving potential for cylindrical structures for launch vehicle systems, such as cryogenic tank barrels and intertanks, produced using the ISC method. Cost-benefit studies conducted at NASA and MTA estimated up to 50% reduction in manufacturing costs is possible compared to conventional multi-piece, welded fabrication methods. Associated reductions in manufacturing lead times have shown that a 40% faster production rate is possible by the ISC process. Following proof of concept and process development efforts, a cylindrical structure manufactured using the ISC method had its maiden flight in 2015 on a Black Brant sounding rocket at NASA Wallops Flight Facility. Further collaboration between ESA, NASA, MTA, LMCO, and Scot Forge led to manufacturing test campaigns at MT Aerospace in 2017 and 2018 that successfully demonstrated ISC fabrication at a scale representative of commercial launch vehicles. This paper summarizes the status and relevant results of the development sequence of the ISC cylinder manufacturing method. Results of cost / benefit and structural performance analyses are also summarized. Moreover, an outlook will be given on further planned activities and development steps.

Martin Stachulla↗

Progress of Coating Stress Compensation of Silicon Mirrors for Lynx X-ray Telescope Mission Concept Using Thermal Oxide Patterning Method

A thermal oxide patterning method has proven to be effective for correcting coating-stress-induced distortion on flat silicon wafers. We report progress on developing this method for correcting curved silicon mirrors distorted by front-side iridium coatings. Owing to the difference in geometry, a finite element model has been established to calculate the appropriate duty cycle maps in thermal oxide hexagon patterns used for compensation. In addition, a photolithographic process, along with three-dimensional printed equipment, has been developed for creating patterns precisely on the back side of curved mirrors. The developed method has been used to recover the original surface shape of two silicon mirrors which are 100-mm long, 0.5-mm thick, having 312-mm radius of curvature, and 30 deg in azimuthal span (Wolter-I geometry). These mirrors’ front sides are sputter-coated by 20-nm iridium layers with ∼-70 N∕m integrated stress. Measurement results show that the developed method can mitigate coating-induced distortion by a factor of ∼5 in RMS height and ∼4 in RMS slope error, corresponding to ∼0.5 arc sec RMS slope error. Residual errors after correction are dominated by midfrequency ripples created during the annealing process, which will be resolved in the future. The presented method is precise and inexpensive and a potential candidate for resolving the coating stress issue for Lynx optics in the future.

x-ray↗

Realistic Thermographic Simulation of Impact Damage with Quadrupole Method

Flash thermography has been shown to be an effective method for detection of delaminations in carbon fiber reinforced polymer (CFRP) composite structures. Improved understanding of the limitations of the technique can be obtained by simulating the inspection process. Time domain finite difference and finite element methods are well suited for such simulations but can be computationally intensive. An alternate method is to solve the Laplace transform of the heat equation, then invert the Laplace transform to produce a time domain response. Often this is referred to as the thermal quadrupole method and is well suited for simulating flash thermography. The quadrupole method has been used extensively for simulating one-dimensional heat flow in multilayer systems with contact resistances at the interfaces. It is also applicable for three dimensional configurations, in particular for simulation of changes in thermographic response due to delaminations in composites. This presentation discusses three-dimensional simulations with realistic shapes and contact resistances at the interfaces between layers. The results of these simulations are compared to thermography measurements on composites specimens with impact induced delaminations.

thermography↗

Realistic Thermographic Simulation of Impact Damage with Quadrupole Method

Flash thermography has been shown to be an effective method for detection of delaminations in carbon fiber reinforced polymer (CFRP) composite structures. Improved understanding of thelimitations of the technique can be obtained by simulating the inspection process. Time domainfinite difference and finite element methods are well suited for such simulations but can be computationally intensive. An alternate method is to solve the Laplace transform of the heat equation, then invert the Laplace transform to produce a time domain response. Often this is referred to as the thermal quadrupole method and is well suited for simulating flash thermography. The quadrupole method has been used extensively for simulating one-dimensional heat flow in multilayer systems with contact resistances at the interfaces. It is also applicable for three dimensional configurations, in particular for simulation of changes in thermographic response due to delaminations in composites. This presentation discusses three-dimensional simulations with realistic shapes and contact resistances at the interfaces between layers. The results of these simulations are compared to thermography measurements on composites specimens with impact induced delaminations.

Thermography↗

Combining RNA-SEQ Datasets from NASA GENELAB: An Evaluation of Correction Methods

Background: Conducting space biology experiments aboard the International Space Station, particularly those utilizing complex model organisms like mice, is expensive and difficult due to limited crew availability, hardware, and space. As a result, sample numbers from these studies are low, reducing the statistical power of any one experiment. Aggregating spaceflight datasets serves as a method to increase sample numbers, allowing for novel insights through bioinformatic analysis of ‘omics data from merged datasets. However, aggregating datasets can introduce unwanted variation including 1) differences in sample handling, processing, and sequencing platforms between datasets (technical variation) as well as 2) differences in experimental design between datasets. Methods: In the present study, NASA GeneLab-hosted RNAseq datasets from mouse liver tissues were used to evaluate several statistical methods to correct for this unwanted variation through two approaches, reference-based and standard. The following correction algorithms were applied with (reference-based) and/or without (standard) considering Universal Mouse RNA Reference samples: ComBat and ComBat_seq from the SVA package, the median polish, empirical Bayes, and ANOVA-based algorithms from the MBatch package, and negative binomial regression normalization in the DESeq2 package. For each approach, after the correction algorithm was applied, differential gene expression (DGE) analysis of flight and ground control samples was performed with the combined data. The robustness of each tool was evaluated using BatchQC to determine statistical differences between datasets before and after correction, Principal Component Analysis to evaluate global gene expression in samples before and after correction, and by comparing DGE analysis of individual datasets and combined datasets before and after correction. Results: The results showed that the reference-based approach introduced several additional (and likely artificial) differentially expressed genes when compared with the respective standard approach. Conclusions: Of the methods tested, standard ComBat_seq and DESeq2 were identified as the most robust correction methods for combining spaceflight mouse liver RNAseq datasets hosted on GeneLab.

Finsam Samson↗

Assessment of Next Generation Airframe Noise Prediction Methods with PAA and ASN Flight Test Data

This paper presents an assessment of the airframe noise prediction capability in the research version of the NASA Aircraft Noise Prediction Program (ANOPP) with the objective of validating the accuracy of the current prediction methods, and more importantly, identifying the potential improvements to update the methods. The assessment compares predictions of the individual noise components currently implemented in ANOPP, including the landing gear, the slat, the flap, and the trailing edge noise, with the data from the Propulsion Airframe Aeroacoustics and Aircraft System Noise flight test conducted as part of the Boeing 2020 ecoDemonstrator program on a Boeing 787-10 Etihad Airways aircraft. For each component, prediction errors and discussions are presented for the noise characteristics, covering not only the source features but also the installation effects. It is shown that the predictions by the research level methods in ANOPP capture the major noise features, but improvement potential is identified in two main categories. The first includes minor features that have not been sufficiently modeled in the current methods, the noise due to interactions of the landing gear wake with the flaps and the slat bracket noise, for example. The second is related to changes in modern aircraft designs, which may alter the noise source mechanisms, and thus, calls for major improvements and updates on the current prediction models, examples including the flap side edge noise for modern flap systems and the trailing edge noise for realistic wings with taper and sweep. To illustrate the progress of the method development, previous generations of airframe noise models in ANOPP are also compared with the flight test data.

Yueping Guo↗

Comparisons of Performance Metrics and Machine Learning Methods on an Entry Descent and Landing Database

This work focuses on evaluating machine learning methods and their applicability to the generation of an aerodynamic database, particularly for trajectory analysis of a capsule during entry, descent, and landing with a focus on uncertainty quantification. The source data to be used is the wind tunnel and computational data for the Integrated Design Assessment Team (IDAT) configuration of the Orion project, which has been publicly released. The methods used to generate the proposed databases are designed to naturally include a prediction interval, which will be evaluated both for their mean response as well as how well the prediction interval performs. These machine learning methods are compared to a traditionally generated database used by the Orion team as a baseline. It is found that while these machine learning methods perform well, the Orion database tends to still outperform them showing that engineering experience is still needed to make the best database possible. However, these methods still provide comparable results with significantly less effort.

Orion↗

Comparisons of Performance Metrics and Machine Learning Methods on an Entry, Descent, and Landing Database

This work focuses on evaluating machine learning methods and their applicability to the generation of an aerodynamic database, particularly for trajectory analysis of a capsule during entry, descent, and landing with a focus on uncertainty quantification. The source data to be used is the wind tunnel and computational data for the Integrated Design Assessment Team (IDAT) configuration of the Orion project, which has been publicly released. The methods used to generate the proposed databases are designed to naturally include a prediction interval, which will be evaluated both for their mean response as well as how well the prediction interval performs. These machine learning methods are compared to a traditionally generated database used by the Orion team as a baseline. It is found that while these machine learning methods perform well, the Orion database tends to still outperform them showing that engineering experience is still needed to make the best database possible. However, these methods still provide comparable results with significantly less effort.

Orion↗

The Relationship Between Two Methods for Estimating Uncertainties in Data Assimilation

This note examines the relationship between two apparently unrelated methods for estimating error statistics or uncertainties of relevance to data assimilation. The first method is due to (Desroziers et al., 2005, Q. J. R. Meteorol. Soc., 131, 3385–3396; referred to as DBCP hereafter) and relies on residual statistics readily available from data assimilation applications. The second method, the three- cornered hat (3CH) developed by Gray and Allan (1974, IEEE 28th Annual Symp. Freq. Control, 243–246), only recently applied to atmospheric sciences, uses three data sets and can derive estimates of relevant error uncertainties as well. The usefulness of both methods lies in them not requiring knowledge of the true value of the quantities at play. DBCP derives its results by relying explicitly on the constraints as- sociated with the data assimilation minimization problem; 3CH is general and its estimates hold as long as errors in the three data sets of choice are uncorrelated. Establishing the relationship between the methods requires applying the 3CH approach to the same observation, background, and analysis data sets used by DBCP. In this case, the same assumptions of DBCP on residual errors allow for cancellation of error cross–covariance terms in 3CH such that two of its corners derive identical estimates for observation and background error covariances as those of DBCP. The error cross–covariance terms associated with the third corner are shown to add up to twice the analysis error covariance so that the 3CH estimate for the third corner recovers the negative of the analysis error covariance. Illustrations of these findings are provided by deriving uncertainties for radio occultation bending angles.

Ricardo Todling↗

A Moving Discontinuous Galerkin Method with Interface Condition Enforcement for Reacting Hypersonic Flows

The necessity of enforcing conservation in computational elements or cells (element conservation) for discontinuous solutions is well understood and respected for solving conservation laws in computational fluid dynamics (CFD). In contrast, interface conservation, where the conservation across cell interfaces is enforced, is long ignored, and yet is also ruled and required by the underlying physics just like element conservation. Violation of the interface conservation across discontinuities is the root cause why an exact discontinuous solution can never be achieved in shock capturing methods. The interface conservation is examined and explored in this talk. Moving discontinuous Galerkin (MDG) finite element method with interface conservation enforcement (MDG-ICE) [1],[2] are then presented for solving compressible flow problems with discontinuities based on the observation that the interface conservation can only be satisfied, only when mesh interfaces are aligned with discontinuities. In the MDG-ICE formulation, both conservative quantities and grid geometry are considered as independent variables. A space-time DG formulation is used to solve the multi-material compressible Euler equations in the standard discontinuous solution space and the discrete grid geometry is solved using a variational formulation in a continuous space. A self-adaptive Levenberg-Marquardt method is utilized to solve the resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation. A number of numerical experiments for a variety of flow problems are conducted to assess the accuracy and performance of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to deliver the designed order of both h- and p-convergence even for discontinuous solutions, and detect all types of interfaces, via interface condition enforcement and satisfy, via grid movement, the compressible Euler equations and the associated interface condition.

Discontinuous Galekin↗

Is There A Relationship Between Cornered-Hat Methods and A Residual Approach to Estimate System Uncertainty?

Recently a relationship has been established between a now traditional residual diagnostic used to estimate observation, background and analysis error covariances of interest to data assimilation and the three-cornered hat (3CH) method. An existing extension of the traditional residual diagnostic uses residuals from a fixed lag-1 Kalman smoother to retrieve system (model) error covariance. It is thus natural to ask if an additional relationship can be established between this extended residual method and some form of cornered--hat method. The answer might seem straightforward. Unlike in the standard residual estimation case, attempting to estimate model error amounts to estimating the statistics of a residual quantity itself. As shown in this presentation, in such cases the 3CH method becomes trivial: the sought out uncertainty can be derived directly from the covariance of the residual quantity at hand. In the particular case of estimating model error, one of the corners would have to be composed of vectors providing estimates of model error; these are not typically available in practice. Therefore, at first glance the present work finds no relationship between the lag-1 smoother residual diagnostic for system error estimation and cornered--hat methods. However, the work suggests that a two-tiered 3CH might be all that is necessary for system uncertainty to be obtained with 3CH.

Ricardo Todling↗