Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational Fluid Dynamic (CFD)”

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 847 records · Page 47

Numerical Prediction of Heat Transfer Coefficients during Xenon Tank Fill in Microgravity

NASA’s Gateway will serve as an orbital outpost to enable sustained human lunar exploration and facilitate access to destinations beyond the Moon. To maintain sustainable deep space missions, on-orbit propellant transfer and refueling must be realized. One of Gateway’s critical components, the Power and Propulsion Element (PPE), will use solar electric propulsion to perform attitude control maneuvers and orbit transfers. Xenon will be used as the propellant for the PPE and stored in a composite overwrapped pressure vessel (COPV). During the propellant refueling process, xenon will compress causing the fluid and tank walls to warm. There is concern that the bond between the tank liner and composite may degrade at elevated temperatures, along with the possibility of the tank exceeding design pressure. While on-orbit testing is being proposed, numerical models are being used to predict the thermal response and reduce risk during a refueling operation in microgravity where natural convection is diminished. Accurate numerical models may be used to inform decisions for follow-on testing, or design and operation of refueling architectures. Challenges of modeling xenon tank fill in microgravity include non-ideal gas behavior, operating near the fluid’s critical point, and lack of experimental flight data. This study presents a computational fluid dynamics (CFD) model with conjugate heat transfer that is used to predict averaged heat transfer coefficients during tank filling. At the anticipated injection flow rates, the CFD results showed that forced convection dominates in microgravity. Initial results from a multi-node Thermal Desktop model that assumed a purely conducting fluid (Nu=1) over-predicted the CFD fluid temperature and pressure, while a similar model allowing for natural convection under-predicted. The heat transfer coefficients predicted by CFD were implemented into the multi-node model, and the COPV thermal response from the two computational approaches were compared with good agreement.

Heat Transfer↗

Investigation of Damping Physics and CFD Tool Validation for Simulation of Baffled Tanks at Variable Slosh Amplitude

Determination of slosh damping is a very challenging task as there is no analytical solution. The damping physics involves the vorticity dissipation which requires the full solution of the nonlinear Navier-Stokes equations. As a result, previous investigations were mainly carried out by extensive experiments. A systematical study is needed to understand the damping physics of baffled tanks, to identify the difference between the empirical Miles equation and experimental measurements, and to develop new semi-empirical relations to better represent the real damping physics. The approach of this study is to use Computational Fluid Dynamics (CFD) technology to shed light on the damping mechanisms of a baffled tank. First, a 1-D Navier-Stokes equation representing different length scales and time scales in the baffle damping physics is developed and analyzed. Loci-STREAM-VOF, a well validated CFD solver developed at NASA MSFC, is applied to study the vorticity field around a baffle and around the fluid-gas interface to highlight the dissipation mechanisms at different slosh amplitudes. Previous measurement data is then used to validate the CFD damping results. The study found several critical parameters controlling fluid damping from a baffle: local slosh amplitude to baffle thickness (A/t), surface liquid depth to tank radius (d/R), local slosh amplitude to baffle width (A/W); and non-dimensional slosh frequency. The simulation highlights three significant damping regimes where different mechanisms dominate. The study proves that the previously found discrepancies between Miles equation and experimental measurement are not due to the measurement scatter, but rather due to different damping mechanisms at various slosh amplitudes. The limitations on the use of Miles equation are discussed based on the flow regime.

Yang, H. Q.↗

Simplified Applied Loads of a Transonic Truss-Braced Wing Derived from Computational Fluid Dynamics

