Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CALCULATION 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.

At least 793 records · Page 44

Polynomial interpolation methods for viscous flow calculations

Higher-order collocation procedures which result in block-tridiagonal matrix systems are derived from (1) Taylor series expansions and from (2) polynomial interpolation, and the relationships between the two formulations, called respectively Hermite and spline collocation, are investigated. A Hermite block-tridiagonal system for a nonuniform mesh is derived, and the Hermite approach is extended in order to develop a variable-mesh sixth-order block-tridiagonal procedure. It is shown that all results obtained by Hermite development can be recovered by appropriate spline polynomial interpolation. The additional boundary conditions required for these higher-order procedures are also given. Comparative solutions using second-order accurate finite difference and spline and Hermite formulations are presented for the boundary layer on a flat plate, boundary layers with uniform and variable mass transfer, and the viscous incompressible Navier-Stokes equations describing flow in a driven cavity.

Rubin, S. G.↗

Upgraded viscous flow analysis of multielement airfoils

A description of an improved version of the NASA/Lockheed multielement airfoil analysis computer program is presented. The improvements include several major modifications of the aerodynamic model as well as substantial changes of the computer code. The modifications of the aerodynamic model comprise the representation of the boundary layer and wake displacement effects with an equivalent source distribution, the prediction of wake parameters with Green's lag-entrainment method, the calculation of turbulent boundary layer separation with the method of Nash and Hicks, the estimation of the onset of confluent boundary layer separation with a modified form of Goradia's method, and the prediction of profile drag with the formula of Squire and Young. The paper further describes the modifications of the computer program for which the structured approach to computer software development was employed. Important aspects of the structured program development such as the functional decomposition of the aerodynamic theory and its numerical implementation, the analysis of the data flow within the code, and the application of a pseudo code are discussed.

Brune, G. W.↗

Upgraded viscous flow analysis of multi-element airfoils

A description of an improved version of the NASA/Lockheed multi-element airfoil analysis computer program is presented. The improvements include several major modifications of the aerodynamic model as well as substantial changes of the computer code. The modifications of the aerodynamic model comprise the representation of the boundary layer and wake displacement effects with an equivalent source distribution, the prediction of wake parameters with Green's lag-entrainment method, the calculation of turbulent boundary layer separation with the method of Nash and Hicks, the estimation of the onset of confluent boundary layer separation with a modified form of Goradia's method, and the prediction of profile drag with the formula of Squire and Young. The modifications of the computer program for which the structured approach to computer software development was employed are also described. Important aspects of the structured program development such as the functional decomposition of the aerodynamic theory and its numerical implementation, the analysis of the data flow within the code, and the application of a pseudo code are discussed.

Brune, G. W.↗

Efficient computation of volume in flow predictions

An efficient method for calculating cell volumes for time-dependent three-dimensional flow predictions by finite volume calculations is presented. Eight arbitrary corner points are considered and the shape face is divided into two planar triangles. The volume is then dependent on the orientation of the partitioning. In the case of a hexahedron, it is noted that any open surface with a boundary that is a closed curve possesses a surface vector independent of the surface shape. Expressions are defined for the surface vector, which is independent of the partitioning surface diagonal used to quantify the volume. Using a decomposition of the cell volume involving two corners, with each the vertex of three diagonals and six corners which are vertices of one diagonal, gives portions which are tetrahedra. The resultant mesh is can be used for time-dependent finite volume calculations one requires less computer time than previous methods.

Vinokur, M.↗

Compressible flow calculations employing the Galerkin/least-squares method

A multielement group, domain decomposition algorithm is presented for solving linear nonsymmetric systems arising in the finite-element analysis of compressible flows employing the Galerkin/least-squares method. The iterative strategy employed is based on the generalized minimum residual (GMRES) procedure originally proposed by Saad and Shultz. Two levels of preconditioning are investigated. Applications to problems of high-speed compressible flow illustrate the effectiveness of the scheme.

Shakib, F.↗

A method of quadrature for calculation of the laminar and turbulent boundary layer in case of plane and rotationally symmetrical flow

For calculation of the characteristic parameters of the boundary layer (momentum-loss thickness and form parameter for the velocity profile), two quadrature formulas are given which are valid for the laminar as well as for the turbulent state of flow. These formulas cover both the two-dimensional and the rotationally symmetrical case. The calculation of the momentum-loss thickness is carried out by a simple integration of the energy theorem. The equation for the form parameter is obtained by coupling of the momentum theorem with the energy theorem. Knowledge of the derivatives of the velocity distribution and of the radius of the body along the length x is not necessary.

