Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “high-order”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Effect of Spatial Filtering in Implicit Large-Eddy Simulations of Separated Flows

The relatively high Reynolds number of turbulent flows encountered in various applications puts these problems well beyond the reach of direct numerical simulation (DNS) at present. Meanwhile, lower-fidelity Reynolds-averaged Navier-Stokes (RANS) calculations are known to be not accurate enough in complex problems, such as smooth-body flow separation and other flows involving highly-unsteady phenomena. Hence, given the current infeasibility of DNS and the unsatisfactory performance of RANS, intermediate techniques such as large-eddy simulation (LES) and hybrid RANS-LES, whose fidelity lie between RANS and DNS, have received much attention for application to various problems of practical importance. Modeling of the effect of missing scales on resolved scales, also known as subgrid-scale (SGS) modeling, is an important subject for LES. SGS models can be broadly categorized as explicit or implicit approaches. The explicit approach is based on an SGS model that explicitly appears in the governing equations expressed in the form of so-called “filtered Navier-Stokes equations”, which describe the evolution of the turbulence scales resolved by the LES grid. The effect of the scales unresolved by the grid is represented by the SGS model. The implicit modeling approach, on the other hand, does not employ an explicit model but instead treats the intrinsic dissipation of the numerical discretization scheme as an implicit SGS model. An LES without an explicit SGS model is commonly termed as an implicit LES (ILES). The relative merits of one SGS modeling approach over another is a subject of ongoing debate. We have opted to employ an ILES methodology, based on high-order compact finite-difference and spatial filtering schemes, in our recent investigations of separated flow problems [1, 2]. Further discussion of our preference of ILES over explicit LES is provided in Uzun and Malik [2]. The spatial filtering operation, described in the next section, is treated as an implicit SGS model for the ILES. Some observations made during the course of our recent investigations, which pointed out to excessive numerical dissipation in certain parts of the flowfield, prompted us to take a closer look at the potential effect of the spatial filter on ILES predictions. This technical note is therefore devoted to spatial filter effects in the context of a high Reynolds number, transonic shock-induced separated flow.

Ali Uzun↗
NASA NTRS2019

eddy Users Manual

eddy is a collection of tools - nonlinear solvers, meshing, post-processing, visualization, optimization, etc. - for performing scale-resolving simulations of multi-physics applications. The framework is designed to enable advanced R&D on a variety of topics by leveraging a mature capability for scale resolving simulations, and simultaneously be an appropriate tool for application analysis and support. Currently, eddy is at a relatively low technical readiness level (TRL), and users and developers should maintain appropriate expectations. The technical details behind eddy are outlined in several publications which can be consulted for more information [1–10]. The solvers are built around an unstructured high-order capability, and heavily utilize the tensor product sum-factorization approach for efficiency. The unsteady formulation utilizes a fully implicit space-time approach with a matrix-free Newton- Krylov method. A primitive steady-state solver is available for testing purposes, but is not expected to converge for all but simple verification cases. The Navier-Stokes fluid solvers do not support either RANS or hybrid-RANS capability, only LES and wall-modeled LES approaches. All of the solvers within eddy support three modes of operation: a primal solve of the full nonlinear problem, and two linearization approaches of the primal solve - the ad joint and the tangent solution. Details on how to select and use these three modes are outlined in Sec. 3.

Murman, Scott M.↗
NASA NTRS2020

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↗
NASA NTRSDate not supplied

TPSAS-NF1676L-12494-DND

Outline: NASA Missions and their requirements for Computational Fluid Dynamics (CFD) methodology, CFD state-of-the-art employed in NASA applications, examples of applications challenging NASA’s CFD capability, desired characteristics of future CFD solvers and the role of high-order methods in their development, and a vision for future production CFD methods.

David M Schuster↗
NASA NTRSDate not supplied

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↗
NASA NTRS2020

GeneLab

GeneLab collects and enables analysis of spaceflight and ground-based spaceflight simulation genomic data, RNA and protein expression, and metabolic profiles. It interfaces with other existing databases containing spaceflight omic data. The 2011 National Research Council (NRC) Decadal Survey on NASA Life and Physical Sciences called for increased opportunities for multi-investigator spaceflight opportunities and greater use of genomic approaches to meet the needs of NASA researchers. To address these recommendations of the NRC Decadal Survey, the Space Life and Physical Sciences Research and Applications Division of NASA's Human Exploration and Operations Mission Directorate has initiated a transition to an Open Science architecture to increase research opportunities, and has developed the GeneLab Platform based on highly leveraged and integrated bioinformatics analytics. GeneLab is an interactive, open-access resource where scientists can upload, download, store, search, share, transfer, and analyze omics data from spaceflight and corresponding analogue experiments. Users can explore GeneLab datasets in the Data Repository, analyze data using the Analysis Platform, visualize high-order data and create collaborative projects using the Collaborative Workspace. Our primary goal is to maximize the utilization of the valuable biological research conducted aboard the International Space Station (ISS) by collecting genomic, transcriptomic, proteomic, and metabolomics data known as “omics”. By providing a portal linking processed data to flight parameters, GeneLab enables exploration of the molecular network responses of terrestrial biology to the space environment. This allows researchers to understand the complex responses of biological systems to the space environment. This technology development activity was transferred from the Human Exploration and Operations Mission Directorate to the Science Mission Directorate Division of Biological and Physical Sciences (BPS) in October 2020.