The transonic truss-braced wing (TTBW) concept is under consideration by the NASA Advanced Air Transport Technology (AATT) project due to its potential to reduce fuel consumption of transport aircraft. To properly evaluate the TTBW concept, detailed structural modeling and subsequent test planning is underway, which requires a simplified set of discrete applied loads that simulate the distributed aerodynamic pressures. Given the unique architecture of the TTBW and its intended operation in the transonic regime, computational fluid dynamics (CFD) is the ideal approach to employ to ensure the aerodynamics of this configuration are sufficiently captured. To convert aerodynamic results into input for structural models, application of a data transfer tool, Discrete Data Transfer Between Dissimilar Meshes (DDTBDM), is presented to resolve aerodynamic pressures from CFD into single load vectors at hardpoints along both the wing and the strut. Evaluation of mapping the aerodynamic forces to hardpoints at rib and spar intersections in a structural model is performed, and the results are also compared with a linear static aeroelastic analysis. The structural responses indicate that distributing forces along only the wing and strut forward spars is insufficient for simulating the distributed aerodynamic loads. Adding moments along the forward spars improved the global response of the wing, but the point moments created unrealistic local rotations within the structure. When distributing loads along both the forward and aft spars of the wing and strut, the global response was less sensitive to the inclusion of point moments, which indicated the arrangement of point forces alone was sufficient input for structural design work and testing and provided a more continuous loading profile.

Finite Element↗

Simplified Applied Loads of a Transonic Truss-Braced Wing Derived from Computational Fluid Dynamic

The transonic truss-braced wing (TTBW) concept is under consideration by the NASA Advanced Air Transport Technology (AATT) project due to its potential to reduce fuel consumption of transport aircraft. To properly evaluate the TTBW concept, detailed structural modeling and subsequent test planning is underway, which requires a simplified set of discrete applied loads that simulate the distributed aerodynamic pressures. Given the unique architecture of the TTBW and its intended operation in the transonic regime, computational fluid dynamics (CFD) is the ideal approach to employ to ensure the aerodynamics of this configuration are sufficiently captured. To convert aerodynamic results into input for structural models, application of a data transfer tool, Discrete Data Transfer Between Dissimilar Meshes (DDTBDM), is presented to resolve aerodynamic pressures from CFD into single load vectors at hardpoints along both the wing and the strut. Evaluation of mapping the aerodynamic forces to hardpoints at rib and spar intersections in a structural model is performed, and the results are also compared with a linear static aeroelastic analysis. The structural responses indicate that distributing forces along only the wing and strut forward spars is insufficient for simulating the distributed aerodynamic loads. Adding moments along the forward spars improved the global response of the wing, but the point moments created unrealistic local rotations within the structure. When distributing loads along both the forward and aft spars of the wing and strut, the global response was less sensitive to the inclusion of point moments, which indicated the arrangement of point forces alone was sufficient input for structural design work and testing and provided a more continuous loading profile.

Finite Element↗

Dynamics of Numerics and CFD

This lecture attempts to illustrate the basic ideas of how the recent advances in nonlinear dynamical systems theory (dynamics) can provide new insights into the understanding of numerical algorithms used in solving nonlinear differential equations (DEs). Examples will be given of the use of dynamics to explain unusual phenomena that occur in numerics. The inadequacy of the use of linearized analysis for the understanding of long time behavior of nonlinear problems will be illustrated, and the role of dynamics in studying the nonlinear stability, accuracy, convergence property and efficiency of using time- dependent approaches to obtaining steady-state numerical solutions in computational fluid dynamics (CFD) will briefly be explained.

Yee, H. C.↗

NASP and SDI Spearhead CFD Developments

The National Aerospace Plane (NASP) program's purpose, as stated by the National Space Council, is to "develop and demonstrate hypersonic technologies with the ultimate goal of single stage to orbit." The council has also directed that "performance of the experimental flight vehicle will be constrained to the minimum necessary to meet the highest priority research, as opposed to operational objectives .... The program will be conducted in such a way as to minimize technical and cost uncertainty associated with the experimental vehicle." The purpose of the Strategic Defense Initiative (SDI), as defined by President Bush, is "...protection from limited ballistic missile strikes, whatever their source." Computational fluid dynamics (CFD) plays a vital role in both endeavors.

Mehta, Unmeel B.↗

Development and application of unified algorithms for problems in computational science

