Search NASA⌕ Search

SEARCH · Search NASA

Results for “grid interactive”

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 361 records · Page 20

Numerical Methodology for Coupled Time-Accurate Simulations of Primary and Secondary Flowpaths in Gas Turbines

Detailed information of the flow-fields in the secondary flowpaths and their interaction with the primary flows in gas turbine engines is necessary for successful designs with optimized secondary flow streams. Present work is focused on the development of a simulation methodology for coupled time-accurate solutions of the two flowpaths. The secondary flowstream is treated using SCISEAL, an unstructured adaptive Cartesian grid code developed for secondary flows and seals, while the mainpath flow is solved using TURBO, a density based code with capability of resolving rotor-stator interaction in multi-stage machines. An interface is being tested that links the two codes at the rim seal to allow data exchange between the two codes for parallel, coupled execution. A description of the coupling methodology and the current status of the interface development is presented. Representative steady-state solutions of the secondary flow in the UTRC HP Rig disc cavity are also presented.

Przekwas, A. J.↗

The Importance of the Cathode Plume and Its Interactions with the Ion Beam in Numerical Simulations of Hall Thrusters

Hall2De is a first-principles, 2-D axisymmetric code that solves the equations of motion for ions, electrons, and neutrals on a magnetic-field-aligned grid. The computational domain downstream of the acceleration channel exit plane is large enough to include self-consistently the cathode boundary. In this paper, we present results from numerical simulations of the H6 laboratory thruster with an internally mounted cathode, with the aim of highlighting the importance of properly accounting for the interactions between the ion beam and cathode plume. The anomalous transport of electrons across magnetic field lines in Hall2De is modelled using an anomalous collision frequency, νanom, yielding νanom approximately equal to omega ce (i.e., the electron cyclotron frequency) in the plume. We first show that restricting the anomalous collision frequency to only regions where the current density of ions is large does not alter the plasma discharge in the Hall thruster as long as the interaction between the ion beam and the cathode plume is captured properly in the computational domain. This implies that the boundary conditions must be placed sufficiently far as to not interfere with the electron transport in this region. These simulation results suggest that electron transport across magnetic field lines occurs largely inside the beam and may be driven by the interactions between beam ions and electrons. A second finding that puts in relevance the importance of including the cathode plume in numerical simulations is on the significance of accounting for the ion acoustic turbulence (IAT), now known to occur in the vicinity of the cathode exit. We have included in the Hall2De simulations a model of the IAT-driven anomalous collision frequency based on Sagdeev's model for saturation of the ion-acoustic instability. This implementation has allowed us to achieve excellent agreement with experimental measurements in the near plume obtained during the operation of the H6 thruster at nominal conditions (300V, 20A) and chamber background pressure of approximately 1.5 x 10(exp -5) Torr. In addition, the numerical results obtained with the latter approach exhibit less sensitivity to background pressure than previous attempts at explaining the features of the plasma properties in the near plume.

Lopez Ortega, Alejandro↗

Using XR for Improving Scientific Discovery With Numerical Weather Models

