Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RANS”

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,009 records · Page 56

CFD Predictions of Soot & CO Emissions Generated by a Partially-Fueled 9-Element Lean-Direct Injection Combustor

A study was undertaken to investigate the CO & soot emissions generated by a partially-fueled 9- element LDI (Lean-Direct Injection) combustor configuration operating in the idle range of jet engine conditions. In order to perform the CFD analysis, several existing soot/chemistry models were implemented into the OpenNCC (Open National Combustion Code). The calculations were based on a Reynolds-Averaged Navier Stokes (RANS) simulation with standard k-epsilon turbulence model, a 62- species jet-a/air chemistry, a 2-equation soot model, & a Lagrangian spray solver. A separate transport equation was solved for all individual species involved in jet-a/air combustion. In the test LDI configuration we examined, only five of the nine injectors were fueled with the major pilot injector operating at an equivalence ratio of near one and the other four main injectors operating at an equivalence ratio near 0.55. The calculations helped to identify several reasons behind the soot & CO formation in different regions of the combustor. The predicted results were compared with the reported experimental data on soot mass concentration (SMC) & emissions index of CO (EICO). The experimental results showed that an increase in either T3 and/or F/A ratio lead to a reduction in both EICO & SMC. The predicted results were found to be in reasonable agreement. However, the predicted EICO differed substantially in one test condition associated with higher F/A ratio.

Raju, M. S.↗

Multigrid Preconditioning for a Space-Time Spectral-Element Discontinuous-Galerkin Solver

In this work we examine a multigrid preconditioning approach in the context of a high- order tensor-product discontinuous-Galerkin spectral-element solver. We couple multigrid ideas together with memory lean and efficient tensor-product preconditioned matrix-free smoothers. Block ILU(0)-preconditioned GMRES smoothers are employed on the coarsest spaces. The performance is evaluated on nonlinear problems arising from unsteady scale- resolving solutions of the Navier-Stokes equations: separated low-Mach unsteady ow over an airfoil from laminar to turbulent ow. A reduction in the number of ne space iterations is observed, which proves the efficiency of the approach in terms of preconditioning the linear systems, however this gain was not reflected in the CPU time. Finally, the preconditioner is successfully applied to problems characterized by stiff source terms such as the set of RANS equations, where the simple tensor product preconditioner fails. Theoretical justification about the findings is reported and future work is outlined.

Franciolini, Matteo↗

Evaluating Combinations of Sentinel-2 Data and Machine-Learning Algorithms for Mangrove Mapping in West Africa

Creating a national baseline for natural resources, such as mangrove forests, and monitoring them regularly often requires a consistent and robust methodology. With freely available satellite data archives and cloud computing resources, it is now more accessible to conduct such large-scale monitoring and assessment. Yet, few studies examine the reproducibility of such mangrove monitoring frameworks, especially in terms of generating consistent spatial extent. Our objective was to evaluate a combination of image processing approaches to classify mangrove forests along the coast of Senegal and The Gambia. We used freely available global satellite data (Sentinel-2), and cloud computing platform (Google Earth Engine) to run two machine learning algorithms, random forest (RF), and classification and regression trees (CART). We calibrated and validated the algorithms using 800 reference points collected using high-resolution images. We further re-ran 10 iterations for each algorithm, utilizing unique subsets of the initial training data. While all iterations resulted in thematic mangrove maps with over 90% accuracy, the mangrove extent ranges between 827-2807 km2 for Senegal and 245-1271 km2 for The Gambia with one outlier for each country. We further report "Places of Agreement" (PoA) to identify areas where all iterations for both methods agree (506.6 km2 and 129.6 km2 for Senegal and The Gambia, respectively), thus have a high confidence in predicting mangrove extent. While we acknowledge the time- and cost-effectiveness of such methods for the landscape managers, we recommend utilizing them with utmost caution, as well as post-classification on-the-ground checks, especially for decision making.

Mondal, Pinki↗

Dynamic Smagorinsky Modeled Large-Eddy Simulations of Turbulence Using Tetrahedral Meshes