A framework is presented for developing computationally unified numerical algorithms for solving nonlinear equations that arise in modeling various problems in mathematical physics. The concept of computational unification is an attempt to encompass efficient solution procedures for computing various nonlinear phenomena that may occur in a given problem. For example, in Computational Fluid Dynamics (CFD), a unified algorithm will be one that allows for solutions to subsonic (elliptic), transonic (mixed elliptic-hyperbolic), and supersonic (hyperbolic) flows for both steady and unsteady problems. The objectives are: development of superior unified algorithms emphasizing accuracy and efficiency aspects; development of codes based on selected algorithms leading to validation; application of mature codes to realistic problems; and extension/application of CFD-based algorithms to problems in other areas of mathematical physics. The ultimate objective is to achieve integration of multidisciplinary technologies to enhance synergism in the design process through computational simulation. Specific unified algorithms for a hierarchy of gas dynamics equations and their applications to two other areas: electromagnetic scattering, and laser-materials interaction accounting for melting.

Shankar, Vijaya↗

Design of Experiments for Both Experimental and Analytical Study of Exhaust Plume Effects on Sonic Boom

Computational fluid dynamics (CFD) analysis has been performed to study the plume effects on sonic boom signature for isolated nozzle configurations. The objectives of these analyses were to provide comparison to past work using modern CFD analysis tools, to investigate the differences of high aspect ratio nozzles to circular (axisymmetric) nozzles, and to report the effects of under expanded nozzle operation on boom signature. CFD analysis was used to address the plume effects on sonic boom signature from a baseline exhaust nozzle. Nearfield pressure signatures were collected for nozzle pressure ratios (NPRs) between 6 and 10. A computer code was used to extrapolate these signatures to a ground-observed sonic boom N-wave. Trends show that there is a reduction in sonic boom N-wave signature as NPR is increased from 6 to 10. As low boom designs are developed and improved, there will be a need for understanding the interaction between the aircraft boat tail shocks and the exhaust nozzle plume. These CFD analyses will provide a baseline study for future analysis efforts. For further study, a design of experiments has been conducted to develop a hybrid method where both CFD and small scale wind tunnel testing will validate the observed trends. The CFD and testing will be used to screen a number of factors which are important to low boom propulsion integration, including boat tail angle, nozzle geometry, and the effect of spacing and stagger on nozzle pairs. To design the wind tunnel experiment, CFD was instrumental in developing a model which would provide adequate space to observe the nozzle and boat tail shock structure without interference from the wind tunnel walls.

Castner, Raymond S.↗

Assessment of Hybrid RANS/LES Turbulence Models for Aeroacoustics Applications

Predicting the noise from aircraft with exposed landing gear remains a challenging problem for the aeroacoustics community. Although computational fluid dynamics (CFD) has shown promise as a technique that could produce high-fidelity flow solutions, generating grids that can resolve the pertinent physics around complex configurations can be very challenging. Structured grids are often impractical for such configurations. Unstructured grids offer a path forward for simulating complex configurations. However, few unstructured grid codes have been thoroughly tested for unsteady flow problems in the manner needed for aeroacoustic prediction. A widely used unstructured grid code, FUN3D, is examined for resolving the near field in unsteady flow problems. Although the ultimate goal is to compute the flow around complex geometries such as the landing gear, simpler problems that include some of the relevant physics, and are easily amenable to the structured grid approaches are used for testing the unstructured grid approach. The test cases chosen for this study correspond to the experimental work on single and tandem cylinders conducted in the Basic Aerodynamic Research Tunnel (BART) and the Quiet Flow Facility (QFF) at NASA Langley Research Center. These configurations offer an excellent opportunity to assess the performance of hybrid RANS/LES turbulence models that transition from RANS in unresolved regions near solid bodies to LES in the outer flow field. Several of these models have been implemented and tested in both structured and unstructured grid codes to evaluate their dependence on the solver and mesh type. Comparison of FUN3D solutions with experimental data and numerical solutions from a structured grid flow solver are found to be encouraging.

Vatsa, Veer N.↗

Analysis of Plume Effects on Sonic Boom Signature for Isolated Nozzle Configurations