Earth science (ES) digital twins will help us understand the complex interactions and interrelationships that make up our Earth system and the impacts of earth science phenomena on it. Our work addresses two underdeveloped areas in current ES digital twin work: improving the understanding and interaction with ES model outputs by using Virtual and Mixed Reality (XR) tools and improving the non-intuitive mapping of continuous ES natural phenomena to gridded reference frames in current numerical models. Traditionally, scientists working on ES view and analyze the results of calculated or measured observables with static 1-dimensional (1D), 2D or 3D plots displayed on flat computer screens or paper. Using such limited mediums, it can be very difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. In addition, numerical models, such as the NASA Goddard Earth Observing System (GEOS) ES model, are almost exclusively formulated, visualized and analyzed in an Eulerian reference frame with fixed grid points in space and time. However, ES phenomena such as convective clouds, hurricanes and wildfire smoke plumes are visualized and analyzed in a Lagrangian reference frame: therefore it is often difficult and unnatural to understand these phenomena in relation to each other, visualized either in an Eulerian or Lagrangian context. In 3D visualizations, data generally takes one of three forms: gridded (e.g., voxelized) data, where space is divided into regions; point clouds, where data is represented as a set of points; and meshes, where objects are rendered as surfaces composed of small polygons (usually triangles). A gridded, Eulerian reference frame has been the default representation for the 2D visual analysis of atmospheric data in part because the numerical methods used to generate atmospheric model data in the first place use a gridded approach, with equations defining the relationships between the physical variables in each of a grid's cells across successive timesteps. In our work, we are particularly interested in data from GEOS. Another reason why gridded representations tend to be used for visualizing data from such models is because trajectories are difficult to interpret from representations on 2D surfaces, due to line-of-sight ambiguity. Instead of a fixed grid from GEOS, we embed a trajectory model to simulate particles' movement throughout a GEOS run. We then ingest these particle trajectories as animated point clouds with a NASA open source XR toolkit, the Mixed Reality Exploration Toolkit (MRET), and merge GEOS data with ES phenomena data onto one combined visualization that the user can intuitively interact with. Efficient rendering of arbitrarily large point clouds is an ongoing challenge being addressed by the computer science community, with the GPU-based optimizations and efficient GPU memory utilization a common theme of recent advances, especially for XR, where sustained high frame rate is mandatory to save the user from suffering due to simulation sickness. In this work, we describe and evaluate our progress in choosing and implementing appropriate methods for rendering arbitrarily large point clouds within MRET for XR. While tracking the XR headset enables the immersion of a user within a 3D scene of a data visualization, tracking of XR handheld controllers or user’s hands enables us to implement intuitive user interactions with the visualized datasets. Conventional tools require a user working with an ES visualization to conduct many interactions to commit their intended selections or manipulations with a visualized dataset; for example to specify a set of points in 3D space. Doing so in a 2D flat screen interface has traditionally required specifying a set of points in three distinct 2D coordinate systems (XY, XZ, and YZ), which is cumbersome. In other scientific domains, it has been shown that specifying or selecting a location or volume in XR using handheld controllers or tracked hands allows for greater speed and accuracy. We anticipate the same will hold true for atmospheric data, and we will share initial results of measuring the utility of such an interface. Notably, as the data being visualized is generated by GEOS as a prediction based on initial conditions, an intended application of our tool is to serve as part of an iterative feedback loop. Through XR, a scientist will review and manipulate a GEOS model run, modifying the conditions as needed to do subsequent runs of GEOS. Thereby, XR-based improvements to speed and accuracy of 3D tagging of points minimizes the effort required by both the scientist and the computer cluster conducting the necessary calculations.

Thomas Grubb↗

Motion and interaction of decaying trailing vortices in spanwise shear wind

A simulation is presented of the drift of trailing vortices in a cross-wind near the ground by an unsteady, two-dimensional, rotational flow field with a concentration of large vorticity in vortical spots (having a finite but small effective size and finite total strength). The problem is analyzed by a combination of the method of matched asymptotic analyses for the decay of the vortical spots and the Euler solution for the unsteady rotational flow. Using the method of averaging, a special numerical method is developed in which the grid size and time step depend only on the length and velocity scales of the background flow and are independent of the effective core size of a vortical spot. The core size can be much smaller than the grid size, whereas the peak velocity in the core is inversely propertional to the spot size. Numerical results are presented to demonstrate the strong interaction between the trajectories of the vortical spots and the change of the vorticity distribution in the background flow field.

Liu, C. H.↗

Diagnosing coupled jet-streak circulations for a northern plains snow band from the operational nested-grid model