GeneLab↗
NASA NTRSDate not supplied

A Framework for Mesh-Geometry Associativity during Mesh Adaptation

A framework has been developed for describing how a computational mesh is associated to the geometry model it discretizes. The target application of this framework is surface mesh adaptation in a CFD flow solver. The framework, called MeshLink, consists of two components. First, a schema has been defined for describing the one-to-many associativity of a surface mesh to the geometry model entities to which it is attached. Second, a high-level library provides a kernel-agnostic wrapper for providing the necessary geometry queries to an application (e.g., mesher, flow solver). Both the schema and library are provided freely and openly. MeshLink’s ability to support solution mesh adaptation on linear and curved meshes for a high-order flow solution is demonstrated on several test cases relevant to the aerospace and automotive industries.

mesh adaptation↗
NASA NTRSDate not supplied

Dust Particle Aeroheating Calculations for Mars Entry Hypersonic Flows

The purpose of the current research effort is to predict particle trajectories in a hypersonic flowfield while accounting for all relevant particle-fluid and particle-particle interactions. The resulting particle solver is titled DUst Simulation & Tracking (DUST) and works in conjunction with the US3D CFD flow solver. The key elements underpinning the current work are an efficient mesh-localiztion algorithm, time-driven hard-sphere model, point-to-point MPI framework, coarse-graining using computational parcels, and high-order Adams-Bashforth time stepping. The coupled US3D-DUST framework will be applied to conduct a multi-physics examination of dust-laden flows around the the Schiarparelli capsule.

Dusty Flows↗
NASA NTRSDate not supplied

Characterization of Collection Efficiency of the Common Research Model Midspan Wing Section in the IRT

This paper presents a preliminary study for the characterization of collection efficiency from icing tests conducted in the Icing Research Tunnel at NASA Glenn Research Center. A test method previously developed for measuring the attachment line maximum collection efficiency of a swept NACA 0012 airfoil model at zero angle of attack was applied to the leading-edge region of a 65%-scale version of the Common Research Model midspan wing section. A correlation for the stagnation line maximum collection efficiency as a function of the modified inertia parameter was obtained with LEWICE3D simulations utilizing a discrete number of drop diameters. It was then compared with the collection efficiency measurement data obtained in the IRT. For the experimental collection efficiency, two ice shape digitization procedures were utilized to extract 2-D chord-wise ice shape profiles, i.e., the Maximum Combined Cross Section or MCCS, and the Minimum Combined Cross Section or MinCCS, at selected span-wise locations from the 3-D scanned ice shapes. The preliminary result showed that the rime ice thickness method can be used to characterize the collection efficiency distribution, for conditions free of ice erosion, in the main ice shape region. However, a high-order statistical processing of the ice scan is needed for estimating the mean ice thickness in areas where feathers are prevalent. From the limited comparison of the experimental and LEWICE3D collection efficiency data in the CRM65 MS model leading edge area, it was shown that the attachment line maximum collection efficiency is reasonably estimated by the best curve-fit correlation in the range of modified inertia parameter tested. A tighter alignment control of the iced and cold clean model scans as well as the model angle of attack are needed to improve the comparison. Further evaluation of this correlation is recommended to assess its applicability for Common Research Model type swept wing icing scaling analysis.

Icing Scaling↗
DOE OSTI2025

Coupling Approaches with Non-matching Grids for Classical Linear Elasticity and Bond-based Peridynamic Models in 1D

Local-nonlocal coupling approaches provide a means to combine the computational efficiency of local models and the accuracy of nonlocal models. To facilitate the coupling of the two models, non-matching grids are often desirable as nonlocal grids usually require a finer resolution than local grids. In that case, it is often convenient to resort to interpolation operators so that models can exchange information in the overlap regions when nodes from the two grids do not coincide. This paper studies three existing coupling approaches, namely 1) a method that enforces matching displacements in an overlap region, 2) a variant that enforces a constraint on the stresses instead, and 3) a method that considers a variable horizon in the vicinity of the interfaces. Further, the effect of the interpolation order and of the grid ratio on the performance of the three coupling methods with non-matching grids is carefully studied on one-dimensional examples using polynomial manufactured solutions. The numerical results show that the degree of the interpolants should be chosen with care to avoid introducing additional modeling errors, or simply minimize these errors, in the coupling approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

An extension of the localized artificial diffusivity method for immiscible and high density ratio flows

