Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “software packages”

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,189 records · Page 66

Comparison of Design Tools for Stress Analysis of Adhesively Bonded Joints

Analytical or semi-analytical models for stress analysis have long been a part of initial adhesive bond design and sizing. Even with the rise of general finite element software and methods, these design models have still remained a preferred method for fast and simple joint analysis. While these methods can yield fairly representative results, the models are constructed on the foundation of geometrical and material simplifications or assumptions that allow closed form or semi-closed form solutions. This study outlines the major differences in basic assumptions for three common design software packages under use at NASA, and shows the ramifications these assumptions in a few exemplar bonded joints.

Composites↗

Statistical Hypothesis Testing in Wavelet Analysis: Theoretical Developments and Applications to Indian Rainfall

Statistical hypothesis tests in wavelet analysis are methods that assess the degree to which a wavelet quantity(e.g., power and coherence) exceeds background noise. Commonly, a point-wise approach is adopted in which a wavelet quantity at every point in a wavelet spectrum is individually compared to the critical level of the point-wise test. However, because adjacent wavelet coefficients are correlated and wavelet spectra often contain many wavelet quantities, the point-wise test can produce many false positive results that occur in clusters or patches. To circumvent the point-wise test drawbacks, it is necessary to implement the recently developed area-wise, geometric, cumulative area-wise, and topological significance tests, which are reviewed and developed in this paper. To improve the computational efficiency of the cumulative area-wise test, a simplified version of the testing procedure is created based on the idea that its output is the mean of individual estimates of statistical significance calculated from the geometric test applied at a set of point-wise significance levels. Ideal examples are used to show that the geometric and cumulative area-wise tests are unable to differentiate wavelet spectral features arising from singularity-like structures from those associated with periodicities. A cumulative arc-wise test is therefore developed to strictly test for periodicities by using normalized arc length, which is defined as the number of points composing a cross section of a patch divided by the wavelet scale in question. A previously proposed topological significance test is formalized using persistent homology profiles (PHPs) measuring the number of patches and holes corresponding to the set of all point-wise significance values. Ideal examples show that the PHPs can be used to distinguish time series containing signal components from those that are purely noise. To demonstrate the practical uses of the existing and newly developed statistical methodologies, a first comprehensive wavelet analysis of Indian rainfall is also provided. An R software package has been written by the author to implement the various testing procedures.

Wind speed↗

Molecular Dynamics Simulations of Silicon Carbide, Boron Nitride and Silicon for Ceramic Matrix Composite Applications

A comprehensive computational molecular dynamics study is presented for crystalline α-SiC (6H, 4H, and 2H SiC), β-SiC (3C SiC), layered boron nitride, amorphous boron nitride and silicon, the constituent materials for high-temperature SiC/SiC compositions. Large-scale Atomic/Molecular Parallel Simulator software package was used. The Tersoff Potential force field was utilized to evaluate their mechanical characteristics of most of the materials, and the Reax force field was used to model silicon when the Tersoff Potential did not provide accurate results. Their mechanical behaviors were evaluated at a strain rate of 10(exp 7)/s and the results agree with the experimental data in the literature. The results are foundational for linking constituent behavior to composite performance, particularly when test data is unavailable or suspect.

Aluko, Olanrewaju↗

Analysis of Mercury Laser Altimeter Crossovers with Improved Mercury and MESSENGER Ephemerides

Based on previous applications of laser altimetry to planetary geodesy at GSFC [Mazarico et al. (2014),(2016)] and taking advantage of new accurate Mercury and MESSENGER orbits by [Genova et al. (2019)], we analyze altimetric crossovers from the MESSENGER (Mercury Surface, Space Environment, Geochemistry and Ranging) Laser Altimeter (MLA) to solve for orbital and geodetic parameters (e.g., rotation and orientation). We present our results based on a new Python software package recently developed at GSFC that can simulate and process altimetry data in a closed-loop. Realistic simulations of MLA data, including an appropriate range noise from the instrument and realistic terrain roughness, are performed in order to fully characterize the robustness of the solution. The simulation results are then applied to our analysis of the full dataset acquired by the MLA instrument.

