Search NASA⌕ Search

SEARCH · Search NASA

Results for “flow solver”

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,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effects of inter-pulse coupling on nanosecond pulsed high frequency discharge ignition in a flowing mixture

This work numerically investigates the effects of non-equilibrium nanosecond plasma discharge pulse rep- etition frequency, pulse number, and flow velocity on the critical ignition volume, minimum ignition energy, and chemistry in a plasma-assisted H 2 /air flow at 300 K and 1 atm using a multi-scale adaptive reduced chemistry solver for plasma assisted combustion (MARCS-PAC). The interactions between discharges/ignition kernels spanning decoupled, partially-coupled and fully-coupled regimes in a pulse train are studied. For a single pulse discharge, increased flow velocity increases the minimum ignition energy required due to the increase of convective heat loss and flame stretch. The results show that the minimum ignition kernel prop- agation speed at the critical ignition kernel volume increases with the flow velocity. The minimum critical ignition volume decreases with the increase of plasma discharge energy. For sequential two-pulse discharges, ignition fails at both decoupled and partially-coupled regimes even when the total discharge energy is above the minimum ignition energy, but succeeds only in the fully-coupled regime at a shorter inter-pulse time. Overlap of the OH radical pool between the sequential two-pulse discharges and the increase of the chemistry effect due to the increase of reduced electric field in the fully-coupled regime contribute to the ignition enhancement. In addition, for two-pulse discharges in the fully-coupled discharge regime, the mixture can be ignited at a total energy below the minimum ignition energy of a single pulse with the same flow conditions. Moreover, for a given total discharge energy with multiple pulsed discharges, the enhancement of the ignition kernel volume has a non-monotonic dependence on discharge frequency and pulse number. The effective ignition enhancement can be achieved with an optimal pulse repetition frequency and pulse number. Furthermore, this work provides a new understanding of the mechanism for repetitive plasma ignition and insights for the optimization of plasma ignition in a reactive flow.

42 ENGINEERING↗
NASA NTRS1992

Unsteady, one-dimensional gas dynamics computations using a TVD type sequential solver

The efficacy of high resolution convection schemes to resolve sharp gradient in unsteady, 1D flows is examined using the TVD concept based on a sequential solution algorithm. Two unsteady flow problems are considered which include the problem involving the interaction of the various waves in a shock tube with closed reflecting ends and the problem involving the unsteady gas dynamics in a tube with closed ends subject to an initial pressure perturbation. It is concluded that high accuracy convection schemes in a sequential solution framework are capable of resolving discontinuities in unsteady flows involving complex gas dynamics. However, a sufficient amount of dissipation is required to suppress oscillations near discontinuities in the sequential approach, which leads to smearing of the solution profiles.

Thakur, Siddharth↗
NASA NTRS1988

CSCM Navier-Stokes thermal/aerodynamic analysis of hypersonic nozzle flows with slot injection and wall cooling

The Conservative Supra-Characteristic Method (CSCM) Navier-Stokes solver is applied to ascertain the problems inherent in the design of a nominal Mach 14 nozzle for NASA-Ames' 3.5-ft Hypersonic Wind Tunnel; attention is given to the effects of boundary layer cooling systems on the aerodynamic redesign of the nozzle throat region. Complete nozzle flowfields are calculated with and without slot injection of either hot or cold fluid into the boundary layer just upstream of the throat, as well as with alternatively adiabatic and cold walls. The CSCM method is capable of resolving subtle differences in the flows.

Codding, William H.↗
NASA NTRS2014

Simulation of Fluid Flow and Collection Efficiency for an SEA Multi-element Probe

