Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity formulations”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Sensitivity of the ACCP Value Framework Science Benefit Score

NASA’s priorities for Earth Science are informed by the 2017-2027 Decadal Survey for Earth Science and Applications from Space of the National Academies of Sciences, Engineering and Medicine. In that document, five Designated Observables are identified as priorities for implementation: Aerosols, Cloud-Convection-Precipitation, Surface Biology, Mass Change, and Surface Deformation Change. In 2019, NASA initiated four multi-year studies to formulate, assess, and recommend potential Observing System concepts that are responsive to the recommendations of the Earth Science Decadal Survey. The Aerosols, Clouds, Convection, and Precipitation (ACCP) study investigates opportunities for synergistic science resulting from implementations that combine both Aerosols science objectives and Clouds, Convection, and Precipitation science objectives. For the ACCP study, a Value Framework was designed to introduce structure, transparency, and traceability in the mission down-select decision process. This paper presents a sensitivity analysis that was conducted to characterize the behavior of the science benefit score of the ACCP Value Framework. it includes derivations of the sensitivity analysis equations, numerical analyses of the behavior of the combined benefit score, and the conclusions drawn from the results of the work that was conducted. The results of this analysis were key to ensuring that the aggregated science benefit scores would enable the team to differentiate and contrast across concepts.

John McLaughlan↗
DOE OSTI2021

Thermal-Mechanical Elastic-Plastic and Ductile Failure Model Calibrations for 6061-T651 Aluminum Alloy from Plate

Numerical simulations of metallic structures undergoing rapid loading into the plastic range require material models that accurately represent the response. In general, the material response can be seen as having four interrelated parts: the baseline response under slow loading, the effect of strain rate, the conversion of plastic work into heat and the effect of temperature. In essence, the material behaves in a thermal-mechanical manner if the loading is fast enough so when heat is generated by plastic deformation it raises the temperature and therefore influences the mechanical response. In these cases, appropriate models that can capture the aspects listed above are necessary. The matters of interest here are the elastic-plastic response and ductile failure behavior of 6061-T651 aluminum alloy under the conditions described above. The work was accomplished by first designing and conducting a material test program to provide data for the calibration of a modular $J_2$ plasticity model with isotropic hardening as well as a ductile failure model. Both included modules that accounted for temperature and strain rate dependence. The models were coupled with an adiabatic heating module to calculate the temperature rise due to the conversion of plastic work to heat. The test program included uniaxial tension tests conducted at room temperature, 150 and 300 C and at strain rates between 10 –4 and 10 3 1/s as well as four geometries of notched tension specimens and two tests on specimens with shear-dominated deformations. The test data collected allowed the calibration of both the plasticity and the ductile failure models. Most test specimens were extracted from a single piece of plate to maintain consistency. Notched tension tests came from a possibly different plate, but from the same lot. When using the model in structural finite element calculations, element formulations and sizes different from those used to model the test specimens in the calibration are likely to be used. A brief investigation demonstrated that the failure model can be particularly sensitive to the element selection and provided an initial guide to compensate in a specific example.

36 MATERIALS SCIENCE↗
NASA NTRS1977

Application of the endochronic theory of viscoplasticity to solid propellants and sandasphalt concrete

Solid propellants, sand-asphalt concrete and hard plastics showed rate sensitive mechanical behavior which, in addition, indicated that these materials have a permanent memory of the strain (or loading) path by which their present state was attained. A constitutive equation was formulated in general three dimensional tensorial form by means of irreversible thermodynamics. By using a very simple analytical form, it was shown that the mechanical behavior of solid propellants and sand-asphalt concrete can be readily described.

Peng, S. T. J.↗
NASA NTRS1979

Recent progress in finite-volume calculations for wing-fuselage combinations

Progress in the application of finite-volume methods to the calculation of transonic potential flows past general wing-body combinations is reviewed. Two different methods of generating boundary-conforming grids are investigated, and the results compared to provide an estimate of solution sensitivity to grid geometry. Both conservative and quasi-conservative difference schemes are used in one of the coordinate systems. Results show that the error introduced by the quasi-conservative formulation seems to be small, although a one-dimensional analysis suggests that schemes of this type do not necessarily produce mass-conserving shocks. Comparison of calculated results with experimental data for realistic fuselage geometries clearly shows the importance of modelling the effect of fuselage geometry upon the wing pressure distribution.

