Advanced Stirling Radioisotope Generator (ASRG) Thermal Power Model in MATLAB

This paper presents a one-dimensional steady-state mathematical thermal power model of the ASRG. It aims to provide a guideline of understanding how the ASRG works and what can change its performance. The thermal dynamics and energy balance of the generator is explained using the thermal circuit of the ASRG. The Stirling convertor performance map is used to represent the convertor. How the convertor performance map is coupled in the thermal circuit is explained. The ASRG performance characteristics under i) different sink temperatures and ii) over the years of mission (YOM) are predicted using the one-dimensional model. Two Stirling converter control strategies, i) fixing the hot-end of temperature of the convertor by adjusting piston amplitude and ii) fixing the piston amplitude, were tested in the model. Numerical results show that the first control strategy can result in a higher system efficiency than the second control strategy when the ambient gets warmer or the general-purpose heat source (GPHS) fuel load decays over the YOM. The ASRG performance data presented in this paper doesn't pertain to the ASRG flight unit. Some data of the ASRG engineering unit (EU) and flight unit that are available in public domain are used in this paper for the purpose of numerical studies

A Summary of the Space-Time Conservation Element and Solution Element (CESE) Method

The space-time Conservation Element and Solution Element (CESE) method for solving conservation laws is examined for its development motivation and design requirements. The characteristics of the resulting scheme are discussed. The discretization of the Euler equations is presented to show readers how to construct a scheme based on the CESE method. The differences and similarities between the CESE method and other traditional methods are discussed. The strengths and weaknesses of the method are also addressed

The Space-Time Conservation Element and Solution Element Method-A New High-Resolution and Genuinely Multidimensional Paradigm for Solving Conservation Laws

Without resorting to special treatment for each individual test case, the 1D and 2D CE/SE shock-capturing schemes described previously (in Part I) are used to simulate flows involving phenomena such as shock waves, contact discontinuities, expansion waves and their interactions. Five 1D and six 2D problems are considered to examine the capability and robustness of these schemes. Despite their simple logical structures and low computational cost (for the 2D CE/SE shock-capturing scheme, the CPU time is about 2 micro-secs per mesh point per marching step on a Cray C90 machine), the numerical results, when compared with experimental data, exact solutions or numerical solutions by other methods, indicate that these schemes can accurately resolve shock and contact discontinuities consistently

Application of the space-time conservation element and solution element method to two-dimensional advection-diffusion problems

The existing 2-D alpha-mu scheme and alpha-epsilon scheme based on the method of space-time conservation element and solution element, which were constructed for solving the linear 2-D unsteady advection-diffusion equation and unsteady advection equation, respectively, are tested. Also, the alpha-epsilon scheme is modified to become the V-E scheme for solving the nonlinear 2-D inviscid Burgers equation. Numerical solutions of six test problems are presented in comparison with their exact solutions or numerical solutions obtained by traditional finite-difference or finite-element methods. It is demonstrated that the 2-D alpha-mu, alpha-epsilon, and nu-epsilon schemes can be used to obtain numerical results which are more accurate than those based on some of the traditional methods but without using any artificial tuning in the computation. Similar to the previous 1-D test problems, the high accuracy and simplicity features of the space-time conservation element and solution element method have been revealed again in the present 2-D test results

Orion Active Thermal Control System Dynamic Modeling Using Simulink/MATLAB

This paper presents dynamic modeling of the crew exploration vehicle (Orion) active thermal control system (ATCS) using Simulink (Simulink, developed by The MathWorks). The model includes major components in ATCS, such as heat exchangers and radiator panels. The mathematical models of the heat exchanger and radiator are described first. Four different orbits were used to validate the radiator model. The current model results were compared with an independent Thermal Desktop (TD) (Thermal Desktop, PC/CAD-based thermal model builder, developed in Cullimore & Ring (C&R) Technologies) model results and showed good agreement for all orbits. In addition, the Orion ATCS performance was presented for three orbits and the current model results were compared with three sets of solutions- FloCAD (FloCAD, PC/CAD-based thermal/fluid model builder, developed in C&R Technologies) model results, SINDA/FLUINT (SINDA/FLUINT, a generalized thermal/fluid network-style solver ) model results, and independent Simulink model results. For each case, the fluid temperatures at every component on both the crew module and service module sides were plotted and compared. The overall agreement is reasonable for all orbits, with similar behavior and trends for the system. Some discrepancies exist because the control algorithm might vary from model to model. Finally, the ATCS performance for a 45-hr nominal mission timeline was simulated to demonstrate the capability of the model. The results show that the ATCS performs as expected and approximately 2.3 lb water was consumed in the sublimator within the 45 hr timeline before Orion docked at the International Space Station

