Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Trajectory Generation”

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,639 records · Page 91

Orbit Determination for the Lunar Reconnaissance Orbiter Using an Extended Kalman Filter

Orbit determination (OD) analysis results are presented for the Lunar Reconnaissance Orbiter (LRO) using a commercially available Extended Kalman Filter, Analytical Graphics' Orbit Determination Tool Kit (ODTK). Process noise models for lunar gravity and solar radiation pressure (SRP) are described and OD results employing the models are presented. Definitive accuracy using ODTK meets mission requirements and is better than that achieved using the operational LRO OD tool, the Goddard Trajectory Determination System (GTDS). Results demonstrate that a Vasicek stochastic model produces better estimates of the coefficient of solar radiation pressure than a Gauss-Markov model, and prediction accuracy using a Vasicek model meets mission requirements over the analysis span. Modeling the effect of antenna motion on range-rate tracking considerably improves residuals and filter-smoother consistency. Inclusion of off-axis SRP process noise and generalized process noise improves filter performance for both definitive and predicted accuracy. Definitive accuracy from the smoother is better than achieved using GTDS and is close to that achieved by precision OD methods used to generate definitive science orbits. Use of a multi-plate dynamic spacecraft area model with ODTK's force model plugin capability provides additional improvements in predicted accuracy.

Orbit Determination↗

Numerical Investigation of Film Coefficient Approximation for Chemically Reacting Boundary-Layer Flows

Aerothermal analysis of spacecraft planetary entry is heavily dependent on heritage engineering models. The film coefficient heat transfer model examined in this paper estimates the convective heating to the vehicle for a laminar, dissociated, chemically reacting boundary layer for an Earth atmosphere. This model requires information about the vehicle and flowfield for a given trajectory point and estimates a proportional relationship between enthalpy potential and convective heat flux. In practice it is the aerothermal engineer who must decide which assumptions are appropriate for his/her application. This work looks at numerous CFD simulations for an arbitrary, axisymmetric flight vehicle to analyze the relative importance of both the mass and energy constraints imposed at the wall boundary, as well as the effect of various diffusion models. Within the subset of tested energy boundary conditions, it is found that the most desirable energy boundary condition is the radiative equilibrium boundary condition, which permits conservative estimates of convective heat flux, but also generates flowfield-dependent spatial thermal distributions along the surface. Other key findings are presented in an effort to make the film coefficient engineering model readily available to design engineers across industry.

Film coefficient↗

Inference offers a metric to constrain dynamical models of neutrino flavor transformation

The multimessenger astrophysics of compact objects presents a vast range of environments where neutrino flavor transformation may occur and may be important for nucleosynthesis, dynamics, and a detected neutrino signal. Development of efficient techniques for surveying flavor evolution solution spaces in these diverse environments, which augment and complement existing sophisticated computational tools, could leverage progress in this field. To this end we continue our exploration of statistical data assimilation (SDA) to identify solutions to a small-scale model of neutrino flavor transformation. SDA is a machine learning formula wherein a dynamical model is assumed to generate any measured quantities. Specifically, we use an optimization formulation of SDA wherein a cost function is extremized via the variational method. Regions of state space in which the extremization identifies the global minimum of the cost function will correspond to parameter regimes in which a model solution can exist. Our example study seeks to infer the flavor transformation histories of two monoenergetic neutrino beams coherently interacting with each other and with a matter background. We require that the solution be consistent with measured neutrino flavor fluxes at the point of detection, and with constraints placed upon the flavor content at various locations along their trajectories, such as the point of emission, and the locations of the Mikheyev-Smirnov-Wolfenstein resonances. We show how the procedure efficiently identifies solution regimes and rules out regimes where solutions are infeasible. Overall, results in this work intimate the promise of this “variational annealing” methodology to efficiently probe an array of fundamental questions that traditional numerical simulation codes render difficult to access.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The generation and morphology of single-crystal silicon carbide wear particles under adhesive conditions

