Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Cases”

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 1,477 records · Page 82

Algorithms for Autonomous GPS Orbit Determination and Formation Flying: Investigation of Initialization Approaches and Orbit Determination for HEO

This report summarizes the efforts to date in processing GPS measurements in High Earth Orbit (HEO) applications by the Colorado Center for Astrodynamics Research (CCAR). Two specific projects were conducted; initialization of the orbit propagation software, GEODE, using nominal orbital elements for the IMEX orbit, and processing of actual and simulated GPS data from the AMSAT satellite using a Doppler-only batch filter. CCAR has investigated a number of approaches for initialization of the GEODE orbit estimator with little a priori information. This document describes a batch solution approach that uses pseudorange or Doppler measurements collected over an orbital arc to compute an epoch state estimate. The algorithm is based on limited orbital element knowledge from which a coarse estimate of satellite position and velocity can be determined and used to initialize GEODE. This algorithm assumes knowledge of nominal orbital elements, (a, e, i, omega, omega) and uses a search on time of perigee passage (tau(sub p)) to estimate the host satellite position within the orbit and the approximate receiver clock bias. Results of the method are shown for a simulation including large orbital uncertainties and measurement errors. In addition, CCAR has attempted to process GPS data from the AMSAT satellite to obtain an initial estimation of the orbit. Limited GPS data have been received to date, with few satellites tracked and no computed point solutions. Unknown variables in the received data have made computations of a precise orbit using the recovered pseudorange difficult. This document describes the Doppler-only batch approach used to compute the AMSAT orbit. Both actual flight data from AMSAT, and simulated data generated using the Satellite Tool Kit and Goddard Space Flight Center's Flight Simulator, were processed. Results for each case and conclusion are presented.

Axelrad, Penina↗

Vertical Field of View Reference Point Study for Flight Path Control and Hazard Avoidance

Researchers within the eXternal Visibility System (XVS) element of the High-Speed Research (HSR) program developed and evaluated display concepts that will provide the flight crew of the proposed High-Speed Civil Transport (HSCT) with integrated imagery and symbology to permit path control and hazard avoidance functions while maintaining required situation awareness. The challenge of the XVS program is to develop concepts that would permit a no-nose-droop configuration of an HSCT and expanded low visibility HSCT operational capabilities. This study was one of a series of experiments exploring the 'design space' restrictions for physical placement of an XVS display. The primary experimental issues here was 'conformality' of the forward display vertical position with respect to the side window in simulated flight. 'Conformality' refers to the case such that the horizon and objects appear in the same relative positions when viewed through the forward windows or display and the side windows. This study quantified the effects of visual conformality on pilot flight path control and hazard avoidance performance. Here, conformality related to the positioning and relationship of the artificial horizon line and associated symbology presented on the forward display and the horizon and associated ground, horizon, and sky textures as they would appear in the real view through a window presented in the side window display. No significant performance consequences were found for the non-conformal conditions.

Comstock, J. Raymond, Jr.↗

Overview of Transonic to Hypersonic Stage Separation Tool Development for Multi-Stage-to-Orbit Concepts

An overview of research efforts at NASA in support of the stage separation and ascent aerothermodynamics research program is presented. The objective of this work is to develop a synergistic suite of experimental, computational, and engineering tools and methods to apply to vehicle separation across the transonic to hypersonic speed regimes. Proximity testing of a generic bimese wing-body configuration is on-going in the transonic (Mach numbers 0.6, 1.05, and 1.1), supersonic (Mach numbers 2.3, 3.0, and 4.5) and hypersonic (Mach numbers 6 and 10) speed regimes in four wind tunnel facilities at the NASA Langley Research Center. An overset grid, Navier-Stokes flow solver has been enhanced and demonstrated on a matrix of proximity cases and on a dynamic separation simulation of the bimese configuration. Steady-state predictions with this solver were in excellent agreement with wind tunnel data at Mach 3 as were predictions via a Cartesian-grid Euler solver. Experimental and computational data have been used to evaluate multi-body enhancements to the widely-used Aerodynamic Preliminary Analysis System, an engineering methodology, and to develop a new software package, SepSim, for the simulation and visualization of vehicle motions in a stage separation scenario. Web-based software will be used for archiving information generated from this research program into a database accessible to the user community. Thus, a framework has been established to study stage separation problems using coordinated experimental, computational, and engineering tools.

