Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Code Analysis”

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 523 records · Page 29

Application of boundary layer control to HSCT low speed configuration

The feasibility of using boundary layer control (BLC) on a high speed civil transport (HSCT) high lift configuration for low speed performance improvement is studied. The possibility of using a part of the suction system previously designed for cruise laminar flow control (LFC) for low speed BLC mode was explored. A suction scheme for BLC was devised for the HSCT high lift configuration. The extent of the suction regions and levels were determined via panel method based inviscid flow analysis coupled with three-dimensional boundary layer analyses. The compatibility of the BLC suction requirements was assessed, and the modifications necessary to operate the system in both modes of operation, were determined. An assessment of the potential aerodynamic performance improvement resulting from an implementation of the BLC concept compared to an optimized simple flap configuration without BLC was made by an Euler code analysis of the simple flap configuration.

Parikh, P. G.↗

Implementation of Stochastic Gradient Descent in an Automated Glow Peak Identification Software for Multiple Thermoluminescent Dosimeter Types

A glow-curve analysis code was previously developed in C++ to analyze thermoluminescent dosimeter glow curves using automated peak detection while applying a first-order kinetics model. A newer version of this code was implemented to improve the automated peak detection and curve fitting models. The Stochastic Gradient Descent Algorithm was introduced to replace the prior approach of taking first and second-order derivatives for peak detection. Additionally, early stopping mechanisms were invoked to improve the previously used Levenberg-Marquardt Algorithm employed for curve fitting. The two software versions were compared through glow curve analysis of different thermoluminescent dosimeter materials and calculation of the corresponding figures of merit. Altogether improvements were shown, namely an increase in the number of peaks detected and a reduction of the mean figure of merit by approximately 46%.

137Cs↗

Effects of Imperfections on the Buckling Response of Compression-Loaded Composite Shells

The results of an experimental and numerical study of the effects of imperfections on the buckling response of unstiffened thin-walled composite cylindrical shells are presented. Results that identify the individual and combined effects of traditional initial geometric shell-wall imperfections and non-traditional shell-wall thickness variations, shell-end geometric imperfections and variations in loads applied to the ends of the shells on the shell buckling response are included. In addition, results illustrating the effects of manufacturing flaws in the form of gaps between adjacent pieces of graphite-epoxy tape in some of the laminate plies are presented in detail. The shells have been analyzed with a nonlinear finite-element analysis code that accurately accounts for these effects on the buckling and nonlinear responses of the shells. The numerical results indicate that traditional and nontraditional initial imperfections can cause a significant reduction in the buckling load of a compression-loaded composite shell. Furthermore, the results indicate that the imperfections couple in a nonlinear manner. The numerical results correlate well with the experimental results. The nonlinear analysis results are also compared to the results from a traditional linear bifurcation buckling analysis. The results suggest that the nonlinear analysis procedure can be used for determining accurate, high-fidelity design knockdown factors for shell buckling and collapse. The results can also be used to determine the effects of manufacturing tolerances on the buckling response of composite shells.

Hilburger, Mark W.↗

Applications of a damage tolerance analysis methodology in aircraft design and production

Objectives of customer mandated aircraft structural integrity initiatives in design are to guide material selection, to incorporate fracture resistant concepts in the design, to utilize damage tolerance based allowables and planned inspection procedures necessary to enhance the safety and reliability of manned flight vehicles. However, validated fracture analysis tools for composite structures are needed to accomplish these objectives in a timely and economical manner. This paper briefly describes the development, validation, and application of a damage tolerance methodology for composite airframe structures. A closed-form analysis code, entitled SUBLAM was developed to predict the critical biaxial strain state necessary to cause sublaminate buckling-induced delamination extension in an impact damaged composite laminate. An embedded elliptical delamination separating a thin sublaminate from a thick parent laminate is modelled. Predicted failure strains were correlated against a variety of experimental data that included results from compression after impact coupon and element tests. An integrated analysis package was developed to predict damage tolerance based margin-of-safety (MS) using NASTRAN generated loads and element information. Damage tolerance aspects of new concepts are quickly and cost-effectively determined without the need for excessive testing.

Woodward, M. R.↗

Thermal stratification potential in rocket engine coolant channels