On 17 March 1989, moderate to heavy snow developed in a 100- to 200-km-wide band extending from South Dakota to northern Michigan. The 4- to 8-inch snowfall within this band was not associated with major cyclogenesis, and developed 500 to 600 km north of a stationary surface front. A diagnostic analysis based on an application of the General Meteorological Package (GEMPAK 5.0) to a numerical simulation from the operational nested-grid model (NGM) is utilized to show that the development of this snow band is related to the interaction of two upper-tropospheric jet streaks and their associated transverse circulation patterns. The eastward propagation of a jet streak from the West Coast toward the middle United States and to the south of a slower-moving jet along the U.S.-Canadian border led to a merger of the ascent maxima associated with the direct and indirect circulations of the northern and southern jets, respectively. The snow band developed as the ascending branches of the jet-streak circulation patterns merged, with the eastward propagation of the heaviest snow linked to the motion of the coupled circulation pattern. The study also demonstrates the usefulness of the operational NGM for providing the higher-resolution datasets required to relate the evolution of jet-streak circulation patterns to the development of mesoscale precipitation bands.

Hakim, Gregory J.↗

Adaptive remeshing for three-dimensional compressible flow computations

A practical solution algorithm for steady 3D Euler flows is presented. This algorithm employs coupling of a surface triangulator, an automatic tetrahedral mesh generator, an unstructured grid flow solver, and an error estimation procedure. The performance of the method is illustrated using a shock interaction problem in high Mach number flow over a swept circular cylinder.

Peraire, J.↗

Tracking of raindrops in flow over an airfoil

The splashback that occurs when raindrops impact an airfoil results in an 'ejecta fog' of small droplets around the leading edge. Acceleration of these droplets by the air flow field is a momentum sink for the air flow and has been hypothesized to contribute to the degradation of airfoil performance in heavy rain. Presented here is a one-way coupled Eulerian-Lagrangian particle tracking scheme to evaluate droplet number densities and momentum sink terms around a NACA 64-210 airfoil section for three rainfall rates. A laminar air flow field is determined with a standard CFD code and input to the particle tracking algorithm. Raindrops are assumed to be non-interacting, non-deforming, non-evaporating, and non-spinning spheres and are tracked through the curvilinear grid used by the air flow code. A simple model is used to simulate impacts and the resulting splashback on the airfoil surface.

Valentine, James R.↗

Modal element method for potential flow in non-uniform ducts: Combining closed form analysis with CFD

An analytical procedure is presented, called the modal element method, that combines numerical grid based algorithms with eigenfunction expansions developed by separation of variables. A modal element method is presented for solving potential flow in a channel with two-dimensional cylindrical like obstacles. The infinite computational region is divided into three subdomains; the bounded finite element domain, which is characterized by the cylindrical obstacle and the surrounding unbounded uniform channel entrance and exit domains. The velocity potential is represented approximately in the grid based domain by a finite element solution and is represented analytically by an eigenfunction expansion in the uniform semi-infinite entrance and exit domains. The calculated flow fields are in excellent agreement with exact analytical solutions. By eliminating the grid surrounding the obstacle, the modal element method reduces the numerical grid size, employs a more precise far field boundary condition, as well as giving theoretical insight to the interaction of the obstacle with the mean flow. Although the analysis focuses on a specific geometry, the formulation is general and can be applied to a variety of problems as seen by a comparison to companion theories in aeroacoustics and electromagnetics.

Baumeister, Kenneth J.↗

Low Dissipative High Order Shock-Capturing Methods using Characteristic-Based Filters

An approach which closely maintains the non-dissipative nature of classical fourth or higher- order spatial differencing away from shock waves and steep gradient regions while being capable of accurately capturing discontinuities, steep gradient and fine scale turbulent structures in a stable and efficient manner is described. The approach is a generalization of the method of Gustafsson and Olsson and the artificial compression method (ACM) of Harten. Spatially non-dissipative fourth or higher-order compact and non-compact spatial differencings are used as the base schemes. Instead of applying a scalar filter as in Gustafsson and Olsson, an ACM like term is used to signal the appropriate amount of second or third-order TVD or ENO types of characteristic based numerical dissipation. This term acts as a characteristic filter to minimize numerical dissipation for the overall scheme. For time-accurate computations, time discretizations with low dissipation are used. Numerical experiments on 2-D vortical flows, vortex-shock interactions and compressible spatially and temporally evolving mixing layers showed that the proposed schemes have the desired property with only a 10% increase in operations count over standard second-order TVD schemes. Aside from the ability to accurately capture shock-turbulence interaction flows, this approach is also capable of accurately preserving vortex convection. Higher accuracy is achieved with fewer grid points when compared to that of standard second-order TVD or ENO schemes. To demonstrate the applicability of these schemes in sustaining turbulence where shock waves are absent, a simulation of 3-D compressible turbulent channel flow in a small domain is conducted.

