Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Numerical Study of Buoyancy and Different Diffusion Effects on the Structure and Dynamics of Triple Flames

Numerical simulations of 2-D triple flames under gravity force have been implemented to identify the effects of gravity on triple flame structure and propagation properties and to understand the mechanisms of instabilities resulting from both heat release and buoyancy effects. A wide range of gravity conditions, heat release, and mixing widths for a scalar mixing layer are computed for downward-propagating (in the same direction with the gravity vector) and upward-propagating (in the opposite direction of the gravity vector) triple flames. Results of numerical simulations show that gravity strongly affects the triple flame speed through its contribution to the overall flow field. A simple analytical model for the triple flame speed, which accounts for both buoyancy and heat release, is developed. Comparisons of the proposed model with the numerical results for a wide range of gravity, heat release and mixing width conditions, yield very good agreement. The analysis shows that under neutral diffusion, downward propagation reduces the triple flame speed, while upward propagation enhances it. For the former condition, a critical Froude number may be evaluated, which corresponds to a vanishing triple flame speed. Downward-propagating triple flames at relatively strong gravity effects have exhibited instabilities. These instabilities are generated without any artificial forcing of the flow. Instead disturbances are initiated by minute round-off errors in the numerical simulations, and subsequently amplified by instabilities. A linear stability analysis on mean profiles of stable triple flame configurations have been performed to identify the most amplified frequency in spatially developed flows. The eigenfunction equations obtained from the linearized disturbance equations are solved using the shooting method. The linear stability analysis yields reasonably good agreements with the observed frequencies of the unstable triple flames. The frequencies and amplitudes of disturbances increase with the magnitude of the gravity vector. Moreover, disturbances appear to be most amplified just downstream of the premixed branches. The effects of mixing width and differential diffusion are investigated and their roles on the flame stability are studied.

Chen, Jyh-Yuan↗
NASA NTRS1999

Numerical Predictions of Wind Turbine Power and Aerodynamic Loads for the NREL Phase II and IV Combined Experiment Rotor

Accurate, reliable and robust numerical predictions of wind turbine rotor power remain a challenge to the wind energy industry. The literature reports various methods that compare predictions to experiments. The methods vary from Blade Element Momentum Theory (BEM), Vortex Lattice (VL), to variants of Reynolds-averaged Navier-Stokes (RaNS). The BEM and VL methods consistently show discrepancies in predicting rotor power at higher wind speeds mainly due to inadequacies with inboard stall and stall delay models. The RaNS methodologies show promise in predicting blade stall. However, inaccurate rotor vortex wake convection, boundary layer turbulence modeling and grid resolution has limited their accuracy. In addition, the inherently unsteady stalled flow conditions become computationally expensive for even the best endowed research labs. Although numerical power predictions have been compared to experiment. The availability of good wind turbine data sufficient for code validation experimental data that has been extracted from the IEA Annex XIV download site for the NREL Combined Experiment phase II and phase IV rotor. In addition, the comparisons will show data that has been further reduced into steady wind and zero yaw conditions suitable for comparisons to "steady wind" rotor power predictions. In summary, the paper will present and discuss the capabilities and limitations of the three numerical methods and make available a database of experimental data suitable to help other numerical methods practitioners validate their own work.

Duque, Earl P. N.↗
NASA NTRS2004

Adaptive Numerical Dissipative Control in High Order Schemes for Multi-D Non-Ideal MHD

