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,135 records · Page 63

Target Tracking with Distributed Sensing and Optimal Data Migration

The paper presents an Extended Kalman Filter based framework for airborne target tracking using adaptive information fusion from multi-modal multi-rate distributed sensors network. First, the tracking algorithm execution location is determined using an optimal data migration strategy, which also computes the associated delays for each sensor data to arrive at the computing location. Next, the fast (zero-delay) sensors information is dynamically fused in the filter correction procedure at the arrival instance of each valid sensor reading. Finally, the target estimation is updated based on the valid slow (delayed) data, which are grouped according to the delay-time steps before application of the Larsen's method. This approach is applied to the synthetic sensor data generated by means of the ground based radar and camera models for the simulated target flight in Reflection simulation environment.

Distributed sensing↗

Target Tracking with Distributed Sensing and Optimal Data Migration

The paper presents an Extended Kalman Filter based framework for airborne target tracking using adaptive information fusion from multi-modal multi-rate distributed sensors network. First, the tracking algorithm execution location is determined using an optimal data migration strategy, which also computes the associated delays for each sensor data to arrive at the computing location. Next, the fast (zero-delay) sensors information is dynamically fused in the filter correction procedure at the arrival instance of each valid sensor reading. Finally, the target estimation is updated based on the valid slow (delayed) data, which are grouped according to the delay-time steps before application of the Larsen's method. This approach is applied to the synthetic sensor data generated by means of the ground based radar and camera models for the simulated target flight in Reflection simulation environment.

Distributed sensing↗

High-Fidelity Simulations of a Tiltwing Vehicle for Urban Air Mobility

NASA is playing an important role in the development of the urban air mobility ecosystem by identifying key research areas and establishing a fleet of conceptual reference vehicles for urban air mobility applications. This paper delves deeper into the study of NASA’s six-passenger tiltwing air taxi aerodynamics using high-fidelity computational fluid dynamics simulations. NASA’s OVERFLOW flow solver has been loosely-coupled with the rotorcraft comprehensive analysis code CAMRAD II. Dual-time stepping, high-order accurate schemes, and the delayed detached-eddy simulation are employed. The tiltwing configuration is analyzed in airplane and helicopter mode, with flight conditions representative of an air taxi mission. Airloads, wake geometry, and performance are studied to identify the unique flow features of this configuration. NASA’s tiltwing vehicle is one of the conceptual designs intended to focus and guide NASA’s research activities in support of aircraft development for vertical take-off and landing air taxi operations.

RVLT↗

Fan Noise Predictions of the NASA Source Diagnostic Test Using Unsteady Simulations with LAVA Part I: Near-Field Aerodynamics and Turbulence

A sliding mesh technique within the Launch, Ascent, and Vehicle Aerodynamics (LAVA) computational framework is validated using the experimental dataset collected as part of the NASA Source Diagnostic Test (SDT) campaign. Two modeling approaches are explored: the unsteady Reynolds-Averaged Navier Stokes (URANS) with Spalart-Allmaras (SA) turbulence model closure, and a hybrid Reynolds-Averaged Navier Stokes/Large Eddy Simulation (RANS/LES) paradigm employing a Zonal Detached Eddy Simulation (ZDES) closure with enhanced shielding protection. Fan stage performance metrics, aerodynamic quantities and turbulent flow structures are analyzed in this work. Initial studies focusing on grid and time-step sensitivity are presented. Sensitivity to different variants of the SA turbulence model is analyzed, supporting the use of the baseline SA model in the production runs. Two conditions are analyzed in detail using URANS and hybrid RANS/LES (HRLES). Mean flow quantities are well-captured by both methods in the low-speed (approach) regime. While URANS misses all the upstream-propagating noise in the inlet due to the rotor-locked tones being evanescent in nature at subsonic fan tip speeds, HRLES captures this broadband component in its pressure field. At the high-speed (sideline) condition, URANS shows better agreement with the SDT data than HRLES in the interstage flow-field. In this regime, URANS captures the tonal content propagating through the inlet, since the tones are now cut-on. Both methods are suitable to capture fan stage performance metrics and mean flow quantities, but only HRLES is able to resolve the fine turbulent structures responsible for broadband noise. The results support the use of the sliding mesh technique implemented in this work for future turbomachinery applications within the LAVA solver framework.