Computational fluid dynamics (CFD) analysis has been performed to study the plume effects on sonic boom signature for isolated nozzle configurations. The objectives of these analyses were to provide comparison to past work using modern CFD analysis tools, to investigate the differences of high aspect ratio nozzles to circular (axisymmetric) nozzles, and to report the effects of underexpanded nozzle operation on boom signature. CFD analysis was used to address the plume effects on sonic boom signature from a baseline exhaust nozzle. Near-field pressure signatures were collected for nozzle pressure ratios (NPRs) between 6 and 10. A computer code was used to extrapolate these signatures to a ground-observed sonic boom N-wave. Trends show that there is a reduction in sonic boom N-wave signature as NPR is increased from 6 to 10. The performance curve for this supersonic nozzle is flat, so there is not a significant loss in thrust coefficient as the NPR is increased. As a result, this benefit could be realized without significant loss of performance. Analyses were also collected for a high aspect ratio nozzle based on the baseline design for comparison. Pressure signatures were collected for nozzle pressure ratios from 8 to 12. Signatures were nearly twice as strong for the two-dimensional case, and trends also show a reduction in sonic boom signature as NPR is increased from 8 to 12. As low boom designs are developed and improved, there will be a need for understanding the interaction between the aircraft boat tail shocks and the exhaust nozzle plume. These CFD analyses will provide a baseline study for future analysis efforts.

Castner, Raymond S.↗

Analysis of Plume Effects on Sonic Boom Signature for Isolated Nozzle Configurations

Computational fluid dynamics (CFD) analysis has been performed to study the plume effects on sonic boom signature for isolated nozzle configurations. The objectives of these analyses were to provide comparison to past work using modern CFD analysis tools, to investigate the differences of high aspect ratio nozzles to circular (axisymmetric) nozzles, and to report the effects of underexpanded nozzle operation on boom signature. CFD analysis was used to address the plume effects on sonic boom signature from a baseline exhaust nozzle. Near-field pressure signatures were collected for nozzle pressure ratios (NPRs) between 6 and 10. A computer code was used to extrapolate these signatures to a ground-observed sonic boom N-wave. Trends show that there is a reduction in sonic boom N-wave signature as NPR is increased from 6 to 10. The performance curve for this supersonic nozzle is flat, so there is not a significant loss in thrust coefficient as the NPR is increased. As a result, this benefit could be realized without significant loss of performance. Analyses were also collected for a high aspect ratio nozzle based on the baseline design for comparison. Pressure signatures were collected for nozzle pressure ratios from 8 to 12. Signatures were nearly twice as strong for the two-dimensional case, and trends also show a reduction in sonic boom signature as NPR is increased from 8 to 12. As low boom designs are developed and improved, there will be a need for understanding the interaction between the aircraft boat tail shocks and the exhaust nozzle plume. These CFD analyses will provide a baseline study for future analysis efforts.

Castner, Raymond S.↗

CFD Simulations of the Supersonic Inflatable Aerodynamic Decelerator (SIAD) Ballistic Range Tests

A series of ballistic range tests were performed on a scaled model of the Supersonic Flight Demonstration Test (SFDT) intended to test the Supersonic Inflatable Aerodynamic Decelerator (SIAD) geometry. The purpose of these experiments were to provide aerodynamic coefficients of the vehicle to aid in mission and vehicle design. The experimental data spans the moderate Mach number range, $3.8-2.0$, with a total angle of attack ($alpha_T$) range, $10o-20o$. These conditions are intended to span the Mach-$alpha$ space for the majority of the SFDT experiment. In an effort to validate the predictive capabilities of Computational Fluid Dynamics (CFD) for free-flight aerodynamic behavior, numerical simulations of the ballistic range experiment are performed using the unstructured finite volume Navier-Stokes solver, US3D. Comparisons to raw vehicle attitude, and post-processed aerodynamic coefficients are made between simulated results and experimental data. The resulting comparisons for both raw model attitude and derived aerodynamic coefficients show good agreement with experimental results. Additionally, near body pressure field values for each trajectory simulated are investigated. Extracted surface and wake pressure data gives further insights into dynamic flow coupling leading to a potential mechanism for dynamic instability.

CFD modeling↗

Higher Order Chemistry Models in the CFD Simulation of Laser-Ablated Carbon Plumes