Sliding friction experiments were performed in vacuum at room temperature on a plane-type SiC surface in contact with iron-based binary alloys. Multiangular and spherical wear particles were found to form as a result of multipass sliding. The multiangular particles were produced by primary and secondary cracking of the 0001, 10(-)10, and 11(-)20 plane-type cleavage planes under the Hertzian stress field or local inelastic deformation zone. When alloy surfaces are in contact with silicon carbide under a load of 0.2 N, the alloy around the contact area is subjected to stresses that are close to the elastic limit in the elastic deformation region and/or exceed it. It was also found that spherical wear particles may be produced by two mechanisms: a penny-shaped fracture along the circular stress trajectories under the local inelastic deformation zone, and the attrition and fatigue of wear particles.

Miyoshi, K.↗

Flow-Field Measurement of Device-Induced Embedded Streamwise Vortex on a Flat Plate

Detailed flow-field measurements were performed downstream of a single vortex generator (VG) using an advanced Stereo Digital Particle Image Velocimetry system. Thc passive flow-control devices examined consisted of a low-profile VG with a device height, h, approximately equal to 20 percent of the boundary-layer thickness, sigma, and a conventional VG with h is approximately sigma. Flow-field data were taken at twelve cross-flow planes downstream of the VG to document and quantify the evolution of embedded streamwise vortex. The effects of device angle of attack on vortex development downstream were compared between the low-profile VG and the conventional VG. Key parameters including vorticity, circulation, trajectory, and half-life radius - describing concentration, strength, path, and size, respectively--of the device-induced streamwise vortex were extracted from the flow-field data. The magnitude of maximum vorticity increases as angle of attack increases for the low-profile VG, but the trend is reversed for the conventional VG, probably due to flow stalling around the larger device at higher angles of attack. Peak vorticity and circulation for the low-profile VG decays exponentially and inversely proportional to the distance downstream from the device. The device-height normalized vortex trajectories for the low-profile VG, especially in the lateral direction, follow the general trends of the conventional VG. The experimental database was used to validate the predictive capability of computational fluid dynamics (CFD). CFD accurately predicts the vortex circulation and path; however, improvements are needed for predicting the vorticity strength and vortex size.

Yao, Chung-Sheng↗

Time-Accurate Computational Fluid Dynamics Simulation of a Pair of Moving Solid Rocket Boosters

Since the Columbia accident, the threat to the Shuttle launch vehicle from debris during the liftoff timeframe has been assessed by the Liftoff Debris Team at NASA/MSFC. In addition to engineering methods of analysis, CFD-generated flow fields during the liftoff timeframe have been used in conjunction with 3-DOF debris transport methods to predict the motion of liftoff debris. Early models made use of a quasi-steady flow field approximation with the vehicle positioned at a fixed location relative to the ground; however, a moving overset mesh capability has recently been developed for the Loci/CHEM CFD software which enables higher-fidelity simulation of the Shuttle transient plume startup and liftoff environment. The present work details the simulation of the launch pad and mobile launch platform (MLP) with truncated solid rocket boosters (SRBs) moving in a prescribed liftoff trajectory derived from Shuttle flight measurements. Using Loci/CHEM, time-accurate RANS and hybrid RANS/LES simulations were performed for the timeframe T0+0 to T0+3.5 seconds, which consists of SRB startup to a vehicle altitude of approximately 90 feet above the MLP. Analysis of the transient flowfield focuses on the evolution of the SRB plumes in the MLP plume holes and the flame trench, impingement on the flame deflector, and especially impingment on the MLP deck resulting in upward flow which is a transport mechanism for debris. The results show excellent qualitative agreement with the visual record from past Shuttle flights, and comparisons to pressure measurements in the flame trench and on the MLP provide confidence in these simulation capabilities.

Strutzenberg, Louise L.↗

Functional Volume Assessment of an Early Version of the Mars Transit Habitat