Numerical simulations of fluid flow and collection efficiency for a Science Engineering Associates (SEA) multi-element probe are presented. Simulation of the flow field was produced using the Glenn-HT Navier-Stokes solver. Three-dimensional unsteady results were produced and then time averaged for the heat transfer and collection efficiency results. Three grid densities were investigated to enable an assessment of grid dependence. Simulations were completed for free stream velocities ranging from 85-135 meters per second, and free stream total pressure of 44.8 and 93.1 kilopascals (6.5 and 13.5 pounds per square inch absolute). In addition, the effect of angle of attack and yaw were investigated by including 5 degree deviations from straight for one of the flow conditions. All but one of the cases simulated a probe in isolation (i.e. in a very large domain without any support strut). One case is included which represents a probe mounted on a support strut within a finite sized wind tunnel. Collection efficiencies were generated, using the LEWICE3D code, for four spherical particle sizes, 100, 50, 20, and 5 micron in diameter. It was observed that a reduction in velocity of about 20% occurred, for all cases, as the flow entered the shroud of the probe. The reduction in velocity within the shroud is not indicative of any error in the probe measurement accuracy. Heat transfer results are presented which agree quite well with a correlation for the circular cross section heated elements. Collection efficiency results indicate a reduction in collection efficiency as particle size is reduced. The reduction with particle size is expected, however, the results tended to be lower than the previous results generated for isolated two-dimensional elements. The deviation from the two-dimensional results is more pronounced for the smaller particles and is likely due to the reduced flow within the protective shroud. As particle size increases differences between the two-dimensional and three dimensional results become negligible. Taken as a group, the total collection efficiency of the elements including the effects of the shroud has been shown to be in the range of 0.93 to 0.99 for particles above 20 microns. The 3D model has improved the estimated collection efficiency for smaller particles where errors in previous estimates were more significant.

Water Content↗
NASA NTRS1990

A review of reduced Navier-Stokes computations for compressible viscous flows

A reduced form of the Navier-Stokes equations, defined by a single composite of the Euler, boundary layer, and triple deck approximations, is considered for the computation of viscous interacting flows. Global pressure or pseudopotential relaxation methods are integrated with coupled sparse matrix direct solvers or coupled strongly implicit ILU algorithms to efficiently capture sharp shocks and regions of recirculation. Solutions are obtained for a variety of 2D and 3D geometries and for Mach numbers (M) spanning the range from incompressible (M = 0) to supersonic (M = 6) flow.

Rubin, S. G.↗
NASA NTRS1998

Numerical Simulations of Unsteady Low-Reynolds-Number Flows Over the APEX Airfoil

Laminar and transitional separation bubbles are an important feature of low-Reynolds-number flows over airfoils. The separation bubbles are unsteady and have a significant impact on the aerodynamic properties of the airfoils. In this paper unsteady low-Reynolds-number separated flows over the APEX airfoil are calculated using a Navier-Stokes solver. The numerical results show the presence of unsteady separation bubbles in the flowfield. An analysis of the numerical results shows that flowfield disturbances are amplified significantly in the separation bubble, leading to periodic vortex shedding. A linear stability analysis of the separated boundary layer is performed and the results show that the dominant wavenumber and frequency in the numerical simulations agree with the most unstable wavenumber and frequency from the linear stability analysis. The numerical results also show the growth and interaction of disturbance waves in the separation bubble. For transonic flows over the APEX airfoil, the calculations show that the presence of shocks causes significant changes in the separation location and consequently, the overall flowfield.

Tatineni, Mahidhar↗
NASA NTRS2006

Validation of a Pressure-Based Combustion Simulation Tool Using a Single Element Injector Test Problem

