Modeling supersonic combustion using a fully implicit numerical method
Explore the source record for details and available documents.
Wilson, Gregory J.↗
SEARCH · Search NASA
Results for “Implicit modeling”
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.
About the coupling of turbulence closure models with averaged Navier-Stokes equations
The MacCormack implicit predictor-corrector model (1981) for numerical solution of the coupled Navier-Stokes equations for turbulent flows is extended to nonconservative multiequation turbulence models, as well as the inclusion of second-order Reynolds stress turbulence closure. A scalar effective pressure turbulent contribution to the pressure field is defined to approximate the effects of the Reynolds stress in strongly sheared flows. The Jacobian matrices of the transport equations are diagonalized to reduce the required computer memory and run time. Techniques are defined for including turbulence in the diagonalization. Application of the method is demonstrated with solutions generated for transonic nozzle flow and for the interaction between a supersonic flat plate boundary layer and a 12 deg compression-expansion ramp.
Vandromme, D.↗
Artificial dissipation models for the Euler equations
Various artificial dissipation models which are used with central difference algorithms for the Euler equations are analyzed for their effect on accuracy, stability and convergence rates. In particular, linear and nonlinear models are investigated using an implicit approximate factorization code (ARC2D) for transonic airfoils. Fully implicit application of the dissipation models is shown to improve robustness and convergence rates. The treatment of dissipation models at boundaries will be examined. It will be shown that accurate, error free solutions with sharp shocks can be obtained using a central difference algorithm coupled with an appropriate nonlinear artificial dissipation model.
Pulliam, T. H.↗
Implicit Large-Eddy Simulations of Hot and Cold Supersonic Jets in Loci-CHEM
Turbulence is a critical factor in most aero- and propulsion-related flows. For many years, the industry-standard approach to incorporating the effect of turbulence in computational fluid dynamics (CFD)calculations has been the Reynolds Averaged Navier Stokes (RANS) modeling approach. RANS implicitly time-averages turbulent motion, and models its effect on the mixing of species, momentum, and energy in a flow field. The approach is economical and often allows a CFD calculation to proceed to steady-state. However, it makes several simplifying assumptions, and the results will only be as good as the RANS turbulence model. Meanwhile, the last two decades have witnessed rapid and explosive growth in computational power, the development of CFD codes which scale well over thousands of processors, and the refinement of unstructured grid-generation tools which facilitate rapid surface and volume gridding of complex geometries. CFD engineering calculations of 108– 109 finite-volume cells have become routine. Thus, large eddy simulation (LES), which attempts to directly resolve the unsteady motion of the largest scales of turbulence, is increasingly of interest for many fluids engineering problems. Additionally, hybrid RANS-LES approaches such as Spalart’s Detached Eddy Simulation (DES), which attempt to marry the best strengths of both methods, are in active use and development. Because it seeks to accurately resolve the unsteady motion of the largest scales of turbulence, LES places certain requirements on the numerical method used to discretize and advance the Navier Stokes equations. In particular, LES requires high time accuracy and minimal intrinsic numerical dispersion and dissipation over a wide range of length scales.
C I Morris↗
Explicit and Implicit Resources: A Simplified Approach to User Requirements Modeling
The International Space Station (ISS) program is dependent on the concept of distributed planning, which places a large responsibility in the hands of the science users. The users assume the responsibility in submitting the operations scheduling requirements in the format required by the scheduling software. The concept of explicit and implicit resources, developed to facilitate the operations scheduling of ISS by non-scheduling experts, is presented. The proposed concepts are defined and examples of scheduling requirements, implemented using explicit and implicit resources, are described together with the requirements modeling and scheduler formatting functions. The software developed to support the implicit and explicit resources concept is discussed.
Hagopian, Jeff↗
An implicit Navier-Stokes code for turbulent flow modeling
This paper presents a numerical approach to calculating turbulent flows employing advanced turbulence models. The main features include a line-by-line Gauss-Seidel algorithm using Roe's approximate Riemann solver, TVD numerical schemes, implicit boundary conditions and a decoupled turbulence-model solver. Based on the problems tested so far, the method has consistently demonstrated its ability in offering accuracy, boundedness and a fast rate of convergence to steady-state solution.
Huang, P. G.↗
Time integration algorithms for the two-dimensional Euler equations on unstructured meshes
Explicit and implicit time integration algorithms for the two-dimensional Euler equations on unstructured grids are presented. Both cell-centered and cell-vertex finite volume upwind schemes utilizing Roe's approximate Riemann solver are developed. For the cell-vertex scheme, a four-stage Runge-Kutta time integration, a fourstage Runge-Kutta time integration with implicit residual averaging, a point Jacobi method, a symmetric point Gauss-Seidel method and two methods utilizing preconditioned sparse matrix solvers are presented. For the cell-centered scheme, a Runge-Kutta scheme, an implicit tridiagonal relaxation scheme modeled after line Gauss-Seidel, a fully implicit lower-upper (LU) decomposition, and a hybrid scheme utilizing both Runge-Kutta and LU methods are presented. A reverse Cuthill-McKee renumbering scheme is employed for the direct solver to decrease CPU time by reducing the fill of the Jacobian matrix. A comparison of the various time integration schemes is made for both first-order and higher order accurate solutions using several mesh sizes, higher order accuracy is achieved by using multidimensional monotone linear reconstruction procedures. The results obtained for a transonic flow over a circular arc suggest that the preconditioned sparse matrix solvers perform better than the other methods as the number of elements in the mesh increases.
Slack, David C.↗
HMI Data Driven Magnetohydrodynamic Model Predicted Active Region Photospheric Heating Rates: Their Scale Invariant, Flare Like Power Law Distributions, and Their Possible Association With Flares
There are many flare forecasting models. For an excellent review and comparison of some of them see Barnes et al. (2016). All these models are successful to some degree, but there is a need for better models. We claim the most successful models explicitly or implicitly base their forecasts on various estimates of components of the photospheric current density J, based on observations of the photospheric magnetic field B. However, none of the models we are aware of compute the complete J. We seek to develop a better model based on computing the complete photospheric J. Initial results from this model are presented in this talk. We present a data driven, near photospheric, 3 D, non-force free magnetohydrodynamic (MHD) model that computes time series of the total J, and associated resistive heating rate in each pixel at the photosphere in the neutral line regions (NLRs) of 14 active regions (ARs). The model is driven by time series of B measured by the Helioseismic & Magnetic Imager (HMI) on the Solar Dynamics Observatory (SDO) satellite. Spurious Doppler periods due to SDO orbital motion are filtered out of the time series of B in every AR pixel. Errors in B due to these periods can be significant.
Goodman, Michael L.↗
Modeling supersonic combustion using a fully-implicit numerical method
A fully-implicit finite-volume algorithm for two-dimensional axisymmetric flows has been coupled to a detailed hydrogen-air reaction mechanism (13 species and 33 reactions) so that supersonic combustion phenomena may be investigated. Numerical computations are compared with ballistic-range shadowgraphs of Lehr (1972) that exhibit two discontinuities caused by a blunt body as it passes through a premixed stoichiometric hydrogen-air mixture. The suitability of the numerical procedure for simulating these double-front flows is shown. The requirements for the physical formulation and the numerical modeling of these flowfields are discussed. Finally, the sensitivity of these external flowfields to changes in certain key reaction rate constants is examined.
Maccormack, Robert W.↗
Extended Magnetohydrodynamics with Embedded Particle-in-Cell Simulation of Ganymede's Magnetosphere
We have recently developed a new modeling capability to embed the implicit particle-in-cell (PIC) model iPIC3D into the Block-Adaptive-Tree-Solarwind-Roe-Upwind-Scheme magnetohydrodynamic (MHD) model. The MHD with embedded PIC domains (MHO-EPIC) algorithm Is a two-way coupled kinetic-fluid model. As one of the very first applications of the MHD-EPIC algorithm, we simulate the Interaction between Jupiter's magnetospherlc plasma and Ganymede's magnetosphere. We compare the MHO-EPIC simulations with pure Hall MHD simulations and compare both model results with Galileo observations to assess the Importance of kinetic effects In controlling the configuration and dynamics of Ganymede's magnetosphere. We find that the Hall MHD and MHO-EPIC solutions are qualitatively similar, but there are significant quantitative differences. In particular. the density and pressure inside the magnetosphere show different distributions. For our baseline grid resolution the PIC solution is more dynamic than the Hall MHD simulation and it compares significantly better with the Galileo magnetic measurements than the Hall MHD solution. The power spectra of the observed and simulated magnetic field fluctuations agree extremely well for the MHD-EPIC model. The MHO-EPIC simulation also produced a few flux transfer events (FTEs) that have magnetic signatures very similar to an observed event. The simulation shows that the FTEs often exhibit complex 3-0 structures with their orientations changing substantially between the equatorial plane and the Galileo trajectory, which explains the magnetic signatures observed during the magnetopause crossings. The computational cost of the MHO-EPIC simulation was only about 4 times more than that of the Hall MHD simulation.
Toth, Gabor↗
Three-dimensional calculations of rotor-airframe interaction in forward flight
A method for analyzing the mutual aerodynamic interaction between a rotor and an airframe model has been developed. This technique models the rotor implicitly through the source terms of the momentum equations. A three-dimensional, incompressible, laminar, Navier-Stokes solver in cylindrical coordinates was developed for analyzing the rotor/airframe problem. The calculations are performed on a simplified model at an advance ratio of 0.1. The airframe surface pressure predictions are found to be in good agreement with wind tunnel test data. Results are presented for velocity and pressure field distributions in the wake of the rotor.
Zori, Laith A. J.↗
Navier-Stokes calculations of rotor-airframe interaction in forward flight
A method for analyzing the mutual aerodynamic interaction between a rotor and an airframe model has been developed. This technique models the rotor implicitly through the source terms of the momentum equations. A three-dimensional, incompressible, laminar, Navier-Stokes solver in cylindrical coordinates was developed for analyzing the rotor-airframe problem. The calculations are performed on a simplified rotor-airframe model at an advance ratio of 0.1. The airframe surface pressure predictions are found to be in good agreement with wind tunnel test data. Results are also presented for velocity and pressure field distributions in the wake of the rotor.
Zori, Laith A. J.↗
The effect of compact implicit differencing in a baroclinic primitive equations model
In an effort to improve forecast accuracy, the horizontal accuracy of gridpoint forecast models at a number of numerical weather prediction (NWP) centers was increased from second order to fourth order. The current GLAS model uses second-order explicit finite-difference formulas in the vertical and fourth-order explicit finite-difference formulas in the horizontal. Experiments are described which indicate the increased forecast accuracy gained by use of compact fourth-order finite differences in a simple baroclinic model. Tables show an increase in forecast accuracy by a factor of 100-150 by use of fourth-order instead of second-order vertical discretization.
Augenbaum, J. M.↗
A tectonic resurfacing model for Venus
Two remarkable aspects of the population of impact craters on Venus are that craters at all sizes are indistinguishable from a random population and that the vast majority of craters have not been significantly modified by tectonic strain or by volcanic flows external to the crater rim, despite evidence from Magellan images that volcanic and tectonic features are widespread on Venus. One interpretation of these observations is that most of the surface dates from the end of a catastrophic global resurfacing event that ceased about 500 My ago, and that the small fraction of craters volcanically embayed or modified by deformation indicates that volcanic and tectonic activity subsequent to that time has been at much lower levels. An alternative model, in which resurfacing occurs episodically in patches a few hundred kilometers in extent and there is a wider spectrum of surface ages, also appears to be consistent with the characteristics of impact craters on Venus. A number of potential mechanisms for catastrophic resurfacing of Venus have been proposed, ranging from geologically sudden convective destabilization of the global lithosphere to strongly time-dependent heat flux and melt generation in the underlying mantle. In most of these geophysical models, resurfacing occurs implicitly or explicitly by volcanism. We explore the hypothesis that, at least in the geologically recent history of Venus, the primary resurfacing mechanism has been tectonic deformation rather than volcanism. We show how such a hypothesis provides at least as good an explanation of a wide range of observations as do volcanic resurfacing models. Finally, we explore the implications of tectonic resurfacing hypothesis for the controversy over the recent resurfacing history of the planet.
Solomon, Sean C.↗
A tectonic resurfacing model for Venus
Two remarkable aspects of the population of impact craters on Venus are as follows: that craters at all sizes are indistinguishable from a random population; and that the vast majority of craters have not been significantly modified by tectonic strain or by volcanic flows external to the crater rim, despite evidence from Magellan images that volcanic and tectonic features are widespread on Venus. One interpretation of these observations is that most of the surface dates from the end of a catastrophic global resurfacing event that ceased about 500 My ago, and that the small fraction of craters volcanically embayed or modified by deformation indicates that volcanic and tectonic activity subsequent to that time has been at much lower levels. An alternative model, in which resurfacing occurs episodically in patches a few hundred kilometers in extent and there is a wider spectrum of surface ages, also appears to be consistent with the characteristics of impact craters on Venus. A number of potential mechanisms for catastrophic resurfacing of Venus have been proposed ranging from geologically sudden convective destabilization of the global lithosphere to strongly time-dependent heat flux and melt generation in the underlying mantle. In most of these geophysical models, resurfacing occurs implicitly or explicitly by volcanism. We explore here the hypothesis that, at least in the geologically recent history of Venus, the primary resurfacing mechanism has been tectonic deformation rather than volcanism. We show how such a hypothesis provides at least as good an explanation of a wide range of observations as do volcanic resurfacing models. Finally, we explore the implications of the tectonic resurfacing hypothesis for the controversy over the recent resurfacing history of the planet.
Solomon, Sean C.↗
Near-Nadiral Normalized Radar Cross Section of the SEA Surface at Ku, Ka, and W-Bands: Comparison of Measurements and Models
Understanding the relationship between wind speed and direction and the near-nadiral normalized radar cross section (NRCS) of the sea surface is important in many oceanographic and atmospheric remote sensing applications: (1) wind speed retrievals in traditional altimeter systems (2) assistance in calibration and path integrated attenuation processing for atmospheric profiling radars The desired wind speed (and direction in some cases) retrieval requires a clear understanding of the relationship between the relevant geophysical quantities and the observed NRCS Such understanding is available from existing electromagnetic models, but the presence of many such models, as well as implicit descriptions of the sea surface, motivates continued evaluation of model performance.
CloudSat↗
Magma oceanography. II - Chemical evolution and crustal formation
A description is presented of an empirical model of fractional crystallization which predicts that slightly modified versions of certain of the proposed whole moon compositions can reproduce the major-element chemistry and mineralogy of most of the primitive highland rocks through equilibrium and fractional crystallization processes combined with accumulation of crystals and trapping of residual liquids. These compositions contain sufficient Al to form a plagioclase-rich crust 60 km thick on top of a magma ocean that was initially no deeper than about 300 km. Implicit in the model are the assumptions that all cooling and crystallization take place at low pressure and that there are no compositional or thermal gradients in the liquid. Discussions of the cooling and crystallization of the proposed magma ocean show these assumptions to be disturbingly naive when applied to the ocean as a whole. However, the model need not be applied to the whole ocean, but only to layers of cooling liquid near the surface.
Longhi, J.↗
Spectral methods for modeling supersonic chemically reacting flow fields
A partial implicit numerical algorithm has been developed for solving the equations describing chemically reacting supersonic flows. The algorithm employs a two-stage Runge-Kutta method for integrating the equations in time and a Chebyshev spectral method for integrating the equations in space. The accuracy and efficiency of the new technique have been assessed by comparison with an existing implicit finite-difference procedure for modeling chemically reacting flows. The comparison showed that the new procedure yielded equivalent accuracy on much coarser grids as compared to the finite-difference procedure with resultant significant gains in computational efficiency.
Drummond, J. P.↗