Eddy-resolving numerical computations of turbulent flows are emerging as viable alternatives to Reynolds Averaged Navier-Stokes (RANS) calculations for flows with an intrinsically steady mean state due to the advances in large-scale parallel computing. In these computations, medium to large turbulent eddies are resolved by the numerics while the smaller or subgrid scales are either modeled or taken care of by the inherent numerical dissipation. To advance the state of the art of unstructured-mesh turbulence simulation capabilities, large eddy simulations (LES) using the dynamic Smagorinsky model (DSM) on tetrahedral meshes are carried out with the space-time conservation element, solution element (CESE) method. In contrast to what has been reported in the literature, the present implementation of dynamic models allows for active backscattering without any ad-hoc limiting of the eddy viscosity calculated from the subgrid-scale model. For the benchmark problems involving compressible isotropic turbulence decay as well as the shock/turbulent boundary layer interaction benchmark problems, no numerical instability associated with kinetic energy growth is observed and the volume percentage of the backscattering portion accounts for about 38-40% of the simulation domain. A slip-wall model in conjunction with the implemented DSM is used to simulate a relatively high Reynolds number Mach 2.85 turbulent boundary layer over a 30° ramp with several tetrahedral meshes and a wall-normal spacing of either Δγ& = 10 or Δγ& = 20. The computed mean wall pressure distribution, separation region size, mean velocity profiles, and Reynolds stress agree reasonably well with experimental data.

CESE method↗

Kestrel Results at Liftoff Conditions for a Space Launch System Configuration in Proximity to the Launch Tower

Aerodynamic data books for Space Launch System vehicles require databases for the integrated forces and moments and section loads during liftoff and transition to the ascent phase of flight. While the force and moment database can be generated from wind tunnel results, computational analyses are necessary to provide the extensive surface information required to generate proper lineloads. Of the two flight regimes, the liftoff problem is the more costly and complex situation to simulate, as it requires modeling of the vehicle in proximity to the launch tower. The effects of massive separation on the leeward pressure fields of both the tower and vehicle are not well captured with RANS methods, necessitating the use of more advanced methods, such as Delayed Detached Eddy Simulation, in conjunction with computational grids sufficiently refined to resolve the wakes. Details on the computational setup for employing the Kestrel flow solver to address the liftoff problem are presented. The methodology involves the use of independent unstructured near-body grids for the vehicle and the tower, overset by a solution adaptive Cartesian off-body grid. Results from the simulations are compared to experimental results from a test in the NASA Langley Research Center 14- by 22-Foot Subsonic Tunnel.

Computational fluid dynamics↗

Comparison of Boundary Layer Similarity Transformations for High Mach Number Flows

Reynolds averaged Navier-Stokes (RANS) simulations of the turbulent boundary layer flow on a flat plate are performed across a Mach number range of 0.3 up to 16 using the computational fluid dynamics (CFD) software, VULCAN-CFD. The simulation results are used to evaluate the ability of the Van Driest transformation and the transformation of Trettel and Larsson to collapse boundary layer velocity profiles under various flow conditions across the range of Mach numbers. The flat platewall boundary conditions and leading edge geometry are varied from adiabatic to isothermal and sharp to blunted to reveal the physics of how these effects impact the performance of the chosen flow transformations. Results indicate that the transformation of Trettel and Larsson produces a better collapse of velocity profile data than the Van Driest transformation regardless of Mach number, wall boundary condition, or leading edge geometry. In addition, the velocity gradients obtained from the transformation of Trettel and Larsson match the normalized untransformed velocity gradients much more closely than those obtained using the Van Driest transformation.

DiGregorio, Nicholas J.↗

Wall-modeled LES of the NASA Juncture Flow Experiment

We perform wall-modeled large eddy simulation (WMLES) of the NASA Juncture Flow experiment (Kegerise & Neuhart, NASA/TM–2019–20286) using an equilibrium wall model and unstructured finite volume solver, to assess its predictive ability for complex geometries. The flow condition simulated corresponds to 5 degrees angle of attack with a freestream Mach number of 0.189, and a Reynolds number based on the mean crank chord of 2.4 million. All simulations neglect the effects of wind tunnel walls, sting and mast present in the experiment. To assess the performance of WMLES on simpler flows with very coarse grids, we simulate flow in a turbulent channel at a friction Reynolds number, Reτ≈2000, and flow through a square duct at anReτ≈1000. The duct flow results using the equilibrium wall model indicate that the skin friction is inaccurate near the corner for grid topologies and resolutions typically used in WMLES, having implications for more complicated flows such as the juncture flow. For the juncture flow simulation, two different approaches are investigated. The first approach consists of a truncated-domain simulation wherein the inflow plane of the computational domain is placed at≈0.2chord length of the wing, and we prescribe the mean flow from a separate Reynolds-averaged Navier-Stokes (RANS) solution along with synthetic turbulence to initiate realistic unsteadiness in the domain. The second approach involves simulating the entire geometry with trips to trigger transition to turbulence. The truncated- and full-domain simulations contain about 62 and 90 million cells, respectively; with 8-10 points per boundary layer thickness in the wall-parallel directions, 14-20 points in the wall-normal direction and a near-wall viscous spacing (∆n+1) of≈100. Preliminary results are encouraging overall in terms of the prediction of wall pressure, wall skin friction, velocity and stresses; but indicate that further work is required to improve the predictions in the separation bubble and wing-fuselage corner regions. The grid resolution used in this study is still fairly coarse, and the results should be interpreted as work in progress.