During the summer of 2020, the NASA Mars Architecture Team (MAT) conducted a functional volume assessment of the Transit Habitat (TH). The TH has evolved substantially since that time and the current TH does not share the same internal configuration, but the insights gained from the assessment remain relevant. In virtually all architectures involving chemical, electric, or nuclear propulsion, the majority of the crew mission is spent aboard the TH. The Earth to Mars transit durations may vary from architecture to architecture, but all are on the order of hundreds of days, regardless of whether an opposition or conjunction class trajectory is selected, and regardless of whether the propulsion system is chemical, electric, or nuclear. Thus, the TH must provide capabilities appropriate to a very long duration. So, despite significant changes in the current Mars architecture over the past few years, this evaluation still provides useful recommendations in the form of habitability guidance that can be applied to current and future TH concepts. This assessment had six primary objectives: provide a sanity check to the BOC-derived TH layout; understand if we can fit the hardware and functional tasks in the volume; provide a high-level assessment of how aggressive the layout is; generate a list of challenges or assumptions necessary to make it work; generate a list of future work to refine understanding; and identify proposed requirements. All of this information is critical to drive habitat sizing studies and key architecture decisions. It was clear that a human-in-the-loop (HITL) evaluation of some kind would be necessary, but several challenges were immediately identified. The most information could be gleaned from a Desert Research and Technology Studies (DRATS) type of mission operations test (MOT), but the TH concept was too low fidelity to construct the type of prototype necessary to conduct a MOT, nor were their financial resources to do so. The next best option – and really the one most appropriate for the BOC’s stage of maturity – is a Virtual Reality (VR) walk-through evaluation. However, NASA was in a shutdown state due to the COVID-19 pandemic and the VR labs were inaccessible. As a result, a tabletop evaluation was created, using .jpg imagery from the ECM CAD model along with Excel-based datasheets. High-level crew living and working functions within the habitat were identified for evaluation. Questionnaires using Likert scales assessed acceptability of habitat functions and sim quality – the degree to which the function was represented in CAD. This enabled the evaluation to be conducted by personnel working remotely. Each function was evaluated individually along with several overarching habitability parameters and vehicle subsystems. The results of this evaluation are discussed, including methodological challenges and rating challenges. Acceptability results are discussed for functions that the participants were able to rate and participant comments for functions that could not be related are also described. Final conclusions are described, including challenges or assumptions needed to make the TH design acceptable, future work needed to refine understanding of the TH, and proposed habitat requirements based on test data.

Transit Habitat↗

Dephasingless laser wakefield acceleration in the bubble regime

Abstract Laser wakefield accelerators (LWFAs) have electric fields that are orders of magnitude larger than those of conventional accelerators, promising an attractive, small-scale alternative for next-generation light sources and lepton colliders. The maximum energy gain in a single-stage LWFA is limited by dephasing, which occurs when the trapped particles outrun the accelerating phase of the wakefield. Here, we demonstrate that a single space–time structured laser pulse can be used for ionization injection and electron acceleration over many dephasing lengths in the bubble regime. Simulations of a dephasingless laser wakefield accelerator driven by a 6.2-J laser pulse show 25 pC of injected charge accelerated over 20 dephasing lengths (1.3 cm) to a maximum energy of 2.1 GeV. The space–time structured laser pulse features an ultrashort, programmable-trajectory focus. Accelerating the focus, reducing the focused spot-size variation, and mitigating unwanted self-focusing stabilize the electron acceleration, which improves beam quality and leads to projected energy gains of 125 GeV in a single, sub-meter stage driven by a 500-J pulse.

43 PARTICLE ACCELERATORS↗

Aerothermodynamic Design of the Mars Science Laboratory Heatshield

Aerothermodynamic design environments are presented for the Mars Science Laboratory entry capsule heatshield. The design conditions are based on Navier-Stokes flowfield simulations on shallow (maximum total heat load) and steep (maximum heat flux, shear stress, and pressure) entry trajectories from a 2009 launch. Boundary layer transition is expected prior to peak heat flux, a first for Mars entry, and the heatshield environments were defined for a fully-turbulent heat pulse. The effects of distributed surface roughness on turbulent heat flux and shear stress peaks are included using empirical correlations. Additional biases and uncertainties are based on computational model comparisons with experimental data and sensitivity studies. The peak design conditions are 197 W/sq cm for heat flux, 471 Pa for shear stress, 0.371 Earth atm for pressure, and 5477 J/sq cm for total heat load. Time-varying conditions at fixed heatshield locations were generated for thermal protection system analysis and flight instrumentation development. Finally, the aerothermodynamic effects of delaying launch until 2011 are previewed.