Yee, H. C.↗

Low Dissipative High Order Shock-Capturing Methods Using Characteristic-Based Filters

An approach which closely maintains the non-dissipative nature of classical fourth or higher- order spatial differencing away from shock waves and steep gradient regions while being capable of accurately capturing discontinuities, steep gradient and fine scale turbulent structures in a stable and efficient manner is described. The approach is a generalization of the method of Gustafsson and Oisson and the artificial compression method (ACM) of Harten. Spatially non-dissipative fourth or higher-order compact and non-compact spatial differencings are used as the base schemes. Instead of applying a scalar filter as in Gustafsson and Olsson, an ACM like term is used to signal the appropriate amount of second or third-order TVD or ENO types of characteristic based numerical dissipation. This term acts as a characteristic filter to minimize numerical dissipation for the overall scheme. For time-accurate computations, time discretizations with low dissipation are used. Numerical experiments on 2-D vortical flows, vortex-shock interactions and compressible spatially and temporally evolving mixing layers showed that the proposed schemes have the desired property with only a 10% increase in operations count over standard second-order TVD schemes. Aside from the ability to accurately capture shock-turbulence interaction flows, this approach is also capable of accurately preserving vortex convection. Higher accuracy is achieved with fewer grid points when compared to that of standard second-order TVD or ENO schemes. To demonstrate the applicability of these schemes in sustaining turbulence where shock waves are absent, a simulation of 3-D compressible turbulent channel flow in a small domain is conducted.

Yee, H. C.↗

Applications of the Space-Time Conservation Element and Solution Element (CE/SE) Method to Computational Aeroacoustic Benchmark Problems

The Internal Propagation problems, Fan Noise problem, and Turbomachinery Noise problems are solved using the space-time conservation element and solution element (CE/SE) method. The problems in internal propagation problems address the propagation of sound waves through a nozzle. Both the nonlinear and linear quasi 1D Euler equations are solved. Numerical solutions are presented and compared with the analytical solution. The fan noise problem concerns the effect of the sweep angle on the acoustic field generated by the interaction of a convected gust with a cascade of 3D flat plates. A parallel version of the 3D CE/SE Euler solver is developed and employed to obtain numerical solutions for a family of swept flat plates. Numerical solutions for sweep angles of 0, 5, 10, and 15 deg are presented. The turbomachinery problems describe the interaction of a 2D vortical gust with a cascade of flat-plate airfoils with/without a downstream moving grid. The 2D nonlinear Euler Equations are solved and the converged numerical solutions are presented and compared with the corresponding analytical solution. All the comparisons demonstrate that the CE/SE method is capable of solving aeroacoustic problems with/without shock waves in a simple and efficient manner. Furthermore, the simple non-reflecting boundary condition used in the CE/SE method which is not based on the characteristic theory works very well in 1D, 2D and 3D problems.

Wang, Xiao-Yen↗

ICEG2D (v2.0) - An Integrated Software Package for Automated Prediction of Flow Fields for Single-Element Airfoils With Ice Accretion

An integrated geometry/grid/simulation software package, ICEG2D, is being developed to automate computational fluid dynamics (CFD) simulations for single- and multi-element airfoils with ice accretions. The current version, ICEG213 (v2.0), was designed to automatically perform four primary functions: (1) generate a grid-ready surface definition based on the geometrical characteristics of the iced airfoil surface, (2) generate high-quality structured and generalized grids starting from a defined surface definition, (3) generate the input and restart files needed to run the structured grid CFD solver NPARC or the generalized grid CFD solver HYBFL2D, and (4) using the flow solutions, generate solution-adaptive grids. ICEG2D (v2.0) can be operated in either a batch mode using a script file or in an interactive mode by entering directives from a command line within a Unix shell. This report summarizes activities completed in the first two years of a three-year research and development program to address automation issues related to CFD simulations for airfoils with ice accretions. As well as describing the technology employed in the software, this document serves as a users manual providing installation and operating instructions. An evaluation of the software is also presented.