The traditional design and analysis practice for advanced propulsion systems, particularly chemical rocket engines, relies heavily on expensive full-scale prototype development and testing. Over the past decade, use of high-fidelity analysis and design tools such as CFD early in the product development cycle has been identified as one way to alleviate testing costs and to develop these devices better, faster and cheaper. Increased emphasis is being placed on developing and applying CFD models to simulate the flow field environments and performance of advanced propulsion systems. This necessitates the development of next generation computational tools which can be used effectively and reliably in a design environment by non-CFD specialists. A computational tool, called Loci-STREAM is being developed for this purpose. It is a pressure-based, Reynolds-averaged Navier-Stokes (RANS) solver for generalized unstructured grids, which is designed to handle all-speed flows (incompressible to hypersonic) and is particularly suitable for solving multi-species flow in fixed-frame combustion devices. Loci-STREAM integrates proven numerical methods for generalized grids and state-of-the-art physical models in a novel rule-based programming framework called Loci which allows: (a) seamless integration of multidisciplinary physics in a unified manner, and (b) automatic handling of massively parallel computing. The objective of the ongoing work is to develop a robust simulation capability for combustion problems in rocket engines. As an initial step towards validating this capability, a model problem is investigated in the present study which involves a gaseous oxygen/gaseous hydrogen (GO2/GH2) shear coaxial single element injector, for which experimental data are available. The sensitivity of the computed solutions to grid density, grid distribution, different turbulence models, and different near-wall treatments is investigated. A refined grid, which is clustered in the vicinity of the solid walls as well as the flame, is used to obtain a steady state solution which may be considered as the best solution attainable with the steady-state RANS methodology. From a design point of view, quick turnaround times are desirable; with this in mind, coarser grids are also employed and the resulting solutions are evaluated with respect to the fine grid solution.

Thakur, Siddarth↗
DOE OSTI2020

Machine learning coupled multi-scale modeling for redox flow batteries

The reaction distribution in macro or device-scale has been studied for redox flow batteries. The reaction distribution on electrode pore-scale structure however is not well understood, lacking especially on how the reaction distribution on the pore-scale may impact the overall performance of a flow battery. This study introduces for the first time a framework of a multi-scale model that provides understanding of the relationship between the pore-scale electrode structure reaction and the device-scale electrochemical reaction uniformity within the flow battery. A reduced order model is constructed based on 128 pore-scale simulations, which provide a quantitative relationship between the battery operation conditions (inlet velocity, current density, inlet concentration) and the surface reaction uniformity for the pore-scale sample. The multi-scale framework upscales this pore-scale surface reaction uniformity to device-scale combined uniformity. Based on the multi-scale model, a time-varying optimization of the inlet velocity is established, leading to significant reduction on pump power consumption with targeted surface reaction uniformity. The multi-scale model establishes the critical link between the micro-structure of a flow battery component and its performance at the macro-scale, therefore providing rationale for further operational or material optimization.

flow batteries, machine learning, multi-scale mode↗
NASA NTRS1981

Numerical simulation of steady transonic flow about airfoils

A computer code has been developed that couples a fast transonic full-potential AF2 solver with both an efficient integral boundary-layer method and a viscous wedge approximation of the shock/boundary-layer interaction. The efficiency of the coupled analysis methods and the method of coupling has resulted in a uniquely efficient analysis tool. The airfoil geometry is modified by the displacement thickness before the shock and the displacement thickness plus the viscous wedge thickness after the shock by considering the viscous effects as an equivalent transpiration boundary condition. The flow about conventional and supercritical airfoils under moderately strong shock situations has been calculated. Comparisons with experimental data indicate that this viscous correction method has improved the accuracy of the full-potential analysis. Furthermore, the computer time required to obtain a converged solution has been reduced.

Lee, S. C.↗
DOE OSTI2020

Operator evolution from the similarity renormalization group and the Magnus expansion

The Magnus expansion is an efficient alternative to solving similarity renormalization group (SRG) flow equations with high-order, memory-intensive ordinary differential equation solvers. The numerical simplifications it offers for operator evolution are particularly valuable for in-medium SRG calculations, though challenges remain for difficult problems involving intruder states. Here we test the Magnus approach in an analogous but more accessible situation, which is the free-space SRG treatment of the spurious bound states arising from a leading-order chiral effective field theory (EFT) potential with very high cutoffs. In this work, we show that the Magnus expansion passes these tests and then use the investigations as a springboard to address various aspects of operator evolution that have renewed relevance in the context of the scale and scheme dependence of nuclear processes. These aspects include SRG operator flow with band- versus block-diagonal generators, universality for chiral EFT Hamiltonians and associated operators with different regularization schemes, and the impact of factorization arising from scale separation. Implications for short-range correlations physics and the possibilities for reconciling high- and low-resolution treatments of nuclear structure and reactions are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1985

