SEARCH · Engineering Papers

Results for “FLOW - CONFORMAL TRANSFORMATION - HODOGRAPH METHOD”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

NASA NTRS2007

Anisotropic Solution Adaptive Unstructured Grid Generation Using AFLR

An existing volume grid generation procedure, AFLR3, was successfully modified to generate anisotropic tetrahedral elements using a directional metric transformation defined at source nodes. The procedure can be coupled with a solver and an error estimator as part of an overall anisotropic solution adaptation methodology. It is suitable for use with an error estimator based on an adjoint, optimization, sensitivity derivative, or related approach. This offers many advantages, including more efficient point placement along with robust and efficient error estimation. It also serves as a framework for true grid optimization wherein error estimation and computational resources can be used as cost functions to determine the optimal point distribution. Within AFLR3 the metric transformation is implemented using a set of transformation vectors and associated aspect ratios. The modified overall procedure is presented along with details of the anisotropic transformation implementation. Multiple two-and three-dimensional examples are also presented that demonstrate the capability of the modified AFLR procedure to generate anisotropic elements using a set of source nodes with anisotropic transformation metrics. The example cases presented use moderate levels of anisotropy and result in usable element quality. Future testing with various flow solvers and methods for obtaining transformation metric information is needed to determine practical limits and evaluate the efficacy of the overall approach.

Marcum, David L.↗

DOE OSTI2025

Nonadiabatic Force Matching for Alchemical Free-Energy Estimation

We propose a method to compute free-energy differences from nonadiabatic alchemical transformations by using flow-based generative models. The method, nonadiabatic force matching, hinges on estimating the dissipation along an alchemical switching process in terms of a nonadiabatic force field that can be learned through stochastic flow matching. The learned field can be used in conjunction with short-time trajectory data to evaluate upper and lower bounds on the alchemical free energy that variationally converge to the exact value if the field is optimal. Applying the method to evaluate the alchemical free energy of atomistic models shows that it can substantially reduce the simulation cost of a free-energy estimate at a negligible loss of accuracy when compared with thermodynamic integration.

Computational chemistry↗

NASA NTRS1982

Analysis and design of steady transonic flow over airfoils by the method of parametric differentiation

The results reported here are based on applying the method of parametric differentiation (MPD) to transform the nonlinear differential equation governing small-disturbance transonic flow to a linear equation. Implicit approximate factorization and monotone methods were used to accelerate convergence of the linear problem by an order of magnitude over successive line over-relaxation. The relative merits of using MPD are discussed in comparison to conventional small-disturbance applications. Several MPD analyses are performed on an array of airfoils. A design procedure utilizing MPD is discussed and demonstrated for two nonlifting cases.

Halt, D. W.↗

NASA NTRS1980

Numerical solutions of transonic flows by parametric differentiation and integral equation techniques

The paper discusses results of an exploratory study of the advantages obtained by combining integral equation methods with the method of parametric differentiation in the treatment of transonic flow problems. In the proposed method, the nonlinear unsteady transonic flow equation for small perturbations is transformed into a linear equation by the use of the method of parametric differentiation. The linear equation is split into a pair of weakly coupled partial differential equations by writing the transformed perturbation potential as the sum of a steady component and an unsteady component. The solution of the steady equation as an integral equation is based on Ogana's treatment (1978). As a test case, the formulation is applied to predict the steady transonic flow over a nonlifting parabolic-arc airfoil.

Sivaneri, N. T.↗

NASA NTRS1991

Growth Of Instabilities In Two Types Of Mixing Layers

Report discusses issues involved in computation of growth of instabilities and vorticity in mixing (shear) layer between two parallel flows having different speeds and separated initially by flat surface. Study focusses upon degree to which these growing flow structures hasten or retard their own further growth, and upon related question of differences between flows computed according to mathematical models of S-layer and T-layer mixing. Transformation method used offers suitable tool for study of such flows.

Kaul, Upender K.↗

NASA NTRS1979

A two dimensional model of viscous slot flow

The effects of viscosity on flows through slots are investigated. Using Briley's formulation of the Alternating Direction Implicit method together with the Thompson transformation, numerical solution of viscous incompressible flow through a slot is simulated, for Reynolds number 1000 and openness ratio one tenth. The problem is formulated in terms of stream function and vorticity. A wall iteration of Israeli type is employed to iteratively induce the no slip condition on velocity.

