Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time step”

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,729 records · Page 96

A Physics-Based Approach to Urban Air Mobility

High-fidelity Computational Fluid Dynamics (CFD) simulations for multi-rotor vehicles have been carried out. The three-dimensional unsteady Navier-Stokes equations are solved on overset grids employing high order accurate schemes, dual-time stepping, and a hybrid turbulence model using NASA's CFD code Over- flow. The vehicles studied consist of small to medium sized drones, and bigger vehicles for future Urban Air Mobility (UAM) applications. The performances for different configurations and rotor mounting are calculated in hover and in forward flight. Understanding the complex flows and the interactions between rotors and with other elements will help design the future multi-rotor vehicles to be quieter, safer, and more efficient.

Urban↗

A Versatile Simulation Framework for Elastodynamic Modeling of Structural Health Monitoring

Structural health monitoring (SHM) has the capacity to reduce failure by detecting damage during service life, by periodic, automated monitoring. Guided Wave (GW) Ultrasound is a common SHM approach for aerospace structures. Modelling the physics of GW SHM systems provides a route for understanding system dependencies, capabilities and limitations as damage evolves during service life. Such a toolset can strengthen the understanding of the connection between GW SHM results and the true material state. The most useful modelling tools are those that provide versatile solutions with respect to the simulated component geometry and computational grid connectivity. This work details a versatile application programming interface (API) for the elastodynamic finite integration technique for modelling GW SHM of metals. The custom code implementation, EFIT-CompCell, allows for the modelling of diverse geometries by automatically balancing the message passing interface parallelization layout. The user provides the basic parameters of the simulation and the software automatically performs an initial balancing based on anticipated computational loads, and establishes the CPU communication patterns for any geometry. This work describes the programming philosophy and code structure used to create EFIT-CompCell and compares its performance and capacity to simulation tools that are more specialized for specific architectures. Results are presented for a simulation of GW SHM of an aluminum fuselage section being tested by the FAA. The simulation consists of 733M voxels which took approximately 70 hours to complete 25000 time steps using 40 Intel Xeon E5-4650v2 Ivy Bridge processor cores.

Gregory, Elizabeth D.↗

Diurnal and Seasonal Solar Induced Chlorophyll Fluorescence and Photosynthesis in a Boreal Scots Pine Canopy

Solar induced chlorophyll fluorescence has been shown to be increasingly an useful proxy for the estimation of gross primary productivity (GPP), at a range of spatial scales. Here, we explore the seasonality in a continuous time series of canopy solar induced fluorescence (hereafter SiF) and its relation to canopy gross primary production (GPP), canopy light use efficiency (LUE), and direct estimates of leaf level photochemical efficiency in an evergreen canopy. SiF was calculated using infilling in two bands from the incoming and reflected radiance using a pair of Ocean Optics USB2000+ spectrometers operated in a dual field of view mode, sampling at a 30 min time step using custom written automated software, from early spring through until autumn in 2011. The optical system was mounted on a tower of 18 m height adjacent to an eddy covariance system, to observe a boreal forest ecosystem dominated by Scots pine. (Pinus sylvestris) A Walz MONITORING-PAM, multi fluorimeter system, was simultaneously mounted within the canopy adjacent to the footprint sampled by the optical system. Following correction of the SiF data for O2 and structural effects, SiF, SiF yield, LUE, the photochemicsl reflectance index (PRI), and the normalized difference vegetation index (NDVI) exhibited a seasonal pattern that followed GPP sampled by the eddy covariance system. Due to the complexities of solar azimuth and zenith angle (SZA) over the season on the SiF signal, correlations between SiF, SiF yield, GPP, and LUE were assessed on SZA <50° and under strictly clear sky conditions. Correlations found, even under these screened scenarios, resulted around ~r2 = 0.3. The diurnal responses of SiF, SiF yield, PAM estimates of effective quantum yield (ΔF/Fm′), and meteorological parameters demonstrated some agreement over the diurnal cycle. The challenges inherent in SiF retrievals in boreal evergreen ecosystems are discussed.

photosynthetic↗

Computational Study of the Side-By-Side Urban Air Taxi Concept

High-fidelity computational fluid dynamics simulations of NASA’s Side-by-side air taxi concept have been carried out. The three-dimensional unsteady Navier-Stokes equations are solved on overset grids using high-order accurate schemes, dual-time stepping, and a hybrid turbulence model. The flow solver has been loosely coupled with a helicopter comprehensive analysis code in order to get the trimmed flight solution. The vehicle simulated is a six-passenger side-by-side intermeshing rotor helicopter with hybrid propulsion for air taxi operations, also known as urban air mobility applications. This concept vehicle is intended to focus and guide NASA research activities in support of aircraft development for emerging aviation markets, in particular vertical take-off and landing air taxi operations.

