Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reference call set”

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.

115 records · Page 7

Thermal Considerations for 2039 Opposition Class Nuclear Electric Propulsion/Chemical Propulsion Crewed Mars Mission

The high specific impulse (Isp) of Nuclear Electric Propulsion (NEP) technology offers the potential for advanced space mission capabilities. However, the five critical technology elements of NEP vehicles have yet to prove technical maturity levels for consideration into mission design. In response to the critical reviews by the NASA Engineering and Safety Center (NESC) and the National Academies of Sciences, Engineering, and Medicine (NASEM), NASA’s Space Nuclear Propulsion (SNP) project created an NEP Technology Maturation Plan (TMP) for focused development of NEP technology. The TMP called for a coordinated set of technology development efforts to meet this objective. The Modular Assembled Radiators for NEP VehicLes (MARVL) Early Career Initiative (ECI) project was initiated to develop a portion of the fifth Critical Technology Element (CTE) of the NEP vehicle: the Primary Heat Rejection Subsystem (PHRS). A target application of a 2039 human-rated Mars mission was outlined in the TMP. For the outlined mission, a NEP vehicle will experience several thermal environments which will impact the design and operation of the PHRS. To maintain radiator temperatures within the required effective temperature range, the effect of natural, induced, and NEP internally generated heat loads on the radiator panel must be well understood. Furthermore, this analysis is critical for analyzing the influence of various orientations and positions of the NEP vehicle relative to nearby celestial bodies throughout the mission. This study conducted a complete enveloping analysis of the thermal environments influencing the NEP vehicle throughout the mission. Thermal analysis was conducted for the radiator panels based on the defined mission environments. This thermal analysis concludes with the selection of ideal radiator orientations for the NEP vehicle, and the identification of worst case hot and cold environmental sink temperatures throughout the mission. For the target application, the environmental sink temperature while the reactor is powered OFF or powered ON ranges from 30 K to 353 K and 2.7 K to 243 K respectively. When considering interplanetary space, the minimum environmental sink temperature when the reactor is powered OFF and the radiators are oriented “edge to Sun” is 2.7 K. The environmental thermal models generated in this study will be used for future studies with the full vehicle system model. The environmental sink temperature curves generated will be used for future radiator and component analysis to predict transient performance in the space environment. The environmental sink temperature and heat rejection capability curves will inform the trade between commissioning orbits that are in consideration. The model may also serve as a useful tool as reference for future crewed space missions, missions involving radiators or temperature sensitive equipment, or other missions requiring analysis of natural orbital thermal environments.

Nuclear Electric Propulsion↗

Geoscience Laser Altimeter System (GLAS) on the ICESat Mission: Initial Science Measurement Performance

The Geoscience Laser Altimeter System is the space lidar on the NASA ICESat mission. Its design combines an altimeter with 5 cm precision with a laser pointing angle determination system and a dual wavelength cloud and aerosol lidar. GLAS measures the range to the Earth s surface with 1064 nm laser pulses. Each laser pulse produces a precision pointing measurement from the stellar reference system (SRS) and an echo pulse waveform, which permits range determination and waveform spreading analysis. The single shot ranging accuracy is < 10 cm for ice surfaces with slopes < 2 degrees. GLAS also measures atmospheric backscatter profiles at both 1064 and 532 nm. The 1064 nm measurements use an analog Si APD detector and measure the height and profile the backscatter signal from thicker clouds. The measurements at 532 nm use photon counting detectors, and will measure the vertical height distributions of optically thin clouds and aerosol layers Before launch, the measurement performance of GLAS was evaluated using a lidar test instrument called the Bench Check Equipment (BCE). The BCE was developed in parallel with GLAS and served as an inverse altimeter, inverse lidar and a stellar source simulator. It was used to simulate the range of expected optical inputs to the GLAS receiver by illuminating its telescope with simulated background light as well as laser echoes with known powers, energy levels, widths and delay times. The BCE also allowed monitoring of the transmitted laser energy, the angle measurements of the SRS, the co-alignment of the transmitted laser beam to the receiver line of sight, and performance of the flight science algorithms. Performance was evaluated during the GLAS development, before and after environmental tests, and after delivery to the spacecraft. The ICESat observatory was launched into a 94 degree inclination, 590 km altitude circular polar orbit on January 12,2003. Beginning in early February, GLAS was powered on tested in stages. Its 1064 nm optical receiver was evaluated in a several tests using both solar background light and an internal test source. Laser 1 was activated on February 20,2003. GLAS operated with Laser 1 for 38 continuous days on orbit using its 1064 nm receiver channel, producing over 130 million individual laser measurements of the Earth s surface and atmosphere. These nadir-pointed measurements fell along the ICESat s ground track, and spanned more than 4 cycles of the initial 8-day ICESat repeat orbit. The initial GLAS measurement set shows strong echo pulses from ranging to the surface topography, oceans, ice sheets and cloud tops, as well as profiles of clouds and aerosols. The GLAS measurements have unprecedented vertical and angular resolution, and show nearly continuous height profiles of ice, land and ocean surfaces or cloud tops, as well profiles of backscatter from thin clouds and aerosol layers. Examples of these GLAS measurements and an initial assessment of its science measurement performance will be presented.