Caughey, D. A.↗
NASA NTRS1994

Monte Carlo prediction of radiative heat transfer in inhomogeneous, anisotropic, nongray media

A Monte Carlo solution technique has been formulated to predict the radiative heat transfer in three-dimensional, inhomogeneous participating media which exhibit spectrally dependent emission and absorption and anisotropic scattering. Details of the technique and selected numerical sensitivities are discussed. The technique was applied to a problem involving a medium composed of a gas mixture of carbon dioxide and nitrogen and suspended carbon particles. A homogeneous medium was modeled to examine the effect of total pressure and carbon-particle concentration on radiative heat transfer. Variation in total pressure, over the range studied, had minimal effect on the amount of heat radiated to the enclosure walls and on the radiative-flux distribution within the medium. Increases in the carbon particle concentration produced significantly higher heat fluxes at the boundaries and altered the radiative flux distribution. The technique was then applied to an inhomogeneous medium to examine effects of specific temperature and carbon particle concentration distributions on radiative heat transfer. For the inhomogeneous conditions examined, the largest radiative flux divergence occurs near the center of the medium and the regions near some enclosure walls act as energy sinks.

HEAT FLUXES↗
NASA NTRSDate not supplied

Open-source Numerical Modeling of Solidification Cracking Susceptibility: Application to Refractory Alloy Systems