Murphy, Kelly J.↗

Advanced Combustion Modeling for Complex Turbulent Flows

The next generation of aircraft engines will need to pass stricter efficiency and emission tests. NASA's Ultra-Efficient Engine Technology (UEET) program has set an ambitious goal of 70% reduction of NO(x) emissions and a 15% increase in fuel efficiency of aircraft engines. We will demonstrate the state-of-the-art combustion tools developed a t Stanford's Center for Turbulence Research (CTR) as part of this program. In the last decade, CTR has spear-headed a multi-physics-based combustion modeling program. Key technologies have been transferred to the aerospace industry and are currently being used for engine simulations. In this demo, we will showcase the next-generation combustion modeling tools that integrate a very high level of detailed physics into advanced flow simulation codes. Combustor flows involve multi-phase physics with liquid fuel jet breakup, evaporation, and eventual combustion. Individual components of the simulation are verified against complex test cases and show excellent agreement with experimental data.

Ham, Frank Stanford↗

Single-Shot Scalar-Triplet Measurements in High-Pressure Swirl-Stabilized Flames for Combustion Code Validation

In support of NASA ARMD's code validation project, we have made significant progress by providing the first quantitative single-shot multi-scalar data from a turbulent elevated-pressure (5 atm), swirl-stabilized, lean direct injection (LDI) type research burner operating on CH4-air using a spatially-resolved pulsed-laser spontaneous Raman diagnostic technique. The Raman diagnostics apparatus and data analysis that we present here were developed over the past 6 years at Glenn Research Center. From the Raman scattering data, we produce spatially-mapped probability density functions (PDFs) of the instantaneous temperature, determined using a newly developed low-resolution effective rotational bandwidth (ERB) technique. The measured 3-scalar (triplet) correlations, between temperature, CH4, and O2 concentrations, as well as their PDF s, also provide a high-level of detail into the nature and extent of the turbulent mixing process and its impact on chemical reactions in a realistic gas turbine injector flame at elevated pressures. The multi-scalar triplet data set presented here provides a good validation case for CFD combustion codes to simulate by providing both average and statistical values for the 3 measured scalars.

Kojima, Jun↗

Comparison of PLIF and CFD Results for the Orion CEV RCS Jets

Nitric-oxide planar laser-induced fluorescence (NO PLIF) was used to visualize and measure centerline streamwise velocity of the Orion Crew Exploration Vehicle (CEV) Reaction Control System (RCS) Jets at NASA Langley Research Center's 31-Inch Mach 10 Air wind tunnel. Fluorescence flow visualizations of pitch, roll, and yaw RCS jets were obtained using different plenum pressures and wind tunnel operating stagnation pressures. For two yaw RCS jet test cases, the PLIF visualizations were compared to computational flow imaging (CFI) images based on Langley Aerothermal Upwind Relaxation Algorithm (LAURA) computational fluid dynamics (CFD) simulations of the flowfield. For the same test cases, the streamwise velocity measurements were compared to CFD. The CFD solution, while showing some unphysical artifacts, generally agree with the experimental measurements.

Ivey, Christopher B.↗

HEMT Frequency Doubler with Output at 300 GHz