Bertone, Stefano↗

Sensor Uncertainty Mitigation and Dynamic Well Clear Volumes in DAIDALUS

This paper presents recent updates to DAIDALUS (Detect & Avoid Alerting Logic for Unmanned Systems), a detect and avoid (DAA) software package for the integration of civil UAS into the airspace. DAIDALUS is the reference implementation of detect and avoid for unmanned aircraft systems chosen by RTCA Special Committee 228 (SC-228), and it is included in its corresponding Minimum Operational Performance Standards (MOPS) document, DO-365. This paper reports on the integration into DAIDALUS of two new capabilities, namely dynamic well clear volumes and sensor uncertainty mitigation.'

Narkawicz, Anthony↗

Quasi-Static and Dynamic Analysis of Composite Panels: Characterization and Validation

The predictive capability of a progressive damage model to simulate the quasi-static and dynamic crushing of composite specimens is evaluated in this article. The material model is commercially available within the ABAQUS finite element software package, i.e., “Damage for Fiber reinforced Composites”, and is often utilized to simulate damage progression and failure in static and dynamic failure scenarios such as crushing of a composite plate. The composite specimens used in this study consist of quasi-isotropic C-shaped and corrugated composite panels comprised of carbon fiber/epoxy braided fabric. The C-shaped panel specimen is used to calibrate the model and optimize material/model properties through comparison with experimental results. Then, using the same model parameters, the quasi-static crush response of the corrugated panel is predicted. The predictive capability of the model is further demonstrated by simulating the dynamic crushing response of both C-shaped and corrugated panels using the same material properties as in the quasi-static crush case with only minor adjustments for the dynamic elastic properties (i.e., rate dependent modulus). Results show that the in-built material model available in ABAQUS can successfully reproduce experimental results for both panels (C-shaped and corrugated) subjected to both quasi-static and dynamic loading scenarios. A small-scale parametric study on physical (experimentally measurable) and nonphysical (purely mathematical) model parameters was also conducted. It can be concluded that, achieving successful simulation results requires a more in-depth understanding of the influence of model parameter variation on the mechanical response of the composite together with the strategies and challenges of the utilized modeling methodology (e.g., FEA mesh density, etc.).

Naghipour Ghezeljeh, Paria↗

Numerical Investigation and Optimization of a Flushwall Injector for Scramjet Applications at Hypervelocity Flow Conditions

An investigation utilizing Reynolds-averaged simulations (RAS) was performed in order to demonstrate the use of design and analysis of computer experiments (DACE) methods in Sandia’s DAKOTA software package for surrogate modeling and optimization. These methods were applied to a flow- path fueled with an interdigitated flushwall injector suitable for scramjet applications at hyper- velocity conditions and ascending along a constant dynamic pressure flight trajectory. The flight Mach number, duct height, spanwise width, and injection angle were the design variables selected to maximize two objective functions: the thrust potential and combustion efficiency. Because the RAS of this case are computationally expensive, surrogate models are used for optimization. To build a surrogate model a RAS database is created. The sequence of the design variables comprising the database were generated using a Latin hypercube sampling (LHS) method. A methodology was also developed to automatically build geometries and generate structured grids for each design point. The ensuing RAS analysis generated the simulation database from which the two objective functions were computed using a one-dimensionalization (1D) of the three-dimensional simulation data. The data were fitted using four surrogate models: an artificial neural network (ANN), a cubic polynomial, a quadratic polynomial, and a Kriging model. Variance-based decomposition showed that both objective functions were primarily driven by changes in the duct height. Multiobjective design optimization was performed for all four surrogate models via a genetic algorithm method. Optimal solutions were obtained at the upper and lower bounds of the flight Mach number range. The Kriging model predicted an optimal solution set that exhibited high values for both objective functions. Additionally, three challenge points were selected to assess the designs on the Pareto fronts. Further sampling among the designs of the Pareto fronts may be required to lower the surrogate model errors and perform more accurate surrogate-model-based optimization.