Edquist, Karl T.↗

Development of a Climatology of Vertically Complete Wind Profiles from Doppler Radar Wind Profiler Systems

This paper describes in detail the QC and splicing methodology for KSC's 50- and 915-MHz DRWP measurements that generates an extensive archive of vertically complete profiles from 0.20-18.45 km. The concurrent POR from each archive extends from April 2000 to December 2009. MSFC NE applies separate but similar QC processes to each of the 50- and 915-MHz DRWP archives. DRWP literature and data examination provide the basis for developing and applying the automated and manual QC processes on both archives. Depending on the month, the QC'ed 50- and 915-MHz DRWP archives retain 52-65% and 16-30% of the possible data, respectively. The 50- and 915-MHz DRWP QC archives retain 84-91% and 85-95%, respectively, of all the available data provided that data exist in the non- QC'ed archives. Next, MSFC NE applies an algorithm to splice concurrent measurements from both DRWP sources. Last, MSFC NE generates a composite profile from the (up to) five available spliced profiles to effectively characterize boundary layer winds and to utilize all possible 915-MHz DRWP measurements at each timestamp. During a given month, roughly 23,000-32,000 complete profiles exist from 0.25-18.45 km from the composite profiles' archive, and approximately 5,000- 27,000 complete profiles exist from an archive utilizing an individual 915-MHz DRWP. One can extract a variety of profile combinations (pairs, triplets, etc.) from this sample for a given application. The sample of vertically complete DRWP wind measurements not only gives launch vehicle customers greater confidence in loads and trajectory assessments versus using balloon output, but also provides flexibility to simulate different DOL situations across applicable altitudes. In addition to increasing sample size and providing more flexibility for DOL simulations in the vehicle design phase, the spliced DRWP database provides any upcoming launch vehicle program with the capability to utilize DRWP profiles on DOL to compute vehicle steering commands, provided the program applies the procedures that this report describes to new DRWP data on DOL. Decker et al. (2015) details how SLS is proposing to use DRWP data and splicing techniques on DOL. Although automation could enhance the current DOL 50-MHz DRWP QC process and could streamline any future DOL 915-MHz DRWP QC and splicing process, the DOL community would still require manual intervention to ensure that the vehicle only uses valid profiles. If a program desires to use high spatial resolution profiles, then the algorithm could randomly add high-frequency components to the DRWP profiles. The spliced DRWP database provides lots of flexibility in how one performs DOL simulations, and the algorithms that this report provides will assist the aerospace and atmospheric communities that are interested in utilizing the DRWP.

Barbre, Robert E., Jr.↗

State predictive information bottleneck

We report the ability to make sense of the massive amounts of high-dimensional data generated from molecular dynamics simulations is heavily dependent on the knowledge of a low-dimensional manifold (parameterized by a reaction coordinate or RC) that typically distinguishes between relevant metastable states, and which captures the relevant slow dynamics of interest. Methods based on machine learning and artificial intelligence have been proposed over the years to deal with learning such low-dimensional manifolds, but they are often criticized for a disconnect from more traditional and physically interpretable approaches. To deal with such concerns, in this work we propose a deep learning based state predictive information bottleneck approach to learn the RC from high-dimensional molecular simulation trajectories. We demonstrate analytically and numerically how the RC learnt in this approach is connected to the committor in chemical physics and can be used to accurately identify transition states. A crucial hyperparameter in this approach is the time delay or how far into the future the algorithm should make predictions about. Through careful comparisons for benchmark systems, we demonstrate that this hyperparameter choice gives useful control over how coarse-grained we want the metastable state classification of the system to be. We thus believe that this work represents a step forward in systematic application of deep learning based ideas to molecular simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Initial Analysis of and Predictive Model Development for Weather Reroute Advisory Use