An active frequency doubler in the form of an InP-based monolithic microwave integrated circuit (MMIC) containing a high-electron-mobility transistor (HEMT) has been demonstrated in operation at output frequencies in the vicinity of 300 GHz. This is the highestfrequency HEMT doubler reported to date, the next-highest-frequency active HEMT doubler having been previously reported to operate at 180 GHz. While the output power of this frequency doubler is less than that of a typical Schottky diode, this frequency doubler is considered an intermediate product of a continuing effort to realize the potential of active HEMT frequency doublers to operate with conversion efficiencies greater than those of passive diode frequency doublers. An additional incentive for developing active HEMT frequency doublers lies in the fact that they can be integrated with amplifiers, oscillators, and other circuitry on MMIC chips. The circuitry of the doubler MMIC (see Figure 1) features grounded coplanar waveguides. Air bridges and vias are used to make contact with the ground plane. The HEMT is biased for Class-A operation (in which current is conducted throughout each cycle of oscillation), which would ordinarily be better suited to linear amplification than to frequency doubling. Ordinarily, class-B operation (in which current is conducted during about half of each cycle of oscillation) would be more suitable for frequency doubling because of the essential nonlinearity of partial-cycle conduction. The reason for the unusual choice of class A was that computational simulations had shown that in this case, the efficiency in class B would be less than in class A. The input matching circuit of this doubler includes transmission lines that afford a good impedance match at the fundamental frequency, plus an open stub to prevent leakage of the second harmonic through the input terminals. The output circuit was designed to suppress the fundamental while providing a good match for the second harmonic.

Samoska, Lorene↗

Mixed Layer Sub-Mesoscale Parameterization - Part 1: Derivation and Assessment

Several studies have shown that sub-mesoscales (SM 1km horizontal scale) play an important role in mixed layer dynamics. In particular, high resolution simulations have shown that in the case of strong down-front wind, the re-stratification induced by the SM is of the same order of the de-stratification induced by small scale turbulence, as well as of that induced by the Ekman velocity. These studies have further concluded that it has become necessary to include SM in ocean global circulation models (OGCMs), especially those used in climate studies. The goal of our work is to derive and assess an analytic parameterization of the vertical tracer flux under baroclinic instabilities and wind of arbitrary directions and strength. To achieve this goal, we have divided the problem into two parts: first, in this work we derive and assess a parameterization of the SM vertical flux of an arbitrary tracer for ocean codes that resolve mesoscales, M, but not sub-mesoscales, SM. In Part 2, presented elsewhere, we have used the results of this work to derive a parameterization of SM fluxes for ocean codes that do not resolve either M or SM. To carry out the first part of our work, we solve the SM dynamic equations including the non-linear terms for which we employ a closure developed and assessed in previous work. We present a detailed analysis for down-front and up-front winds with the following results: (a) down-front wind (blowing in the direction of the surface geostrophic velocity) is the most favorable condition for generating vigorous SM eddies; the de-stratifying effect of the mean flow and re-stratifying effect of SM almost cancel each other out,

Canuto, V. M.↗

An Efficient Deterministic Approach to Model-based Prediction Uncertainty Estimation

Prognostics deals with the prediction of the end of life (EOL) of a system. EOL is a random variable, due to the presence of process noise and uncertainty in the future inputs to the system. Prognostics algorithm must account for this inherent uncertainty. In addition, these algorithms never know exactly the state of the system at the desired time of prediction, or the exact model describing the future evolution of the system, accumulating additional uncertainty into the predicted EOL. Prediction algorithms that do not account for these sources of uncertainty are misrepresenting the EOL and can lead to poor decisions based on their results. In this paper, we explore the impact of uncertainty in the prediction problem. We develop a general model-based prediction algorithm that incorporates these sources of uncertainty, and propose a novel approach to efficiently handle uncertainty in the future input trajectories of a system by using the unscented transformation. Using this approach, we are not only able to reduce the computational load but also estimate the bounds of uncertainty in a deterministic manner, which can be useful to consider during decision-making. Using a lithium-ion battery as a case study, we perform several simulation-based experiments to explore these issues, and validate the overall approach using experimental data from a battery testbed.