Shenoy, Rajiv R.↗

Optimization of a Solver for Computational Materials and Structures Problems on NVIDIA Volta and AMD Instinct GPUs

The Scalable Implementation of Finite Elements by NASA (ScIFEN) is a software package developed to solve complex computational materials and structures problems using the finite element method (FEM). In this paper, we describe optimization techniques to speed up the linear solver computation that occurs within the ScIFEN application. We consider GPUs from two different vendors, NVIDIA and AMD as our target platforms for optimization and highlight differences in performance and optimization techniques. The NVIDIA GPU Volta V100 is used in the Summit system deployed at Oak Ridge National Laboratory, and the new exascale system, Frontier, will be using AMD Radeon Instinct GPU. We evaluated the performance of various optimization techniques on test matrices, ranging in size from100K to 4M, that are representative of ScIFEN applications. The linear solver computation is memory-bound on both GPUs. Our experiments show that on the NVIDIA GPU we obtained up to79%of the theoretical peak bandwidth, while the AMD GPU achieved 59%. Overall, the NVIDIA V100 GPU outperforms the AMD MI 25 GPU1. We observed an overall speedup of up to37X on an NVIDIA V100 compared to an Intel Skylake 12-coremachine. The solver for a 4M degree of freedom system took under 2.5 seconds.

Mohammad Zubair↗

Aircraft System Identification from Multisine Inputs and Frequency Responses

The identification of aircraft flight dynamics is often performed using frequency responses, which are nonparametric models that quantify the steady-state magnitude and phase of a dynamic system response to sinusoidal inputs, as a function of frequency. Frequency responses are computed from measured input and output data, and then model parameters, such as stability and control derivatives, are estimated to best fit a parametric model to the empirical frequency response data. The utility of this approach is due to the familiarity of engineers with frequency responses, a number of theoretical and practical advantages under specific conditions, the availability of software packages, and many other reasons.

System identification↗

Atomic Data for Plasma Spectroscopy: The CHIANTI Database, Improvements and Challenges

CHIANTI is an atomic database and software package for modeling emission lines and continua from hot astrophysical plasmas. It is freely available to all researchers and has been widely used in the Heliophysics and Astrophysics communities for almost 25 years. In this review, we summarize the properties of the current version of the database and give an overview of the relevant atomic processes. We also discuss progress towards a complete implementation of collisional-radiative modeling, simultaneously solving for atomic level and ion populations for individual elements.

A&M databases↗

Providing Long-term Solar and Meteorological Resource Information from NASA Research through the GIS-Enabled POWER Web Services Portal