The localized artificial diffusivity (LAD) method is widely regarded as the preferred multi-material regularization scheme for the compact finite difference method, because it is conservative, easy to implement, and generally robust for a wide range of multi-material problems. However, traditional LAD methods face significant challenges when applied to flows with large density ratios and when maintaining thermodynamic equilibrium across material interfaces. These limitations arise from the formulation of the artificial diffusivity flux and the reliance on enthalpy diffusion for interface regularization. Additionally, traditional LAD methods struggle to ensure stability under large density ratio conditions, fail to maintain a finite interface thickness, and are therefore unsuitable for modeling immiscible interfaces. Here, in this work, we discuss the origins of these issues in traditional LAD methods and propose modifications which enable the simulation of large density ratio and immiscible flows. The proposed method targets the artificial diffusion fluxes at gradients and ringing in the volume fraction, rather than the mass fraction in traditional methods, to consistently regularize large density ratio interfaces. Furthermore, the proposed method introduces an artificial bulk density diffusion term to enforce equilibrium conditions across interfaces. To address the challenge of modeling immiscible flows, a conservative diffuse interface term is incorporated into the formulation to ensure a finite interface thickness. Specific consideration is taken in the design of the method to ensure that these crucial properties are maintained for N -material flows. The effectiveness of the proposed method is demonstrated through a series of canonical test cases, and its accuracy is validated by comparison with experimental data on micro-bubble collapse in water. These results highlight the method’s robustness and its ability to overcome the limitations of traditional LAD approaches.

Artificial diffusivity↗
NASA NTRS2004

Nonlinear optics and crystalline whispering gallery mode resonators

We report on our recent results concerning fabrication of high-Q whispering gallery mode (WGM) crystalline resonators, and discuss some possible applications of lithium niobate WGM resonators in nonlinear optics and photonics. In particular, we demonstrate experimentally a tunable third-order optical filter fabricated from the three metalized resonators; and report observation of parametric frequency dobuling in a WGM resonator made of periodically poled lithium niobate (PPLN).

whispering gallery mode↗
NASA NTRS2020

Entropy Stable h/p-Nonconforming Discretization with the Summation-by-Parts Property for the Compressible Euler and Navier–Stokes Equations

In this paper, we extend the entropy conservative/stable algorithms presented by Del Rey Fernandez and coauthors for the compressible Euler and Navier-Stokes equations on nonconforming p-refined/coarsened curvilinear grids to h/p refinement/coarsening. The main difficulty in developing nonconforming algorithms is the construction of appropriate coupling procedures across nonconforming interfaces. Here, we utilize a computationally simple and efficient approach based upon using decoupled interpolation operators. The resulting scheme is entropy conservative/stable and elementwise conservative. Numerical simulations of the isentropic vortex and viscous shock propagation con firm the entropy conservation/stability and accuracy properties of the method (achieving ~ p + 1 convergence), which are comparable to those of the original conforming scheme. Simulations of the Taylor{Green vortex at R(e) = 1,600 and turbulent flow past a sphere at R(e(infinity)) = 2,000 show the robustness and stability properties of the overall spatial discretization for unstructured grids. Finally, to demonstrate the entropy conservation property of a fully-discrete explicit entropy stable algorithm with h=p refinement/coarsening, we present the time evolution of the entropy function obtained by simulating the propagation of the isentropic vortex using a relaxation Runge-Kutta scheme.

Nonconforming interfaces↗
NASA NTRSDate not supplied

Explicit Discontinuous Galerkin Methods for Conservation Laws

The two explicit DG methods in this study are based on a ‘predictor-corrector’ formulation, the first introduced by Lörcher, Gassner, and Munz (2007, 2008) called space–time expansion discontinuous Galerkin or STE-DG scheme, and the second, introduced independently by the author (Huynh 2006, 2013) called the upwind moment scheme. The predictor step of the two methods is essentially identical using a Cauchy-Kovalevsky (CK) procedure, which involves no interaction of the data among neighboring cells. The corrector step also shares the same space-time integration formulation and is where interaction of the data among neighboring cells takes place; the difference, however, is in how the resulting space-time volume integral is estimated. As a consequence of the different estimates, for the case of advection in one spatial dimension (1D), the moment scheme has a CFL (Courant-Friedrichs-Lewy) condition of 1 for all p and is accurate to order 2p+1, i.e., it possesses the super accuracy property, whereas the STE-DG method has a more restrictive CFL condition and is accurate to the expected order of p+1. For 1D advection, compared with the CFL conditions of 1/(2p+1) of standard RK-DG (Runge-Kutta) scheme where space and time discretization are of the same order, the moment scheme allows a significantly larger time step size. It also turns out that the scheme yields a result identical to Van Leer’s scheme III (1977), which amounts to shifting the data a distance of advection corresponding to the time step and projecting the result onto the space of polynomial solutions. Contrary to Van Leer’s approach, however, the space-time ‘predictor-corrector’ formulation facilitates extensions to the case of systems of equations. Concerning 2D extensions, in the case of advection, when the flow is along the diagonal direction, the CFL conditions for the moment schemes become restrictive as will be shown by Fourier (Von Neumann) stability and accuracy analyses. Since the moment scheme employs the right Radau points as collocation points in time, the method is closely related to the implicit Radau IIA scheme, which is stable for any time step size. The role of Radau IIA in relieving stability restriction for these explicit DG schemes remains to be explored

Discontinuous Galerkin↗