The goal is to extend our adaptive numerical dissipation control in high order filter schemes and our new divergence-free methods for ideal MHD to non-ideal MHD that include viscosity and resistivity. The key idea consists of automatic detection of different flow features as distinct sensors to signal the appropriate type and amount of numerical dissipation/filter where needed and leave the rest of the region free of numerical dissipation contamination. These scheme-independent detectors are capable of distinguishing shocks/shears, flame sheets, turbulent fluctuations and spurious high-frequency oscillations. The detection algorithm is based on an artificial compression method (ACM) (for shocks/shears), and redundant multi-resolution wavelets (WAV) (for the above types of flow feature). These filter approaches also provide a natural and efficient way for the minimization of Div(B) numerical error. The filter scheme consists of spatially sixth order or higher non-dissipative spatial difference operators as the base scheme for the inviscid flux derivatives. If necessary, a small amount of high order linear dissipation is used to remove spurious high frequency oscillations. For example, an eighth-order centered linear dissipation (AD8) might be included in conjunction with a spatially sixth-order base scheme. The inviscid difference operator is applied twice for the viscous flux derivatives. After the completion of a full time step of the base scheme step, the solution is adaptively filtered by the product of a 'flow detector' and the 'nonlinear dissipative portion' of a high-resolution shock-capturing scheme. In addition, the scheme independent wavelet flow detector can be used in conjunction with spatially compact, spectral or spectral element type of base schemes. The ACM and wavelet filter schemes using the dissipative portion of a second-order shock-capturing scheme with sixth-order spatial central base scheme for both the inviscid and viscous MHD flux derivatives and a fourth-order Runge-Kutta method are denoted.

Yee, H. C.↗
NASA NTRS1997

Steady-State Cycle Deck Launcher Developed for Numerical Propulsion System Simulation

One of the objectives of NASA's High Performance Computing and Communications Program's (HPCCP) Numerical Propulsion System Simulation (NPSS) is to reduce the time and cost of generating aerothermal numerical representations of engines, called customer decks. These customer decks, which are delivered to airframe companies by various U.S. engine companies, numerically characterize an engine's performance as defined by the particular U.S. airframe manufacturer. Until recently, all numerical models were provided with a Fortran-compatible interface in compliance with the Society of Automotive Engineers (SAE) document AS681F, and data communication was performed via a standard, labeled common structure in compliance with AS681F. Recently, the SAE committee began to develop a new standard: AS681G. AS681G addresses multiple language requirements for customer decks along with alternative data communication techniques. Along with the SAE committee, the NPSS Steady-State Cycle Deck project team developed a standard Application Program Interface (API) supported by a graphical user interface. This work will result in Aerospace Recommended Practice 4868 (ARP4868). The Steady-State Cycle Deck work was validated against the Energy Efficient Engine customer deck, which is publicly available. The Energy Efficient Engine wrapper was used not only to validate ARP4868 but also to demonstrate how to wrap an existing customer deck. The graphical user interface for the Steady-State Cycle Deck facilitates the use of the new standard and makes it easier to design and analyze a customer deck. This software was developed following I. Jacobson's Object-Oriented Design methodology and is implemented in C++. The AS681G standard will establish a common generic interface for U.S. engine companies and airframe manufacturers. This will lead to more accurate cycle models, quicker model generation, and faster validation leading to specifications. The standard will facilitate cooperative work between industry and NASA. The NPSS Steady-State Cycle Deck team released a batch version of the Steady-State Cycle Deck in March 1996. Version 1.1 was released in June 1996. During fiscal 1997, NPSS accepted enhancements and modifications to the Steady-State Cycle Deck launcher. Consistent with NPSS' commercialization plan, these modifications will be done by a third party that can provide long-term software support.

VanDrei, Donald E.↗
NASA NTRS1997

Numerical Schemes for the Hamilton-Jacobi and Level Set Equations on Triangulated Domains

Borrowing from techniques developed for conservation law equations, numerical schemes which discretize the Hamilton-Jacobi (H-J), level set, and Eikonal equations on triangulated domains are presented. The first scheme is a provably monotone discretization for certain forms of the H-J equations. Unfortunately, the basic scheme lacks proper Lipschitz continuity of the numerical Hamiltonian. By employing a virtual edge flipping technique, Lipschitz continuity of the numerical flux is restored on acute triangulations. Next, schemes are introduced and developed based on the weaker concept of positive coefficient approximations for homogeneous Hamiltonians. These schemes possess a discrete maximum principle on arbitrary triangulations and naturally exhibit proper Lipschitz continuity of the numerical Hamiltonian. Finally, a class of Petrov-Galerkin approximations are considered. These schemes are stabilized via a least-squares bilinear form. The Petrov-Galerkin schemes do not possess a discrete maximum principle but generalize to high order accuracy.