Thompson David S.↗

Visual Computing Environment

The Visual Computing Environment (VCE) is a NASA Lewis Research Center project to develop a framework for intercomponent and multidisciplinary computational simulations. Many current engineering analysis codes simulate various aspects of aircraft engine operation. For example, existing computational fluid dynamics (CFD) codes can model the airflow through individual engine components such as the inlet, compressor, combustor, turbine, or nozzle. Currently, these codes are run in isolation, making intercomponent and complete system simulations very difficult to perform. In addition, management and utilization of these engineering codes for coupled component simulations is a complex, laborious task, requiring substantial experience and effort. To facilitate multicomponent aircraft engine analysis, the CFD Research Corporation (CFDRC) is developing the VCE system. This system, which is part of NASA's Numerical Propulsion Simulation System (NPSS) program, can couple various engineering disciplines, such as CFD, structural analysis, and thermal analysis. The objectives of VCE are to (1) develop a visual computing environment for controlling the execution of individual simulation codes that are running in parallel and are distributed on heterogeneous host machines in a networked environment, (2) develop numerical coupling algorithms for interchanging boundary conditions between codes with arbitrary grid matching and different levels of dimensionality, (3) provide a graphical interface for simulation setup and control, and (4) provide tools for online visualization and plotting. VCE was designed to provide a distributed, object-oriented environment. Mechanisms are provided for creating and manipulating objects, such as grids, boundary conditions, and solution data. This environment includes parallel virtual machine (PVM) for distributed processing. Users can interactively select and couple any set of codes that have been modified to run in a parallel distributed fashion on a cluster of heterogeneous workstations. A scripting facility allows users to dictate the sequence of events that make up the particular simulation.

Lawrence, Charles↗

Forebody and Inlet Design for the HIFiRE 2 Flight Test

A forebody and inlet have been designed for the HIFiRE 2 scramjet flight test. The test will explore the operating, performance, and stability characteristics of a simple hydrocarbon-fueled scramjet combustor as it transitions from dual-mode to scramjet-mode operation and during supersonic combustion at Mach 8+ flight conditions. Requirements for the compression system were derived from inlet starting and combustor inflow requirements as well as physical size constraints. The design process is described. A planar, fixed geometry, mixed compression concept was used to produce laterally uniform flow at the inlet entrance and a conservative amount of internal contraction with respect to inlet starting. A grid sensitivity study was performed so that important flow physics caused by three-dimensional shock boundary layer interactions could be captured with confidence. Results from low Mach number operability studies, nominal trajectory cases, and high dynamic pressure heat load cases are discussed. The forebody and inlet solutions provide information for on-going combustor calculations, mass capture across the trajectory for fuel system design, and surface heating rates for thermal/structural analysis. The design has a one freestream Mach number margin for inlet starting, exceeds the high Mach number combustor entrance pressure requirement, produces high quality flow at the inlet exit for all Mach numbers and vehicle attitudes in the design space, and fits inside the booster shroud.

Ferlemann, Paul G.↗

Experimental Database with Baseline CFD Solutions: 2-D and Axisymmetric Hypersonic Shock-Wave/Turbulent-Boundary-Layer Interactions