Cooke, C. H.↗

NASA NTRS1985

Unsteady Heat Transfer Due to Time-dependent Free Stream Velocity

The present work was undertaken in order to study the unsteady combined convection from a horizontal circular cylinder to a transverse flow. A coordinate perturbation method is used to transform the governing set of partial differential equations into a system of ordinary differential equations. The free stream time-dependent velocity was assumed to be sinusoidal and the boundary layer response due to both low as well as high frequencies of oscillation will be studied. Currently numerical solutions are being obtained for the distribution of the unsteady Nusselt number and the friction factor.

Gorla, R. S. R.↗

NASA NTRS1977

Unsteady linearized transonic flow analysis for slender bodies

An unsteady linearized formulation based on Oswatitsch-Keune's parabolic method is developed to analyze transonic flow past oscillating slender bodies. In contrast to the widely used integral transform method, it is shown that all solutions can be derived by a simpler method directly in the physical plane. By various expansion procedures, low-frequency solutions then are derived according to two clearly defined frequency ranges. Adams-Sears' iteration is employed to account for the second-order effects. Stability derivatives are compared with available theories and data. It is found that the derivatives depend more sensitively on thickness than on the reduced frequency. Finally, a critical assessment of the present method is given.

Liu, D. D.↗

NASA NTRS1993

Boundary conditions and the simulation of low Mach number flows

The problem of accurately computing low Mach number flows, with the specific intent of studying the interaction of sound waves with incompressible flow structures, such as concentrations of vorticity is considered. This is a multiple time (and/or space) scales problem, leading to various difficulties in the design of numerical methods. Concentration is on one of these difficulties - the development of boundary conditions at artificial boundaries which allow sound waves and vortices to radiate to the far field. Nonlinear model equations are derived based on assumptions about the scaling of the variables. Then these are linearized about a uniform flow and exact boundary conditions are systematically derived using transform methods. Finally, useful approximations to the exact conditions which are valid for small Mach number and small viscosity are computed.

Hagstrom, Thomas↗

NASA NTRS2019

Characterization of Thermal Protection Systems: An Analysis of Optical & Thermal Properties

Thermal Protection System Materials & thermal control coatings were characterized by a variety of methodologies including heat flow meter method (HFM), dynamic scanning calorimetry (DSC), fourier transform infrared spectroscopy (FTIR), and Ultra-Violet & Visible light spectroscopy (UV-VIS). Physical properties were measured to aid in the construction of thermal response models for a number a projects including MEDLI and MEDLI2. Comparing pre-flight predictions with actual flight data can allow for a reduction of margins and improve designs. The models enable design of TPS for flight missions and analysis of flight data to understand aerothermal environments experienced on the mission. This presentation reveals critical information: heat capacity (cp vs T), thermal conductivity ( vs T), and emissivity (E vs T) for Super Lightweight Ablative (SLA), Phenolic Impregnated Carbon Ablative (PICA-D), & Heatshield for Extreme Entry Environment Technology (HEEET). These properties were also measured for various thermal control coatings manufactured by AZ Technology. A portion of this work is incomplete and highlights questions to investigate in the future.

Cole, Brian↗

NASA NTRS1977

Application of the method of integral relations to laminar boundary layers in three dimensions

The method of integral relations is extended to general three-dimensional compressible laminar boundary layer flows. The transformation employed to transform the basic three-dimensional compressible boundary layer equations into quasi-incompressible form is an extension of the Howarth transformation. The resulting system of differential equations is integrated numerically by the method of integral relations as proposed by Dorodnitsyn. To demonstrate the accuracy of the method, it is applied to calculation of the parabolic flow over a flat plate and the boundary flow over an infinite yawed cylinder, for which solutions are known. It is then applied to the flow over a flat plate disturbed by a cylinder normal to the plate, for which a finite-difference solution is available for comparison. It is finally applied to calculating the crossflow velocity variation for supersonic flow over a swept wedge.

Holt, M.↗

NASA NTRS1987

An engineering study of hybrid adaptation of wind tunnel walls for three-dimensional testing

