Check-Cases for Lunar Six-Degree-of-Freedom (6DOF) Simulations
Explore the source record for details and available documents.
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.
Explore the source record for details and available documents.
We present an implementation of a propulsor model based on body force method into the OVERFLOW computational fluid dynamics code to model turbofan engines and propulsors of similar type. The model estimates the forces imparted on the fluid by the blade camber surfaces as a body force source terms in the momentum and energy equations over grids that represent the rotor and stator of a fan stage. We tested the implementation on three cases: 1) Source Diagnostics Test (SDT) fan with R4 rotors, which has extensive test data on aerodynamic performance and rotor wake surveys, 2) A stand-alone Aeronaut TF8000 propulsor which is an off-the-shelf propulsor used on model aircraft and 3) The D8 aircraft model with TF8000 propulsors placed in a wind tunnel. Despite missing some of the features in the velocity profiles near the endwalls (i.e. hub and casing surfaces), the OVERFLOW simulations with body force model predicted area-averaged flow speed and total pressure rise through the SDT propulsor within a few percent of the LDV (Laser Doppler Velocimetry) measurements. In the case of TF8000 propulsor on the D8 airframe, the model under-predicted mechanical power coefficient by several percent of the wind tunnel test results when the horizontal force balance condition over the airframe is targeted by tuning rotor speed. By investigating the upstream influence of the rotor swirl, it was found out that the induced swirl velocity effects upstream were relatively small and they rapidly vanished before reaching one fan diameter upstream of the fan face. The body force model provided insights on aerothermodynamics and aeromechanics of boundary layer ingesting propulsor; these insights could not be obtained by using the uniform pressure jump model.
On 8 July 1989, a very strong microburst was detected by the Low-Level Windshear Alert system (LLWAS), within the approach corridor just north of Denver Stapleton Airport. The microburst was encountered by a Boeing 737-200 in a 'go-around' configuration which was reported to have lost considerable air speed and altitude during penetration. Data from LLWAS revealed a pulsating microburst with an estimated peak velocity change of 48 m/s. Wilson et al. reported that the microburst was accompanied by no apparent visible clues such as rain or virga, although blowing dust was present. Weather service hourly reports indicated virga in all quadrants near the time of the event. A National Center for Atmospheric Research (NCAR) research Doppler radar was operating; but according to Wilson et al., meaningful velocity could not be measured within the microburst due to low radar-reflectivity factor and poor siting for windshear detection at Stapleton. This paper presents results from the three-dimensional numerical simulation of this event, using the Terminal Area Simulation System (TASS) model. The TASS model is a three-dimensional nonhydrostatic cloud model that includes parameterizations for both liquid and ice phase microphysics, and has been used in investigations of both wet and dry microburst case studies. The focus of this paper is the pulsating characteristic and the very-low radar reflectivity of this event. Most of the surface outflow contained no precipitation. Such an event may be difficult to detect by radar.
A thermal test of an enclosure housing the AMPS electronics cards was performed. The enclosure was fully populated with 17 nonfunctioning AMPS electronics cards. Heaters were applied to the cards to simulate their operation. A cold plate was installed on the top of the case. The case was then insulated to simulate its installation within an electronics rack with the face and rear of the enclosure left open. The cold plate was connected to a chiller which provided cooling water at a fixed input temperature to the cold plate. Flow rates from the chiller of 0.0 L/min, 1.4 L/min and 1.8 L/min were tested. The test was operated at each flow rate over a range of heater power levels until the steady state temperature of the electronics cards was reached. Heater power levels from 2.5 W/card to 45 W/card were tested. The data collected was used to determine the steady state operating temperature of the cards, the temperature distribution between the cards, the time to reach steady-state temperature, the change in steady state temperature with change in heater power, the temperature distribution between the cards at the steady state condition and the heat lost to the surrounds by the enclosure. One of the main goals of the testing was to determine how much waste heat the cards could generate and still maintain a touch temperature at or below 40°C. Based on this requirement, for the 3 cooling rates tested, the maximum card waste heat that could be generated by each individual card was 3.9 W, 16.3 W and 17.6 W respectively. An analysis was also set up to provide an analytical model to estimate the card temperature under different cold plate flow rates and card heater power levels. The results off the analysis was compared to that of the experiment of the range of flow rates and power levels tested.
Detailed design issues associated with liquid rocket engine injectors and combustion chamber operation require CFD methodology which simulates highly three-dimensional, turbulent, vaporizing, and combusting flows. The primary utility of such simulations involves predicting multi-dimensional effects caused by specific injector configurations. SECA, Inc. and Engineering Sciences, Inc. have been developing appropriate computational methodology for NASA/MSFC for the past decade. CFD tools and computers have improved dramatically during this time period; however, the physical submodels used in these analyses must still remain relatively simple in order to produce useful results. Simulations of clustered coaxial and impinger injector elements for hydrogen and hydrocarbon fuels, which account for real fluid properties, is the immediate goal of this research. The spray combustion codes are based on the FDNS CFD code' and are structured to represent homogeneous and heterogeneous spray combustion. The homogeneous spray model treats the flow as a continuum of multi-phase, multicomponent fluids which move without thermal or velocity lags between the phases. Two heterogeneous models were developed: (1) a volume-of-fluid (VOF) model which represents the liquid core of coaxial or impinger jets and their atomization and vaporization, and (2) a Blob model which represents the injected streams as a cloud of droplets the size of the injector orifice which subsequently exhibit particle interaction, vaporization, and combustion. All of these spray models are computationally intensive, but this is unavoidable to accurately account for the complex physics and combustion which is to be predicted, Work is currently in progress to parallelize these codes to improve their computational efficiency. These spray combustion codes were used to simulate the three test cases which are the subject of the 2nd International Workshop on-Rocket Combustion Modeling. Such test cases are considered by these investigators to be very valuable for code validation because combustion kinetics, turbulence models and atomization models based on low pressure experiments of hydrogen air combustion do not adequately verify analytical or CFD submodels which are necessary to simulate rocket engine combustion. We wish to emphasize that the simulations which we prepared for this meeting are meant to test the accuracy of the approximations used in our general purpose spray combustion models, rather than represent a definitive analysis of each of the experiments which were conducted. Our goal is to accurately predict local temperatures and mixture ratios in rocket engines; hence predicting individual experiments is used only for code validation. To replace the conventional JANNAF standard axisymmetric finite-rate (TDK) computer code 2 for performance prediction with CFD cases, such codes must posses two features. Firstly, they must be as easy to use and of comparable run times for conventional performance predictions. Secondly, they must provide more detailed predictions of the flowfields near the injector face. Specifically, they must accurately predict the convective mixing of injected liquid propellants in terms of the injector element configurations.
Direct numerical simulation of turbulent reacting flows places extreme demands on computational resources. At the present time, simulations can be performed only for greatly simplified reaction systems and for very low Reynolds numbers. Direct simulation of more realistic cases occurring at higher Reynolds number and including multiple species and numerous chemical reactions will exceed available computational resources far into the future. Because of this, there is a clear need to develop the technique of large eddy simulation for reacting flows. Unfortunately this task is complicated by the fact that combustion arises from chemical reactions that occur at the smallest scales of the flow. Capturing the large-scale behavior without resolving the small-scale details is extremely difficult in combustion problems. Thus LES modeling for turbulent combustion encounters difficulties not present in modeling momentum transport, in which the main effect of the small scales is to provide dissipation. The difficulty is more pronounced in premixed combustion, where detailed chemistry plays an essential role in determining the flame speed (or overall burning rate); in nonpremixed combustion infinite rate chemistry can be assumed, eliminating the small scale features to a first approximation.
Accurately modeling wind farm performance in complex atmospheric flows remains a challenge. This paper presents the modeling results of the American WAKE experimeNt (AWAKEN) wind farm benchmark, a collaborative effort involving 16 research groups from academia and industry within the International Energy Agency Wind Technology Collaboration Programme Task 57. The study evaluates a diverse suite of simulation tools, ranging from fast-running engineering wake models to high-fidelity large-eddy simulations, against a diurnal case study observed during the AWAKEN campaign. The benchmark utilized a three-phase structure to progressively assess model performance as observational data availability increased. Initial blind predictions showed that higher-fidelity models did not uniformly outperform simpler simulation tools. A distinct spatial bias was observed where models struggled to resolve the interplay between a low-level jet, wakes, and terrain-induced flow acceleration. In subsequent phases, leveraging additional measurements for model improvement led to a reduction in mean absolute error across the model ensemble; however, this effect was most pronounced in engineering wake models, where targeted calibration reduced error by up to 40~\%. Overall, the study demonstrates that inflow characterization remains a primary prerequisite for accuracy, particularly for models relying on coarse forcing datasets. While the limited ability to resolve local terrain-flow interactions under single-day conditions represent a recognized constraint, the overall findings on wake modeling and real-world validation still provide valuable guidance for model application and for mitigating this limitation.
Microphysical processes in deep convective clouds are sensitive to the number concentrations of cloud condensation nuclei (CCN) and ice nucleating particles (INPs), but the effects of INPs are less studied. Modeling studies investigating the effects of INPs and/or CCN on deep convection typically retain a volume‐dependent raindrop freezing relation. The resulting neglect of aerosol accumulation in raindrops via drop collisions has likely produced unrealistic storm responses to INPs in past studies. To address this deficiency, a new immersion freezing algorithm was developed and embedded in a bulk microphysics scheme that freezes both cloud drops and raindrops using the same immersion freezing INP (IF‐INP) activity spectrum based on measurements. Multiple idealized simulations of a single case of deep convection observed during the Clouds, Aerosols, and Complex Terrain Interactions (CACTI) field campaign were conducted, with microphysical differences produced by independently altering IF‐INP temperature dependencies and CCN number concentrations from their observed values. Surface precipitation in all simulations resulted almost exclusively from riming graupel that melted upon descending to the surface. Rainfall and cold pools were substantially and systematically weakened with increased CCN due to decreased graupel riming rates but were relatively insensitive to variations in the magnitude and slope of IF‐INP spectra due to compensating depletion of supercooled liquid water. These compensating processes were a consequence of the accumulation of IF‐INPs in raindrops, encouraging caution in studying IF‐INP effects upon thunderstorms using traditional volume‐dependent drop freezing relationships.
This paper describes the power-hardware-in-the- loop (PHIL) evaluation of the feasibility of service restoration solutions determined by the PowerModelsONM.jl tool. This tool incorporates microgrids and the networking of microgrids into its determination of an optimal service restoration solution. The paper presents PHIL simulation results for a case study based on a real distribution feeder with multiple microgrids, showcas- ing the effectiveness of networked microgrids in aiding system restoration after an outage. The study leverages high-fidelity, real- time electromagnetic transient models to ensure accuracy in the simulation results. This work is the final output from the Resilient Operation of Networked Microgrids (RONM) project funded by the U.S. Department of Energy Office of Electricity Microgrid Program and led by Los Alamos National Laboratory. RONM focused on the application of PowerModelsONM.jl to enhance the resilience of distribution systems.
Automobiles approaching a controlled intersection need to brake and come to a full stop when the signal transitions from green to red, and the vehicle must later accelerate to normal speeds after the signal changes back to green. These stops and starts associated with normal signal changes lead to unnecessary energy consumption and vehicle emissions. Previous studies have revealed that the optimization of traffic intersections' signal controls and coordination facilitates smoother traffic flows with reduced stop-and-go driving, which can significantly reduce traffic congestion and unnecessary fuel waste. This paper presents an interactive visual analytics dashboard that allows transportation planners to explore and analyze energy consumption patterns resulting from temporally varying traffic signal phases at multiple intersections along a major transportation corridor using traffic simulation outputs. The visual dashboard is implemented as an accessible and responsive web application and employs a combination of visualization techniques to cover multiple aspects of vehicle acceleration and braking at multiple adjacent intersections along a corridor. The paper presents a case study of a simulated traffic scenario on the Shallowford Road traffic corridor located in Chattanooga, Tennessee to demonstrate the capability of the visual dashboard.
The DSS Application is a computer program comprising a Windows version of the UNIX-based Decelerator System Simulation (DSS) coupled with an Excel front end. The DSS is an executable code that simulates the dynamics of airdropped cargo from first motion in an aircraft through landing. The bare DSS is difficult to use; the front end makes it easy to use. All inputs to the DSS, control of execution of the DSS, and postprocessing and plotting of outputs are handled in the front end. The front end is graphics-intensive. The Excel software provides the graphical elements without need for additional programming. Categories of input parameters are divided into separate tabbed windows. Pop-up comments describe each parameter. An error-checking software component evaluates combinations of parameters and alerts the user if an error results. Case files can be created from inputs, making it possible to build cases from previous ones. Simulation output is plotted in 16 charts displayed on a separate worksheet, enabling plotting of multiple DSS cases with flight-test data. Variables assigned to each plot can be changed. Selected input parameters can be edited from the plot sheet for quick sensitivity studies.
ABSTRACT It has been recently shown numerically that there exists an inverse transfer of magnetic energy in decaying, non-helical, magnetically dominated, magnetohydrodynamic turbulence in three dimensions (3D). We suggest that magnetic reconnection is the underlying physical mechanism responsible for this inverse transfer. In the two-dimensional (2D) case, the inverse transfer is easily inferred to be due to smaller magnetic islands merging to form larger ones via reconnection. We find that the scaling behaviour is similar between the 2D and 3D cases, i.e. the magnetic energy evolves as t−1, and the magnetic power spectrum follows a slope of k−2. We show that on normalizing time by the magnetic reconnection time-scale, the evolution curves of the magnetic field in systems with different Lundquist numbers collapse on to one another. Furthermore, transfer function plots show signatures of magnetic reconnection driving the inverse transfer. We also discuss the conserved quantities in the system and show that the behaviour of these quantities is similar between the 2D and 3D simulations, thus making the case that the dynamics in 3D could be approximately explained by what we understand in 2D. Lastly, we also conduct simulations where the magnetic field is subdominant to the flow. Here, too, we find an inverse transfer of magnetic energy in 3D. In these simulations, the magnetic energy evolves as t−1.4 and, interestingly, a dynamo effect is observed.
The Hackett method (a Wall Pressure Signature Method) was selected to be adapted for the 12 ft Wind Tunnel WIAC system. This method uses limited measurements of the static pressure at the wall, in conjunction with the solid wall boundary condition, to determine the strength and distribution of singularities representing the test article. The singularities are used in term for estimating wall interference at the model location. Hackett's method will have to be formulated for application to the unique geometry of the 12 ft tunnel. The WIAC code will be validated by conducting numerically simulated experiments rather than actual wind tunnel experiments. The simulations will be used to generate both free air and confined wind tunnel flow fields for each of the test articles over a range of test configurations. Specifically the pressure signature at the test section wall will be computed for the confined case to provide the simulated 'measured' data. These data will serve as the input for the WIAC method. The performance of the WIAC method then may be evaluated by comparing the corrected parameters with those for the free air simulation.
Efforts to make increasingly use of suitable computer programs in the design of hardware have the potential to reduce expenditures. In this context, NASA has evaluated the benefits provided by software tools through an application to the Environmental Control and Life Support (ECLS) system. The present paper is concerned with the benefits obtained by an employment of simulation tools in the case of the Air Revitalization System (ARS) of a Space Station life support system. Attention is given to the ARS functions and components, a computer program overview, a SAND (solid amine water desorbed) bed model description, a model validation, and details regarding the simulation benefits.
A Wall Signature method, the Hackett method, has been selected to be adapted for the 12-ft Wind Tunnel wall interference assessment/correction (WIAC) system in the present phase. This method uses limited measurements of the static pressure at the wall, in conjunction with the solid wall boundary condition, to determine the strength and distribution of singularities representing the test article. The singularities are used in turn for estimating wall interferences at the model location. The Wall Signature method will be formulated for application to the unique geometry of the 12-ft Tunnel. The development and implementation of a working prototype will be completed, delivered and documented with a software manual. The WIAC code will be validated by conducting numerically simulated experiments rather than actual wind tunnel experiments. The simulations will be used to generate both free-air and confined wind-tunnel flow fields for each of the test articles over a range of test configurations. Specifically, the pressure signature at the test section wall will be computed for the tunnel case to provide the simulated 'measured' data. These data will serve as the input for the WIAC method-Wall Signature method. The performance of the WIAC method then may be evaluated by comparing the corrected parameters with those for the free-air simulation. Each set of wind tunnel/test article numerical simulations provides data to validate the WIAC method. A numerical wind tunnel test simulation is initiated to validate the WIAC methods developed in the project. In the present reported period, the blockage correction has been developed and implemented for a rectangular tunnel as well as the 12-ft Pressure Tunnel. An improved wall interference assessment and correction method for three-dimensional wind tunnel testing is presented in the appendix.
A Wall Signature method originally developed by Hackett has been selected to be adapted for the Ames 12-ft Wind Tunnel WIAC system in the project. This method uses limited measurements of the static pressure at the wall, in conjunction with the solid wall boundary condition, to determine the strength and distribution of singularities representing the test article. The singularities are used in turn for estimating blockage wall interference. The lifting interference will be treated separately by representing in a horseshoe vortex system for the model's lifting effects. The development and implementation of a working prototype will be completed, delivered and documented with a software manual. The WIAC code will be validated by conducting numerically simulated experiments rather than actual wind tunnel experiments. The simulations will be used to generate both free-air and confined wind-tunnel flow fields for each of the test articles over a range of test configurations. Specifically, the pressure signature at the test section wall will be computed for the tunnel case to provide the simulated 'measured' data. These data will serve as the input for the WIAC method--Wall Signature method. The performance of the WIAC method then may be evaluated by comparing the corrected data with those of the free-air simulation.
A Wall Signature method originally developed by Hackett was selected to be adapted for the Ames 12-ft Wind Tunnel Wall Interference Assessment/Correction (WIAC) System in the project. This method uses limited measurements of the static pressure at the wall, in conjunction with the solid wall boundary condition, to determine the strength and distribution of singularities representing the test article. The singularities are used in turn for estimating wall interference at the model location. The lifting interference will be treated separately by representing in a horseshoe vortex system for the model's lifting effects. The development and implementation of a working prototype will be completed, delivered, and documented with a software manual. The WIAC code will be validated by conducting numerically simulated experiments rather than actual wind tunnel experiments. The simulations will be used to generate both free-air and confined wind-tunnel flow fields for each of the test articles over a range of test configurations. Specifically, the pressure signature at the test section wall will be computed for the tunnel case to provide the simulated 'measured' data. These data will serve as the input for the WIAC method - Wall Signature method. The performance of the WIAC method then may be evaluated by comparing the corrected parameters with those for the free-air simulation. The following two additional tasks are included in the supplement No. 1 to the basic Grant. On-line wall interference calculation: The developed wall signature method (modified Hackett's method) for Ames 12-ft Tunnel will be the pre-computed coefficients which facilitate the on-line calculation of wall interference; and support system effects estimation: The effects on the wall pressure measurements due to the presence of the model support systems will be evaluated.
A Wall Signature method originally developed by Hackett has been selected to be adapted for the Ames 12-ft Wind Tunnel WIAC system in the project. This method uses limited measurements of the static pressure at the wall, in conjunction with the solid wall boundary condition, to determine the strength and distribution of singularities representing the test article. The singularities are used in turn for estimating wall interference at the model location. The development and implementation of a working prototype will be completed, delivered and documented with a software manual. The WIAC code will be validated by conducting numerically simulated experiments rather than actual wind tunnel experiments. The simulations will be used to generate both free-air and confined wind-tunnel flow fields for each of the test articles over a range of test configurations. Specifically, the pressure signature at the test section wall will be computed for the tunnel case to provide the simulated 'measured' data. These data will serve as the input for the WIAC method--Wall Signature method. The performance of the WIAC method then may be evaluated by comparing the corrected parameters with those for the free-air simulation. The following two additional tasks are included: (1) On-line wall interference calculation: The developed wall signature method (modified Hackett's method) for Ames 12-ft Tunnel will be the pre-computed coefficients which facilitate the on-line calculation of wall interference, and (2) Support system effects estimation: The effects on the wall pressure measurements due to the presence of the model support systems will be evaluated.