AATT↗

Fan Noise Predictions of the NASA Source Diagnostic Test Using Unsteady Simulations with LAVA Part I: Near-Field Aerodynamics and Turbulence

A sliding mesh technique within the Launch, Ascent, and Vehicle Aerodynamics (LAVA) computational framework is validated using the experimental dataset collected as part of the NASA Source Diagnostic Test (SDT) campaign. Two modeling approaches are explored: the unsteady Reynolds-Averaged Navier Stokes (URANS) with Spalart-Allmaras (SA) turbulence model closure, and a hybrid Reynolds-Averaged Navier Stokes/Large Eddy Simulation (RANS/LES) paradigm employing a Zonal Detached Eddy Simulation (ZDES) closure with enhanced shielding protection. Fan stage performance metrics, aerodynamic quantities and turbulent flow structures are analyzed in this work. Initial studies focusing on grid and time-step sensitivity are presented. Sensitivity to different variants of the SA turbulence model is analyzed, supporting the use of the baseline SA model in the production runs. Two conditions are analyzed in detail using URANS and hybrid RANS/LES (HRLES). Mean flow quantities are well-captured by both methods in the low-speed (approach) regime. While URANS misses all the upstream-propagating noise in the inlet due to the rotor-locked tones being evanescent in nature at subsonic fan tip speeds, HRLES captures this broadband component in its pressure field. At the high-speed (sideline) condition, URANS shows better agreement with the SDT data than HRLES in the interstage flow-field. In this regime, URANS captures the tonal content propagating through the inlet, since the tones are now cut-on. Both methods are suitable to capture fan stage performance metrics and mean flow quantities, but only HRLES is able to resolve the fine turbulent structures responsible for broadband noise. The results support the use of the sliding mesh technique implemented in this work for future turbomachinery applications within the LAVA solver framework.

AATT↗

Model Predictive Control in the Three-body Problem Using Invariant Funnels As Terminal Sets

This paper describes a method for augmenting Model Predictive Control techniques using invariant funnels computed in the Circular Restricted Three-Body Problem. We use ellipsoids that roughly approximate the boundary of the invariant funnel as convex terminal sets for a short look-ahead optimization problem at each time step. We apply this method to a hypothetical low-thrust mission to land on Jupiter’s mooon Europa and show that including the invariant funnels as terminal sets reduces the amount of control effort required by almost an order of magnitude.

Close, Sigrid↗

Multi-Fidelity Computational Analysis of a Quiet Single-Main Rotor Helicopter for Air Taxi Operations

A multi-fidelity computational study has been conducted on NASA’s quiet single-main rotor helicopter concept for urban air mobility using computational fluid dynamics and rotorcraft comprehensive analysis tools. Two approaches have been followed in order to simulate the rotor flows: (1) a rotor disk model, which models the rotors as source terms, and (2) a high-fidelity approach using high-order accurate schemes and dual-time stepping, which simulates the rotor with its individual rotating blade grids. The delayed detached-eddy simulation model has been employed. A loose-coupling approach between the flow solver and a rotorcraft comprehensive code is utilized to include vehicle trim and motion variables. The quiet single-main rotor helicopter is one of the conceptual designs intended to focus and guide NASA’s research activities in support of aircraft development for vertical take-off and landing air taxi operations.

RVLT↗

Inferred Sea Level Prediction in the NASA GMAO Seasonal Forecasting System