Abshire, James B.↗

Flight motor set 360L001 (STS-26R). (Reconstructed dynamic loads analysis)

A transient analysis was performed to correlate the predicted versus measured behavior of the Redesigned Solid Rocket Booster (RSRB) during Flight 360L001 (STS-26R) liftoff. Approximately 9 accelerometers, 152 strain gages, and 104 girth gages were bonded to the motors during this event. Prior to Flight 360L001, a finite element model of the RSRB was analyzed to predict the accelerations, strains, and displacements measured by this developmental flight instrumentation (DFI) within an order of magnitude. Subsequently, an analysis has been performed which uses actual Flight 360L001 liftoff loading conditions, and makes more precise predictions for the RSRB structural behavior. Essential information describing the analytical model, analytical techniques used, correlation of the predicted versus measured RSRB behavior, and conclusions, are presented. A detailed model of the RSRB was developed and correlated for use in analyzing the motor behavior during liftoff loading conditions. This finite element model, referred to as the RSRB global model, uses super-element techniques to model all components of the RSRB. The objective of the RSRB global model is to accurately predict deflections and gap openings in the field joints to an accuracy of approximately 0.001 inch. The model of the field joint component was correlated to Referee and Joint Environment Simulation (JES) tests. The accuracy of the assembled RSRB global model was validated by correlation to static-fire tests such DM-8, DM-9, QM-7, and QM-8. This validated RSRB global model was used to predict RSRB structural behavior and joint gap opening during Flight 360L001 liftoff. The results of a transient analysis of the RSRB global model with imposed liftoff loading conditions are presented. Rockwell used many gage measurements to reconstruct the load parameters which were imposed on the RSRB during the Flight 360L001 liftoff. Each load parameter, and its application, is described. Also presented are conclusions and recommendations based on the analysis of this load case and the resulting correlation between predicted and measured RSRB structural behavior.

Call, V. B.↗

Automated Knowledge Discovery From Simulators

A computational method, SimLearn, has been devised to facilitate efficient knowledge discovery from simulators. Simulators are complex computer programs used in science and engineering to model diverse phenomena such as fluid flow, gravitational interactions, coupled mechanical systems, and nuclear, chemical, and biological processes. SimLearn uses active-learning techniques to efficiently address the "landscape characterization problem." In particular, SimLearn tries to determine which regions in "input space" lead to a given output from the simulator, where "input space" refers to an abstraction of all the variables going into the simulator, e.g., initial conditions, parameters, and interaction equations. Landscape characterization can be viewed as an attempt to invert the forward mapping of the simulator and recover the inputs that produce a particular output. Given that a single simulation run can take days or weeks to complete even on a large computing cluster, SimLearn attempts to reduce costs by reducing the number of simulations needed to effect discoveries. Unlike conventional data-mining methods that are applied to static predefined datasets, SimLearn involves an iterative process in which a most informative dataset is constructed dynamically by using the simulator as an oracle. On each iteration, the algorithm models the knowledge it has gained through previous simulation trials and then chooses which simulation trials to run next. Running these trials through the simulator produces new data in the form of input-output pairs. The overall process is embodied in an algorithm that combines support vector machines (SVMs) with active learning. SVMs use learning from examples (the examples are the input-output pairs generated by running the simulator) and a principle called maximum margin to derive predictors that generalize well to new inputs. In SimLearn, the SVM plays the role of modeling the knowledge that has been gained through previous simulation trials. Active learning is used to determine which new input points would be most informative if their output were known. The selected input points are run through the simulator to generate new information that can be used to refine the SVM. The process is then repeated. SimLearn carefully balances exploration (semi-randomly searching around the input space) versus exploitation (using the current state of knowledge to conduct a tightly focused search). During each iteration, SimLearn uses not one, but an ensemble of SVMs. Each SVM in the ensemble is characterized by different hyper-parameters that control various aspects of the learned predictor - for example, whether the predictor is constrained to be very smooth (nearby points in input space lead to similar output predictions) or whether the predictor is allowed to be "bumpy." The various SVMs will have different preferences about which input points they would like to run through the simulator next. SimLearn includes a formal mechanism for balancing the ensemble SVM preferences so that a single choice can be made for the next set of trials.

Burl, Michael↗

Statistical Mechanics of Turbulent Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V.↗

MHD Turbulence and Magnetic Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V↗

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↗