Solid wall tunnels having only upper and lower walls flexing are described. An algorithm for selecting the wall contours for both 2 and 3 dimensional wall flexure is presented and numerical experiments are used to validate its applicability to the general test case of 3 dimensional lifting aircraft models in rectangular cross section wind tunnels. The method requires an initial approximate representation of the model flow field at a given lift with wallls absent. The numerical methods utilized are derived by use of Green's source solutions obtained using the method of images; first order linearized flow theory is employed with Prandtl-Glauert compressibility transformations. Equations are derived for the flexed shape of a simple constant thickness plate wall under the influence of a finite number of jacks in an axial row along the plate centerline. The Green's source methods are developed to provide estimations of residual flow distortion (interferences) with measured wall pressures and wall flow inclinations as inputs.

Brown, Clinton↗

DOE OSTI2024

Ensemble variational Fokker-Planck methods for data assimilation

Particle flow filters solve Bayesian inference problems by smoothly transforming a set of particles into samples from the posterior distribution. Particles move in state space under the flow of an McKean-Vlasov-Itˆo process. This work introduces the Variational Fokker-Planck (VFP) framework for data assimilation, a general approach that includes previously known particle flow filters as special cases. The McKean-Vlasov-Itˆo process that transforms particles is defined via an optimal drift that depends on the selected diffusion term. It is established that the underlying probability density - sampled by the ensemble of particles - converges to the Bayesian posterior probability density. For a finite number of particles the optimal drift contains a regularization term that nudges particles toward becoming independent random variables. Based on this analysis, we derive computationally-feasible approximate regularization approaches that penalize the mutual information between pairs of particles, and avoid particle collapse. Moreover, the diffusion plays a role akin to a particle rejuvenation approach that aims to alleviate particle collapse. The VFP framework is very flexible. Different assumptions on prior and intermediate probability distributions can be used to implement the optimal drift, and localization and covariance shrinkage can be applied to alleviate the curse of dimensionality. A robust implicit-explicit method is discussed for the efficient integration of stiff McKean- Vlasov-Itˆo processes. Here, the effectiveness of the VFP framework is demonstrated on three progressively more challenging test problems, namely the Lorenz ’63, Lorenz ’96 and the quasi-geostrophic equations.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2015

In-Flight Aeroelastic Stability of the Thermal Protection System on the NASA HIAD, Part II: Nonlinear Theory and Extended Aerodynamics

Conical shell theory and a supersonic potential flow aerodynamic theory are used to study the nonlinear pressure buckling and aeroelastic limit cycle behavior of the thermal protection system for NASA's Hypersonic Inflatable Aerodynamic Decelerator. The structural model of the thermal protection system consists of an orthotropic conical shell of the Donnell type, resting on several circumferential elastic supports. Classical Piston Theory is used initially for the aerodynamic pressure, but was found to be insufficient at low supersonic Mach numbers. Transform methods are applied to the convected wave equation for potential flow, and a time-dependent aerodynamic pressure correction factor is obtained. The Lagrangian of the shell system is formulated in terms of the generalized coordinates for all displacements and the Rayleigh-Ritz method is used to derive the governing differential-algebraic equations of motion. Aeroelastic limit cycle oscillations and buckling deformations are calculated in the time domain using a Runge-Kutta method in MATLAB. Three conical shell geometries were considered in the present analysis: a 3-meter diameter 70 deg. cone, a 3.7-meter 70 deg. cone, and a 6-meter diameter 70 deg. cone. The 6-meter configuration was loaded statically and the results were compared with an experimental load test of a 6-meter HIAD. Though agreement between theoretical and experimental strains was poor, the circumferential wrinkling phenomena observed during the experiments was captured by the theory and axial deformations were qualitatively similar in shape. With Piston Theory aerodynamics, the nonlinear flutter dynamic pressures of the 3-meter configuration were in agreement with the values calculated using linear theory, and the limit cycle amplitudes were generally on the order of the shell thickness. The effect of axial tension was studied for this configuration, and increasing tension was found to decrease the limit cycle amplitudes when the circumferential elastic supports were neglected, but resulted in more complex behavior when the supports were included. The nominal flutter dynamic pressure of the 3.7-meter configuration was significantly lower than that of the 3-meter, and it was found that two sets of natural modes coalesce to flutter modes near the same dynamic pressure. This resulted in a significant drop in the limit cycle frequencies at higher dynamic pressures, where the flutter mode with the lower frequency becomes more critical. Pre-buckling pressure loads and the aerodynamic pressure correction factor were studied for all geometries, and these effects resulted in significantly lower flutter boundaries compared with Piston Theory alone. The maximum dynamic pressure predicted by aerodynamic simulations of a proposed 3.7-meter HIAD vehicle was still lower than any of the calculated flutter dynamic pressures, suggesting that aeroelastic effects for this vehicle are of little concern.