Barth, Timothy J.↗
NASA NTRS2007

Development of Numerical Tools for the Investigation of Plasma Detachment from Magnetic Nozzles

A multidimensional numerical simulation framework aimed at investigating the process of plasma detachment from a magnetic nozzle is introduced. An existing numerical code based on a magnetohydrodynamic formulation of the plasma flow equations that accounts for various dispersive and dissipative processes in plasmas was significantly enhanced to allow for the modeling of axisymmetric domains containing three.dimensiunai momentum and magnetic flux vectors. A separate magnetostatic solver was used to simulate the applied magnetic field topologies found in various nozzle experiments. Numerical results from a magnetic diffusion test problem in which all three components of the magnetic field were present exhibit excellent quantitative agreement with the analytical solution, and the lack of numerical instabilities due to fluctuations in the value of del(raised dot)B indicate that the conservative MHD framework with dissipative effects is well-suited for multi-dimensional analysis of magnetic nozzles. Further studies will focus on modeling literature experiments both for the purpose of code validation and to extract physical insight regarding the mechanisms driving detachment.

Sankaran, Kamesh↗
NASA NTRS2008

Numerical Determination of Critical Conditions for Thermal Ignition

The determination of ignition or thermal explosion in an oxidizing porous body of material, as described by a dimensionless reaction-diffusion equation of the form .tu = .2u + .e-1/u over the bounded region O, is critically reexamined from a modern perspective using numerical methodologies. First, the classic stationary model is revisited to establish the proper reference frame for the steady-state solution space, and it is demonstrated how the resulting nonlinear two-point boundary value problem can be reexpressed as an initial value problem for a system of first-order differential equations, which may be readily solved using standard algorithms. Then, the numerical procedure is implemented and thoroughly validated against previous computational results based on sophisticated path-following techniques. Next, the transient nonstationary model is attacked, and the full nonlinear form of the reaction-diffusion equation, including a generalized convective boundary condition, is discretized and expressed as a system of linear algebraic equations. The numerical methodology is implemented as a computer algorithm, and validation computations are carried out as a prelude to a broad-ranging evaluation of the assembly problem and identification of the watershed critical initial temperature conditions for thermal ignition. This numerical methodology is then used as the basis for studying the relationship between the shape of the critical initial temperature distribution and the corresponding spatial moments of its energy content integral and an attempt to forge a fundamental conjecture governing this relation. Finally, the effects of dynamic boundary conditions on the classic storage problem are investigated and the groundwork is laid for the development of an approximate solution methodology based on adaptation of the standard stationary model.

Luo, W.↗
NASA NTRS2010

Experimental and Numerical Investigation of Flow Properties of Supersonic Helium-Air Jets

Heated high speed subsonic and supersonic jets operating on- or off-design are a source of noise that is not yet fully understood. Helium-air mixtures can be used in the correct ratio to simulate the total temperature ratio of heated air jets and hence have the potential to provide inexpensive and reliable flow and acoustic measurements. This study presents a combination of flow measurements of helium-air high speed jets and numerical simulations of similar helium-air mixture and heated air jets. Jets issuing from axisymmetric convergent and convergent-divergent nozzles are investigated, and the results show very strong similarity with heated air jet measurements found in the literature. This demonstrates the validity of simulating heated high speed jets with helium-air in the laboratory, together with the excellent agreement obtained in the presented data between the numerical predictions and the experiments. The very close match between the numerical and experimental data also validates the frozen chemistry model used in the numerical simulation.

Miller, Steven A. E.↗
NASA NTRS2010

Black-Hole Binaries, Gravitational Waves, and Numerical Relativity