Ventura Diaz, Patricia↗

Computational Study of NASA's Quadrotor Urban Air Taxi Concept

High-fidelity computational fluid dynamics simulations have been carried out in order to analyze NASA's quadrotor urban air taxi concept for urban air mobility, also know as on-demand mobility applications. High-order accurate schemes, dual-time stepping, and the delayed detached-eddy simulation model have been employed. The ow solver has been loosely coupled with a rotorcraft comprehensive analysis code. The vehicle simulated is a six-passenger quadrotor for air taxi operations. A study of power reduction as a function of the rear-rotor to front-rotors vertical separation has been performed, for a quad-rotor without the airframe, in cruise flight conditions. Then, the quadrotor without the airframe has been simulated in hover. The airloads and wake geometries are analyzed. To finish the study the complete quadrotor vehicle is presented. NASA's quadrotor air taxi concept is one of the many concepts being developed by NASA in support of aircraft development for vertical take-o and landing air taxi operations.

Ventura Diaz, Patricia↗

Enabling Dynamic Vehicle Analyses With Improved Atmospheric Attenuation Models in Glenn Research Center Communication Analysis Suite

To aid in meeting the NASA objective of returning humans to the Moon, the Glenn Research Center’s Communication Analysis Suite was augmented with two distinct capabilities. The first capability added was the vehicle propagator. This allows the addition of dynamic aircraft and ground vehicles around any celestial body within the solar system during an analysis. This functionality interpolates the position and velocity of the vehicle relative to a celestial body at the time steps analyzed using the type of path and either a series of waypoints or a direction and duration of travel. The implications of this new capability include lunar rovers and/or drones, such as Dragonfly, where the vehicle propagator will analyze the communications architecture. The newly created vehicle propagator is now in use in communications studies for the 2024 lunar missions, simulating the movement of lunar rovers across the Moon’s southern pole. The second capability added was the augmentation of the atmospheric attenuation model. The previous model did not have a uniform low-elevation attenuation model due to the trigonometric approximation for path length and the exponential nature of low-elevation scintillation. User-defined weather parameters were also added to the updated atmospheric attenuation model. The previous model solely used tabular data based upon the season and location of the transmitting antenna. Multiple simulations of the same configuration now return different results based on the differing weather parameters. Cognitive communications analysis efforts can use this second capability to generate neural network training data based on differing weather conditions at utilized ground stations, a critical step in allowing neural networks to learn how weather parameters impact communications performance.

vehicle propagation↗

TPSAS-NF1676L-18428-DND

The Challenges of Unsteady Adjoint-Based Design: - The adjoint approach still provides all of the sensitivities at the same cost as analysis, and the 20x estimate still applies for the expense of an optimization - But every simulation is now an unsteady problem - Where the steady adjoint solver linearized about a single solution (the steady-state), the unsteady adjoint solver must essentially do this at every physical time step

Eric Nielsen↗

TPSAS-NF1676L-10767-DND

Version 2.1ae of the NASA CFD code OVERFLOW has been used to simulate the flow past tandem cylinders. The multiblock grid system consisted of 15 overset meshes with a total grid count of 15.7 million points. The grid extended 3 cylinder diameters in the spanwise direction with periodic boundary conditions at the ends. Riemann conditions were used at the outer boundaries, and a no-slip condition on the cylinder surfaces. The 3rd order HLLC upwind spatial discretization scheme was used with a dual-time stepping algorithm. A modified version of the Delayed Detached Eddy Simulation hybrid RANS/LES turbulence model proposed by Spalart was used in a fully turbulent mode. The code was run on 53 Core 2 Duo nodes using MPI across nodes and OPENMP within each node. Simulation results will be compared with the available experimental data.

David P. Lockard↗

Coupling a Particle Trajectory Capability to the US3D CFD Flow Solver

The purpose of the current research effort is to incorporate particle trajectory equations into the US3D CFD flow solver. The integrated code will allow higher-fidelity modeling of particle-fluid interactions and allow for the simulation of physical effects such as surface heating augmentation due to particle impacts. The key elements underpinning the current work are a robust tracking algorithm that can identify the encompassing mesh elements for individual particles and a point to point MPI frame-work. These developments combined with accurate property models and high-order time-stepping deliver realistic simulations on larger computational grids for a wide range of particle mass loadings. The coupled US3D-particle code will be applied to the Schiarparelli capsule and Earth-based experimental test cases.