FLOW, TURBULENT↗

Analyses of Fatigue Crack Growth and Closure Near Threshold Conditions for Large-Crack Behavior

A plasticity-induced crack-closure model was used to study fatigue crack growth and closure in thin 2024-T3 aluminum alloy under constant-R and constant-K(sub max) threshold testing procedures. Two methods of calculating crack-opening stresses were compared. One method was based on a contact-K analyses and the other on crack-opening-displacement (COD) analyses. These methods gave nearly identical results under constant-amplitude loading, but under threshold simulations the contact-K analyses gave lower opening stresses than the contact COD method. Crack-growth predictions tend to support the use of contact-K analyses. Crack-growth simulations showed that remote closure can cause a rapid rise in opening stresses in the near threshold regime for low-constraint and high applied stress levels. Under low applied stress levels and high constraint, a rise in opening stresses was not observed near threshold conditions. But crack-tip-opening displacement (CTOD) were of the order of measured oxide thicknesses in the 2024 alloy under constant-R simulations. In contrast, under constant-K(sub max) testing the CTOD near threshold conditions were an order-of-magnitude larger than measured oxide thicknesses. Residual-plastic deformations under both constant-R and constant-K(sub max) threshold simulations were several times larger than the expected oxide thicknesses. Thus, residual-plastic deformations, in addition to oxide and roughness, play an integral part in threshold development.

Newman, J. C., Jr.↗

Efficient Approximation for Structural Optimization Under Multiple Constraints

The cooperative agreement covered work between August 1995 and August 1997. The focus of the work was efficient approximations of structural response and sensitivity. The effort proceeded in three directions as follows: (1) Development of an approximation extended to efficient sensitivity approximations and demonstrated for structural models for the High Speed Civil Transport; (2) Preliminary development of the adjoint method for calculating sensitivity derivatives; and (3) A review of method for fast exact reanalysis. Attachments of papers which were submitted during this period are included.

Haftka, Raphael T.↗

Finite Volume Numerical Methods for Aeroheating Rate Calculations from Infrared Thermographic Data

The use of multi-dimensional finite volume numerical techniques with finite thickness models for calculating aeroheating rates from measured global surface temperatures on hypersonic wind tunnel models was investigated. Both direct and inverse finite volume techniques were investigated and compared with the one-dimensional semi -infinite technique. Global transient surface temperatures were measured using an infrared thermographic technique on a 0.333-scale model of the Hyper-X forebody in the Langley Research Center 20-Inch Mach 6 Air tunnel. In these tests the effectiveness of vortices generated via gas injection for initiating hypersonic transition on the Hyper-X forebody were investigated. An array of streamwise orientated heating striations were generated and visualized downstream of the gas injection sites. In regions without significant spatial temperature gradients, one-dimensional techniques provided accurate aeroheating rates. In regions with sharp temperature gradients due to the striation patterns two-dimensional heat transfer techniques were necessary to obtain accurate heating rates. The use of the one-dimensional technique resulted in differences of 20% in the calculated heating rates because it did not account for lateral heat conduction in the model.

Daryabeigi, Kamran↗

Improved Optimization Method for Array Feed Calculations

Reflector antenna optimization schemes using array feeds have been used to recover antenna losses resulting from antenna distortions and aberrations and to generate contour coverage patterns.

focal plane array feed antennas gain optimization ↗

Finite Volume Numerical Methods for Aeroheating Rate Calculations from Infrared Thermographic Data

The use of multi-dimensional finite volume heat conduction techniques for calculating aeroheating rates from measured global surface temperatures on hypersonic wind tunnel models was investigated. Both direct and inverse finite volume techniques were investigated and compared with the standard one-dimensional semi-infinite technique. Global transient surface temperatures were measured using an infrared thermographic technique on a 0.333-scale model of the Hyper-X forebody in the NASA Langley Research Center 20-Inch Mach 6 Air tunnel. In these tests the effectiveness of vortices generated via gas injection for initiating hypersonic transition on the Hyper-X forebody was investigated. An array of streamwise-orientated heating striations was generated and visualized downstream of the gas injection sites. In regions without significant spatial temperature gradients, one-dimensional techniques provided accurate aeroheating rates. In regions with sharp temperature gradients caused by striation patterns multi-dimensional heat transfer techniques were necessary to obtain more accurate heating rates. The use of the one-dimensional technique resulted in differences of 20% in the calculated heating rates compared to 2-D analysis because it did not account for lateral heat conduction in the model.