Understanding the predictions of general relativity for the dynamical interactions of two black holes has been a long-standing unsolved problem in theoretical physics. Black-hole mergers are monumental astrophysical events ' releasing tremendous amounts of energy in the form of gravitational radiation ' and are key sources for both ground- and spacebased gravitational wave detectors. The black-hole merger dynamics and the resulting gravitational waveforms can only he calculated through numerical simulations of Einstein's equations of general relativity. For many years, numerical relativists attempting to model these mergers encountered a host of problems, causing their codes to crash after just a fraction of a binary orbit cnuld be simulated. Recently ' however, a series of dramatic advances in numerical relativity has ' for the first time, allowed stable / robust black hole merger simulations. We chronicle this remarkable progress in the rapidly maturing field of numerical relativity, and the new understanding of black-hole binary dynamics that is emerging. We also discuss important applications of these fundamental physics results to astrophysics, to gravitationalwave astronomy, and in other areas.

Kelly, Bernard J.↗
NASA NTRS2010

Error and Uncertainty Quantification in the Numerical Simulation of Complex Fluid Flows

The failure of numerical simulation to predict physical reality is often a direct consequence of the compounding effects of numerical error arising from finite-dimensional approximation and physical model uncertainty resulting from inexact knowledge and/or statistical representation. In this topical lecture, we briefly review systematic theories for quantifying numerical errors and restricted forms of model uncertainty occurring in simulations of fluid flow. A goal of this lecture is to elucidate both positive and negative aspects of applying these theories to practical fluid flow problems. Finite-element and finite-volume calculations of subsonic and hypersonic fluid flow are presented to contrast the differing roles of numerical error and model uncertainty. for these problems.

Barth, Timothy J.↗
NASA NTRS2013

Experimental Validation of Numerical Simulations for an Acoustic Liner in Grazing Flow

A coordinated experimental and numerical simulation effort is carried out to improve our understanding of the physics of acoustic liners in a grazing flow as well our computational aeroacoustics (CAA) method prediction capability. A numerical simulation code based on advanced CAA methods is developed. In a parallel effort, experiments are performed using the Grazing Flow Impedance Tube at the NASA Langley Research Center. In the experiment, a liner is installed in the upper wall of a rectangular flow duct with a 2 inch by 2.5 inch cross section. Spatial distribution of sound pressure levels and relative phases are measured on the wall opposite the liner in the presence of a Mach 0.3 grazing flow. The computer code is validated by comparing computed results with experimental measurements. Good agreements are found. The numerical simulation code is then used to investigate the physical properties of the acoustic liner. It is shown that an acoustic liner can produce self-noise in the presence of a grazing flow and that a feedback acoustic resonance mechanism is responsible for the generation of this liner self-noise. In addition, the same mechanism also creates additional liner drag. An estimate, based on numerical simulation data, indicates that for a resonant liner with a 10% open area ratio, the drag increase would be about 4% of the turbulent boundary layer drag over a flat wall.

Tam, Christopher K. W.↗
NASA NTRS2013

Numerical Investigation of a Fuselage Boundary Layer Ingestion Propulsion Concept

In the present study, a numerical assessment of the performance of fuselage boundary layer ingestion (BLI) propulsion techniques was conducted. This study is an initial investigation into coupling the aerodynamics of the fuselage with a BLI propulsion system to determine if there is sufficient potential to warrant further investigation of this concept. Numerical simulations of flow around baseline, Boundary Layer Controlled (BLC), and propelled boundary layer controlled airships were performed. Computed results showed good agreement with wind tunnel data and previous numerical studies. Numerical simulations and sensitivity analysis were then conducted on four BLI configurations. The two design variables selected for the parametric study of the new configurations were the inlet area and the inlet to exit area ratio. Current results show that BLI propulsors may offer power savings of up to 85% over the baseline configuration. These interim results include the simplifying assumption that inlet ram drag is negligible and therefore likely overstate the reduction in power. It has been found that inlet ram drag is not negligible and should be included in future analysis.

Elmiligui, Alaa A.↗
NASA NTRS2018

Numerical Investigation of Roughness Effects on Transition on Spherical Capsules