The potential for rocket engine coolant channel flow stratification was computationally studied. A conjugate, 3-D, conduction/advection analysis code (SINDA/FLUINT) was used. Core fluid temperatures were predicted to vary by over 360 K across the coolant channel, at the throat section, indicating that the conventional assumption of a fully mixed fluid may be extremely inaccurate. Because of the thermal stratification of the fluid, the walls exposed to the rocket engine exhaust gases will be hotter than an assumption of full mixing would imply. In this analysis, wall temperatures were 160 K hotter in the turbulent mixing case than in the full mixing case. The discrepancy between the full mixing and turbulent mixing analyses increased with increasing heat transfer. Both analysis methods predicted identical channel resistances at the coolant inlet, but in the stratified analysis the thermal resistance was negligible. The implications are significant. Neglect of thermal stratification could lead to underpredictions in nozzle wall temperatures. Even worse, testing at subscale conditions may be inadequate for modeling conditions that would exist in a full scale engine.

Kacynski, Kenneth J.↗

An Analysis of the Macroscopic Tensile Behavior of a Nonlinear Nylon Reinforced Elastomeric Composite System Using MAC/GMC

A special class of composite laminates composed of soft rubbery matrices and stiff reinforcements made of steel wires or synthetic fibers is examined, where each constituent behaves in a nonlinear fashion even in the small strain domain. Composite laminates made of piles stacked at alternating small orientation angles with respect to the applied axial strain are primarily dominated by the nonlinear behavior of the reinforcing fibers. However; composites with large ply orientations or those perpendicular to the loading axis, will approximate the behavior of the matrix phase and respond in even a more complex fashion for arbitrarily stacked piles. The geometric nonlinearity due to small cord rotations during loading was deemed here to have a second order effect and consequently dropped from any consideration. The user subroutine USRMAT within the Micromechanics Analysis Code with the Generalized Method of Cells (MAC/GMC), was utilized to introduce the constituent material nonlinear behavior. Stress-strain behavior at the macro level was experimentally generated for single and multi ply composites comprised of continuous Nylon-66 reinforcements embedded in a carbon black loaded rubbery matrix. Comparisons between the predicted macro composite behavior and experimental results are excellent when material nonlinearity is included in the analysis. In this paper, a brief review of GMC is provided, along with a description of the nonlinear behavior of the constituents and associated constituent constitutive relations, and the improved macro (or composite) behavior predictions are documented and illustrated.

Assaad, Mahmoud↗

Analytical Round Robin for Elastic-Plastic Analysis of Surface Cracked Plates: Phase I Results

An analytical round robin for the elastic-plastic analysis of surface cracks in flat plates was conducted with 15 participants. Experimental results from a surface crack tension test in 2219-T8 aluminum plate provided the basis for the inter-laboratory study (ILS). The study proceeded in a blind fashion given that the analysis methodology was not specified to the participants, and key experimental results were withheld. This approach allowed the ILS to serve as a current measure of the state of the art for elastic-plastic fracture mechanics analysis. The analytical results and the associated methodologies were collected for comparison, and sources of variability were studied and isolated. The results of the study revealed that the J-integral analysis methodology using the domain integral method is robust, providing reliable J-integral values without being overly sensitive to modeling details. General modeling choices such as analysis code, model size (mesh density), crack tip meshing, or boundary conditions, were not found to be sources of significant variability. For analyses controlled only by far-field boundary conditions, the greatest source of variability in the J-integral assessment is introduced through the constitutive model. This variability can be substantially reduced by using crack mouth opening displacements to anchor the assessment. Conclusions provide recommendations for analysis standardization.

Wells, D. N.↗

Boundary Layer Stability and Laminar-Turbulent Transition Analysis with Thermochemical Nonequilibrium Applied to Martian Atmospheric Entry

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.

Transition↗

Some useful innovations with TRASYS and SINDA-85

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.

Amundsen, Ruth M.↗

Nontangent, Developed Contour Bulkheads for a Single-Stage Launch Vehicle

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.

Wu, K. Chauncey↗

Incompressible Navier-Stokes Computations for the Development of a Ventricular Assist Device

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.

Kiris, Cetin↗

Some Useful Innovations with Trasys and Sinda-85

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.

Amundsen, Ruth M.↗

Integrated Security-Safety Consequence Study of a Heat Pipe Reactor

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.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Nontangent, Developed Contour Bulkheads for a Wing-Body Single Stage Launch Vehicle

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.

Wu, K. Chauncey↗

TPSAS-NF1676L-13638-DND

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.

K Chauncey Wu↗

Deciphering Degradation: Machine Learning on Real-World Performance Data

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.

14 SOLAR ENERGY↗

Environmentally induced discharges on satellites

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.

Stevens, N. J.↗

Griffin: A MOOSE-based reactor physics application for multiphysics simulation of advanced nuclear reactors

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.

97 MATHEMATICS AND COMPUTING↗