Daryabeigi, Kamran↗

J. E. Nichols and D. M. Peteet reply

The quantification of northern peatland carbon stocks is vital to our understanding of the global carbon cycle, and we welcome the opportunity to respond to the comments submitted in the preceding “Matters Arising” (Yu et al., Ratcliffe et al.). In response to Yu et al., we provide additional clarification of how data from the Neotoma Paleoecological Database (NPD) are stored and used. We also repeated our quantification of northern peat carbon stocks using only radiocarbon data vetted by other publications. However, our conclusions are unchanged by the elimination of the additional data. We respond to Ratcliffe et al. by addressing the conflict between the “time-history” and “inventory methods.” We calculated northern peatland carbon stocks by the inventory method using our dataset and carefully quantify the uncertainties. We find that while the mean values for carbon stock differ by twofold, the uncertainties around those means are larger than previously acknowledged and are overlapping. This result emphasizes the need for continued work by the carbon cycle science community to understand and quantify uncertainty in terrestrial carbon stock estimates.

Peatlands↗

Recent Developments in Tolerancing Methods for Imaging Spectrometers

We discuss detailed tolerancing methods developed for imaging spectrometers at NASA Jet Propulsion Laboratory, California Institute of Technology using the Earth Surface Mineral Dust Source Investigation (EMIT) imaging spectrometer as an illustrative example. We tolerance five metrics simultaneously: along-track response function, cross-track response function, spectral response function, spectral centroid uniformity, and spatial centroid uniformity. A method to calculate tolerancing sensitivities for each metric directly, a method to statistically combine Monte Carlo files from multiple tolerancing runs, and an example summary error budgets that communicate the key and driving tolerances for each metric are discussed. These methods facilitate rapid and semi-automated assessment of the predicted performance of imaging spectrometer systems from design through to assembly and launch life cycle, using metrics that are directly relevant to the extraction of accurate spectroscopic data from these instruments.

Mouroulis, Pantazis↗

Thermal Management System Modeling in the Heat Transport System Simulation (HeaTSSPy) Package

This paper describes the development of a thermal management system (TMS) concept design and analysis software package called Heat Transport System Simulation (HeaTSSPy). Built within Python using the OpenMDAO framework, HeaTSSPy can be used to size and optimize an active (using liquid/air heat exchangers) or passive (using finned heat sinks) TMS. The package makes use of modular TMS elements that allow for the creation of different system architectures and includes components such as heat sinks, heat exchangers, liquid pumps, fans, ducts, air inlets, air nozzles, and liquid pipes. Modeling methods for these components include a combination of physics-based analytical and empirical equations that relate component sizing criteria to system performance. The HeatSSPy heat sink methods are fully detailed within this paper, while the heat exchanger methods are described and referenced from previous work. This paper also uses high-fidelity simulation to validate two different methods of calculating thermal resistance using CFD results. Once the methods are fully described, the code is exercised to compare an active TMS with that of a passive TMS. Design criteria for this study include rejected heat, system altitude, Mach number, and ambient temperature. These criteria are used to develop a TMS system with estimated performance metrics such as weight, drag, and operational power. Results of this paper show the crossover point when a passive system begins to weigh more than an active system for a given heat rejection.

Thermal management system↗

Sensitivity analysis of discrete structural systems: A survey

Methods for calculating sensitivity derivatives for discrete structural systems are surveyed, primarily covering literature published during the past two decades. Methods are described for calculating derivatives of static displacements and stresses, eigenvalues and eigenvectors, transient structural response, and derivatives of optimum structural designs with respect to problem parameters. The survey is focused on publications addressed to structural analysis, but also includes a number of methods developed in nonstructural fields such as electronics, controls, and physical chemistry which are directly applicable to structural problems. Most notable among the nonstructural-based methods are the adjoint variable technique from control theory, and the Green's function and FAST methods from physical chemistry.

Adelman, H. M.↗

Sensitivity analysis of discrete structural systems

Methods for calculating sensitivity derivatives for discrete structural systems are surveyed, primarily covering literature published during the past two decades. Methods are described for calculating derivatives of static displacements and stresses, eigenvalues and eigenvectors, transient structural response, and derivatives of optimum structural designs with respect to problem parameters. The survey is focused on publications developed in nonstructural fields such as electronics, controls, and physical chemistry which are directly applicable to structural problems. Most notable among the nonstructural-based methods are the adjoint variable technique from control theory, and the Green's function and FAST methods from physical chemistry.

Adelman, H. M.↗