Introduction. Alloys such as aluminum, nickel-base, and austenitic stainless steels are susceptible to solidification cracking during welding and 3D printing. Compositional optimization is one method used to effectively mitigate solidification cracking of those alloy systems. With the surge in hypersonic and in-space propulsion activities, refractory metals (Nb, Mo, Ta, W, and Re) and their alloy derivatives are increasing in importance due to their extreme high melting point and retention of high-temperature strength; however, their chemistry was most typically optimized to promote ductility during mechanical operations such as drawing and forming. Welding of such alloys has been a challenge due to a number of issues including solidification cracking, atmospheric contamination (O, C, and N), as well as a shift in ductile-to-brittle transition to higher temperature following grain growth induced by welding. Compositional optimization of refractory alloys for solidification cracking resistance in particular is desirable as their usage increases with the advent of advanced manufacturing methods such as 3D printing. This work evaluates the effect of compositional variation in refractory metal systems on the solidification cracking susceptibility with the goals of optimizing existing alloys and joining process techniques, and formulating new alloys with increased solidification cracking resistance. Experimental Procedures. A python code was developed in a Jupyter notebook environment (Michael and Sowards, 2023) to facilitate the calculation of crack susceptibility index proposed by Kou (2015). Composition is entered as a single point, or as a 1-D or 2-D array. The notebook calls pycalphad (Otis and Liu, 2017 and Bocklund et al, 2020) to calculate the evolution of fraction solid as a function of temperature (under either Scheil or equilibrium assumptions) and then evaluates steepness of the fraction solid curve near the terminal stage of solidification to predict solidification cracking resistance. Open source thermodynamic databases available at online repositories are used (van de Walle). The process is setup in an automated fashion to generate plots that show variation in solidification cracking susceptibility according to composition on 1-D line plots or 2-D contour plots. The Jupyter notebook and crack susceptibility algorithm was also integrated with a widely used commercial CALPHAD code for validation and alloy exploration. Results and Discussion. The crack susceptibility model was first validated against a series of refractory alloy compositions evaluated in past work which utilized a specialized Varestraint test built inside a vacuum chamber environment (Lessman and Gold, 1971). The alloys tested in the Varestraint apparatus included T-111 (Ta-8W-2Hf), ASTAR-811C (Ta-8W-1Re-0.7Hf-0.025C), FS-85 (Nb-27Ta-10W-1Zr), T-222 (Ta-9.6W-2.4Hf-0.01C), Ta-10W, B-66 (Nb-5Mo-5V-1Zr), and SCb-291 (Nb-10W-10Ta). The initial test of the model showed a strong correlation with empirical Varestraint data, i.e., a Spearman rank correlation between model predictions and hot cracking measurements was observed to be greater than 0.8. Following the validation, a set of refractory metal binary mixtures was investigated to evaluate sensitivity of Nb, Mo, W, and Ta to C, N, and O content. A series of plots were produced that suggest ppmw ranges of C, N, and O where solidification cracking increases significantly and reaches a maximum. Also comparative ranking of each primary refractory metal to each interstitial was produced. For example C produces greater cracking response in Mo whereas O produces greater cracking response in Ta and Nb. Such compositional values have utility in setting limits on pickup of these interstitial elements during welding and printing rather than using a one-size-fits-all approach. Furthermore, the results have use in determining additive powder recycling requirements, which is especially pertinent for refractory metal powders due to their high cost compared to conventional alloys. Another application created thousands of hypothetical alloys within the nominal specified composition range of two widely used refractory alloys C103 (Nb-10Hf-1Ti) and TZM (Mo-0.5Ti-0.1Zr). The cracking index was calculated for the alloys and results were fed into machine learning regression techniques including Multiple Linear Regression, Ridge Regression, and Lasso Regression to determine relative potency each alloying element had on computed solidification cracking index. A series of linear equations were produced that relate composition of C103 and TZM to solidification cracking index. The crack susceptibility of C103 for example is described by an equation of the form: cracking index ~ O + 0.667*C + 0.635*N + 0.00037*Ta – 0.0008*Hf (in wt.%) From that equation, it is clear that O has strong propensity to induce solidification cracking. Interestingly, Hf is shown to reduce calculated cracking response. Finally, realizing the potential of this method to discover new refractory alloy formulations across the period table that have low solidification cracking sensitivity, the code was applied to new untested alloy systems including W-Zr-C, W-Ta-C, and others. Conclusions. In summary, an open source numerical method has been developed using Python code to calculate Kou’s crack susceptibility index. The method was applied to refractory metals which are inherently difficult to study from a weldability testing standpoint since inert shielding gas is not sufficient and welding is typically done in vacuum, especially in light of findings presented here where oxygen has profound influence on solidification cracking. This work revealed the effect of compositional variations on a series of refractory metals and showed the framework defined here will be useful in 1) the development of new alloys that have improved weldability and 3D printability, 2) placing compositional limits on existing alloys, and 3) ensuring adequate controls of manufacturing processes such as 3D printing where powder reuse is critical. Keywords. pycalphad; Python; refractory metals; solidification cracking. References. B. Bocklund et. al. (2020) http://doi.org/10.5281/zenodo.3630657. S. Kou. (2015) https://doi.org/10.1016/j.actamat.2015.01.034. G.G. Lessmann and R.E. Gold. Welding Journal, issue 1, pp. 1-s – 8-s (1971). F.N. Michael and J.W. Sowards. NASA/TM-20230002218 (2023). R. Otis and Z.-K. Liu. (2017) http://doi.org/10.5334/jors.140. A. Van de Wallle et. al. (2018) https://doi.org/10.1016/j.calphad.2018.04.003.

pycalphad↗
DOE OSTI2026

Benchmarking Cylindrical Blast Wave Theory Against the OSIRIS-REx Sample Return Capsule Reentry

Weak shock theory based on cylindrical blast waves has been used to interpret meteor infrasound, but it has not been systematically benchmarked against a non-ablating hypersonic source with independently known parameters. The objective of this study is not to propose a new theoretical framework, but to evaluate the operational validity of the existing suite of blast radius formulations against a high-fidelity ground truth dataset. The OSIRIS-REx Sample Return Capsule reentry on 24 September 2023 provides such a benchmark because the capsule geometry, trajectory, and infrasound emission points are constrained from mission data and ray tracing, reducing source-side uncertainty associated with ablation. Using observations from 39 infrasound stations, this benchmarking study evaluates six published blast radius (${R}_{0}$ ) formulations and three weak-shock transition coefficients ( C ) within a stratified atmospheric propagation model to predict signal period and peak overpressure. The benchmarking identifies the Sakurai formulation as the best-performing formulation for non-ablating bodies, with the Jones/Plooster formulation performing comparably when a physically appropriate C is adopted. Sakurai and Jones/Plooster yield linear-period median absolute percentage residuals of 9% and 11%, respectively. The period predictions show only weak sensitivity to C at these propagation distances. The Mach-diameter approximation commonly used in meteor studies overestimates ${R}_{0}$ by more than a factor of 3 in the absence of ablation. Finally, these results establish a performance baseline for applying cylindrical blast wave theory to effectively non-ablating hypersonic bodies and demonstrate that the signal period is a robust observable for constraining ${R}_{0}$.