A Sahai↗

Investigating the Numerical Stability of Using an Impedance Boundary Condition to Model Broadband Noise Scattering With Acoustic Liners

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for achieving aircraft noise reduction and are characterized by a frequency-dependent impedance value. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic wave scattering by geometric bodies treated with acoustic liners. A Broadband Impedance Model will be discussed in which the liner impedance is specified along a wide range of frequencies. The liner impedance boundary condition will be derived and coupled with a time-domain boundary integral equation to model acoustic scattering by a flat plate consisting of both un-lined and lined surfaces. The scattering solution will be obtained iteratively using both spatial and temporal basis functions and the stability will be demonstrated through eigenvalue analysis. Stability will be assessed for its dependence on time step, spatial discretization, as well as temporal basis function order. Both second- and third-order Lagrange temporal basis functions are considered.

acoustics↗

Investigating the Numerical Stability of Using an Impedance Boundary Condition to Model Broadband Noise Scattering With Acoustic Liners

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for achieving aircraft noise reduction and are characterized by a frequency-dependent impedance value. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic wave scattering by geometric bodies treated with acoustic liners. A Broadband Impedance Model will be discussed in which the liner impedance is specified along a wide range of frequencies. The liner impedance boundary condition will be derived and coupled with a time-domain boundary integral equation to model acoustic scattering by a flat plate consisting of both un-lined and lined surfaces. The scattering solution will be obtained iteratively using both spatial and temporal basis functions and the stability will be demonstrated through eigenvalue analysis. Stability will be assessed for its dependence on time step, spatial discretization, as well as temporal basis function order. Both second- and third-order Lagrange temporal basis functions are considered.

acoustics↗

Scale resolving simulations of the NASA Juncture Flow Model using the LAVA solver

Two distinct simulation methodologies: Delayed Detached Eddy Simulation (DDES) and stress based Wall-modelled Large Eddy Simulation (WMLES) are evaluated using structured overset curvilinear grids for the NASA juncture flow model. While both methodologies are shown to mitigate the primary shortcomings of steady state Reynolds Averaged Navier Stokes(RANS) simulations, several unresolved aspects are identified. Strong sensitivity to RANS-typegrid refinement is observed in the DDES with a substantial deterioration of the solution quality with increasing spatial resolution associated with deficiencies in the shielding function. Mean profiles for attached boundary layers on the fuselage show spurious inflections suggesting modelled stress depletion on finer grids. Lower numerical dissipation in terms of spatial discretization and time step size is seen to improve the solution quality, and the advantage of using RC and QCR2000 corrections in the underlying RANS closure is demonstrated for DDES on RANS-type meshes. Equilibrium Wall-Modelled LES used grids that resolved the tripping dots over the fuselage nose and the wing leading edge, consistent with the experiment setup. These simulations resulted in a cost-competitive approach compared to DDES on RANS-type grids. Although the agreement between WMLES predictions of first and second order single point statistics with experimental measurements is promising, a fundamental shortcoming is noted in terms of an overshoot in streamwise momentum in corner regions of the wing-fuselage juncture. This excess momentum subsequently delays the onset of separation,thereby resulting in an underprediction in length of the separation bubble. Some quantitative sensitivity to numerical discretization is observed; lowering of numerical dissipation shows better agreement with the experiment.

T^3↗

A Heat and Mass Transfer Model of the Orion European Service Module Propulsion Sub-System

A tool for simulating the thermodynamics of the fluid components of the Orion European Service Module propulsion subsystem has been developed. The 1-dimensional,time-stepping method tracks energy and mass movement through the propellant and pressurant tanks during burn and coast phases of a simulated mission. The code is fast and flexible, allowing easy modifications and incorporation into other computational tools. The results from the simulation have been used to support mission planning, development of flight rules and to guide experimental testing of components and systems.

Spacecraft propulsion↗

AN INTRODUCTION TO THE GEONEX LEVEL-1G PRODUCTS: TOP-OF-ATMOSPHERE REFLECTANCE AND BRIGHTNESS TEMPERATURE

This paper introduces the GeoNEX (Geostationary-NASA Earth eXchange) Level-1G products of top-of-atmosphere (TOA) reflectance and brightness temperature. The products use data streams from the latest geostationary (GEO) sensors including the GOES-16/17 ABI and the Himawari-8/9 AHI. The GeoNEX processing pipeline starts by converting digital numbers to physical quantities with the latest radiometric calibration information. It integrates algorithms to automatically detect and remove residual geolocation errors, to estimate the pixel-wise data-acquisition time, and to accurately calculate the solar illumination angles for each pixel in the domain at every time step. The outputs are reprojected to a globally tiled common grid in geographic coordinates designed to facilitate inter-comparisons and/or synergies between the GeoNEX products and existing Earth observation datasets from polar-orbiting satellites. Therefore, the GeoNEX L1G products provide accurate and consistent TOA reflectance and brightness temperature datasets for scientific analyses and downstream product development.