Goldman, Benjamin D.↗

NASA NTRS2000

Description of Panel Method Code ANTARES

Panel method code ANTARES was developed to compute wall interference corrections in a rectangular wind tunnel. The code uses point doublets to represent blockage effects and line doublets to represent lifting effects of a wind tunnel model. Subsonic compressibility effects are modeled by applying the Prandtl-Glauert transformation. The closed wall, open jet, or perforated wall boundary condition may be assigned to a wall panel centroid. The tunnel walls can be represented by using up to 8000 panels. The accuracy of panel method code ANTARES was successfully investigated by comparing solutions for the closed wall and open jet boundary condition with corresponding Method of Images solutions. Fourier transform solutions of a two-dimensional wind tunnel flow field were used to check the application of the perforated wall boundary condition. Studies showed that the accuracy of panel method code ANTARES can be improved by increasing the total number of wall panels in the circumferential direction. It was also shown that the accuracy decreases with increasing free-stream Mach number of the wind tunnel flow field.

Ulbrich, Norbert↗

NASA NTRS2001

Stable Algorithm For Estimating Airdata From Flush Surface Pressure Measurements

An airdata estimation and evaluation system and method, including a stable algorithm for estimating airdata from nonintrusive surface pressure measurements. The airdata estimation and evaluation system is preferably implemented in a flush airdata sensing (FADS) system. The system and method of the present invention take a flow model equation and transform it into a triples formulation equation. The triples formulation equation eliminates the pressure related states from the flow model equation by strategically taking the differences of three surface pressures, known as triples. This triples formulation equation is then used to accurately estimate and compute vital airdata from nonintrusive surface pressure measurements.

Whitmore, Stephen, A.↗

NASA NTRS1975

Shear flow aerodynamics - Lifting surface theory

A lifting surface theory based on a parallel shear flow model is presented for steady, incompressible flows. The theory is intended to account approximately for the presence of a boundary layer. The method of Fourier transforms is used to calculate the pressure on a surface of infinite extent and arbitrary contour. Immediately above the surface is a region of sheared flow (the boundary layer), outside of which the flow velocity is constant. The Fourier transform of the pressure on this surface is used to derive the shear flow equivalent to the kernel function of classical potential flow lifting surface theory. The kernel function provides an integral relation between the upwash at a given point on the surface and the pressure everywhere on the surface. This relation is treated as an integral equation for the pressure, and is solved numerically. Computations are presented for the lift and pitching moment on a flat plate in two-dimensional flow, and for flat, rectangular wings of aspect ratio 1, 2, and 5. As expected, the shear layer decreases the lift curve slope; however, the shear layer (whose thickness is constant along the wing chord) has little effect on the center of pressure.

Ventres, C. S.↗

NASA NTRS1979

A new approach to the solution of large, full matrix equations: A two-dimensional potential flow feasibility study

An approach to the solution of matrix problems resulting from integral equations of mathematical physics is presented. Based on the inherent smoothness in such equations, the problem is reformulated using a set of orthogonal basis vectors, leading to an equivalent coefficient problem which can be of lower order without significantly impairing the accuracy of the solution. This approach was evaluated using a two-dimensional Neumann problem describing the inviscid, incompressible flow over an airfoil. Two different kinds of mode functions were investigated, namely eigenfunction series and Fourier series. The method using Fourier series was found preferable. It uses all of the coefficients from a Fast Fourier Transform algorithm in an approximate method which exploits the known structure of the transformed coefficient matrix and very promising results for the flow over a realistic airfoil are obtained. On the basis of the results presented here, an order of magnitude reduction in this computer time can be expected for such problems as compared with the time for a direct matrix solution.

James, R. M.↗