atmospheric entry↗
NASA NTRS1991

Near-wall modeling of compressible turbulent flow

A near-wall two-equation model for compressible flows is proposed. The model is formulated by relaxing the assumption of dynamic field similarity between compressible and incompressible flows. A postulate is made to justify the extension of incompressible models to ammount for compressibility effects. This requires formulation the turbulent kinetic energy equation in a form similar to its incompressible counterpart. As a result, the compressible dissipation function has to be split into a solenoidal part, which is not sensitive to changes of compressibility indicators, and a dilatational part, which is directly affected by these changes. A model with an explicit dependence on the turbulent Mach number is proposed for the dilatational dissipation rate.

So, Ronald M. C.↗
NASA NTRS1995

Driving Force Analysis in an Infinite Anisotropic Plate with Multiple Crack Interactions

The methodology and a rigorous solution formulation are presented for stress intensity factors (SIF's, k) and total strain energy release rates (SERR, G(sub T)) of a multicracked plate, that has fully interacting cracks and is subjected to a far-field arbitrary stress state. The fundamental perturbation problem is derived, and the steps needed to formulate the system of singular integral equations whose solution gives rise to the evaluation of the SIF's are identified. Parametric studies are conducted for two, three, and four crack problems. The sensitivity and characteristics of the model is demonstrated.

Binienda, Wieslaw K.↗
NASA NTRS1995

Thermomechanical postbuckling of multilayered composite panels with cutouts

The results of a study of the detailed thermomechanical postbuckling response characteristics of flat unstiffened composite panels with central circular cutouts are presented. The panels are subjected to combined temperature changes and applied edge loading (or edge displacements). The analysis is based on a first-order shear deformation plate theory. A mixed formulation is used with the fundamental unknowns consisting of the generalized displacements and the stress resultants of the plate. The postbuckling displacements, transverse shear stresses, transverse shear strain energy density, and their sensitivity coefficients are evaluated. The sensitivity coefficients measure the sensitivity of the post-buckling response to variations in the different lamination and material parameters of the panel. Numerical results are presented showing the effects of the variations in the hole diameter, laminate stacking sequence, fiber orientation, and aspect ratio of the panel on the thermomechanical postbuckling response and its sensitivity to changes in panel parameters.

Noor, Ahmed K.↗
NASA NTRS1997

Evolutionary Design of Controlled Structures

Basic physical concepts of structural delay and transmissibility are provided for simple rod and beam structures. Investigations show the sensitivity of these concepts to differing controlled-structures variables, and to rational system modeling effects. An evolutionary controls/structures design method is developed. The basis of the method is an accurate model formulation for dynamic compensator optimization and Genetic Algorithm based updating of sensor/actuator placement and structural attributes. One and three dimensional examples from the literature are used to validate the method. Frequency domain interpretation of these controlled structure systems provide physical insight as to how the objective is optimized and consequently what is important in the objective. Several disturbance rejection type controls-structures systems are optimized for a stellar interferometer spacecraft application. The interferometric designs include closed loop tracking optics. Designs are generated for differing structural aspect ratios, differing disturbance attributes, and differing sensor selections. Physical limitations in achieving performance are given in terms of average system transfer function gains and system phase loss. A spacecraft-like optical interferometry system is investigated experimentally over several different optimized controlled structures configurations. Configurations represent common and not-so-common approaches to mitigating pathlength errors induced by disturbances of two different spectra. Results show that an optimized controlled structure for low frequency broadband disturbances achieves modest performance gains over a mass equivalent regular structure, while an optimized structure for high frequency narrow band disturbances is four times better in terms of root-mean-square pathlength. These results are predictable given the nature of the physical system and the optimization design variables. Fundamental limits on controlled performance are discussed based on the measured and fit average system transfer function gains and system phase loss.

Masters, Brett P.↗
NASA NTRS2003

Estimating Terrestrial Carbon Exchange from Space: How Often and How Well?

Data from a new space mission measuring integrated light-use efficiency could provide a breakthrough in understanding of global carbon, water, and energy dynamics, and greatly improve the accuracy of model predictions for terrestrial carbon cycles and climate. Over the past decade, Gamon and others have shown that changes in photo-protective pigments are sensitive indicators of declines in light-use efficiency of plants and plant canopies. The requirements for integrated diurnal measurements from space need to be defined, before a space mission can be formulated successfully using this concept. We used towerbased CO 2 flux data as idealized proxies for remote measurements, examining their sampling properties. Thousands of half-hourly CO 2 flux measurements are needed before their average begins to converge on an average annual net CO 2 exchange. Estimates of daily integrated fluxes (i.e., diurnal curves) are more statistically efficient, especially if the spacing between measured days is quasiregular, rather than random. Using a few measurements per day one can distinguish among days with different net CO 2 exchanges. Fluxes sampled between mid-morning to mid-afternoon are more diagnostic than early morning or late afternoon measurements. Similar results (correlation >0.935) were obtained using 2 measurements per day with high accuracy ([:plusmn:]5%), 3 measurements per day with medium accuracy ([:plusmn:] 10%), or 5 measurements per day at lower accuracy ([:plusmn:]20%). An observatory in a geosynchronous or near-geosynchronous orbit could provide appropriate observations, as could a multi-satellite constellation in polar orbits, but there is a potential trade-off between the required number of observations per day and quality of each observation.

Knox, Robert G.↗
NASA NTRS1996

A Comparison of Computational Aeroacoustic Prediction Methods for Transonic Rotor Noise

This paper compares two methods for predicting transonic rotor noise for helicopters in hover and forward flight. Both methods rely on a computational fluid dynamics (CFD) solution as input to predict the acoustic near and far fields. For this work, the same full-potential rotor code has been used to compute the CFD solution for both acoustic methods. The first method employs the acoustic analogy as embodied in the Ffowcs Williams-Hawkings (FW-H) equation, including the quadrupole term. The second method uses a rotating Kirchhoff formulation. Computed results from both methods are compared with one other and with experimental data for both hover and advancing rotor cases. The results are quite good for all cases tested. The sensitivity of both methods to CFD grid resolution and to the choice of the integration surface/volume is investigated. The computational requirements of both methods are comparable; in both cases these requirements are much less than the requirements for the CFD solution.

Brentner, Kenneth S.↗
DOE OSTI2024

Diffeomorphism invariance and quantum mechanical paradoxes

Paradoxes in gravitational physics, such as grandfather paradoxes with closed timelike curves or the AMPS paradox, are often constructed in a weak gravity regime but upon further examination engender strong gravitational responses such as unstable Cauchy horizons or firewalls. In contrast, some proposed paradoxes in nonrelativistic quantum mechanics ignore gravity completely. Such nonrelativistic proposed paradoxes are often not gauge invariant — they use local operators which should be forbidden by diffeomorphism-invariant quantum gravity. Ignoring this complication is reasonable if there is a consistent weak gravity regime where the paradox can be formulated with gauge invariant observables up to some order in Newton’s constant. Here, we show that this approach remains inconsistent due to a lack of analyticity of the necessary solutions. As a consequence, there is no consistent weak gravity, diffeomorphism invariant embedding of such paradoxes—just as in other gravitational paradoxes they generate a strong gravitational response and are in principle sensitive to quantum gravity.

quantum gravity↗
NASA NTRS1982

On-orbit propulsion requirements and performance assessment of ion propulsion subsystems for future GEO large satellite missions

It is pointed out that a large class of future spacecraft, referred to as large space structures (LSS), will require advanced stationkeeping thrust subsystems. The present investigation is concerned with the performance requirements of such advanced stationkeeping thrust subsystems. An analytical model is employed to evaluate the sensitivity of the total stationkeeping thrust-system mass, for geosynchronous spacecraft, to variations in the mission parameters and thrust-system performance, taking into account ion propulsion thrust systems. The model is formulated for geosynchronous missions and considers only the N-S orbit perturbations due to sun-moon forces and the E-W orbit perturbations due to solar pressure, since these are the dominant orbit perturbations of LSS missions. The thrust system is scaled in relation to the performance capabilities of the existing 8- and 30-cm diameter ion-thruster technology and consists of separate N-S and E-W thrusters.

Wessel, F. J.↗
NASA NTRS1990

Laser Atmospheric Wind Sounder (LAWS) phase 1. Volume 2

This report summarizes and documents the results of the 12-month phase 1 work effort. The objective of phase 1 was to establish the conceptional definition of the laser atmospheric wind sounder (LAWS) sensor system, including accommodations analyses to ensure compatibility with the Space Station Freedom (SSF) and the Earth Observing System (EOS) Polar Orbiting Platform (POP). Various concepts were investigated with trade studies performed to select the configuration to be carried forward to the phase 2 Preliminary Design Definition. A summary of the LAWS system and subsystem trade studies that were performed leading to the baseline design configuration is presented in the appendix. The overall objective of the LAWS Project is to define, design, and implement an operational space based facility, LAWS, for accurate measurement of Earth wind profiles. Phase 1 addressed three major areas: (1) requirements definition; (2) instrument concepts and configurations; and (3) performance analysis. For the LAWS instrument concepts and configurations, the issues which press the technological state of the art are reliable detector lifetime and laser performance and lifetime. Lag angle compensation, pointing accuracy, satellite navigation, and telescope design are significant technical issues, but they are considered to be currently state of the art. The primary issues for performance analysis concern interaction with the atmosphere in terms of backscatter and attenuation, wind variance, and cloud blockage. The phase 1 tasks were formulated to address these significant technical issues and demonstrate the technical feasibility of the LAWS concept. Primary emphasis was placed on analysis/trade and identification of candidate concepts. Promising configurations were evaluated for performance, sensitivities, risks, and budgetary costs. Lockheed's baseline LAWS configuration is presented.

Source record↗
NASA NTRS1992

Software issues in three dimensional continuum shape optimization employing boundary formulations

The paper addresses the issue of how individual computational techniques for the accurate and economical calculation of information required for large-scale 3D continuum structural shape optimization via boundary element analysis (BEA) formulations can be prudently selected and incorporated in large-scale BEA programs, employing either direct or iterative equation solvers. A series of techniques that allows for the incorporation of economical shape design sensitivity analysis capability with a minimal amount of software development is described. The use of reanalysis is shown to be the key ingredient associated with each of these techniques. It is concluded that, from a software engineering perspective, capabilities that facilitate 3D shape optimization can be implemented in BEA programs with only modest investments in program development.

Kane, James H.↗
NASA NTRS2001

Development and Lab-Scale Testing of a Gas Generator Hybrid Fuel in Support of the Hydrogen Peroxide Hybrid Upper Stage Program

As part of a NASA funded contract to develop and demonstrate a gas generator cycle hybrid rocket motor for upper stage space motor applications, the development and demonstration of a low sensitivity, high performance fuel composition was undertaken. The ultimate goal of the development program was to demonstrate successful hybrid operation (start, stop, throttling) of the fuel with high concentration (90+%) hydrogen peroxide. The formulation development and lab-scale testing of a simple DOT Class 1.4c gas generator propellant is described. Both forward injected center perforated and aft injected end burner hybrid combustion behavior were evaluated with gaseous oxygen and catalytically decomposed 90% hydrogen peroxide. Cross flow and static environments were found to yield profoundly different combustion behaviors, which were further governed by binder type, oxidizer level and, significantly, oxidizer particle size. Primary extinguishment was accomplished via manipulation of PDL behavior and oxidizer turndown, which is enhanced with the hydrogen peroxide system. Laboratory scale combustor results compared very well with 11-inch and 24-inch sub-scale test results with 90% hydrogen peroxide.

Lund, Gary K.↗