Production of single-walled carbon nanotubes (SWNT) has taken place for a number of years and by a variety of methods such-as laser ablation, chemical vapor deposition, and arc-jet ablation. Yet, little is actually understood about the exact chemical kinetics and processes that occur in SWNT formation. In recent time, NASA Johnson Space Center has devoted a considerable effort to the experimental evaluation of the laser ablation production process for SWNT originally developed at Rice University. To fully understand the nature of the laser ablation process it is necessary to understand the development of the carbon plume dynamics within the laser ablation oven. The present work is a continuation of previous studies into the efforts to model plume dynamics using computational fluid dynamics (CFD). The ultimate goal of the work is to improve understanding of the laser ablation process, and through that improved understanding, refine the laser ablation production of SWNT. Fig. 1 shows a basic schematic of the laser-ablation oven at NASA-JSC. Construction of the facility is simple in concept. Two concentric quartz tubes of 1.5 mm thickness form the inner and outer tubes with inside diameters of 2.2 and 5.08 cm respectively. At one end of the inner tube are located two 60 Hz pulsed lasers operating at 1064 nm and 532 nm wavelength with beam diameters of 5 mm aligned coaxially with the longitudinal axis of the inner quartz tube. For standard nanotube production runs, a 10 ns 532 nm pulse is followed 50 ns later by a 10 ns 1064 nm pulse. Each pulse is of 300 mJ energy. A target of carbon graphite with approximately 1% nickel and cobalt catalysts is located at the other end of the inner quartz tube. In the ordinary processing of SWNT, a base flow of 100 sccm of argon is maintained from the laser location and exits past the carbon target at a pressure of 66.7 kPa. These conditions yield a baseline mass flow through the chamber of 2.723x10(exp -6)kg/s of argon. The whole oven facility is heated to a temperature of 1473 K prior to nanotube production runs. Upon laser irradiation, part of the carbon target ablates immediately and forms a carbon vapor plume that penetrates into the argon base flow towards the laser initially at supersonic velocities. In the previous studies either a single carbon species, C3, was used to model the plume development, or a simplified 11 species carbon model reduced down from the Krestinin and Moravsky model for full fullerene chemistry was used. While both of these options yielded significant results, it was felt that the actual chemistry occ&g in the carbon plume might have a greater affect on the plume than assumed. Indeed, in the earlier attempts at modeling the carbon plume, several thermophysical characteristics could never be matched to experimental observations of plume development - mainly the propagation distance of the plume itself. In the present study, two additional chemistry models will be used to duplicate the previous studies simulations of the carbon plume. The first chemistry model used in this study is again a reduced form of the Krestinin and Moravsky rates. However the highest order carbon species allowed has been increased from C6 to C30 - therefore allowing the simulation of up to a half of the standard C60 fullerene. The second chemistry model investigated is a reduced form of a full carbon nanotube model developed at NASA-JSC. The C30 studies have already been accomplished at the present time, and the reduced SWNT model studies are currently underway. To pursue the current study, one sacrifice had to be made in that the simulation grid spacing had to be increased from 0.5 mm spacing to 1 mm spacing for the sake of computational efficiency since computational effort is proportional to the square of the number of grid points multiplied by the number of species considered. propagation that is far more in line with the experimental results observed by Puretzky et al as shown in Fig. 2. e C6 studies had yielded a far greater propagation in previous studies. In addition, chemical species development with the C30 model indicates that many higher order carbon species are produced outside of the plume proper (indicated by plotting contours of the background argon concentrations in Fig. 3) - this result was not observed in previous studies. In fact, some species primarily occurred outside of the plume itself - as shown for C27 in Fig. 4 when it is compared to Fig. 3. It could be asserted that this has occurred because all of the C27 in the plume had already been consumed in the formation of C30, but this does not seem to be indicated over time. Several other factors that arose in the previous studies have also been made more clear by the use of the higher order chemical models - one being that the use of c6 as an indicator species was mistaken. C6 is the only carbon species in the previous studies that was not injected into the flowfield as a boundary condition; it was therefore hoped that this species would provide insight into the formation of higher order carbon species for comparison to full SWNT production. But, when the plot of total mass in the plume is examined on a species by species basis in Fig. 5, it is seen that Cg was a fairly insignificant contributor to the total carbon mass in the plume and would not provide information on higher order carbon formation. the thermophysical characteristics of the carbon plume as well as simulate the carbon plume using the reduced SWNT model to provide an even better simulation of full chemistry effects upon plume propagation.