Solar and other renewable energy systems are optimized using geophysical parameters describing information about the solar resource and meteorological environments and how those environments may be changing in time. The Prediction of Worldwide Renewable Energy Resource (POWER) team at the NASA LaRC Research Center deployed its first version of the revitalized POWER geophysical parameter website that employs Esri Geographic Information System (GIS) tools. This web application provides access to both time series and climatological data sets spanning from a few days behind real-time back to the early 1980’s with a spatial resolution of 0.5 x 0.5 degree, thus extending over 30 years. The GIS tools enable, generate and store climatological averages using spatial queries and calculations in a spatial database resulting in greater accessibility to government agencies, industry and individuals. There is an API that provides data through a URL coding and can also invoked inside user software packages. Output data formats now include ASCII, CSV, geoTiff, JSON and .netCDF. Since May 2018, over 10.4 TB of data have been delivered to fulfill over 110 million data requests from 240,000 unique user IPs (https://power.larc.nasa.gov). This presentation provides an overview of this project and the current version of the new POWER web capabilities through to the end usage. Surface solar parameters are adapted from both the NASA/GEWEX (Global and Energy Water Cycle Exchange) SRB (Surface Radiation Budget) project and the CERES (Clouds and Earth’s Radiant Energy System) data products. The meteorological parameters are adapted from NASA’s MERRA-2 (Modern Era Retrospective-analysis for Research and Applications). We review uncertainty of various basic parameters using surface measurements. After the introduction, we specifically discuss various clusters of parameters for equator pointing tilted surfaces that provide for the estimation of optimal tilt angle by month and year. The discussion will emphasize a few examples of how solar and building engineers are using the data products. We will then preview the capabilities of Version 2 featuring examples of hourly solar and meteorological data products, customized user reports and expanded web mapping services from CERES data products. The web services provide a unique, expandable resource for renewable energy systems engineers to design, evaluate and optimize to environments worldwide.

Solar Resource↗

Implementation of Charged Particle Behavior in Discrete Element Method (DEM) Simulations

To understand the behavior of charged lunar regolith when perturbed by lunar landers, it is important to couple the grain dynamics with mechanical and electrical particle interactions. To accomplish this, improvements have been made to a discrete element method (DEM) software package to include both short- and long-range interactions between spherical particles. Short-range interactions rely on contact between the particles, such as electrical conduction and triboelectric charge transfer driven by work functions. Long-range interactions act at a distance between every pairing of particles, such as electrostatic forces and gravitational forces. Results from simulations between a few particles are compared with theory to verify these added behaviors prior to scaling up to more complex scenarios. The radii, initial charges, electrical conductivities, work functions, and separation of the particles are varied and the resultant charges as well as the time required to reach the final state are determined.

Electrostatics↗

Effects of Spin-Orbit Resonances and Tidal Heating on the Inner Edge of the Habitable Zone

Much attention has been given to the climate dynamics and habitable boundaries of synchronously rotating planets around low mass stars. However, other rotational states are possible, including spin–orbit resonant configurations, particularly when higher eccentricity orbits can be maintained in a system. Additionally, the oscillating strain as a planet moves from periastron to apoastron results in friction and tidal heating, which can be an important energy source. Here, we simulate the climate of ocean-covered planets near the inner edge of the habitable zone around M to solar stars with the NASA GISS ROCKE-3D general circulation model, and leverage the planetary evolution software package, VPLanet, to calculate tidal heating rates for Earth-sized planets orbiting 2600 and 3000 K stars. This study is the first to use a 3D general circulation model that implements tidal heating to investigate habitability for multiple resonant states. We find that for reference experiments without tidal heating, the resonant state has little impact on the radial position of the inner edge because for a given stellar flux, higher-order states tend to be warmer than synchronous rotators, but for a given temperature, have drier upper atmospheres. However, when strong tidal heating is present, the rotational component implies a strong dependence of habitable conditions on the system evolution and rotational state. Since tidal and stellar heating both decrease rapidly with orbital distance, this results in a compact orbital width separating temperate and uninhabitable climates. We summarize these results and also compare ROCKE-3D to previously published simulations of the inner edge.

Spin-Orbit Resonances↗

Thermodynamic constraints on the textural evolution of eucrite, EET 90020

Basaltic eucrites, which are possibly parts of the Vestan crust, formed as either lava flows or intrusions, and can offer insights into early crust formation. However most eucrites have experienced thermal metamorphism, which exacerbates the challenges of understanding these samples. A variety of petrogenetic models such as, partial melting of a primitive source, fractional crystallization of magmas emplaced in a primitive crust, and/or partial melting of a eucritic source coupled with melt mixing and assimilation have been proposed in order to explain observed petrologic and chemical characteristics of eucrites. EET 90020 is an unbrecciated eucrite that has experienced significant thermal metamorphism, with temperatures of metamorphic equilibration calculated to be ~840-1040°C. Its petrologic history remains contentious, in part due to differences in trace element analyses [10,17, this work], which has led to multiple metamorphic interpretations [8,10-11]. Here, we combined petrologic observations, chemical analyses, and thermodynamic modeling, to interpret micro-domain textures identified in EET 90020 and better constrain its petrologic history. Ultimately, the development of such textures is a direct result of the geologic processes operating on a young Vesta or similar asteroid. Sample description: We identified three textural domains in EET 90020 (Figure 1); (I) A coarse grain domain with granoblastic plagioclase and pyroxene, (II) a medium grain domain with curved grain boundaries between plagioclase and pyroxene, and (III) a fine grain domain dominated by tridymite that fills in interstitial space between spherical plagioclase. Major element chemistry in phases is consistent across domains. Pyroxene have distinctive high Ca lamellae (Wo40.92En28.1Fs31.0) forming from the low-Ca pigeonite host (Wo3.5En33.7Fs62.8). There are small amounts of fayalitic olivine along pyroxene/oxide grain boundaries, where pyroxenes have reacted with ilmenite/chromite clasts. Plagioclase is anorthitic (Avg ~An88), with little variation (An86-92.5). Trace element bulk rock analyses were collected via ICP-MS and individual mineral analyses with LA-ICP-MS (Figure 2). The bulk sample has an observable enrichment in light rare earth elements, with a notable depletion in Eu. Plagioclase is enriched in LREEs, depleted in HREEs, and has a positive Eu anomaly while pyroxene is enriched in HREEs with a negative Eu anomaly. Variations in trace element analyses is likely influenced by the presence of phosphates because are highly enriched in REEs and have a dramatic effect on the bulk domain compositions. Thermodynamic modeling: The software package, Perple_X, was used to calculate mineral phase equilibria over a range of conditions using a Gibbs free energy minimization approach [12]. Models used all major/minor elements except P2O5. Thermodynamic properties from [13] were used to constrain endmember phase stabilities. Activity models were used to describe mixing in phases with solid solution. Isochemical P-T phase diagrams were constructed for the bulk thin section, and the coarse, medium and fine grain domain compositions. Domain compositions were determined using microprobe analyses of phases identified in each domain and observed modal mineral abundancies. Thermodynamically stable pyroxene compositions were extracted from model results and compared to pyroxene compositions collected via EMPA, allowing us to calculate temperatures for metamorphic equilibrium. Results: Results that use the coarse domain composition (Fig. 3a) demonstrate overlap at T ~ 1020°C indicating high T equilibration. There was no overlap for the medium and fine grain domains, or the bulk composition (e.g., Fig. 3b) indicating that these are not equilibrium assemblages. These results were replicated using data from [8,20]. Discussion: The disequilibrium results from the bulk rock, fine grain and medium grain models imply that EET 90020 was modified during and/or after peak metamorphism (open system behavior). However, this is in contrast to the coarse domain where it appears as if equilibrium was maintained. We suggest that this conflict is explained by highly localized equilibrium occurring at millimeter length-scales, which has been observed in terrestrial metamorphism where aqueous fluid infiltrate [14] or melt loss occurs [15]. It is unlikely that aqueous alteration was significant in modifying the bulk composition of EET 90020, because it lacks hydrated minerals. However, the presence of melt is supported by the presence of an abundance of curved grain boundaries and spherical grains in the medium and fine grain domains This contrasts the coarse domain, where grain boundaries are angular. Given the evidence for disequilibrium in the medium and fine grain domains, it would be inappropriate to draw additional conclusions about textural evolution using thermodynamic models of those domains. However, modeling results of the coarse grain domain, where equilibrium was maintained, can provide useful insights into the geologic evolution of EET 90020. We approximate maximum metamorphic temperatures, using the Ca component in pyroxene, at ~1020°C (Fig. 4), which is consistent with previous two-pyroxene thermometry [11]. At this temperature, up to 10 vol % melt can be produced in the coarse grain domain, which is at the boundary of minimum melt needed in order to segregate and form a melt network [16]. Tridymite and plagioclase would be the first minerals to melt out during heating starting ~990°C. This is consistent with the observed melt textures formed by tridymite and plagioclase in the fine grain domain. We suggest that a partial melting model best explains the textural development of EET 90020. Upon heating, melt generated was either trapped, resulting in no net change in the localized rock composition, or migrated, resulting in disequilibrium between pyroxene and the surrounding matrix. A positive Eu anomaly in plagioclase indicates that all of the plagioclase was melted previously. This is consistent with previous studies that have suggested that EET 90020 represents a residual eucrite that experienced partial melting and subsequent melt loss [8]. We speculate that partial melting could have occurred when EET 90020 was heated in the asteroid’s lower crust or adjacent to an intrusion [7]. Future work will study other, texturally heterogeneous samples to see if EET 90020 presents a unique circumstance or whether crustal partial melting was ubiquitous on the eucrite asteroid.

Eucrites↗

Sensitivity of Multispectral Imager Liquid Water Cloud Microphysical Retrievals to the Index of Refraction

A cloud property retrieved from multispectral imagers having spectral channels in the shortwave infrared (SWIR) and/or midwave infrared (MWIR) is the effective particle radius (CER), a radiatively relevant weighting of the cloud particle size distribution. The physical basis of the CER retrieval is the dependence of SWIR/MWIR cloud reflectance on the cloud particle single scattering absorption, which in turn depends on the complex index of refraction of bulk liquid water (or ice) in addition to the cloud particle size. There is a general consistency in the choice of the liquid water index of refraction by the cloud remote sensing community, largely due to the few available independent datasets and compilations. Here we examine the sensitivity of CER retrievals to the available laboratory index of refraction datasets in the SWIR and MWIR using the retrieval software package that produces NASA’s standard MODIS/VIIRS continuity cloud products. The sensitivity study incorporates two laboratory index of refraction datasets that include measurements at supercooled water temperatures, one in the SWIR [Kou et al., 1993] and one in the MWIR [Wagner et al., 2005]. Neither has been broadly utilized in the cloud remote sensing community. It is shown that these two new datasets can significantly change CER retrievals (e.g., 1-2 μm) relative to common datasets used by the community. Further, index of refraction data for a 265 K water temperature results in more consistent retrievals between the two spectrally distinct 2.2 μm atmospheric window channels on MODIS and VIIRS. As a result, the 265 K values from the SWIR and MWIR index of refraction datasets were adopted for use in the production version of the continuity cloud product. The results indicate the need to better understand temperature-dependent bulk water absorption and uncertainties in these spectral regions.

MODIS↗

Climate Ocean Modeling on a Beowulf Class System

With the growing power and shrinking cost of personal computers, the availability of fast ethernet interconnections, and public domain software packages, it is now possible to combine them to build desktop parallel computers at a fraction of what it would cost to buy systems of comparable power from supercomputer companies.

Beowulf↗

Integrated Modeling Methodology Validation Using the Micro-Precision Interferometer Testbed

This paper validates the integrated modeling methodology used for design and performance evaluation of complex opto-mechanical systems, particularly spaceborne interferometers. The methodology integrates structural modeling, optical modeling, and control system design into a common environment, the Integrated Modeling of Optical Systems (IMOS) software package.

opto-mechanical↗

Modeling of Guided Waves for Aerospace Applications

Advancements in computer hardware has led to new possibilities for rapid modeling and simulation capabilities across many scientific fields. Nondestructive evaluation (NDE) can benefit from increased use of simulation tools to guide optimization of inspection and health monitoring methods, enhance understanding of data, aid in development of defect characterization methods, and generate data sets for use with machine learning and model-assisted probability of detection. Recent work at NASA has entailed development and benchmarking of both custom simulation codes and commercial simulation tools for ultrasonic wave propagation. This paper describes recent work at NASA in modeling of guided waves in composites and other aerospace materials. Results and computational speeds for a composite benchmark case are reported for a custom finite difference Rotated Staggered Grid code and for the commercial finite element software package, Pogo. Recent progress in linking NDE models to parametric analysis tools is also discussed.

Nondestructive evaluation↗