The state of the boundary layer on space re-entry vehicles significantly affects the design of the thermal protection system. However, the physical mechanism that leads to the laminar-turbulent boundary-layer transition on blunt spherical capsules remains an open question in literature. This work numerically assesses the potential of roughness-induced non-modal disturbance growth on re-entry capsules with a spherical-section forebody by optimal transient-growth theory and direct numerical simulation. Two different sets of wind-tunnel experiments are considered. Optimal transient-growth studies have been performed for the blunt capsule experiments at Mach 5.9 in the Hypersonic Ludwieg tube Braunschweig (HLB) of the Technische Universitat Braunschweig. In some of these measurements, the capsule model was equipped with a specifically designed patch of distributed micron-sized surface roughness. The transient-growth results for the HLB capsule are compared to corresponding numerical data for a Mach 6 blunt capsule experiment in the Adjustable Contour Expansion (ACE) facility of the Texas A&M University (TAMU) at lower Reynolds number. Similar trends are observed for both configurations. In particular, a rather low maximum energy gain is noted for the surface temperature conditions of the experiments. It is shown that the surface temperature dependence of the optimal transient-growth results is very similar for both capsule configurations. Moreover, the generation of stationary disturbances by well-defined roughness patches on the capsule surface is studied for the conditions of the HLB experiment using direct numerical simulations (DNS). To help explain the observed laminar-turbulent transition downstream of the roughness patch in some of the HLB capsule experiments, additional simulations were carried out to study the evolution of unsteady perturbations within the steady disturbance flow field due to the roughness patch. However, the DNS did not provide any indication of modal or non-modal disturbance growth in the wake of the roughness patch, and hence, the physical mechanism underlying the observed onset of transition remains unknown.

Hein, Stefan↗
NASA NTRS2018

A Numerical and Experimental Investigation of Flow Separation Control over a Wall-Mounted Hump Model

Numerical simulations and wind tunnel experiments were performed to investigate flow separation and its control over a wall-mounted hump model. Three-dimensional, unsteady fluid dynamic simulations were used to capture the relevant flow physics of unsteady flow separation and its control via active and passive flow control methods. Surface pressure measurements and surface oilflow visualization were obtained in the experiments. Current measurements documented a considerably thinner incoming boundary layer, but similar pressure distribution, compared to the reference data reported in the CFD (Computational Fluid Dynamics) Validation (CFDVAL 2004) Experiments. Consistent with the data available in the literature, numerical simulations predicted a longer separation bubble, which manifested itself as a shift in the pressure distributions. Streamwise vortices generated by passive vortex generators increased the suction pressure and provided substantial pressure recovery. Numerical simulations of vortex generators qualitatively predicted the effects of control, and revealed patterns of attached/separated flow regions. On the other hand, active flow control using steady discrete jets reduced, but was not able to eliminate, flow separation for the tested configuration. Although the major disagreement between the numerical and experimental results was in the pressure recovery region due to the inaccurate prediction of separation bubbles, the relative pressure recoveries from the baselines were comparable. Surface oilflow visualization and simulated surface streamlines were qualitatively in good agreement revealing the key flow features.

Koklu, Mehti↗
NASA NTRS2020

Re-Computation of Numerical Results Contained in NACA Report No. 685

In an engineering note published in the Journal of Aircraft in the year 2000, Thomas A. Zeiler made generally known that some of the early works on aeroelastic flutter by Theodore Theodorsen and I.E. Garrick (NACA Report Nos. 496, 685, and 741) contained numerical errors in some of their numerical examples. Some of the plots containing numerical errors were later reproduced in two classic aeroelasticity texts (BAH and BA). Because these foundational papers and texts are often used in graduate courses on aeroelasticity, Zeiler recommended that an effort be undertaken to employ the computational resources available today (digital computers) to recompute the example problems in these early works and to publish the results to provide a complete and error-free set of numerical examples. This paper presents recomputed theoretical results contained in NACA Report No. 685 (NACA 685), “Mechanism of Flutter, A Theoretical and Experimental Investigation of the Flutter Problem,” by Theodore Theodorsen and I.E. Garrick. The recomputations were performed employing the solution method described in NACA 685, but using modern computational tools. With some exceptions, the magnitudes and trends of the original results were in good-to-excellent agreement with the recomputed results, a surprising but gratifying result considering that the NACA 685 results were computed “by hand” using pencil, paper, slide rules, and mechanical calculators called comptometers. Checks on the recomputations (about 25% were checked) were performed using the so-called 𝑝𝑝-method of flutter solution. In all cases, including those where the original and recomputed results differed significantly, the checks were in excellent agreement with the recomputed results.