Prahladh S Iyer↗

Verification of Anisotropic Mesh Adaptation for Turbulent Simulations over ONERA M6 Wing

Unstructured anisotropic mesh adaptation is known to be an efficient way to control discretization errors in Computational Fluid Dynamics (CFD) simulations. Method verification is required to provide the confidence for routine use in production analysis. The current work aims at verification of anisotropic mesh adaptation for RANS simulations over the ONERA M6 wing. The present verification study is performed using four different flow solvers, three different implementations of the metric field, and three mesh mechanics packages. Two of the flow solvers use stabilized finite-element discretizations (FUN3D-SFE and GGNS), one uses finite-volume discretization (FUN3D-FV), and the last one uses mixed finite-volume and finite element discretizations (Wolf). The mesh adaptation is based on an error estimator that aims to control the quadratic error term in the linear interpolation of Mach number. Two sets of adaptations were performed; the first one controls the interpolation error in L2 norm and the second one controls the interpolation error in L4 norm. Convergence studies were performed on the forces and the pitching moment using all four solvers, and the results are compared with previously verified convergence studies on fixed (nonadapted) meshes. Both forces and pitching moment on adapted meshes are found to be converging to the fine mesh values faster than those on fixed meshes. In addition to forces and moments, convergence of surface pressure and skin friction coefficients at various measurement locations on the wing are also presented. Adapted-mesh surface pressure distributions agree with the fine fixed mesh pressure distributions. Adapted-mesh skin friction distributions contain high frequency noise with mean values approaching the fixed mesh pressure skin friction distributions.

Aravind Balan↗

TPSAS-NF1676L-31264-DND

Nonlinear unsteady aerodynamic models (TSD, Euler/Navier-Stokes (RANS), LES) are complex, expensive solutions. CFD used by industry mostly for steady/static loads cases. Limit Cycle Oscillation (LCO) is an equilibrium point of the nonlinear aeroelastic system induced by nonlinear unsteady aerodynamics (shock dynamics, separated flow) and/or nonlinear structural dynamics (geometric and material nonlinearities). While Euler solutions can generate LCOs, viscosity yields improved accuracy (shock strength, position). Most common aeroelastic CFD codes currently consist of nonlinear unsteady aerodynamics with linear structural dynamics. Van der Pol system is a simple example of LCOs.

Walter A Silva↗

TPSAS-NF1676L-12973-DND

CFD methods are employing semi-empirical models used in analysis of hypersonic airbreathing engine flow paths. RANS (or Favre averaged) codes have models for turbulent stresses, mass and energy transport; turbulence chemistry interactions. LES methods have models for subgrid scale turbulence. Models depend upon experimental data for validation. Information on mean flow and statistics of the turbulent fluctuation in flow properties is useful. Data requirements: simple well defined supersonic combustion flows with well-known boundary conditions, time and spatially resolved, good instrument precision, and converged statistics. Approach: dual-pump CARS and planar laser-induced fluorescence imaging of OH radical (PLIF).

Andrew Cutler↗

TPSAS-NF1676L-14454-DND

FUN3D Core Capabilities - Established as a research code in late 1980’s; now supports numerous internal and external efforts across the speed range - Solves 2D/3D steady and unsteady Euler and RANS equations on node-based mixed element grids for compressible and incompressible flows - General dynamic mesh capability: any combination of rigid / overset / morphing grids, including 6-DOF effects - Aeroelastic modeling using mode shapes, full FEM, CC, etc. - Constrained / multipoint adjoint-based design and mesh adaptation - Distributed development team using agile/extreme software practices including 24/7 regression and performance testing - Capabilities fully integrated, online documentation, training videos, tutorials

Eric J Nielsen↗

OpenNEX: An open collaboration platform for the earth science community