Daigle, Matthew J.↗

Ice Accretion Modeling using an Eulerian Approach for Droplet Impingement

A three-dimensional Eulerian analysis has been developed for modeling droplet impingement on lifting bodes. The Eulerian model solves the conservation equations of mass and momentum to obtain the droplet flow field properties on the same mesh used in CFD simulations. For complex configurations such as a full rotorcraft, the Eulerian approach is more efficient because the Lagrangian approach would require a significant amount of seeding for accurate estimates of collection efficiency. Simulations are done for various benchmark cases such as NACA0012 airfoil, MS317 airfoil and oscillating SC2110 airfoil to illustrate its use. The present results are compared with results from the Lagrangian approach used in an industry standard analysis called LEWICE.

Kim, Joe Woong↗

Progress on the Development of a Simulation Environment for Optical Communications

Deep space laser communications require extremely accurate beam pointing to take advantage of the narrow beams achievable at optical wavelengths. This pointing accuracy must be achieved in the presence of spacecraft basebody motion which may exceed laser pointing requirements by orders of magnitude.In this paper a model of an optical band transceiver pointing control system is developed that can be used to predict performance under various operating scenarios. The transceiver model consists of an electro-mechanical model of the telescope platform and isolator. A novel photon counting detector array is used in the simulation as the focal plane detector. In the simulation we are able to inject various cases of spacecraft basebody motion based on both flight data and future mission spacecraft jitter requirements. Various models of uplink beacon flux levels and atmospheric scintillation are also available for analysis.Using these models, detector processing and control functions are implemented in the simulation. A complete acquisition sequence is demonstrated with blind acquisition and tracking of the modulated uplink beam and positioning of the downlink beam on the focal plane array. These simulations predict that pointing requirements will be met with representative disturbance models and uplink beam scintillation.

pointing control↗

EMI from Spacecraft Docking Systems Spacecraft Charging - Plasma Contact Potentials

The plasma contact potential of a visiting vehicle (VV), such as the Orion Service Module (SM), is determined while docking at the Orion Crew Exploration Vehicle (CEV). Due to spacecraft charging effects on-orbit, the potential difference between the CEV and the VV can be large at docking, and an electrostatic discharge (ESD) could occur at capture, which could degrade, disrupt, damage, or destroy sensitive electronic equipment on the CEV and/or VV. Analytical and numerical models of the CEV are simulated to predict the worst-case potential difference between the CEV and the VV when the CEV is unbiased (solar panels unlit: eclipsed in the dark and inactive) or biased (solar panels sunlit: in the light and active).

Norgard, John D.↗

The Impact of MISR-Derived Injection Height Initialization on Wildfire and Volcanic Plume Dispersion in the HYSPLIT Model

The dispersion of particles from wildfires, volcanic eruptions, dust storms, and other aerosol sources can affect air quality and other environmental factors downwind. Aerosol injection height is one source attribute that mediates downwind dispersion, as wind speed and direction can vary dramatically with elevation. Using plume heights derived from space-based, multi-angle imaging, we examine the impact of initializing plumes with satellite-measured vs. nominal (model-calculated or VAAC observations) injection height on the simulated dispersion of six large aerosol plumes. When there are significant differences in nominal vs. satellite-derived particle injection heights, or if one injection height is within the planetary boundary layer (PBL) and the other is above the PBL, differences in simulation results can arise. In the cases studied with significant nominal vs. satellite-derived injection height differences, the NOAA Air Resources Laboratory's Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model tends to represent plume evolution better if the injection height in the model is constrained by hyper-stereo satellite retrievals.

Plume heights↗

A Geofence Violation Prevention Mechanism for Small UAS

