SEARCH · Engineering Papers
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.
Neuroscience initiatives aim to develop new technologies and tools to measure and manipulate neuronal circuits. In this paper, to deal with the massive amounts of data generated by these tools, the authors envision the co-location of open data repositories in standardized formats together with high-performance computing hardware utilizing open source optimized analysis codes.
As Martian atmospheric entry vehicles increase in size to accommodate larger payloads, transitional ow may need to be taken into account in the design of the heat shield in order to reduce heat shield mass. The mass of the Thermal Protection System (TPS) comprises a significant portion of the vehicle mass, and a reduction of this mass would result in fuel savings. The current techniques used to design entry shields generally assume fully turbulent flow when the vehicle is large enough to expect transitional flow, and while this worst-case scenario provides a greater factor of safety it may also result in overdesigned TPS and unnecessarily high vehicle mass. Greater accuracy in the prediction of transition would also reduce uncertainty in the thermal and aerodynamic loads. Stability analysis, using e(sup N) -based methods including Linear Stability Theory (LST) and the Parabolized Stability Equations (PSE), offers a physics-based method of transition prediction that has been thoroughly studied and applied in perfect gas flows, and to a more limited extent in reacting and nonequilibrium flows. These methods predict the amplification of a known disturbance frequency and allow identification of the most unstable frequency. Transition is predicted to occur at a critical amplification or N Factor, frequently determined through experiment and empirical correlations. The LAngley Stability and TRansition Analysis Code (LASTRAC), with modifications for thermochemically reacting flows and arbitrary gas mixtures, will be presented with LST results on a simulation of a high enthalpy CO2 gas wind tunnel test relevant to Martian atmospheric entry. The results indicate transition caused by modified Tollmien-Schlichting waves on the leeward side, which are predicted to be more stable and cause transition slightly downstream when thermochemical nonequilibrium is included in the stability analysis for the same mean flow solution.
Steam generator tube rupture(SGTR) accident is one characteristic accident of High Temperature Gas-Cooled Reactor (HTGR), which may cause a potential hazard to reactor system. Because the amount of water ingress into reactor core is related closely to the severity of the accident, the research on water-ingress is significant for the evaluation on inherent safety characteristics of HTGR. Both different accident scenario and the system design will have influence on the amount of the ingress-water under SGTR. In this paper, a pebble-bed high temperature gas cooled reactor is modeled by thermal-hydraulic system analysis code and the water-ingress transient under SGTR is simulated to study the influence of the protective action of emergency shutdown on the water ingress. The research indicates that the protective action has no obvious influence on the amount of the ingress-water within the current studied range. The amount of ingress-water is controlled below the design value to satisfy the demand of reactor safety. The analysis results will provide some support for the HTGR research. (authors)
Several innovative methods were used to allow more efficient and accurate thermal analysis using SINDA-85 and TRASYS, including model integration and reduction, planetary surface calculations, and model animation. Integration with other modeling and analysis codes allows an analyst to import a geometry from a solid modeling or computer-aided design (CAD) software package, rather than building the geometry 'by hand.' This is more efficient as well as potentially more accurate. However, the use of solid modeling software often generates large analytical models. The problem of reducing large models was elegantly solved using the response of the transient derivative to a forcing step function. The thermal analysis of a lunar rover implemented two unusual features of the TRASYS/SINDA system. A little-known TRASYS routine SURFP calculates the solar heating of a rover on the lunar surface for several different rover positions and orientations. This is used not only to determine the rover temperatures, but also to automatically determine the power generated by the solar arrays. The animation of transient thermal results is an effective tool, especially in a vivid case such as the 14-day progress of the sun over the lunar rover. An animated color map on the solid model displays the progression of temperatures.
Risk importance measures are indexes that are used to rank systems, structures and components (SSCs) using risk-informed methods. The most used/known measures are: Risk Reduction Worth (RRW), Risk Achievement Worth (RAW), Birnbaum (B) and Fussell-Vesely (FV). Once obtained from classical Probabilistic Risk Analysis (PRA) analyses, these risk measures can be effectively employed to optimize component testing and maintenance. In contrast to classical PRA methods, Dynamic PRA methods couple stochastic methods with safety analysis codes to determine risk associate to complex systems such as nuclear plants. Compared to classical PRA methods, they can evaluate with higher resolution the safety impact of timing and sequencing of events on the accident progression. The objective of this paper is to present a series of algorithms that can be used to determine classical risk importance measures (RRW, RAW, B and FV) along with newly developed ones from a Dynamic PRA analysis.
Dry weights for single-stage launch vehicles that incorporate nontangent, developed contour bulkheads are estimated and compared to a baseline vehicle with 1.414 aspect ratio ellipsoidal bulkheads. Weights, volumes, and heights of optimized bulkhead designs are computed using a preliminary design bulkhead analysis code. The dry weights of vehicles that incorporate the optimized bulkheads are predicted using a vehicle weights and sizing code. Two optimization approaches are employed. A structural-level method, where the vehicle's three major bulkhead regions are optimized separately and then incorporated into a model for computation of the vehicle dry weight, predicts a reduction of4365 lb (2.2 %) from the 200,679-lb baseline vehicle dry weight. In the second, vehicle-level, approach, the vehicle dry weight is the objective function for the optimization. For the vehicle-level analysis, modified bulkhead designs are analyzed and incorporated into the weights model for computation of a dry weight. The optimizer simultaneously manipulates design variables for all three bulkheads to reduce the dry weight. The vehicle-level analysis predicts a dry weight reduction of 5129 lb, a 2.6% reduction from the baseline weight. Based on these results, nontangent, developed contour bulkheads may provide substantial weight savings for single stage vehicles.
An incompressible flow analysis code, INS3D, has been applied to the development of a mechanical heart assist device. The solution method is based on the artificial compressibility approach and uses an implicit-upwind differencing scheme together with a Gauss-Seidel line relaxation method. The equations are solved in steadily rotating reference frames and the centrifugal and the Coriolis force terms are included as source terms. The resulting computational procedure is validated for liquid rocket engine analysis and applied subsequently to analyze a Ventricular Assist Device (VAD). A new design configuration is developed which includes an inducer upstream of the impeller main blades, and substantial improvement is observed in the performance of the VAD.
Several innovative methods have been used to allow more efficient and accurate thermal analysis using SINDA-85 and TRASYS, including model integration and reduction, planetary surface calculations, and model animation. Integration with other modeling and analysis codes allows an analyst to import a geometry from a solid modeling or computer-aided design (CAD) software package, rather than building the geometry "by hand." This is more efficient as well as potentially more accurate. However, the use of solid modeling software often generates large analytical models. The problem of reducing large models has been elegantly solved using the response of the transient derivative to a forcing step function. The thermal analysis of a lunar rover implemented two unusual features of the TRASYS/SINDA system. A little-known TRASYS routine SURFP calculates the solar heating of a rover on the lunar surface for several different rover positions and orientations. This is used not only to determine the rover temperatures, but also to automatically determine the power generated by the solar arrays. The animation of transient thermal results is an effective tool, especially in a vivid case such as the 14-day progress of the sun over the lunar rover. An animated color map on the solid model displays the progression of temperatures.
This report presents a preliminary demonstration of integrated security/safety consequence modeling to support potential physical security exemption justifications under the new Title 10 Code of Federal Regulations Part 53 licensing pathway, which allows exemptions if sabotage-induced radiological consequences do not exceed 25 rem at the site boundary two hours after release without mitigative actions. A hypothetical heat pipe microreactor subjected to direct sabotage by a 150 lbs. TNT-equivalent high-explosive device was analyzed using a simplified radionuclide release estimate based on DOE-HDBK-3010-94, implemented in the MELCOR severe accident analysis code to model release transport and attenuation through reactor/building compartments, with MELCOR outputs coupled to MACCS to compute offsite doses and the distance-to-threshold metric D 25rem . Six scenarios were evaluated to examine sensitivities to release pathway/building configuration and radionuclide form (vapor versus aerosols with different size distributions). Results show a positive correlation between leak path factor and D 25rem and indicate that reactor building compartmentalization can substantially reduce consequences (a single-story configuration reduced D 25rem by roughly an order of magnitude relative to a cavity-only configuration), while vapor versus aerosol assumptions had limited impact for the modeled two-story case on the spatial grid used.
Dry weights for a SSTO vehicle which incorporates nontangent, developed contour bulkheads are estimated and compared to a baseline vehicle with 1.41 4 aspect ratio ellipsoidal bulkheads, Weights, volumes and heights of optimized bulkhead designs are computed using a preliminary design bulkhead analysis code. The dry weight of a vehicle which incorporates the optimized bulkheads is predicted using a vehicle weights and sizing code. Two optimization approaches are employed. A structural-level method, where the vehicle s three major bulkhead regions are optimized separately and then incorporated into a model for computation of the vehicle dry weight, predicts a reduction of 4365 Ib (2.2 percent) from the 200,679 Ib baseline vehicle dry weight. In the second, vehicle-level, approach, the vehicle dry weight is the objective function for the optimization. During the vehicle- level analysis, modified bulkhead designs are first analyzed, then incorporated into the weights model for computation of a dry weight. The optimizer simultaneously manipulates design variables for all three bulkheads to reduce the dry weight. The vehicle-level analysis predicts a dry weight reduction of 5129 Ib, a 2.6 percent reduction from the baseline value. These results suggest that nontangent, developed contour bulkheads may provide substantial weight savings for SSTO vehicles.
Dry weights for single-stage launch vehicles that incorporate nontangent, developed contour bulkheads are estimated and compared to a baseline vehicle with 1.414 aspect ratio ellipsoidal bulkheads. Weights, volumes, and heights of optimized bulkhead designs are computed using a preliminary design bulkhead analysis code. The dry weights of vehicles that incorporate the optimized bulkheads are predicted using a vehicle weights and sizing code. Two optimization approaches are employed. A structural-level method, where the vehicle's three major bulkhead regions are optimized separately and then incorporated into a model for computation of the vehicle dry weight, predicts a reduction of 4365 Ib (2.2%) from the 200,679-lb baseline vehicle dry weight. In the second, vehicle-level, approach, the vehicle dry weight is the objective function for the optimization. For the vehicle-level analysis, modified bulkhead designs are analyzed and incorporated into the weights model for computation of a dry weight. The optimizer simultaneously manipulates design variables for all three bulkheads to reduce the dry weight. The vehicle- level analysis predicts a dry weight reduction of 5129 Ib, a 2.6% reduction from the baseline weight. Based on these results, nontangent, developed contour bulkheads may provide substantial weight savings for single stage vehicles.
The solar investment community is in need of up-to-date and data-derived system degradation rates to use in financial models that calculate the risk and expectation of energy revenue. To date, financial models choose degradation rates from past studies that are limited and may not be representative of solar projects in development today. As the solar industry expands in all sectors (residential, commercial, and utility), there is a growing amount of energy generation data available to analyze. With the help of the open source degradation analysis code (RdTools) it is now possible to derive degradation on an ongoing basis and continuously provide up-to-date durability statistics to the investment community. Ongoing and accurate statistics of fielded photovoltaic systems allows financial stakeholders to constrain energy revenue projections, lower financial risk, and increase the bankability of solar. This presentation will provide a mid-term progress update on a 2 year SETO-funded degradation project that involves applying RdTools analysis on kWh Analytics’ industry database. With results from the RdTools analysis in hand, machine learning methods are employed to quantify how much variance of system performance degradation can be explained by predictor variables such as environmental factors, equipment bill of materials characteristics, and system design information. The goal of this presentation is to promote the usage of RdTools on an ongoing to data owners and to solicit feedback from researchers and stakeholders on how to maximize impact of the SETO-funded kWh Analytics degradation project.
A discharge process whose trigger conditions are a negative exposed metallic surface surrounded by a less negative dielectric, and a large voltage gradient at a dielectric/metal interface is proposed. Analysis of SCATHA data for a discharge substantiates the postulation. Surface discharges cause a small transient charge transfer to space which results in voltage transients. A method of computing these transients, based on the charge lost through the capacitance to space and a fraction of charge stored in the dielectric at the discharge source was developed. It gives an estimate of the discharge transients at the discharge site, which is used as input for coupling code analysis of structure/system response. The transient computations were applied to a three-axis stabilized, geosynchronous satellite for both sunlight and eclipse charging. The energy of the transient pulses are about 1 mJ for sunlight discharge and 8 mJ for eclipse. Changing of selected coatings on the satellite would relieve the stress.
Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor physics application for multiphysics simulations of advanced reactor designs jointly developed by Idaho National Laboratory and Argonne National Laboratory. This paper summarizes the motivation, significance, architecture, design, and features of Griffin. Griffin offers flexible and extensible features to address the challenges associated with advanced reactor designs. These features range from fundamental particle transport to specific reactor physics tasks. The features cover a wide range including on-the-fly and traditional two-step cross-section generation methods, steady-state and transient transport solvers suitable for both heterogeneous and homogeneous models, high-fidelity depletion where thousands of isotopes can be tracked and low-fidelity depletion characterized by burnup, etc. The most fundamental aspect that sets Griffin apart from other reactor analysis codes is that it is developed based on the MOOSE framework. A modular development approach is strongly enforced, with multiphysics being an essential element considered since the beginning of Griffin’s development. Griffin links various MOOSE physics modules and couples to other MOOSE-based applications and non-MOOSE-based applications for multiphyiscs simulations. Griffin includes three modules: ISOXML for preparing and managing multigroup cross sections, radiation transport for solving the neutron transport equation, and reactor analysis for user-oriented reactor physics analysis functionalities. Griffin uses various finite element methods for spatial discretization, multigroup approximation for energy discretization and discrete ordinates method, spherical harmonics expansion method, and diffusion approximation for streaming direction discretization to solve the neutron transport equation. Griffin’s flexibility is evidenced through Griffin’s various applications to fast reactor, high-temperature reactor, pebble bed reactor, molten salt reactor, and microreactor designs. Griffin development follows the software quality assurance procedure for MOOSE-based applications and with software requirements consistent with the ASME NQA-1 standard. Griffin has been adopted into the reactor analysis system for the U.S. NRC and is in use at U.S. companies, universities and national laboratories.
Three helically-grooved seals were tested and the results were compared to the MTI code SPIRALG. A smooth annular seal was tested at six eccentricity ratios from 0 to 0.5. The following are concluded in this viewgraph presentation: (1) Helical-grooved seals provide a substantial reduction in cross-coupled stiffness coefficients. Negative k(sub xy) values are obtained for no-swirl or low-swirl cases. (2) SPIRALG is completely unsuitable for the type of seal tested, namely, turbulent flow, wide grooves and lands, etc. (3) A good analysis code is needed to guide the design of helically-grooved annular seals including groove and smooth sections.
Welding residual stresses in thin plate A516-70 steel and 2219-T87 aluminum butt weldments were measured by the strain-gage hole drilling and X-ray diffraction methods. The residual stress data were used to construct 3D strain fields which were modeled as thermally induced strains. These 3D strain fields were then analyzed with the WARP31) FEM fracture analysis code in order to predict their effect on fatigue and on fracture. For analyses of fatigue crack advance and subsequent verification testing, fatigue crack growth increments were simulated by successive saw-cuts and incremental loading to generate, as a function of crack length, effects on crack growth of the interaction between residual stresses and load induced stresses. The specimen experimental response was characterized and compared to the WARM linear elastic and elastic-plastic fracture mechanics analysis predictions. To perform the fracture analysis, the plate material's crack tearing resistance was determined by tests of thin plate M(T) specimens. Fracture analyses of these specimen were performed using WARP31D to determine the critical Crack Tip Opening Angle [CTOA] of each material. These critical CTOA values were used to predict crack tearing and fracture in the weldments. To verify the fracture predictions, weldment M(T) specimen were tested in monotonic loading to fracture while characterizing the fracture process.
An analytical investigation was conducted to study the influence of various parameters on predicting the aeroelastic loads and stability of a full-scale hingeless rotor in hover and forward flight. The CAMRAD/JA (Comprehensive Analytical Model of Rotorcraft Aerodynamics and Dynamics, Johnson Aeronautics) analysis code is used to obtain the analytical predictions. Data are presented for rotor blade bending and torsional moments as well as inplane damping data obtained for rotor operation in hover at a constant rotor rotational speed of 425 rpm and thrust coefficients between 0.0 and 0.12. Experimental data are presented from a test in the wind tunnel. Validation of the rotor system structural model with experimental rotor blade loads data shows excellent correlation with analytical results. Using this analysis, the influence of different aerodynamic inflow models, the number of generalized blade and body degrees of freedom, and the control-system stiffness at predicted stability levels are shown. Forward flight predictions of the BO-105 rotor system for 1-G thrust conditions at advance ratios of 0.0 to 0.35 are presented. The influence of different aerodynamic inflow models, dynamic inflow models and shaft angle variations on predicted stability levels are shown as a function of advance ratio.
Cyber-physical testing provides a unique platform to enable the design of resilient space structures. This hybrid approach requires the development of a structural model that accounts for various hazards (e.g., micrometeorite and debris impact) and interacts with physical tests and other sub-system models (e.g., thermal) of the space habitat. A two-dimensional finite element analysis code was developed in MATLAB to facilitate the evaluation of potential designs under operating and unexpected loads and prepare the computational framework for eventually performing cyber-physical testing. The code’s efficiency was enhanced by using an object-oriented programming approach that reduced data transfer between functions. In this study, the code is implemented to predict the response of a dome-style structure made of regolith concrete to impact loading and identify the force magnitude that will cause the tensile strength to be exceeded in domes with different thicknesses.