Geostationary satellite, GOES-16, Himawari-8, NASA↗

Adaptive Stress Testing: Using Reinforcement Learning to Find Failures in Safety-Critical Systems

Emerging applications in artificial intelligence, such as driverless cars and autonomous aircraft promise to be more efficient, cheaper to operate, and always available. However, ensuring the safety of these systems remains a major challenge to their certification and adoption. These autonomous systems are expected to routinely make safety-critical decisions where failures can have serious consequences including loss of life and property. Testing and validation techniques aim to identify and diagnose potential failures before the system is deployed. However, finding failure scenarios in autonomous systems can be very challenging due to high-dimensional and continuous state spaces, interaction with large environments over many time steps, and the rarity of failures. This talk presents Adaptive Stress Testing (AST), a simulation-based testing framework for finding the most likely path to a failure event of a safety-critical system. The key idea of AST is that stress testing can be formulated as a Partially Observable Markov Decision Process (POMDP), which enables reinforcement learning techniques to be used for finding failure events. Reinforcement learning algorithms can efficiently explore the search space and have been shown to scale to very large systems. We present applications of AST to find failures in various safety-critical systems including the aircraft collision avoidance systems, autonomous cars, and small unmanned aerial vehicles.

autonomous vehicles↗

Scalability of Cohesive Fatigue Analyses Using Explicit Solvers

A cohesive fatigue law has been integrated into a constitutive material model compatible with an explicit finite element solver. The cohesive fatigue model response is based on engineering approximations of the endurance limit and the Goodman diagram. This approach can predict stress-life diagrams for crack initiation, the Paris law regime, and transient effects of crack initiation and stable tearing. Simplified cyclic loading is utilized so that the applied load(or displacement) corresponds to the peak load of a fatigue cycle. Loads are held constant during fatigue while damage develops with increasing solution increments. An automatically-calculated ratio of fatigue cycles per solution increment controls the rate of damage growth, ensuring that damage growth is modeled with a sufficient minimum number of increments and damage growth advances to a minimum desired extent within the explicit analysis step time. The compatibility with an explicit finite element solver enables the analysis of structures that are computationally intractable for implicit finite element solvers. Scalability studies are conducted for geometrically nonlinear problems involving fiber-reinforced composite structures that exhibit fatigue damage growth of interacting matrix cracks and delaminations

Frank A Leone↗

Multibody Simulation of REMORA CubeSat Docking to and Pushing a Spent Rocket Booster

This paper details the multibody simulation of three phases of critical importance to the feasibility of the REMORA CubeSat space debris mitigation concept: the final approach of the CubeSat to a spent rocket booster; the grappling of the spent rocket booster using a robotic arm; and the pushing of the spent rocket booster to divert its course from another on-orbit asset. The extension of a robotic mobility and manipulation modeling toolkit (M3TK) from multibody dynamics simulation of manipulators and ground vehicles to simulation of orbital robotics is outlined. This includes the identification of the appropriate parameters required to concisely and generically describe thruster loads, thruster mixing, spacecraft control, and spacecraft navigation for the purpose of on-orbit robotics simulation. A high-level spacecraft navigator commands maneuvers to target spacecraft states. A PID spacecraft controller takes the target states and calculates desired forces and torques. A thruster mixer solves a quadratic program to determine the optimal thruster firing times for the propulsion system. Pulse width modulated actuation of eight canted cold gas thrusters is used in the simulated approach to a rocket nozzle from a distance of 200 m. A five degree of freedom robotic arm is controlled to position a pair of pincers to grasp the rocket nozzle. Contact dynamics are used to accurately simulate the grasping of the rocket nozzle by the pincers. A similar, but separate simulation is performed to assess the ability of the REMORA CubeSat to push the large spent rocket booster. This diversion maneuver makes use of an additional, larger thruster, and pushes the rocket booster in excess of 400 m. Appropriate motor control gains on the robotic arm are found to be higher during the pushing phase than those which are appropriate during free motion; this increase promotes rigidity of the arm and allows it to properly direct the pushing force. Challenges encountered in time step selection for numerical stability of the simulation are also discussed.

Mukherjee, Rudranarayan↗