The ability to safely confine the trajectories of small UAS to a specific geographical area is a key enabler for capabilities that require operating in close proximity to populated areas as well as other users of the airspace. These capabilities require highly reliable geofencing algorithms. In particular, these algorithms must promptly alert imminent breaches of keep-in/keep-out geofences by considering factors such as the vehicle speed and uncertainties in the state of the aircraft. This paper presents a novel approach to the prevention of geofence boundary violation based on closure rate constraints. These constraints are incorporated into a control framework to effectively prevent fence breaches. Simulation results illustrating an example use case of this framework are presented.

Balachandran, Swee↗

InSight's Reconstructed Aerothermal Environments

The InSight Mars Lander successfully landed on the surface on November 26, 2018. This poster will describe the methodologies and margins used in developing the aerothermal environments for design of the thermal protection systems (TPS), as well as a prediction of as-flown environments based on the best estimated trajectory. The InSight mission spacecraft design approach included the effects of radiant heat flux to the aft body from the wake for the first time on a US Mars Mission, due to overwhelming evidence in ground testing for the European ExoMars mission (2009/2010) [1] and 2010 tests in the Electric Arc Shock Tube (EAST) facility [2]. The radiant energy on an aftbody was also recently confirmed via measurement on the Schiaparelli mission [3]. In addition, the InSight mission expected to enter the Mars atmosphere during the dust storm season, so the heatshield TPS was designed to accommodate the extra recession due to the potential dust impact. This poster will compare the predicted aerothermal environments using the reconstructed best estimated trajectory to the design environments. Design Approach: The InSight spacecraft was planned to be a near-design-to-print copy of the Phoenix spacecraft. The determination of the heatshield TPS requirements was approached as if it was a new design due to the new requirement of flying through a dust storm. The baseline for aftbody was build-to-print, and all analyses focused on ensuring adequate margin. This proved to be a challenge because the Phoenix aftbody was designed to withstand only convective heating and the InSight aftbody was evaluated for both convective and radiative heating. Aerothermal environments were predicted using the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) and the Data Parallel Line Relaxation (DPLR) CFD codes, and the Nonequilibrium Radiative Transport and Spectra Program (NEQAIR) utilizing bounding design trajectories derived from Monte Carlo analyses from the Program to Optimize Simulated Trajectories II (POST2). In all cases, super-catalytic flowfields were assigned to ensure the most conservative heating results. Two trajectories were evaluated: 1) the trajectory with the maximum heat flux was utilized to determine the flowfield characteristics and the viability of the selection of TPS materials; and 2) the trajectory with the maximum heat load was used to determine the required thicknesses of the TPS materials. Evaluation of the MEDLI data [4], along with ground test data [5] led to the determination of whether or not the flow would transition from laminar to turbulent on the heatshield, which also determined the TPS sizing location for the heatshield. Aerothermal margins were added for the convective heating and developed for the radiative heating. TPS material sizing was determined with the Reaction Kinetic Ablation Program (REKAP) and the Fully Implicit Ablation and Thermal Analysis program (FIAT) using a three-branched approach to account for aerothermal, material response, and material properties uncertainties. In addition, the heatshield recession was augmented by an analysis of the effect of entry through a potential dusty atmosphere using a methodology developed in References [6] and [7]. These analyses resulted in an increase to the Phoenix heatshield TPS thickness. Reconstruction Efforts: Once the best estimated trajectory is reconstructed by the team, the LAURA/HARA (High-Temperature Aerothermo-dynamic Radiation model) and DPLR/NEQAIR code pairs will be used to predict the as-flown aerothermal conditions. In these runs, fully-catalytic flowfields will be assigned because it is a more physically accurate description of the chemistry in the flow. Once again, determination of the onset of turbulence on the heatshield will be evaluated. The as-flown aerothermal environments will then be compared to the design environments.

Beck, R. A.↗

Cartesian Mesh Simulations for the Third AIAA Sonic Boom Prediction Workshop

