Use of a transcendental approximation in transient conduction analysis
Boundary condition solution and transcendental approximation for transient heat conduction problem
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.
Boundary condition solution and transcendental approximation for transient heat conduction problem
An analytical technique was developed to predict the behavior of a rotor system subjected to sudden unbalance. The technique is implemented in the Turbine Engine Transient Rotor Analysis (TETRA) computer program using the component element method. The analysis was particularly aimed toward blade-loss phenomena in gas turbine engines. A dual-rotor, casing, and pylon structure can be modeled by the computer program. Blade tip rubs, Coriolis forces, and mechanical clearances are included. The analytical system was verified by modeling and simulating actual test conditions for a rig test as well as a full-engine, blade-release demonstration.
The development of two advanced reduced-basis methods, the force derivative method and the Lanczos method, and two widely used modal methods, the mode displacement method and the mode acceleration method, for transient structural analysis of unconstrained structures is presented. Two example structural problems are studied: an undamped, unconstrained beam subject to a uniformly distributed load which varies as a sinusoidal function of time and an undamped high-speed civil transport aircraft subject to a normal wing tip load which varies as a sinusoidal function of time. These example problems are used to verify the methods and to compare the relative effectiveness of each of the four reduced-basis methods for performing transient structural analyses on unconstrained structures. The methods are verified with a solution obtained by integrating directly the full system of equations of motion, and they are compared using the number of basis vectors required to obtain a desired level of accuracy and the associated computational times as comparison criteria.
The results of both linear (stability and synchronous response) and transient nonlinear analyses are reported. Dynamic coefficients were developed for the HPFTP interstage seals, and introduced into the rotordynamic model. The influence on HPFTP rotordynamics of a change in interstage seals from the smooth stepped design to a smooth straight configuration was examined. The sensitivity of the stability and synchronous results to changes in bearing stiffnesses and damping was determined. The influence on rotordynamic stability of a change from the stiff symmetric bearing carrier design to an asymmetric bearing carrier configuration was also studied.
NASA’s Dawn spacecraft entered orbit around the dwarf planet Ceres on March 6, 2015 thereby becoming the first mission to orbit two celestial bodies outside of the Earth-moon system. This feat was made possible due to the use of electric propulsion (EP) to provide the 11 km/s of delta-V that was necessary. This type of propulsion is also enabling for many of NASA’s future missions, such as the proposed Asteroid Redirect Robotic Mission (ARRM). Specifically, the Hall Effect Rocket with Magnetic Shielding (HERMeS) is the propulsive design baseline due to its long operational lifetime and high specific impulse. Extensive modeling is required to verify the thruster's predicted thermal performance and survival during the bounding mission phases. Since the thruster would constrain plasma on the order of tens of thousands of hours, components in the thruster must be able to accommodate much higher sustained temperatures than typical spacecraft hardware. Furthermore, the high temperatures necessitate care in analyzing thermal gradients within the thruster. In this paper, we describe a thermal model of the HERMeS thruster and focus upon the thruster's transient performance under different initial firing conditions. The thruster's transient thermal response is critical to verify margin in the thermal stresses inherent in the thruster and to determine safe start-up procedures to protect the thruster's critical components. This model is used to perform transient thermal analysis of the thruster during startup and shutdown, provide temperature maps for structural analysis, and plan for thermal cycle testing conducted at JPL. The thermal cycle testing includes thruster starts at the cold temperature limits and operation at maximum expected external heat flux. This work will support laboratory environmental testing of the HERMeS thruster, which is being conducted at JPL in 2016.
The convergence properties of three modal methods for solving structural transient response (STR) problems are examined by comparing the accuracy of the total number of modes used to represent the responses. The goal is to identify a method which involves as few degrees of freedom as possible. Two mode superposition methods, the mode displacement method (MDM) and the mode acceleration method (MAM), are fitted with equivalent expressions for damped response. An error norm is employed to compare the accuracies of the methods in calculating displacements, moments, and shear forces for a viscously-damped cantilevered beam experiencing several dynamic loading conditions and levels of damping. The MAM converges quicker and uses fewer dof, mainly because it has a pseudo-static response term which accounts for some of the flexibility of higher modes. A third method, which integrates the convolution integral iteratively over time, generates a higher order approximation to the STR with a reduced number of modes. It converges fastest among the three, but requires definition of an appropriate forcing function. Finally, the MAM and MDM models are applied to study the effects of the forcing frequency and the initial conditions.
Analysis and solution are presented for transient thermal stresses in a free heat-generating flat plate and a free, hollow-generating cylinder as a result of sudden environmental changes. The technique used and graphical results obtained are of interest to the heat transfer industry.
The application of adaptive finite element methods to the solution of transient heat conduction problems in two dimensions is investigated. The computational domain is represented by an unstructured assembly of linear triangular elements and the mesh adaptation is achieved by local regeneration of the grid, using an error estimation procedure coupled to an automatic triangular mesh generator. Two alternative solution procedures are considered. In the first procedure, the solution is advanced by explicit timestepping, with domain decomposition being used to improve the computational efficiency of the method. In the second procedure, an algorithm for constructing continuous lines which pass only once through each node of the mesh is employed. The lines are used as the basis of a fully implicit method, in which the equation system is solved by line relaxation using a block tridiagonal equation solver. The numerical performance of the two procedures is compared for the analysis of a problem involving a moving heat source applied to a convectively cooled cylindrical leading edge.
A computer program to calculate transient and steady state temperatures, pressures, and coolant flows in a cooled turbine blade or vane with an impingement insert is described. Input to the program includes a description of the blade geometry, coolant supply conditions, outside thermal boundary conditions and wheel speed. Coolant-side heat transfer coefficients are calculated internally in the program, with the user specifying the mode of heat transfer at each internal flow station. Program output includes the temperature at each node, the coolant pressures and flow rates, and the inside heat transfer coefficients. A sample problem is discussed.
A computer program to calculate transient and steady state temperatures, pressures, and coolant flows in a cooled turbine blade or vane with an impingement insert is described. Input to the program includes a description of the blade geometry, coolant supply conditions, outside thermal boundary conditions and wheel speed. Coolant-side heat transfer coefficients are calculated internally in the program, with the user specifying the mode of heat transfer at each internal flow station. Program output includes the temperature at each node, the coolant pressures and flow rates, and the inside heat transfer coefficients. A sample problem is discussed.
The status of an effort to increase the efficiency of calculating transient temperature fields in complex aerospace vehicle structures is described. The advantages and disadvantages of explicit and implicit algorithms are discussed. A promising set of implicit algorithms, known as the GEAR package is described. Four test problems, used for evaluating and comparing various algorithms, have been selected and finite element models of the configurations are discribed. These problems include a space shuttle frame component, an insulated cylinder, a metallic panel for a thermal protection system and a model of the space shuttle orbiter wing. Calculations were carried out using the SPAR finite element program, the MITAS lumped parameter program and a special purpose finite element program incorporating the GEAR algorithms. Results generally indicate a preference for implicit over explicit algorithms for solution of transient structural heat transfer problems when the governing equations are stiff. Careful attention to modeling detail such as avoiding thin or short high-conducting elements can sometimes reduce the stiffness to the extent that explicit methods become advantageous.
In a fault tolerant control (FTC) system, a parameter varying FTC law is reconfigured based on fault parameters estimated by fault detection and isolation (FDI) modules. FDI modules require some time to detect fault occurrences in aero-vehicle dynamics. This paper illustrates analysis of a FTC system based on estimated fault parameter transient behavior which may include false fault detections during a short time interval. Using Lyapunov function analysis, the upper bound of an induced-L2 norm of the FTC system performance is calculated as a function of a fault detection time and the exponential decay rate of the Lyapunov function.
A reactivity-initiated accident (RIA) is a design-basis accident under which critical heat flux (CHF) is likely to be exceeded. The operational margin for RIAs is currently determined using steady-state CHF lookup tables, which provide conservative estimates relative to transient CHF phenomena. The Transient Reactor Test Loop (TRTL) facility at Oregon State University is capable of performing out-of-pile rapid heating experiments representative of a RIA at conditions representative of a pressurized water reactor (PWR). Here, to further our understanding of and ability to predict transient CHF under PWR conditions, we performed a sensitivity analysis on a RELAP5-3D model of the TRTL facility coupled to the RAVEN code framework to define a proposed experimental test matrix to be performed at the TRTL facility. We then implemented a flow boiling CHF correlation into RELAP5-3D and performed a secondary sensitivity analysis inspecting the impact of the built-in RELAP5-3D CHF and heat transfer multipliers on both the prediction of CHF and key safety parameters, such as peak cladding temperature and heat flux. The results show that the multiplier with the highest influence toward the prediction of CHF occurrence and the safety parameters is the transient CHF multiplier. Operational performance envelopes have been developed for each of the test matrix cases and will be used for validation once the experiments are performed. The TRTL facility is currently performing shakedown testing to verify system performance prior to proceeding with the experimental campaign. Restart testing results include pump curve restart testing, pressure tests, and heater rod thermocouple transients.
This report focuses on analysis and determination of the causes of transient water levels measured between 1964 and 2017 near the Cheshire underground nuclear test, which was detonated on February 14, 1976 in Area 20 (USDOE, 2015) of the Pahute Mesa (PM) testing areas of the Nevada National Security Site (NNSS) located in southern Nevada (Figure 1). Transient water levels near Cheshire indicate influences of a water-production well and yet-unexplained water level rise attributed to testing effects (Jackson and Fenelon, 2018; Jackson et al. 2021). Transient groundwater flow is relevant to the implementation of the hydrologic source term (HST), which prescribes the quantity and release rates of radionuclide source terms to be applied at individual test locations in assessment of regional-scale radionuclide transport (Carle and Tompson, 2021; Carle et al., 2021), as previously investigated for NNSS testing areas at PM (Pawloski et al., 2001), Frenchman Flat (Carle et al, 2007), Yucca Flat (Mcnab, 2009; Tompson, 2009), and Rainier Mesa (Carle, 2011; Tompson et al., 2011).
This report focuses on analysis and determination of the causes of transient water levels measured between 1964 and 2017 near the Cheshire underground nuclear test, which was detonated on February 14, 1976 in Area 20 (USDOE, 2015b) of the Pahute Mesa (PM) testing areas of the Nevada National Security Site (NNSS) located in southern Nevada (Figure 1). Transient water levels near Cheshire indicate influences of a water-production well and yet-unexplained water level rise attributed to testing effects (Jackson and Fenelon, 2018; Jackson et al. 2021). Transient groundwater flow is relevant to the implementation of the hydrologic source term (HST), which prescribes the quantity and release rates of radionuclide source terms to be applied at individual test locations in assessment of regional-scale radionuclide transport (Carle and Tompson, 2021; Carle et al., 2021), as previously investigated for NNSS testing areas at PM (Pawloski et al., 2001), Frenchman Flat (Carle et al, 2007), Yucca Flat (Carle et al., 2009; Mcnab, 2009; Tompson et al., 2009), and Rainier Mesa (Carle, 2011; Tompson et al., 2011).
Studies aimed at an increase in the efficiency of calculating transient temperature fields in complex aerospace vehicle structures are reported. The advantages and disadvantages of explicit and implicit algorithms are discussed and a promising set of implicit algorithms with variable time steps, known as GEARIB, is described. Test problems, used for evaluating and comparing various algorithms, are discussed and finite element models of the configurations are described. These problems include a coarse model of the Space Shuttle wing, an insulated frame tst article, a metallic panel for a thermal protection system, and detailed models of sections of the Space Shuttle wing. Results generally indicate a preference for implicit over explicit algorithms for transient structural heat transfer problems when the governing equations are stiff (typical of many practical problems such as insulated metal structures). The effects on algorithm performance of different models of an insulated cylinder are demonstrated. The stiffness of the problem is highly sensitive to modeling details and careful modeling can reduce the stiffness of the equations to the extent that explicit methods may become the best choice. Preliminary applications of a mixed implicit-explicit algorithm and operator splitting techniques for speeding up the solution of the algebraic equations are also described.
While several data-driven models have been developed for transient stability assessment, how to consider the uncertainties from load and renewable generations and provide interpretation of data-driven assessment results are still open. This paper proposes an explainable Bayesian Neural Network (BNN) for probabilistic transient stability assessment (TSA). By extracting the uncertainties from loads and wind farms, the BNN model can make a reliable prediction and quantify the prediction uncertainties. We also develop the Gradient Shap algorithm to make the global and local explanations for the probabilistic TSA model, a significant advantage over existing black-box data-driven methods. Numerical results on the modified IEEE 39-bus system show that the proposed method outperforms the existing methods in terms of prediction accuracy and uncertainty quantification capabilities. The explainability of the proposed method allows system operators to design preventive controls for enhancing system stability.