Reliable predictions of sea level anomalies on seasonal timescales with lead times of 1 to 9 months may have relevance to stakeholders – for example, in the advance deployment of resources for coastal flood mitigation. Routine prediction and analysis may also highlight physical processes associated with sea level change and modeling capabilities on seasonal and other timescales. These forecasts may represent interannual changes in the seasonal slope of the ocean surface, teleconnection effects such as the El Niño/Southern Oscillation phenomenon, and variations in seasonal hydrology including precipitation and coastal runoff. Coupled atmosphere/ocean models are routinely used in the seasonal prediction of temperature anomalies, precipitation anomalies, sea ice cover, and climate indices such as the Niño3.4 predictions under the North American Multi-Model Ensemble (NMME) protocol. Within the limits of their configuration, these complex Earth-system models have a potential for depicting regional changes in oceanic column properties, including the sea surface height. Seasonal prediction models generally have no representation of long-term mass contributions from melting land ice, or changes in vertical land motion; their output may be more specifically characterized as predictions of the ocean dynamic sea level. In practice however, the sea surface height prognostic variable is substantially compromised by the forecast model response to initial conditions. Imbalances between the initial, observed hydrologic cycle and the forecast model state produce abrupt adjustments in the model sea surface height. As a result, most seasonal prediction systems employ a constraint on the globally-averaged sea surface height that is applied at each time step. This essentially renders the prognostic sea surface height variable as unserviceable. Several approaches have previously been used to retrieve sea level information from seasonal forecasts beyond the use of the sea surface height variable. Here, we extend a method of relating other prognostic values, including ocean circulation and climate indices, to observed sea level variations. We use the merged altimetry record of the NASA MEaSUREs Gridded Sea Surface Height Anomalies data set and monthly revised local reference gauge observations from the National Oceanography Centre Permanent Service for Mean Sea Level (PSMSL) to evaluate derived prognostic variables from the NASA Global Modeling and Assimilation Office subseasonal-to-seasonal system version 2.1 (GMAO S2S v2.1). We focus on results for the midlatitudes with particular emphasis on US gauge locations. As shown in previous studies, prognostic ENSO-related indices in boreal winter are well correlated with gauge observations for the US west coast, but also for other locations in the southeastern US. Other forecast climate indices such as the North Atlantic Oscillation have relations to sea level that are limited both seasonally and spatially. As expected, surface atmospheric pressure (e.g., inverse barometer effect) is found to be particularly well correlated with observed sea level. We provide a characterization of forecast skill for seasonal sea level with this method.

Richard I Cullather↗

LGM Paleoclimate Constraints Inform Cloud Parameterizations and Equilibrium Climate Sensitivity in CESM2

The Community Earth System Model version 2 (CESM2) simulates a high equilibrium climate sensitivity (ECS > 5°C) and a Last Glacial Maximum (LGM) that is substantially colder than proxy temperatures. In this study, we examine the role of cloud parameterizations in simulating the LGM cooling in CESM2. Through substituting different versions of cloud schemes in the atmosphere model, we attribute the excessive LGM cooling to the new CESM2 schemes of cloud microphysics and ice nucleation. Further exploration suggests that removing an inappropriate limiter on cloud ice number (NoNimax) and decreasing the time-step size (substepping) in cloud microphysics largely eliminate the excessive LGM cooling. NoNimax produces a more physically consistent treatment of mixed-phase clouds, which leads to an increase in cloud ice content and a weaker shortwave cloud feedback over mid-to-high latitudes and the Southern Hemisphere subtropics. Microphysical substepping further weakens the shortwave cloud feedback. Based on NoNimax and microphysical substepping, we have developed a paleoclimate-calibrated CESM2 (PaleoCalibr), which simulates well the observed twentieth century warming and spatial characteristics of key cloud and climate variables. PaleoCalibr has a lower ECS (∼4°C) and a 20% weaker aerosol-cloud interaction than CESM2. PaleoCalibr represents a physically more consistent treatment of cloud microphysics than CESM2 and is a valuable tool in climate change studies, especially when a large climate forcing is involved. Our study highlights the unique value of paleoclimate constraints in informing the cloud parameterizations and ultimately the future climate projection.

equilibrium climate sensitivity↗

LS-DYNA User-Defined Internal Ballistic Modeling

- LS-DYNA Explicit Finite Element Analysis software and CADPROG Internal Ballistic Analysis code are integrated as a single modeling tool to analyze pyro-mechanical devices. - Combining these two codes allows the ballistic calculations to be performed without modeling the propellant in an FEA mesh. At each LS-DYNA time step, CADPROG is analytically run using the kinematic data (stroke, velocity, chamber volume) fed back from the LS-DYNA simulation. Calculated pressure is then applied to a piston or any applicable surface mesh in LS-DYNA model interactively. -All external loads such as friction, damping loads, assisting/resisting loads, locking/unlocking loads are modeled in LS-DYNA and coupled with pressure loads from CADPROG to create a fully defined equation of motion -Both codes are written in Fortran 77. CADPROG is inserted in LS-DYNA source code as a user-defined subroutine and re-compiled to generate a custom executable.