The Space-Time Conservation Element and Solution Element Method: A New High-Resolution and Genuinely Multidimensional Paradigm for Solving Conservation Laws

A new high resolution and genuinely multidimensional numerical method for solving conservation laws is being, developed. It was designed to avoid the limitations of the traditional methods. and was built from round zero with extensive physics considerations. Nevertheless, its foundation is mathmatically simple enough that one can build from it a coherent, robust. efficient and accurate numerical framework. Two basic beliefs that set the new method apart from the established methods are at the core of its development. The first belief is that, in order to capture physics more efficiently and realistically, the modeling, focus should be placed on the original integral form of the physical conservation laws, rather than the differential form. The latter form follows from the integral form under the additional assumption that the physical solution is smooth, an assumption that is difficult to realize numerically in a region of rapid chance. such as a boundary layer or a shock. The second belief is that, with proper modeling of the integral and differential forms themselves, the resulting, numerical solution should automatically be consistent with the properties derived front the integral and differential forms, e.g., the jump conditions across a shock and the properties of characteristics. Therefore a much simpler and more robust method can be developed by not using the above derived properties explicitly

The method of space-time conservation element and solution element-applications to one-dimensional and two-dimensional time-marching flow problems

A nontraditional numerical method for solving conservation laws is being developed. The new method is designed from a physicist's perspective, i.e., its development is based more on physics than numerics. Even though it uses only the simplest approximation techniques, a 2D time-marching Euler solver developed recently using the new method is capable of generating nearly perfect solutions for a 2D shock reflection problem used by Helen Yee and others. Moreover, a recent application of this solver to computational aeroacoustics (CAA) problems reveals that: (1) accuracy of its results is comparable to that of a 6th order compact difference scheme even though nominally the current solver is only of 2nd-order accuracy; (2) generally, the non-reflecting boundary condition can be implemented in a simple way without involving characteristic variables; and (3) most importantly, the current solver is capable of handling both continuous and discontinuous flows very well and thus provides a unique numerical tool for solving those flow problems where the interactions between sound waves and shocks are important, such as the noise field around a supersonic over- or under-expansion jet

Application of the space-time conservation element and solution element method to shock-tube problem

An Euler solver based on the method of space-time conservation element and solution element is in this paper to simulate shock-tube flows involving shock waves, contact discontinuities, expansion waves and their intersections. Seven test problems are considered to examine the capability of this method. The numerical results, when compared with exact solutions and/or numerical solutions by other methods, indicate that the present method can accurately resolve strong shock and contact discontinuities without using any ad hoc techniques which are used only at the neighborhood of a discontinuity

Thermal Analysis on Cryogenic Liquid Hydrogen Tank on an Unmanned Aerial Vehicle System

Thermal analyses are performed on the liquid hydrogen (LH2) tank designed for an unmanned aerial vehicle (UAV) powered by solar arrays and a regenerative proton-exchange membrane (PEM) fuel cell. A 14-day cruise mission at a 65,000 ft altitude is considered. Thermal analysis provides the thermal loads on the tank system and the boiling-off rates of LH2. Different approaches are being considered to minimize the boiling-off rates of the LH2. It includes an evacuated multilayer insulation (MLI) versus aerogel insulation on the LH2 tank and aluminum versus stainless steel spacer rings between the inner and outer tank. The resulting boil-off rates of LH2 provided by the one-dimensional model and three-dimensional finite element analysis (FEA) on the tank system are presented and compared to validate the results of the three-dimensional FEA. It concludes that heat flux through penetrations by conduction is as significant as that through insulation around the tank. The tank system with MLI insulation and stainless steel spacer rings result in the lowest boiling-off rate of LH2

CFD Simulation on the J-2X Engine Exhaust in the Center-Body Diffuser and Spray Chamber at the B-2 Facility

A computational fluid dynamic (CFD) code is used to simulate the J-2X engine exhaust in the center-body diffuser and spray chamber at the Spacecraft Propulsion Facility (B-2). The CFD code is named as the space-time conservation element and solution element (CESE) Euler solver and is very robust at shock capturing. The CESE results are compared with independent analysis results obtained by using the National Combustion Code (NCC) and show excellent agreement