Simulation results are presented for all cases from the Third AIAA Sonic Boom Prediction Workshop. An inviscid, embedded-boundary Cartesian-mesh flow solver is used in conjunction with adjoint-based mesh adaptation to compute nearfield pressure signatures. Specialized techniques are applied to maximize accuracy and minimize cost on Cartesian meshes. The Richardson-based error estimate highlights regions of the signatures most sensitive to mesh refinement. Timing results and coarse, medium, and fine mesh sizes for nearfield cases demonstrate that the parallel decomposition approach is efficient in both computational time and wall-clock. Pressure signals are propagated to the ground using an augmented Burgers’ equation solver to predict boom carpets. Ground signatures and loudness metrics are presented for a standard atmosphere as well as more realistic atmospheric profiles, which affect overall noise levels and can significantly widen the boom carpet. Mesh convergence studies show that high sampling frequencies, around 500 kHz, are required for propagation, and the sampling frequency increases at large off-track angles with longer acoustic ray paths and propagation times. The numerical methods yield accurate results for predicting low sonic boom signatures while being among the least computationally expensive of the workshop.

ARMD↗

Computational Aerodynamic Analysis in Support of the CRM Tail Cone Thruster Configuration Wind Tunnel Test

NASA’s Advanced Air Transport Technology (AATT) project is breaching the boundaries of aircraft design in pursuit of eco-friendly solutions that are compatible with urban noise comfort levels. Boundary layer ingesting (BLI) propulsion systems promise to reduce fuel burn with additional potential benefits in noise reduction. Type-II BLI systems of the STARC-ABL type are the subject of a test campaign planned for fiscal year 2022 in the National Transonic Facility (NTF), for which a CRM-based model with a retrofitted tail cone thruster (TCT) has been designed. The present work is a precursor to the NTF test, where the 240 cases planned for the experiment were simulated using the Launch, Ascent and Vehicle Aerodynamics (LAVA)computational framework. Solution sensitivity to angle of attack, engine operating conditions, and the presence of the supporting structure (sting) in the wind tunnel test are analyzed in the extensive dataset. The main flow features contributing to the inlet distortion are identified as the vertical tail wake, wing downwash and fuselage upsweep vortices, with the latter two experiencing the greatest sensitivity to angle of attack. Finally, results from an inlet-guide-vane(IGV) design/integration study are presented. The LAVA team and the turbomachinery design team at NASA’s Ames and Glenn Research Centers are collaborating in an effort to reduce flow distortion upstream of the fan by means of integrating an IGV system into the CRM-TCT model. A significant improvement in flow distortion metrics has been achieved since the initial design iteration. Results employing an actuator zone model with realistic radially-varying thrust profiles to simulate first-order fan effects within LAVA are presented.

AATT↗

The simulation of the impact of the spatial distribution of vegetation on the urban microclimate: A case study in Mostaganem

Increasing vegetation is one of the solutions proposed to mitigate the phenomenon of urban heat islands. The present research aims to study the effect of the spatial distribution of vegetation on the urban microclimate during the summer period. The city of Mostaganem with its Mediterranean climate was chosen as a case study. Based on meteorological data collected on a hot day in July, the study is carried out using the ENVI-met 4.4.5. Scenarios of different vegetation layouts were carried out. Air, mean radiant and surface temperatures are the microclimatic parameters that have been analyzed for each of the scenarios, and in order to study the outdoor thermal comfort, the Predicted Mean Vote (PMV) is chosen as the comfort index, it compares the average value of the votes of a large group of people on the thermal sensation scale. The study revealed that tree alignments can have a cooling effect on the city temperature close to 1.2◦C. However, its impact on reducing the surface temperature can reach up to 4 ◦C. The introduction of vegetation on the site had the overall effect of lowering the air temperature and therefore relatively improving the level of comfort compared to the reference scenario. But, in order to reduce discomfort to an acceptable level during the day, vegetation can be simultaneously accommodated by city parameters such as urban form.

Nora Bachir↗