A database compilation of hypersonic shock-wave/turbulent boundary layer experiments is provided. The experiments selected for the database are either 2D or axisymmetric, and include both compression corner and impinging type SWTBL interactions. The strength of the interactions range from attached to incipient separation to fully separated flows. The experiments were chosen based on criterion to ensure quality of the datasets, to be relevant to NASA's missions and to be useful for validation and uncertainty assessment of CFD Navier-Stokes predictive methods, both now and in the future. An emphasis on datasets selected was on surface pressures and surface heating throughout the interaction, but include some wall shear stress distributions and flowfield profiles. Included, for selected cases, are example CFD grids and setup information, along with surface pressure and wall heating results from simulations using current NASA real-gas Navier-Stokes codes by which future CFD investigators can compare and evaluate physics modeling improvements and validation and uncertainty assessments of future CFD code developments. The experimental database is presented tabulated in the Appendices describing each experiment. The database is also provided in computer-readable ASCII files located on a companion DVD.

Database↗

Simulations of Cryogenic Line Chilldown with Advanced Sub-Grid Wall Boiling Models

A meso-scale model developed at MIT [2][11] for boiling processes in water was adapted for cryogens and demonstrated for chilldown of propellant transfer lines in both liquid nitrogen and hydrogen. The sub-grid boiling model accurately captures the contributions to heat transfer from the generation of bubble nuclei, growth, and interaction of the bubbles in the microlayer as well as quenching of the boiling surface following bubble departure. It was adapted for cryogenic fluids using thermodynamic scaling concepts taking into account non-dimensional pressures and temperatures that are scaled by the corresponding critical values for the fluid. The boiling model was demonstrated for line chilldown in liquid nitrogen (Darr et al. [6]). The predicted wall temperature at which quenching occurs was close to the test data while the slope of the temperature curve after quenching is initiated, showing a steeper variation than the test data. Simulations were also performed for liquid hydrogen by simulating experiments of Hartwig et al. [7]. Chilldown times in liquid hydrogen are much more rapid due to higher heat transfer in the film boiling regime and accounting for the higher turbulence levels were found to be important. Furthermore, at the much lower fluid temperatures in liquid hydrogen flows, accounting for the variable thermal properties of the solid material is critical and has a dramatic impact on the quench times. The need for additional studies to better understand the evolution of the quench front in liquid hydrogen was noted.

Computational Fluid Dynamics↗

Simulations of Cryogenic Line Chilldown with Advanced Sub-Grid Wall Boiling Models

A meso-scale model developed at MIT [2][11] for boiling processes in water was adapted for cryogens and demonstrated for chilldown of propellant transfer lines in both liquid nitrogen and hydrogen. The sub-grid boiling model accurately captures the contributions to heat transfer from the generation of bubble nuclei, growth, and interaction of the bubbles in the microlayer as well as quenching of the boiling surface following bubble departure. It was adapted for cryogenic fluids using thermodynamic scaling concepts taking into account non-dimensional pressures and temperatures that are scaled by the corresponding critical values for the fluid. The boiling model was demonstrated for line chilldown in liquid nitrogen (Darr et al. [6]). The predicted wall temperature at which quenching occurs was close to the test data while the slope of the temperature curve after quenching is initiated, showing a steeper variation than the test data. Simulations were also performed for liquid hydrogen by simulating experiments of Hartwig et al. [7]. Chilldown times in liquid hydrogen are much more rapid due to higher heat transfer in the film boiling regime and accounting for the higher turbulence levels were found to be important. Furthermore, at the much lower fluid temperatures in liquid hydrogen flows, accounting for the variable thermal properties of the solid material is critical and has a dramatic impact on the quench times. The need for additional studies to better understand the evolution of the quench front in liquid hydrogen was noted.

Computational Fluid Dynamics↗

Comparison of solutions of the Navier-Stokes and interacting boundary layer equations for separated turbulent flow

Comparisons are presented between numerical solutions of the Navier-Stokes and interacting boundary-layer equations for the subsonic separated turbulent flow over an axisymmetric boattail configuration for which experimental data was available. Comparisons were made using a standard algebraic turbulence model and its modification using a streamwise relaxation model. Grid size studies were performed to provide an approximate assessment of the sensitivity of each method to mesh size. In general, good agreement was observed between the solutions of the Navier-Stokes and interacting boundary-layer equations; however, the interaction approach responded more strongly to the modified turbulence model thereby giving better agreement with experimental data.

Vatsa, V. N.↗