pyrotechnics↗

Tail-Cone Thruster Fan Rotor Flutter Analysis

The aeroelastic stability of a tail cone thruster fan was analyzed for three structural modes and various nodal diameter patterns. Computational fluid dynamics modeling with blade vibrations was used in the analyses performed at design (peak-efficiency) and two near-stall conditions at the design rotational speed. Aerodynamic damping was used to assess flutter since this integrally bladed rotor design will have minimal structural damping. Aerodynamic damping was found to be low at the design condition. For the first mode, aerodynamic damping was seen to decrease towards stall, vanishing for a zero nodal diameter pattern. Further decrease in mass flow led to negative aerodynamic damping (flutter) for the first and second nodal diameters in the near stall condition. The flowfield at this lowest mass flow rate exhibited some additional unsteadiness at higher frequency. The sensitivity of results to time step and amplitudes of vibration was examined and found to be very satisfactory. The largest uncertainty in the results comes from the observed sensitivity of the aerodynamic damping to details of the variation of grid deformations in the blade passage mesh from the blade passage to the boundaries. For the parametric variations analyzed, the calculated aerodynamic damping is lower if the grid deformations near the blade surface are higher. Further work is needed to understand the cause of this sensitivity and find possible solutions to reduce the effects on flutter predictions.

Fan↗

CFD Modeling & Validation of Single Bubble Growth During Cavitation in Zero Gravity

Cryogenic Fluid Management (CFM) in zero gravity represents a vital area of research as NASA pursues future long-duration space missions. Success of these missions relies heavily on the ability to store and transfer the cryogenic liquid fuels efficiently and safely. In the Zero Boil-Off Tank (ZBOT-1) Experiment in 2018, unexpected and intense phase change occurred at the screened Liquid Acquisition Device (LAD) during tank pressure control that was later attributed to cavitation. The cavitation ultimately resulted in pump failure. Since cryogen cavitation during storage and transfer of propellants in microgravity has been a source of serious concern for the aerospace community, a series of drop tower experiments with liquid methane were also performed at the University of Bremen to carefully study cavitation from a single bubble nucleus in an engineered micron-sized wall cavity. As the test cell is depressurized during the 10-second drop, cavitation causes bubble growth from the nucleus. Here, in a collaboration between the Case Western Reserve University and University of Bremen, we present a two-phase Computational Fluid Dynamics (CFD) model of bubble growth from a single wall cavity during depressurization in microgravity. Both 2D axisymmetric and 3D parametric simulation case studies are performed to show the evolution of the bubble growth. The model uses a Volume of Fluid (VOF) approach with and without adaptive mesh refinement and adaptive time stepping to capture the growth dynamics using an explicit formulation with algebraic reconstruction of the interface. CFD results show good agreement with corresponding experiment. Based on the velocity, temperature, and volume fraction fields generated by the validation case studies, detailed discussion and analysis of fluid flow and heat transfer during bubble growth process are presented to delineate the intricacies of the cavitation process. Strengths and shortcomings of the present numerical approach in capturing this phenomenon are also highlighted.

Two Phase Flow↗

Wake Disturbance-Based Separation Criteria Study for Urban Air Mobility Aircraft

This study continues previous research evaluating the aircraft flight dynamics in the presence of wake disturbances. A model for fixed-wing-generated wake effects is input to the closed-loop system of a Revolutionary Vertical Lift Technology (RVLT) Lift+Cruise aircraft operating in the fixed-wing Cruise configuration. The wake effects are implemented as a time-varying nonlinear aerodynamic model that affect the stability, dynamic, and control derivatives of the aircraft at every time step. The final product of this study is a reduced separation minimum based on further simulation testing with the nonlinear RVLT Lift+Cruise aircraft model, as well as additional wake inputs from rotorcraft and static urban canyons.

Urban Air Mobility↗