In response to severe weather conditions, traffic management coordinators (TMCs) reroute air traffic around regions of airspace affected by the severe weather. Presently, acceptable reroutes are specified by issuing reroute advisories. Advisories, based on reroutes published in the National Severe Weather Playbook, consist of several routes describing routing options. These reroutes are selected by TMCs based on their understanding of weather conditions and their previous experience dealing with similar weather conditions. Providing recommendations and analysis of available reroute options could assist the TMCs in making rerouting decisions. While reroute advisories have historically been based on Playbook plays, plays are frequently modified or combined to generate an advisory. Thus, it is natural to present rerouting suggestions to TMCs as modified Playbook plays. The challenge here is to compare advisory subroutes and play subroutes in a meaningful way in order to map advisories back to the Playbook plays upon which they may have been based. In this work, several subroute comparison metrics are presented and discussed. One of these metrics is then used to match advisories to Playbook plays. With the development and implementation of NextGen technologies, there is currently a trend of moving away from the pre-defined Playbook plays (and related advisories) in favor of a more precise specification of trajectories using GPS based navigation tools, such as area navigation (RNAV). Moving towards ubiquitous use of RNAV routes, the objective of rerouting flights will be met with more precisely specified and accurately flown RNAV routes. In the final version of this paper, significant RNAV routes will be identified by finding RNAV routes that are flown frequently. This will require a trajectory comparison technique to compare actually flown flight tracks to specified RNAV routes. The metrics developed here will be used or modified for this task.

weather required reroutes↗

Stability Analysis of Streaks Induced By Optimized Vortex Generators

Numerical computations are performed to investigate the potential for transition control in an axisymmetric boundary layer via fully realizable, streamwise stationary streaks induced by an azimuthally periodic array of surface mounted vortex generators (VGs). Previous work has shown that suitable streaks of this type can significantly reduce the growth of Mack’s second mode instabilities, but large streak amplitudes can make the flow susceptible to previously absent streak instabilities that can become the leading cause of transition. Here, we use the adjoint capabilities of the SU2 flow solver to optimize the VG shape to maximize the reduction in the growth of second-mode disturbances while also preventing the streak amplitudes from reaching large enough values to precipitate an earlier onset of transition via streak instabilities. The geometry and the freestream flow conditions are selected to match a relevant trajectory lo-cation from the HIFiRE-1 flight experiment. Results show that the optimized VGs can increase the mean streak amplitude by 117% with respect to a manually developed baseline design. The stability of this optimized basic state is analyzed via the plane-marching parabolized stability equations, predicting a fully laminar flow over the entire cone, or equivalently, yielding transition delay of 130% versus the 17% for the baseline VGs.

Boundary layer transition↗

Manipulation strategies for massive space payloads

Control for the bracing strategy is being examined. It was concluded earlier that trajectory planning must be improved to best achieve the bracing motion. Very interesting results were achieved which enable the inverse dynamics of flexible arms to be calculated for linearized motion in a more efficient manner than previously published. The desired motion of the end point beginning at t=0 and ending at t=t sub f is used to calculate the required torque at the joint. The solution is separated into a causal function that is zero for t is less than 0 and an accusal function which is zero for t is greater than t sub f. A number of alternative end point trajectories were explored in terms of the peak torque required, the amount of anticipatory action, and other issues. The single link case is the immediate subject and an experimental verification of that case is being performed. Modeling with experimental verification of closed chain dynamics continues. Modeling effort has pointed out inaccuracies that result from the choice of numerical techniques used to incorporate the closed chain constraints when modeling our experimental prototype RALF (Robotic Arm Large and Flexible). Results were compared to TREETOPS, a multi body code. The experimental verification work is suggesting new ways to make comparisons with systems having structural linearity and joint and geometric nonlinearity. The generation of inertial forces was studied with a small arm that will damp the large arm's vibration.

Book, Wayne J.↗

A gravity field comparison and analysis to support Magellan navigation

