Determination of radiation interchange factors
Method utilizes an analytical model, a method of solution which is compatible with digital computer analysis, and a generalized computer program to carry out the computations.
Bobco, R. P.↗
SEARCH · Engineering Papers
Results for “factorization methods”
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.
Modeling the Effects of a Backward-Facing Step on Boundary-Layer Transition
We model transition to turbulence in a two-dimensional boundary layer downstream of a backward-facing step (BFS) along a flat plate. With the goal of evaluating the available engineering models for predicting the effects of step excrescences on the transition characteristics, two separate methodologies are used to monitor the streamwise shift in the transition onset location as the step height and the flow speed are varied across the range of a previously reported experiment involving step-height-to-local-displacement-thickness ratios of 0 < h/δ* < 1.6. Unlike the variable N -factor method from the previous literature, both of these methods are general in scope and do not involve any empirical correlations that are specific to step excrescences. The first of these techniques involves an N -factor method that directly accounts for the change in boundary-layer instability characteristics due to the step. Stability computations using the harmonic linearized Navier-Stokes equations (HLNSE), which fully account for the nonparallel-mean-flow effects close to the BFS, indicate that the measured transition locations at nearly all test conditions ( h/δ* < 1.3) correlate well with a computed N -factor of N tr = 7.6, demonstrating a successful stability-based transition criterion related to step excrescences. Linear stability theory, which does not account for nonparallel effects, demonstrates reasonable agreement with the HLNSE results, yielding good predictions for the overall trends, but predicts a somewhat earlier onset of transition than HLNSE. The other methodology used in this work involves transport-equation-based transition models. We first show that the Langtry-Menter y - Re θt transition model cannot accurately predict the location of transition onset for moderate BFS heights because it is unable to accurately account for the flow history effects. Along with the Langtry-Menter transition model, we also show the amplification factor transport model does not produce accurate transition locations for subsonic flow over steps even though it accounts for some flow history effects.
Transition↗
Modeling the Effects of a Backward-Facing Step on Boundary-Layer Transition
We model transition to turbulence in a two-dimensional boundary layer downstream of a backward-facing step (BFS) along a flat plate. With the goal of evaluating the available engineering models for predicting the effects of step excrescences on the transition characteristics, two separate methodologies are used to monitor the streamwise shift in the transition onset location as the step height and the flow speed are varied across the range of a previously reported experiment involving step-height-to-local-displacement-thickness ratios of 0 < h/δ* < 1.6. Unlike the variable N -factor method from the previous literature, both of these methods are general in scope and do not involve any empirical correlations that are specific to step excrescences. The first of these techniques involves an N -factor method that directly accounts for the change in boundary-layer instability characteristics due to the step. Stability computations using the harmonic linearized Navier-Stokes equations (HLNSE), which fully account for the nonparallel-mean-flow effects close to the BFS, indicate that the measured transition locations at nearly all test conditions ( h/δ* < 1.3) correlate well with a computed N -factor of N tr = 7.6, demonstrating a successful stability-based transition criterion related to step excrescences. Linear stability theory, which does not account for nonparallel effects, demonstrates reasonable agreement with the HLNSE results, yielding good predictions for the overall trends, but predicts a somewhat earlier onset of transition than HLNSE. The other methodology used in this work involves transport-equation-based transition models. We first show that the Langtry-Menter y - Re θt transition model cannot accurately predict the location of transition onset for moderate BFS heights because it is unable to accurately account for the flow history effects. Along with the Langtry-Menter transition model, we also show the amplification factor transport model does not produce accurate transition locations for subsonic flow over steps even though it accounts for some flow history effects.
Transition↗
Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations
Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework rely upon transition correlations that were derived from a limited number of low-speed experiments. Furthermore, these models often account for only a subset of the relevant transition mechanisms and/or cannot accurately predict the sensitivity of those mechanisms to the changes in significant flow parameters. A preceding investigation had targeted the assessment of the transport-equation-based transition models in NASA's OVERFLOW 2.3b solver, namely, the amplification factor transport (AFT-2017b) equation model coupled with the Spalart-Allmaras RANS model and the Langtry-Menter transition models (LM2009 without crossflow effects and LM2015 including the modeling of crossflow transition) implemented with Menter’s shear-stress transport equation (SST2003) RANS model. Comparisons with recent measurements at transonic freestream conditions on the Common Research Model with Natural Laminar Flow (CRM-NLF) reinforced our earlier finding that all three of the above models significantly underpredict the reported extent of the laminar flow region over the entire span of the wing, regardless of the dominant instability mechanism(s) underlying the onset of the transition process. The underprediction of the laminar flow extent was attributed to the failure of the above models in accounting for the stabilizing effect of compressibility on the amplification of Tollmien-Schlichting instabilities. Based on previous linear stability studies related to compressibility effects, the present work proposes modifications to the two classes of transition models that reduce to the original form of each model at low subsonic speeds and do not require any nonlocal flow information or additional transport equation(s). The modifications are shown to significantly improve the predicted laminar extent of the flow and compare well against the data from the CRM-NLF experiment. Additionally, a previous assessment of transition prediction based on the dual, nonparallel N -factor method in conjunction with linear parabolized stability equations (PSE) is extended to additional angles of attack to provide the first comprehensive assessment of transition models based on nonparallel disturbance amplification over the CRM-NLF. In general, the transition criterion based on the dual, nonparallel N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at R eMAC = 15 million for all angles of attack investigated herein and provides additional validation of the improved predictions from the compressibility-corrected transition models.
CFD modeling↗
Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations
Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework rely upon transition correlations that were derived from a limited number of low-speed experiments. Furthermore, these models often account for only a subset of the relevant transition mechanisms and/or cannot accurately predict the sensitivity of those mechanisms to the changes in significant flow parameters. A preceding investigation had targeted the assessment of the transport-equation-based transition models in NASA's OVERFLOW 2.3b solver, namely, the amplification factor transport (AFT-2017b) equation model coupled with the Spalart-Allmaras RANS model and the Langtry-Menter transition models (LM2009 without crossflow effects and LM2015 including the modeling of crossflow transition) implemented with Menter’s shear-stress transport equation (SST2003) RANS model. Comparisons with recent measurements at transonic freestream conditions on the Common Research Model with Natural Laminar Flow (CRM-NLF) reinforced our earlier finding that all three of the above models significantly underpredict the reported extent of the laminar flow region over the entire span of the wing, regardless of the dominant instability mechanism(s) underlying the onset of the transition process. The underprediction of the laminar flow extent was attributed to the failure of the above models in accounting for the stabilizing effect of compressibility on the amplification of Tollmien-Schlichting instabilities. Based on previous linear stability studies related to compressibility effects, the present work proposes modifications to the two classes of transition models that reduce to the original form of each model at low subsonic speeds and do not require any nonlocal flow information or additional transport equation(s). The modifications are shown to significantly improve the predicted laminar extent of the flow and compare well against the data from the CRM-NLF experiment. Additionally, a previous assessment of transition prediction based on the dual, nonparallel N -factor method in conjunction with linear parabolized stability equations (PSE) is extended to additional angles of attack to provide the first comprehensive assessment of transition models based on nonparallel disturbance amplification over the CRM-NLF. In general, the transition criterion based on the dual, nonparallel N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at R eMAC = 15 million for all angles of attack investigated herein and provides additional validation of the improved predictions from the compressibility-corrected transition models.
CFD modeling↗
Eigenvalue Sensitivity Computations for Linear Stability Theory
To realize the drag reduction benefit of boundary-layer transition control strategies, it is crucial to integrate transition prediction into the vehicle design through an optimization process. The integration of transition prediction based on linear stability analysis into adjoint d design optimization requires coupling an adjoint enabled computational fluid dynamics (CFD) solver with an adjoint enabled linear stability code. In particular, the boundary-layer transition location is often predicted using the N-factor method based on linear stability theory (LST). Thus, sensitivity of the linear-stability eigenvalues constitute an essential building block for optimizing the laminar flow performance. The present paper describes an implementation of LST eigenvalue sensitivity analysis that can be easily coupled with a CFD solver. Specifically, we describe a discrete adjoint formulation for the transition location prediction based on the N-factor method. The verification of this formulation is carried out by comparing the adjoint-based sensitivity of the local growth rate of a given instability mode with respect to the disturbance frequency, and the adjoint-based sensitivity of the transition location with respect to spanwise wavenumber with those sensitivities computed using a finite-difference approximation. Finally, the adjoint LST formulation is applied to flat-plate boundary-layer flows at transonic, supersonic, and hypersonic conditions, to determine the behavior and sensitivities of the transition location with respect to a range of disturbance spanwise wavenumbers.
Boundary Layer Transition↗
Factorization and reduction methods for optimal control of distributed parameter systems
A Chandrasekhar-type factorization method is applied to the linear-quadratic optimal control problem for distributed parameter systems. An aeroelastic control problem is used as a model example to demonstrate that if computationally efficient algorithms, such as those of Chandrasekhar-type, are combined with the special structure often available to a particular problem, then an abstract approximation theory developed for distributed parameter control theory becomes a viable method of solution. A numerical scheme based on averaging approximations is applied to hereditary control problems. Numerical examples are given.
Burns, J. A.↗
Radiative transfer in rectangular enclosures - A discretized exchange factor solution
The discretized exchange factor method is used to analyze radiative exchange in a rectangular enclosure. The results compare excellently with those of other methods, especially the zonal method. Since the direct exchange factors are between nodal points no integration is necessary for evaluation of these factors. It is found that the present approach provides accurate results even when a small number of nodes is used.
Naraghi, M. H. N.↗
Explicit and implicit solution of the Navier-Stokes equations on a massively parallel computer
The design, implementation, and performance of a two-dimensional time-accurate Navier-Stokes solver for the CM2 supercomputer are described. The program uses a single processor for each grid point. Two different time-stepping methods have so far been implemented: an explicit third-order Runge-Kutta method and an implicit approximation-factorization method. The CM2 results are checked against those of a mature well-vectorized Cray 2 program, both for correctness and performance. The code is found to be correct, and the performance in some cases is up to several times that of the Cray 2.
Levit, Creon↗
Extension of the discrete generalized multigroup method using SPH factors
The discrete generalized multigroup (DGM) method provides a way to treat the energy dependence of neutron transport similarly to the standard multigroup approximation. However, DGM uses an orthogonal basis to retain the energy dependence in higher-order terms. Using correction factors similar to traditional Superhomogénéisation (SPH) factors, the DGM method may be extended to produce cross sections that are homogenized over both space and energy. Additionally, since some fine-group energy dependence is retained, the resulting homogenized cross sections are more problem-independent than cross sections homogenized by SPH factors alone. In particular, a 44-group set of cross sections is collapsed to approximately a 1% error in the pincell fission densities for a test problem using three DOF per coarse energy group.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Estimating Eigenenergies from Quantum Dynamics: A Unified Noise-Resilient Measurement-Driven Approach
Ground state energy estimation in physical, chemical, and materials sciences is one of the most promising applications of quantum computing. In this work, we introduce a new hybrid approach that finds the eigenenergies by collecting real-time measurements and post-processing them using the machinery of dynamic mode decomposition (DMD). From the perspective of quantum dynamics, we establish that our approach can be formally understood as a stable variational method on the function space of observables available from a quantum many-body system. We also provide strong theoretical and numerical evidence that our method converges rapidly even in the presence of a large degree of perturbative noise, and show that the method bears an isomorphism to robust matrix factorization methods developed independently across various scientific communities. Our numerical benchmarks on spin and molecular systems demonstrate an accelerated convergence and a favorable resource reduction over state-of-the-art algorithms. The DMD-centric strategy can systematically mitigate noise and stands out as a leading hybrid quantum-classical eigensolver.
Shen, Yizhi↗
Multigrid acceleration of the flux split Euler equations
Multigrid acceleration is applied to a flux-split algorithm for solving the Euler equations in two and three dimensions. The basic algorithm is an implicit spatially-split approximate factorization method. The stability of the scheme in comparison to other factorization is examined. Results are presented for two-dimensional airfoil flows and three-dimensional wing flows which demonstrate substantially improved convergence with the multigrid algorithm. An asymptotic spectral radius of 0.89 and 0.93 is attained for a 97 x 17 x 17 wing solution at subcritical and supercritical conditions, respectively.
Anderson, W. K.↗
Implicit flux-split schemes for the Euler equations
Recent progress in the development of implicit algorithms for the Euler equations using the flux-vector splitting method is described. Comparisons of the relative efficiency of relaxation and spatially-split approximately factored methods on a vector processor for two-dimensional flows are made. For transonic flows, the higher convergence rate per iteration of the Gauss-Seidel relaxation algorithms, which are only partially vectorizable, is amply compensated for by the faster computational rate per iteration of the approximately factored algorithm. For supersonic flows, the fully-upwind line-relaxation method is more efficient since the numerical domain of dependence is more closely matched to the physical domain of dependence. A hybrid three-dimensional algorithm using relaxation in one coordinate direction and approximate factorization in the cross-flow plane is developed and applied to a forebody shape at supersonic speeds and a swept, tapered wing at transonic speeds.
Thomas, J. L.↗
Verification of redundancy management design
Statistical method checks designs by simulating system operating conditions and adding error factors. Method has potential applicability to commercial and industrial situations where redundancy management system is used to detect and isolate failed components.
Gelderloos, H. C.↗
Projection methods for the numerical solution of Markov chain models
Projection methods for computing stationary probability distributions for Markov chain models are presented. A general projection method is a method which seeks an approximation from a subspace of small dimension to the original problem. Thus, the original matrix problem of size N is approximated by one of dimension m, typically much smaller than N. A particularly successful class of methods based on this principle is that of Krylov subspace methods which utilize subspaces of the form span(v,av,...,A(exp m-1)v). These methods are effective in solving linear systems and eigenvalue problems (Lanczos, Arnoldi,...) as well as nonlinear equations. They can be combined with more traditional iterative methods such as successive overrelaxation, symmetric successive overrelaxation, or with incomplete factorization methods to enhance convergence.
Saad, Youcef↗
Aerostructural safety factor criteria
The present modification of the conventional safety factor method for aircraft structures evaluation involves the expression of deterministic safety factors in probabilistic tolerance limit ratios; these are found to involve a total of three factors that control the interference of applied and resistive stress distributions. The deterministic expression is extended so that it may furnish a 'relative ultimate safety' index that encompasses all three distribution factors. Operational reliability is developed on the basis of the applied and the yield stress distribution interferences. Industry standards are suggested to be derivable from factor selections that are based on the consequences of failure.
Verderaime, V.↗
Temperature and pressure effects on capacitance probe cryogenic liquid level measurement accuracy
The inaccuracies of liquid nitrogen and liquid hydrogen level measurements by use of a coaxial capacitance probe were investigated as a function of fluid temperatures and pressures. Significant liquid level measurement errors were found to occur due to the changes in the fluids dielectric constants which develop over the operating temperature and pressure ranges of the cryogenic storage tanks. The level measurement inaccuracies can be reduced by using fluid dielectric correction factors based on measured fluid temperatures and pressures. The errors in the corrected liquid level measurements were estimated based on the reported calibration errors of the temperature and pressure measurement systems. Experimental liquid nitrogen (LN2) and liquid hydrogen (LH2) level measurements were obtained using the calibrated capacitance probe equations and also by the dielectric constant correction factor method. The liquid levels obtained by the capacitance probe for the two methods were compared with the liquid level estimated from the fluid temperature profiles. Results show that the dielectric constant corrected liquid levels agreed within 0.5 percent of the temperature profile estimated liquid level. The uncorrected dielectric constant capacitance liquid level measurements deviated from the temperature profile level by more than 5 percent. This paper identifies the magnitude of liquid level measurement error that can occur for LN2 and LH2 fluids due to temperature and pressure effects on the dielectric constants over the tank storage conditions from 5 to 40 psia. A method of reducing the level measurement errors by using dielectric constant correction factors based on fluid temperature and pressure measurements is derived. The improved accuracy by use of the correction factors is experimentally verified by comparing liquid levels derived from fluid temperature profiles.
Edwards, Lawrence G.↗
Stability theory and transition prediction applied to a general aviation fuselage
The linear stability of a fully three-dimensional boundary layer formed over a general aviation fuselage was investigated. The location of the onset of transition was estimated using the N-factor method. The results were compared with existing experimental data and indicate N-factors of approximately 8.5 on the side of the fuselage and 3.0 near the top. Considerable crossflow existed along the side of the body, which significantly affected the unstable modes present in the boundary layer. Fair agreement was found between the predicted frequency range of linear instability modes and available experimental data concerning the spectral content of the boundary layer.
Spall, R. E.↗