Turbulence Model Validation through Joint Experimental / Computational Studies of Separated Flow Over a Three-Dimensional Tapered Bump: Part II - Computational Investigation

This report serves as documentation of activities related to Computational Fluid Dynamics (CFD) analysis applied to the Smooth Body Separation Experiment designed to produce validation-quality experimental data from wind-tunnel tests performed at the University of Notre Dame. All CFD analysis utilized the Boeing CFD (BCFD) flow solver. BCFD was used to analyze flow over the Boeing Speed Bump with second-order accuracy in space and time. Inviscid fluxes were computed using the HLLE scheme for steady-state simulations, and the Low-Dissipation Roe scheme for time-accurate simulations. The working fluid was modeled as calorically perfect air using Sutherland's law for viscosity. Thermal diffusivity was determined by assuming a constant molecular Prandtl number of 0.72, and a turbulent Prandtl number of 0.9. Reynolds Averaged Navier-Stokes (RANS) simulations were conducted using the one-equation Spalart-Allmaras turbulence model with rotation correction3 and the quadratic constitutive relation (SA-RC-QCR), as well as Menter's two-equation, Shear-Stress Transport (SST) model. Time-accurate computations were conducted using Delayed Detached eddy Simulation (DDES) with the one-equation RANS model, and Large Eddy Stress balance (LESb) with the two-equation RANS model. Implicit dual-time stepping was used to achieve second-order temporal accuracy. The experimental investigation conducted during the course of this study is described in detail in an accompanying report.

Flow separation↗

Insights into Long-Term Global Precipitation from IMERG and GPCP

Multiple satellite-based datasets provide estimates of the long-term record of global precipitation. Each necessarily includes both the real atmospheric behavior and a collection of artifacts driven by the input data sources and design choices in the retrievals and dataset construction. For this presentation, the Integrated Multi-satellitE Retrievals for GPM (IMERG) products from the Global Precipitation Measurement (GPM) mission’s U.S. Science Team are designed as a modern High-Resolution Precipitation Product (HRPP), while the Global Precipitation Climatology Project’s (GPCP) products are designed as a modern Climate Data Record (CDR). Although they share some common input data sources, the computational process for each is rather different. Comparing the precipitation estimates from each algorithm gives us insights into likely artifacts and the natural variations that are common to both. The goal of an HRPP is to give the best estimate of precipitation at each time step in the dataset, generally meaning that “all possible available” data are used. In contrast, the CDR is intended to provide a precipitation record that has relatively homogeneous statistics, necessary for climate analysis. Each analysis strives for both goals, of course, but this means that the HRPP uses the disparate collection of satellites whose statistics don’t quite match the ideal record, while the CDR is computed from fewer satellites that provide a relatively homogeneous set of precipitation statistics. This presentation employs the current versions of the IMERG and GPCP products, V07B and V3.2, respectively. A variety of analyses have demonstrated that IMERG has more artifacts, as expected, which tend to result from changes in the source of intercalibration for the various input satellites, namely the Combined Radar-Radiometer Algorithm using TMI and PR during the TRMM era and using GMI and DPR-Ku during the GPM era. As well, a shift occurs when the altitudes of the TRMM and GPM satellites change as a result of orbit boosts. For the most part the mean precipitation shows good continuity across these boundaries, but there are noticeable changes in the respective histograms. The latter has taken on more importance in recent years due to increased scrutiny on extremes which implicitly focus on the upper end of the precipitation histograms. The GPCP product is more homogeneous, with only one major change in calibrator, transitioning from the SSMI series of satellites to the SSMIS series of satellites in 2009. The GPCP analysis does a good job of minimizing artifacts in the means at this boundary, while the histograms show a smaller, but still noticeable shift in the histograms, a result that is similar to data boundaries in IMERG. The GPCP Daily, which is a month-by-month rescaling of the daily accumulated IMERG Final product, shows more consistency than the equivalent daily IMERG, but inherits the histogram shifts from IMERG. Given this overlay of artifacts, there is enough consistency between the products to illustrate some important long-term variations, including interannual variations, not all of which are easily attributable to ENSO events, and trends that are large regionally, but comparatively small when averaged across the globe.