Scott, C. D.↗

CFD Analysis of Modular Thrusters Performance

The effective performance of modular thrusters in an aerospike configuration is difficult to determine. Standard analytical tools are applicable to conventional nozzle shapes, but are limited when applied to an aerospike nozzle (An aerospike nozzle is an altitude compensating external nozzle). Three baseline nozzle shapes are derived using standard analytical procedures. The baseline nozzle sizes are restricted to fill a volume envelope. The three shapes are an axi-symmetric round nozzle, a two dimensional planar square exit nozzle, and a super elliptic round to nearly square nozzle. The integrated (thruster/aerospike) performance of the three nozzles is determined through the use of three dimensional viscous computational fluid dynamic (CFD) calculation where complex features of the flow field can be accurately captured. The resulting installed performance is then used to evaluate the efficiency of these nozzle shapes for aerospike applications. The determination of effective performance of a thruster nozzle integrated into an aerospike nozzle require the solution of the three dimensional turbulent Navier-Stokes equations. The model used in this study consisted of two zones; one of the upstream thruster cowl surface so freestream conditions can be accurately predicted, and two, the aerospike surface beginning with the thruster outflow and extending to the end of the aerospike surface. The numerical grid consisted of over 120,000 nodes and used symmetry on the thruster centerline and edge. A two species non-reacting chemistry model was used to capture the variation of fluid properties between the hot plume base and freestream air. From the results of the three baseline nozzle aerospike calculations, the effictive performance of the nozzle was determined. The flow fields of these calculations do show some variation between the cases. Recirculation zones on the cowl surface is predicted for the two dimensional planar nozzle and a smaller one for the super elliptic nozzle. The recirculation is caused by the strong pressure gradient between the plume and freestream flows. The axi-symmetric nozzle results indicate recirculation zones on the thruster face. These recirculation zones smooth the pressure gradient between the plume and freestream flow limiting the formation of recirculation on the cowl surface. Thruster to thruster interaction is evident for the axi-symmetric and supper elliptic calculation while the two dimensional planar nozzle did not have any lateral expansion in the nozzle, so thruster to thruster interaction is limited. The integrated performance results, at the altitude choosen, show very little variation between the three thruster shapes. This result allows for nozzle shape determination based on additional considerations (thermal, structural, weight) besides performance.

Ungewitter, Ronald J.↗

X-59 CFD and Sonic Boom Assessments

X-59 is a NASA experimental demonstrator aircraft designed to mitigate the adverse impact of sonic boom through advanced aircraft shaping methodologies. The Lockheed Martin (LM) Corporation was awarded a contract to manufacture the X-59 to demonstrate unique boom shaping technologies and address system integration challenges to be able to limit maximum sonic boom perceived level of loudness to be less than 75 dB across the entire carpet while cruising at a Mach number of 1.4. Initial flights and envelope expansion of the X-59 is expected to begin in 2022, with community tests expected at multiple locations starting from early 2024. This paper discusses the Computational Fluid Dynamics (CFD) analysis and sonic boom assessment of the X-59.

John M. Morgenstern↗

Ignition Transient Calculations in the Space Shuttle Solid Rocket Motor