Boyd Perry, III↗
NASA NTRSDate not supplied

Quantification of Numerical Uncertainty via Nonlinear Dynamical Approach

Motivations (Ensure a Higher Level of Confidence in the Predictability & Reliability of Numerical Simulation for Multiscale Complex Nonlinear Fluid Problems) - The last two decades have been an era when computation is ahead of analysis & when very large scale practical computations are increasingly used in poorly understood multiscale complex nonlinear physical problems & non-traditional fields (Especially when computations offer the ONLY way of generating this type of data limited simulations). - At present some of the numerical uncertainties can be explained and minimized by traditional numerical analysis and standard CFD practices. However, such practices, usually based on linearized analysis, MIGHT NOT be sufficient for strongly nonlinear and/or stiff problems. - We need a good understanding of the nonlinear behavior of numerical schemes being used as an integral part of code verification, validation and certification.

HEC↗
NASA NTRSDate not supplied

Nodal Numerical Modeling of Submerged Helium Injection in a Cryogenic Propellant Tank

Subcooling of cryogenic propellant by helium injection is one of the most effective methods for suppressing bulk boiling and keeping subcooled propellant conditions for pre-launch, launch, and post-launch pressurization applications. For tank pressurization, submerged helium injection can substantially reduce helium consumption by infusing gaseous propellant into the tank ullage. This paper presents a mathematical model of the helium bubbling process in liquid oxygen to estimate the amount of oxygen vapor absorbed by the rising helium bubbles and the amount of subcooling of liquid oxygen due to evaporative heat and mass transfer. This mathematical model was incorporated in a simulation model of tank pressurization built with Generalized Fluid System Simulation Program (GFSSP), a general-purpose flow network code developed at NASA/Marshall Space Flight Center. The numerical predictions of subcooling have been compared with the experimental data of Cho et al. which investigated the propellant subcooling effect as a function of system pressure, helium injection temperature, and flowrate for a non-drained submerged injection system. The numerical predictions of helium consumption have been compared with the test data from a NASA Centaur test vehicle which included both direct and submerged injection with draining of propellants. The numerical model developed with GFSSP has been validated against two sets of experimental data and has been shown to predict both propellant subcooling and helium consumption to within 30% in most cases. The test data used for the model validation were taken in 1-g, but the mass diffusion model was developed to be applied in both 1-g and micro-g environments.

Cryogenics↗
NASA NTRSDate not supplied

Convolution-Based Numerical Solutions of Transient Temperature Fields during Powder Bed Fusion Additive Manufacturing: Theory, Accuracy, and Computational Cost

Powder bed fusion (PBF) additive manufacturing has found numerous applications in the aerospace domain. However, components fabricated via PBF have a complex time-temperature history that significantly impacts subsequent mechanical performance. This study examines convolution-based numerical solutions of transient temperature fields that support simulations involving arbitrary beam shapes and paths during PBF. The convolutional approach is verified through comparisons with analytical solutions of the temperature field. The computational speed and accuracy of the method are assessed through comparisons with other explicit and implicit numerical techniques. In addition, the straightforward translation of the approach from a CPU to a GPU implementation and the resultant performance improvement are presented. The role of the technique in predicting microstructure evolution during PBF (for a greater process-structure-property-performance framework) is also demonstrated. This work supports the development of computational materials methods for understanding and controlling the time-temperature history during PBF.

powder bed fusion↗