A Numerical Method of Calculating Propeller Noise Including Acoustic Nonlinear Effects

Using the transonic flow fields(s) generated by the NASPROP-E computer code for an eight blade SR3-series propeller, a theoretical method is investigated to calculate the total noise values and frequency content in the acoustic near and far field without using the Ffowcs Williams - Hawkings equation. The flow field is numerically generated using an implicit three dimensional Euler equation solver in weak conservation law form. Numerical damping is required by the differencing method for stability in three dimensions, and the influence of the damping on the calculated acoustic values is investigated. The acoustic near field is solved by integrating with respect to time the pressure oscillations induced at a stationary observer location. The acoustic far field is calculated from the near field primitive variables as generated by NASPROP-E computer code using a method involving a perturbation velocity potential as suggested by Hawkings in the calculation of the acoustic pressure time-history at a specified far field observed location. the methodologies described are valid for calculating total noise levels and are applicable to any propeller geometry for which a flow field solution is available.

Korkan, K. D.↗
NASA NTRS1991

Periodic blade loads of a high speed propeller at small angle of attack

A code has been developed to predict the periodic aerodynamic loads of an advanced turboprop propeller. The analytical formulation accounts for flow three-dimensionality and flow periodicity due to the propeller inclination. The flow past the blade sections is computed using a thin layer Navier-Stokes solver. An iterative procedure is used to account for the induced axial and rotational velocities. The viscous periodic results are obtained for an eight-bladed Hamilton Standard SR-7L advanced propeller at a cruise Mach number of 0.813 and 35,000 ft. altitude. The results are shown for flow field quantities and performance parameters during the blade passage in the plane of rotation illustrating the periodic nature of blade flow separation and shocks. The time averaged coefficients of thrust and power are computed and compared with available flight test data. The results obtained show excellent agreement at cruise conditions for small nacelle angles of attack.

Takallu, M. A.↗
NASA NTRS1990

ULTRA-SHARP nonoscillatory convection schemes for high-speed steady multidimensional flow

For convection-dominated flows, classical second-order methods are notoriously oscillatory and often unstable. For this reason, many computational fluid dynamicists have adopted various forms of (inherently stable) first-order upwinding over the past few decades. Although it is now well known that first-order convection schemes suffer from serious inaccuracies attributable to artificial viscosity or numerical diffusion under high convection conditions, these methods continue to enjoy widespread popularity for numerical heat transfer calculations, apparently due to a perceived lack of viable high accuracy alternatives. But alternatives are available. For example, nonoscillatory methods used in gasdynamics, including currently popular TVD schemes, can be easily adapted to multidimensional incompressible flow and convective transport. This, in itself, would be a major advance for numerical convective heat transfer, for example. But, as is shown, second-order TVD schemes form only a small, overly restrictive, subclass of a much more universal, and extremely simple, nonoscillatory flux-limiting strategy which can be applied to convection schemes of arbitrarily high order accuracy, while requiring only a simple tridiagonal ADI line-solver, as used in the majority of general purpose iterative codes for incompressible flow and numerical heat transfer. The new universal limiter and associated solution procedures form the so-called ULTRA-SHARP alternative for high resolution nonoscillatory multidimensional steady state high speed convective modelling.

Leonard, B. P.↗
DOE OSTI2020

Estimating posterior quantity of interest expectations in a multilevel scalable framework