The orbit characteristics of Magellan are such that the spacecraft is extremely sensitive to gravity perturbations; therefore, the precision orbit determination required to construct the radar images is very dependent on accurate gravitational field modeling. As such, a comparison and analysis is presented of four global harmonic Venusian gravity fields generated by two-way S-band Doppler tracking data from the Pioneer Venus Orbiter in order to support the precise navigational requirements of the Magellan spacecraft. The four fields are: (1) one of degree and order 10 (a 10 x 10) from Mottinger et al. (1985), (2) an 18 x 18 from Bills et al. (1987), (3) a 21 x 21 from McNamee et al. (1990), and (4) a 36 x 36 from Smith et al. (1991). Least-square fits are done with each gravity field on 36 8-rev Magellan X-band data arcs, spaced at 10-deg longitude intervals, completely encircling Venus. Trajectory differences are also performed for two distinct cases, one where the gravity field modeling coincides, and the other where the gravity fields are very different. Although all of the above fields are found to model Venus gravity fairly well, the results show that there is still much room for improvement in Venus gravity modeling.

Kronschnabl, G. R.↗

Traffic Aware Planner for Cockpit-Based Trajectory Optimization

The Traffic Aware Planner (TAP) software application is a cockpit-based advisory tool designed to be hosted on an Electronic Flight Bag and to enable and test the NASA concept of Traffic Aware Strategic Aircrew Requests (TASAR). The TASAR concept provides pilots with optimized route changes (including altitude) that reduce fuel burn and/or flight time, avoid interactions with known traffic, weather and restricted airspace, and may be used by the pilots to request a route and/or altitude change from Air Traffic Control. Developed using an iterative process, TAP's latest improvements include human-machine interface design upgrades and added functionality based on the results of human-in-the-loop simulation experiments and flight trials. Architectural improvements have been implemented to prepare the system for operational-use trials with partner commercial airlines. Future iterations will enhance coordination with airline dispatch and add functionality to improve the acceptability of TAP-generated route-change requests to pilots, dispatchers, and air traffic controllers.

Woods, Sharon E.↗

Effects of Dead‐End Fractures on Non‐Fickian Transport in Three‐Dimensional Discrete Fracture Networks

Abstract Understanding mechanistic causes of non‐Fickian transport in fractured media is important for many hydrogeologic processes and subsurface applications. This study elucidates the effects of dead‐end fractures on non‐Fickian transport in three‐dimensional (3D) fracture networks. Although dead‐end fractures have been identified as low‐velocity regions that could delay solute transport, the direct relation between dead‐end fractures and non‐Fickian transport has been elusive. We systematically generate a large number of 3D discrete fracture networks with different fracture length distributions and fracture densities. We then identify dead‐end fractures using a novel graph‐based method. The effect of dead‐end fractures on solute residence time maximizes at the critical fracture density of the percolation threshold, leading to strong late‐time tailing. As fracture density increases beyond the percolation threshold, the network connectivity increases, and dead‐end fractures diminish. Consequently, the increase in network connectivity leads to a reduction in the degree of late‐time tailing. We also show that dead‐end fractures can inform about main transport paths, such as the mean tortuosity of particle trajectories. This study advances our mechanistic understanding of solute transport in 3D fracture networks.

58 GEOSCIENCES↗

radEq Add-On Module for CFD Solver Loci-CHEM

Loci-CHEM to be applied to flow velocities where surface radiation due to heating from compression and friction becomes significant. The module adds a radiation equilibrium boundary condition to the computational fluid dynamics (CFD) code to produce accurate results. The module expanded the upper limit for accurate CFD solutions of Loci-CHEM from Mach 4 to Mach 10 based on Space Shuttle Orbiter Re-Entry trajectories. Loci-CHEM already has a very promising architecture and performance, but absence of radiation equilibrium boundary condition limited the application of Loci-CHEM to below Mach 4. The immediate advantage of the add-on module is that it allows Loci-CHEM to work with supersonic flows up to Mach 10. This transformed Loci-CHEM from a rocket engine- heritage CFD code with general subsonic and low-supersonic applications, to an aeroheating code with hypersonic applications. The follow-on advantage of the module is that it is a building block for additional add-on modules that will solve for the heating generated at Mach numbers higher than 10.

McCloud, Peter↗