The work presented is part of an effort to develop a multidimensional ignition transient model for large solid propellant rocket motors. On the Space Shuttle, the ignition transient in the slot is induced when the igniter, itself a small rocket motor, is fired into the head-end portion of the main rocket motor. The computational results presented in this paper consider two different igniter configurations. The first configuration is a simulated Space Shuttle RSRM igniter which has one central nozzle that is parallel to the centerline of the motor. The second igniter configuration has a nozzle which is canted at an angle of 45 deg from the centerline of the motor. This paper presents a computational fluid dynamic (CFD) analyses of certain flow field characteristics inside the solid propellant star grain slot of the Space Shuttle during the ignition transient period of operation for each igniter configuration. The majority of studies made to date regarding ignition transient performance in solid rocket motors have concluded that the key parameter to be determined is the heat transfer rate to the propellant surface and hence the heat transfer coefficient between the gas and the propellant. In this paper the heat transfer coefficients, pressure and velocity distributions are calculated in the star slot. In order to validate the computational method and to attempt to establish a correlation between the flow field characteristics and the heat transfer rates a series of cold flow experimental investigations were conducted. The results of these experiments show excellent qualitative and quantitative agreement with the pressure and velocity distributions obtained from the CFD analysis. The CFD analysis utilized a classical pipe flow type correlation for the heat transfer rates. The experimental results provide an excellent qualitative comparison with regard to spatial distribution of the heat transfer rates as a function of nozzle configuration and igniter pressure. The results indicate that from a quantitative point of view that the pipe flow correlation gives reasonably good results. Furthermore, there appears to be a direct correlation between igniter pressure and an average Reynolds number in the star grain slot. This may lead to a simple method for modifying the convection heat transfer correlation. Calculated results of pressure-vs-time for the first 200 msec of motor firing of the Space Shuttle RSRM support the trends shown for the heat transfer rate comparisons between the cold flow CFD and experimental data.

Jenkins, Rhonald M.↗

Axisymmetric Plume Simulations with NASA's DSMC Analysis Code

A comparison of axisymmetric Direct Simulation Monte Carlo (DSMC) Analysis Code (DAC) results to analytic and Computational Fluid Dynamics (CFD) solutions in the near continuum regime and to 3D DAC solutions in the rarefied regime for expansion plumes into a vacuum is performed to investigate the validity of the newest DAC axisymmetric implementation. This new implementation, based on the standard DSMC axisymmetric approach where the representative molecules are allowed to move in all three dimensions but are rotated back to the plane of symmetry by the end of the move step, has been fully integrated into the 3D-based DAC code and therefore retains all of DAC s features, such as being able to compute flow over complex geometries and to model chemistry. Axisymmetric DAC results for a spherically symmetric isentropic expansion are in very good agreement with a source flow analytic solution in the continuum regime and show departure from equilibrium downstream of the estimated breakdown location. Axisymmetric density contours also compare favorably against CFD results for the R1E thruster while temperature contours depart from equilibrium very rapidly away from the estimated breakdown surface. Finally, axisymmetric and 3D DAC results are in very good agreement over the entire plume region and, as expected, this new axisymmetric implementation shows a significant reduction in computer resources required to achieve accurate simulations for this problem over the 3D simulations.

Stewart, B. D.↗

Modeling the Effects of Liquid-to-Gas Density Ratio on Slosh Dynamics

Mechanical models are commonly used in Guidance, Navigation and Controls (GN&C) system-level models to represent propellant slosh forces on launch vehicles and spacecraft. The slosh model parameters are typically predicted using semi-empirical analytical methods that are based on experimental data for high liquid-to-gas density ratios (such as water and air at standard sea level conditions). Model parameter calculations typically neglect any contribution to slosh dynamics from the gas phase. However, some cryogenic propellant systems operate in conditions where the density ratio can be orders of magnitude smaller. Analytical and computational modeling was used in this effort to investigate the effects of liquid-to-gas density ratio on slosh dynamics. No experimental data nor previous studies are available at this time that quantify the effect of liquid-to-gas density ratio on slosh dynamics. However, in a recent CFD study, results showed that density ratio can have a significant effect on slosh model parameters. In this study, a dual pendulum slosh model was created that distinctly represents both liquid and gas phase dynamics. Solution of the model, which invokes the Euler-Lagrange equations of motion, demonstrates that slosh frequency is a combination of liquid slosh frequency and gas phase slosh frequency. Additionally, it shows that slosh mass is reduced due to the opposing motion of the gas phase with respect to that of the liquid phase. The slosh model parameter trends were verified using Computational Fluid Dynamics (CFD) analysis results for various liquid-to-gas density ratios, tanks, and fill levels.

Christopher D Moore↗