Scalable approaches for uncertainty quantification are necessary for characterizing prediction confidence in large-scale subsurface flow simulations with uncertain permeability. To this end we explore a multilevel Monte Carlo approach for estimating posterior moments of a particular quantity of interest, where we employ an element-agglomerated algebraic multigrid (AMG) technique to generate the hierarchy of coarse spaces with guaranteed approximation properties for both the generation of spatially correlated random fields and the forward simulation of Darcy's law to model subsurface flow. In both these components (sampling and forward solves), we exploit solvers that rely on state-of-the-art scalable AMG. To illustrate the applicability of this approach, numerical tests are performed on two 3D examples-a unit cube and an egg-shaped domain with an irregular boundary-where the scalability of each simulation as well as the scalability of the overall algorithm are demonstrated.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Accelerating Subsurface Transport Simulation on Heterogeneous Clusters

Reactive transport numerical models simulate chemical and microbiological reactions that occur along a flowpath. These models have to compute reactions for a large number of locations. They solve the set of ordinary differential equations (ODEs) that describes the reaction for each location through the Newton-Raphson technique. This technique involves computing a Jacobian matrix and a residual vector for each set of equation, and then solving iteratively the linearized system by performing Gaussian Elimination and LU decomposition until convergence. STOMP, a well known subsurface flow simulation tool, employs matrices with sizes in the order of 100x100 elements and, for numerical accuracy, LU factorization with full pivoting instead of the faster partial pivoting. Modern high performance computing systems are heterogeneous machines whose nodes integrate both CPUs and GPUs, exposing unprecedented amounts of parallelism. To exploit all their computational power, applications must use both the types of processing elements. For the case of subsurface flow simulation, this mainly requires implementing efficient batched LU-based solvers and identifying efficient solutions for enabling load balancing among the different processors of the system. In this paper we discuss two approaches that allows scaling STOMP's performance on heterogeneous clusters. We initially identify the challenges in implementing batched LU-based solvers for small matrices on GPUs, and propose an implementation that fulfills STOMP's requirements. We compare this implementation to other existing solutions. Then, we combine the batched GPU solver with an OpenMP-based CPU solver, and present an adaptive load balancer that dynamically distributes the linear systems to solve between the two components inside a node. We show how these approaches, integrated into the full application, provide speed ups from 6 to 7 times on large problems, executed on up to 16 nodes of a cluster with two AMD Opteron 6272 and a Tesla M2090 per node.

GPU↗
DOE OSTI2017

Self-Excited, Multi-kHz Dynamics in a Linear, Semi-Bounded Detonation Channel

A linear, semi-bounded detonation channel was studied to characterize the injection, mixing, and ignition process in a two dimensional analog of a rotating detonation wave combustor (RDWC). The linear channel was developed to provide large optical accessibility in support of high-fidelity diagnostics. The canonical geometry was also chosen to support the generation of future numerical validation. A high degree of modularity was designed into the combustor to allow for testing multiple injector con figurations, and for reconfiguring the high-frequency instrumentation and optical window locations. Initial experiments have shown the presence of self-excited, multi-kHz combustion waves that bear the distinct characteristics of detonation waves. High-frequency pressure data details the dynamic response to detonation propagation in the combustor. The role of manifold-chamber coupling is investigated as the mechanism behind the formation of self-excited detonation waves in the combustor. We explore the formation of the waves using Schlieren imaging techniques and examine the structures present within the reaction zone and the interface between reactants and products. Companion computations were carried out using an in-house CFD solver GEMS(General Equations and Mesh Solver) to develop understanding of the structure and dynamics of the flow, flame, and their coupled interactions.

Schwinn, Kyle S.↗
NASA NTRS1990

Euler and Navier-Stokes computations for two-dimensional geometries using unstructured meshes

A general purpose unstructured mesh solver for steady-state two-dimensional inviscid and viscous flows is described. The efficiency and accuracy of the method are enhanced by the simultaneous use of adaptive meshing and an unstructured multigrid technique. A method for generating highly stretched triangulations in regions of viscous flow is outlined, and a procedure for implementing an algebraic turbulence model on unstructured meshes is described. Results are shown for external and internal inviscid flows and for turbulent viscous flow over a multi-element airfoil configuration.

Mavriplis, D. J.↗