Satellite data from the past several decades provide the most consistent record of land-surface processes that form the basis for scientific assessments of the impacts of climate variations and changes on the environment and human social-economical activities. During this time, scientific research on the characterization and assessment of environmental changes had tended to focus on large-scale land-surface changes with significant social-economic impacts. Increasingly, attention is shifting toward changes that occur more locally and that most directly relate to the everyday life of the majority of the population. In addition, there has been needs to develop management and policy decision support systems that are based on local environmental information. Almost at the same time, the advancement in sensor technology has allowed us to collect an unprecedented volume of environmental data. These data must be curated and analyzed to extract useful information for research and decision support purposes. Established in 2013 and funded by NASA, the Open NASA Earth eXchange (OpenNEX; https://opennex.org/ , Jia et al., 2019 ) project partnered with Amazon Web Services (AWS) to make available a large amount of Earth observing data, modeling results, and analysis tools on the AWS. OpenNEX provides researchers, developers, educators, and ordinary users with easy access to an integrated Earth science computational and data platform, enabling citizen scientists and application developers to realize the full value of NASA data assets and software tools. To encourage the public's engagement in this project, NASA ran virtual workshops and prize competitions. The virtual workshops provided online lectures and tutorials about how prominent scientists used the data in their research and the tutorials gave examples how to use the tools to interrogate the data in the Amazon cloud. Finally, the prize competitions allowed much wider participation in the OpenNEX project and enable testing the non-traditional projects and out-of-box ideas. OpenNEX has continued to evolve and mature. Here, we highlight new features and functionalities available to the community.

Jian Zhang↗

Variable Mixing Nozzle Design with Slotted Vortex Generators for Jet Noise Reduction

A new variable geometry turbofan nozzle concept is presented with the dual goals of airport noise reduction and high propulsive efficiency at cruise, while employing only a single moving part. The nozzle utilizes a number of curved vanes to function as vortex generators (VGs) during takeoff and initial climb. Streamwise vortices provide an increase in mixing between the jet and external stream, with a resulting decrease in high frequency noise generation and a noise source distribution which is more amenable to airframe-based shielding. Unlike conventional VGs, the vanes are positioned without transverse incidence, and a pressure difference across the surface of each vane is created via an adjacent slot. During high altitude flight, when nozzle efficiency concerns outweigh any preference for jet noise reduction, all slots may be closed through rotation of a slotted ring inside the nozzle. RANS-based numerical analysis is performed to calculate propulsive metrics, determine noise source characteristics and understand various design parameter sensitivities. A new formulation is developed to approximately correct computed noise intensities for any throttle adjustments required to maintain takeoff thrust, and a rough estimate of shielding effectiveness is proposed by means of integration over a spatially distributed noise source. A total of 19 different nozzle geometries are considered in the present study. Results indicate significant noise benefits at takeoff, including a 2 kHz fly-over noise reduction of roughly 2-3 dB, while allowing for less than 15% of the net thrust reduction at supersonic cruise calculated for a fixed-penetration chevron nozzle.

Jet Noise↗

Human Landing System Storable Propellant Architecture: Mission Design, Guidance, Navigation, and Control

In preparation of the NextSTEP-2 Appendix H awards, NASA’s Human Landing System (HLS) Program (HLSP) went through an abbreviated Design Analysis Cycle (DAC), or miniDAC, to characterize and quantify how a storable propellant design would affect and improve overall system performance of the Government reference design. As part of that miniDAC, HLSP ran multiple trade studies and assessments, including aspects of Mission Design, Flight Mechanics (FM), and Guidance, Navigation, and Control (GN\&C). The storable propellant designs were assessed through the baseline mission targeting the Near Rectilinear Halo Orbit (NRHO) and several alternate staging orbits. This paper focuses on the FM and GN\&C trades and assessments completed as part of the miniDAC as well as the GN\&C elements of the miniDAC Master Equipment List (MEL) developed for the storable propellant lander configurations.

A S Craig↗

Trajectory Simulation Using Multi Model Monte Carlo with Python (MXMCPy)

EDL (Entry, Descent and Landing) is the process from a vehicle approaching a surface to landing on it, such as a Mars rover approaching the planet before landing. POST2 (Program to Optimize Simulated Trajectories 2) is Langley’s primary EDL simulation tool and is used NASA-wide for simulations. POST2 can generate highly accurate results by running a precise, but time consuming, Monte Carlo (MC) simulation hundreds or thousands of times. Though POST2 can produce highly accurate results, it can take unrealistic time spans to generate these results, which has created a need to speed up the simulations. The new NASA software MXMCPy offers various ways to speed up the simulations while getting just as precise results. Instead of running high-precision POST2 simulations many times for traditional MC, MXMCPy can run fewer high-precision POST2 simulations and many less precise POST2 simulations and merge the results. MXMCPy contains 30+ different methods which will each suggest different allocations between model precision levels, which result in results of varying precision based on the POST2 simulation. I created Python and Bash code to automate the 5 steps of MXMCPy’s application to POST2. I also tested the precision of traditional Monte Carlo simulations to MXMCPy aided simulations and found that MXMCPy can achieve substantially more precise solutions at the same computer runtime. I learned Test Driven Development (TDD), a software programming workflow which involves writing computer-automated tests before writing the code which is being tested. These tests are ran every time the code is changed and they can find glitches in the code much quicker than a human can. This programming workflow saved me a lot of time because the automated tests could tell me exactly where the code had stopped working. I plan on using this software development method for future academic and professional software projects. I have greatly enjoyed my work at NASA, so I have been applying to NASA internships and Pathways positions. In addition, I plan on applying what I have learned about Test Driven Development to my computer science courses next semester

James Warner↗

Transition Analysis for the CRM-NLF Wind Tunnel Configuration

This paper reports the results of a comprehensive linear stability analysis of the boundary layer flow over the common research model with natural laminar flow (CRM-NLF) aircraft configuration. The flow conditions match selected test conditions from a recent wind tunnel experiment in the National Transonic Facility at the NASA Langley Research Center. Previous work has shown that the measured onset of laminar-turbulent transition during the experiments can be correlated with the linear amplification of Tollmien-Schlichting (TS) and stationary crossflow (CF) instabilities in the swept wing boundary layer. However, a significant scatter ( N ∈ (4,9)) was observed in the values of the logarithmic amplification factor along the measured transition front. This previous analysis was based on an approximate basic state (based on a boundary layer code with conical flow approximation) and parallel stability computations without surface curvature effects. Here, we examine the effects of these various approximations with the goal of quantifying the resulting changes in the N-factor correlations. Specifically, both linear stability theory (LST) and the parabolized stability equations (PSE)are used in conjunction with an accurate definition of the laminar boundary layer flow as computed with a Navier-Stokes solver with a Reynolds-Averaged-Navier-Stokes (RANS) based turbulence model within the turbulent parts of the flow. Furthermore, the effects of instability wave propagation within a fully three-dimensional boundary layer are also evaluated by integrating the disturbance growth rates along suitably chosen, curvilinear (i.e., nonplanar) propagation trajectories. The results of this analysis are also used in an accompanying paper by Venkatachari et al. to develop improved, physics based transition predictions for the same CRM-NLF configuration.

Boundary layer transition↗

Investigation of Improved Wall Heat Flux Behavior Observed When Computing Hypersonic Flow Using a Face-averaged Nodal-gradient Approach On Tetrahedral Grids

We investigate the origin of a significant reduction in computed numerical “noise” in the spatial variation of computed wall heat flux observed when using a node-centered weighted least squares(WLSQ) gradient approach when solving the Reynolds-averaged Navier-Stokes (RANS) equations for hypersonic boundary layer flows on highly stretched tetrahedral grids. This investigation is conducted using a 2nd-order, cell-centered, finite-volume, discretization which we employ to contrast results obtained using two WLSQ gradient approaches: a conventional 3-D, cell-centered, node-neighbor, WLSQ gradient (NN-CCG) approach and a 3-D, face-averaged, node-centered gradient (F-ANG)approach. The F-ANG approach, is shown to produce significantly less “noise” in the spatial variation of wall heat flux, than the NN-CCG approach and we demonstrate the cause of this “noise” via a series of numerical experiments where we compute hypersonic laminar and turbulent boundary layer flows.We further investigate the level of grid resolution required for a tetrahedral grid to approach the fidelity of a grid resolved hexahedral grid solution obtained using an equivalent structured grid method.

Jeffery A. White↗

Kibo Robot Programming Challenge 5/21/2020 Checkout Activity Report

This a report on JAXA’s Kibo Robot Programing Challenge (RPC) Checkout activity which took place on 2020-05-21. This activity is an on-orbit test for the upcoming competition organized by JAXA and to be ran by September 2020. It aimed to check the multiple aspects of the operation of Bumble in support of this competition. The activity had three primary objectives: a) verification of four different positions to be visited during the competition, namely Start Position, P1, P2, and P3, b) demonstration of Bumble autonomously moving from the Start Position to P3 without keep-in/keep-out zones and c) demonstration of Bumble autonomously moving from the Start Position to P3 with JAXA’s defined keep-in/keep-out zones. The report emphasizes the confirmation of JAXA’s previously calculated and commanded coordinates for each position and the Astrobee Facilities Team (AFT) estimation of new coordinates in case those coordinates did not match JAXA’s desired target positions per ground and crew procedures. Four estimated coordinates are proposed and various suggestions on operation procedures are given to improve the maneuverability of Bumble during the actual day of the competition.

Astrobee↗