precipitation↗

Image Moment-Based Extended Object Tracking for Complex Motions

A novel image moment-based model for shape estimation and tracking of an extended target moving with a complex trajectory is presented. The proposed extended object tracking algorithm is based on multiple noisy measurement points sampled from the target at each time step. The shape of the object, approximated by an ellipse, is estimated using a combination of image moments. Dynamic models of image moments for constant velocity and coordinated turn motions are mathematically derived. An unscented Kalman filter - interacting multiple model (UKF-IMM) method is used to track the object and estimate its shape. A likelihood function based on average log-likelihood is derived for the IMM filter. Simulation results of the proposed UKF-IMM algorithm with the image momentbased models are presented that show the estimation of the shape of the object moving in a complex trajectory. The intersection over union (IoU), and the root mean square errors (RMSEs) of the position and velocity of the centroid of the ellipse are used as metrics. The comparison results of the proposed algorithm with a benchmark algorithm from literature based on the IoU and RMSE metrics are presented.

Extended object tracking↗

GlennICE Simulation of 24, 28 and 36 Inch Diameter eVTOL Propellers in Forward Flight

Flow solutions and water collection results are presented for three different propeller sizes. Flow solutions were generated using the NASA FUN3D software, and water collection was simulated using NASA GlennICE software. The geometries considered represent experimental articles that have been tested in the NASA Icing Research Tunnel. Each configuration has four blades. For the experiment, the same spinner centerbody is used for each propeller size. The different sized blades are similar for radii larger than about 40% of total span. Operating conditions were chosen to match operating points where experimental data was collected. Results are presented for both single bin simulations, as well as 7-bin Langmuir D distributions. For 80 microns MVD the Langmuir D distribution showed small differences compared to the single bin. For 15 microns MVD the Langmuir D distribution had a larger impact, especially on the spinner. In the experiment, the centerbody was unchanged as propeller size was changed. In this paper, additional simulations were completed with the 36 inch case perfectly scaled down to 24 and 28 inch. The perfect geometric scaling allows for an assessment of the effect of constant centerbody size on the results. While slight differences in pressure distribution on the propeller are visible between scaled and non-scaled centerbody, the effect seems quite minimal regarding water collection. All simulations were run very cold to ensure rime ice, since handling runback water in the non-inertial frame is beyond the scope of the present work. Ice shapes were generated using a single time step growth. The total mass on the blade and spinner was generally underpredicted for 15 microns MVD, and overpredicted for 80 microns MVD. When comparing ice shapes at operating conditions the numerical results produce very similar ice shapes.

Aircraft Icing↗

GlennICE Simulation of 24, 28 and 36 Inch Diameter eVTOL Propellers in Forward Flight

Flow solutions and water collection results are presented for three different propeller sizes. Flow solutions were generated using the NASA FUN3D software, and water collection was simulated using NASA GlennICE software. The geometries considered represent experimental articles that have been tested in the NASA Icing Research Tunnel. Each configuration has four blades. For the experiment, the same spinner centerbody is used for each propeller size. The different sized blades are similar for radii larger than about 40% of total span. Operating conditions were chosen to match operating points where experimental data was collected. Results are presented for both single bin simulations, as well as 7-bin Langmuir D distributions. For 80 microns MVD the Langmuir D distribution showed small differences compared to the single bin. For 15 microns MVD the Langmuir D distribution had a larger impact, especially on the spinner. In the experiment, the centerbody was unchanged as propeller size was changed. In this paper, additional simulations were completed with the 36 inch case perfectly scaled down to 24 and 28 inch. The perfect geometric scaling allows for an assessment of the effect of constant centerbody size on the results. While slight differences in pressure distribution on the propeller are visible between scaled and non-scaled centerbody, the effect seems quite minimal regarding water collection. All simulations were run very cold to ensure rime ice, since handling runback water in the non-inertial frame is beyond the scope of the present work. Ice shapes were generated using a single time step growth. The total mass on the blade and spinner was generally underpredicted for 15 microns MVD, and overpredicted for 80 microns MVD. When comparing ice shapes at operating conditions the numerical results produce very similar ice shapes.

Aircraft Icing↗