Engineering topics

Kaul, Upender K.

Publications and source records attributed to Kaul, Upender K..

At least 19 records

NASA NTRS2019

Extending A Correction Method for Unsteady Transonic Aerodynamics to Variable Camber Continuous Trailing Edge Flap

In this study, we have extended our correction method for transonic aerodynamics to baseline and variable camber continuous trailing edge (VCCTEF) geometries of super-critical airfoils at cruise conditions. This correction method modifies Theodorsen theory, which is based on incompressible potential flow over a flat plate. The method is partly based on CFD RANS simulations of oscillating NASA generic transport model (GTM) baseline cruise and VCCTEF airfoil geometries in pitch and shows promise for addressing transonic flutter problems of aeroelasticity. In the earlier study1, the Theodorsen functions were modulated with non-linear functions of Mach number, M, and reduced frequency, k, and thus the correction was demonstrated for the symmetrical NACA0012 airfoil. In the present study, the modulation functions are also allowed to account for the effects of camber and thickness of the airfoil. This is demonstrated by applying the correction method to various symmetric NACA00xx airfoil geometries (NACA001, NACA0002, NACA0012, NACA0015 and NACA0018), the NASA GTM super-critical baseline and VCCTEF airfoil geometries, e.g., the circular airfoil, VCCTEF222. First results indicate that the correction method accounts for a large variation of thickness of airfoils and camber. The present correction method will guide the development of a new state space model for the VCCTEF system and eventually a new transfer function that will be incorporated in a new aeroelastic framework leading to an appropriate transonic flutter model for use in the future aircraft systems in development under the NASA Advanced Air Transportation Technologies (AATT) project.

Kaul, Upender K.↗

NASA NTRS2019

A Correction Method for Unsteady Transonic Aerodynamics

This paper presents a new method for transonic pitching airfoils based on a RANS CFD study and the Theodorsen model of an oscillating pitching flat plate. This study quantifies the deviation of the lift coefficient predictions using CFD from that obtained using the Theodorsen model, which is based on the incompressible potential flow assumption. The present method corrects this theoretical model by modulating the Theodorsen functions by coefficient functions that depend on the reduced frequency and the Mach number. It is demonstrated that the modified theoretical model predicts lift coefficient in good agreement with the CFD results in the Mach number range from incompressible (M =0.2) to transonic (M =0.755) flow for a range of reduced frequencies typical of transonic flutter. The simulations are first validated by comparing pitching NACA0012 airfoil results with experimental results at transonic flight conditions, which establishes the requirements for a grid converged unsteady transonic solution. The hysteresis loop, Cl versus α, attains a grid independent solution that compares well with experiment. The present correction method will guide the development of a new state space model for the Variable Camber Continuous Trailing Edge Flap (VCCTEF) system and eventually a new transfer function that will be incorporated in a new aeroelastic framework leading to an appropriate transonic flutter model for use in the future aircraft systems in development under the NASA Advanced Air Transportation Technologies (AATT) project.

Kaul, Upender K.↗

NASA NTRS2017

Semi-Supervised Learning of Lift Optimization of Multi-Element Three-Segment Variable Camber Airfoil

This chapter describes a new intelligent platform for learning optimal designs of morphing wings based on Variable Camber Continuous Trailing Edge Flaps (VCCTEF) in conjunction with a leading edge flap called the Variable Camber Krueger (VCK). The new platform consists of a Computational Fluid Dynamics (CFD) methodology coupled with a semi-supervised learning methodology. The CFD component of the intelligent platform comprises of a full Navier-Stokes solution capability (NASA OVERFLOW solver with Spalart-Allmaras turbulence model) that computes flow over a tri-element inboard NASA Generic Transport Model (GTM) wing section. Various VCCTEF/VCK settings and configurations were considered to explore optimal design for high-lift flight during take-off and landing. To determine globally optimal design of such a system, an extremely large set of CFD simulations is needed. This is not feasible to achieve in practice. To alleviate this problem, a recourse was taken to a semi-supervised learning (SSL) methodology, which is based on manifold regularization techniques. A reasonable space of CFD solutions was populated and then the SSL methodology was used to fit this manifold in its entirety, including the gaps in the manifold where there were no CFD solutions available. The SSL methodology in conjunction with an elastodynamic solver (FiDDLE) was demonstrated in an earlier study involving structural health monitoring. These CFD-SSL methodologies define the new intelligent platform that forms the basis for our search for optimal design of wings. Although the present platform can be used in various other design and operational problems in engineering, this chapter focuses on the high-lift study of the VCK-VCCTEF system. Top few candidate design configurations were identified by solving the CFD problem in a small subset of the design space. The SSL component was trained on the design space, and was then used in a predictive mode to populate a selected set of test points outside of the given design space. The new design test space thus populated was evaluated by using the CFD component by determining the error between the SSL predictions and the true (CFD) solutions, which was found to be small. This demonstrates the proposed CFD-SSL methodologies for isolating the best design of the VCK-VCCTEF system, and it holds promise for quantitatively identifying best designs of flight systems, in general.

Intelligent Systems↗

NASA NTRS2016

Lift Optimization Study of a Multi-Element Three-Segment Variable Camber Airfoil

This paper reports a detailed computational high-lift study of the Variable Camber Continuous Trailing Edge Flap (VCCTEF) system carried out to explore the best VCCTEF designs, in conjunction with a leading edge flap called the Variable Camber Krueger (VCK), for take-off and landing. For this purpose, a three-segment variable camber airfoil employed as a performance adaptive aeroelastic wing shaping control effector for a NASA Generic Transport Model (GTM) in landing and take-off configurations is considered. The objective of the study is to define optimal high-lift VCCTEF settings and VCK settings/configurations. A total of 224 combinations of VCK settings/configurations and VCCTEF settings are considered for the inboard GTM wing, where the VCCTEFs are configured as a Fowler flap that forms a slot between the VCCTEF and the main wing. For the VCK settings of deflection angles of 55deg, 60deg and 65deg, 18, 19 and 19 vck configurations, respectively, were considered for each of the 4 different VCCTEF deflection settings. Different vck configurations were defined by varying the horizontal and vertical distance of the vck from the main wing. A computational investigation using a Reynolds-Averaged Navier-Stokes (RANS) solver was carried out to complement a wind-tunnel experimental study covering three of these configurations with the goal of identifying the most optimal high-lift configurations. Four most optimal high-lift configurations, corresponding to each of the VCK deflection settings, have been identified out of all the different configurations considered in this study yielding the highest lift performance.

Study↗

NASA NTRS2015

A 3-D Computational Study of a Variable Camber Continuous Trailing Edge Flap (VCCTEF) Spanwise Segment

Results of a computational study carried out to explore the effects of various elastomer configurations joining spanwise contiguous Variable Camber Continuous Trailing Edge Flap (VCCTEF) segments are reported here. This research is carried out as a proof-of-concept study that will seek to push the flight envelope in cruise with drag optimization as the objective. The cruise conditions can be well off design such as caused by environmental conditions, maneuvering, etc. To handle these off-design conditions, flap deflection is used so when the flap is deflected in a given direction, the aircraft angle of attack changes accordingly to maintain a given lift. The angle of attack is also a design parameter along with the flap deflection. In a previous 2D study,1 the effect of camber was investigated and the results revealed some insight into the relative merit of various camber settings of the VCCTEF. The present state of the art has not advanced sufficiently to do a full 3-D viscous analysis of the whole NASA Generic Transport Model (GTM) wing with VCCTEF deployed with elastomers. Therefore, this study seeks to explore the local effects of three contiguous flap segments on lift and drag of a model devised here to determine possible trades among various flap deflections to achieve desired lift and drag results. Although this approach is an approximation, it provides new insights into the "local" effects of the relative deflections of the contiguous spanwise flap systems and various elastomer segment configurations. The present study is a natural extension of the 2-D study to assess these local 3-D effects. Design cruise condition at 36,000 feet at free stream Mach number of 0.797 and a mean aerodynamic chord (MAC) based Reynolds number of 30.734x10(exp 6) is simulated for an angle of attack (AoA) range of 0 to 6 deg. In the previous 2-D study, the calculations revealed that the parabolic arc camber (1x2x3) and circular arc camber (VCCTEF222) offered the best L/D characteristics and minimum drag in cruise. In the present 3-D study, calculations show that for the same C(sub t), the 3-D circular arc camber wing segment produces the largest drag for a given lift, larger than either of the two 2-D configurations, as was also conjectured in the previous study. This study indicates a wing stall around 4.5 deg angle of attack.

Camber↗

NASA NTRS2012

Effect of Inflow Boundary Conditions on Hovering Tilt-Rotor Flows

A study on the effect of inflow turbulence boundary conditions on the local flow on and around a body in flight has been carried out. The study has been carried out using OVERFLOW2 flow solver using the default Spalart-Allmaras turbulence model in OVERFLOW. Many OVERFLOW turbulent flow simulations have been reported using the SA-fv3 model1 over the years. The present study demonstrates that the turbulence levels imposed as a boundary condition (b.c.) in the far-field as implemented in OVERFLOW are not correct. In fact, very low level of turbulence at the far-field boundary as implemented in OVERFLOW results in the SA-fv3 model predicting transition-like profiles on a given body. By choosing sufficiently high levels of the Reynolds number of turbulence, Rt, as a boundary condition, this anomalous behavior of SA-fv3 model is eliminated. Since numerous papers using OVERFLOW have been presented in the literature including that by the author2 using low level of inflow turbulence, it will be beneficial to the CFD community at large and in particular to the OVERFLOW community to understand this effect of the inflow b.c. in the SA-fv3 model. Various results reported over the years in domains such as high-lift applications, drag prediction applications and rotorcraft flow applications can be revisited using the new turbulence inflow boundary condition prescription as suggested in this study. It has been demonstrated in the literature that the effect of inflow turbulence level on the downstream development of flow is too important to be ignored. Demonstrative results are shown and compared with experiment for the case of hovering XV-15 rotor flow at a tip Reynolds number of 4.9 x 10(exp 6) and Mach number of 0.69 in Figs. 1 and 2 below. Fig. 1(a,b) shows skin friction predictions at radial station of r/R = 0.28 corresponding to 5 different inflow boundary conditions obtained with the SA-fv3 model. As the inflow turbulence is increased monotonically, from a value of Rt = 0.1 to 100.0 as shown in Fig. 1, a fully developed turbulent skin friction profile is realized progressively as shown in Fig. 1(b). Similarly skin friction results shown in Fig. 2(a,b) corresponding to radial stations of r/R = 0.72 reinforce the postulate made here that as the inflow turbulence level is increased gradually, fully turbulent on-the-blade simulation is realized uniformly. In the full paper extensive results on skin friction, pressure and velocity profiles as well as the Figure of Merit will be presented and compared with experiment that will illustrate the substantial effect of the inflow boundary conditions on the development of hovering rotor flow downstream.

Kaul, Upender K.↗

NASA NTRS2009

Finite-difference simulation and visualization of elastodynamics in time-evolving generalized curvilinear coordinates

Modeling and simulation of free and forced structural vibrations is essential to an overall structural health monitoring capability. In the various embodiments, a first principles finite-difference approach is adopted in modeling a structural subsystem such as a mechanical gear by solving elastodynamic equations in generalized curvilinear coordinates. Such a capability to generate a dynamic structural response is widely applicable in a variety of structural health monitoring systems. This capability (1) will lead to an understanding of the dynamic behavior of a structural system and hence its improved design, (2) will generate a sufficiently large space of normal and damage solutions that can be used by machine learning algorithms to detect anomalous system behavior and achieve a system design optimization and (3) will lead to an optimal sensor placement strategy, based on the identification of local stress maxima all over the domain.

Kaul, Upender K.↗

NASA NTRS2007

Enhanced Elliptic Grid Generation

An enhanced method of elliptic grid generation has been invented. Whereas prior methods require user input of certain grid parameters, this method provides for these parameters to be determined automatically. "Elliptic grid generation" signifies generation of generalized curvilinear coordinate grids through solution of elliptic partial differential equations (PDEs). Usually, such grids are fitted to bounding bodies and used in numerical solution of other PDEs like those of fluid flow, heat flow, and electromagnetics. Such a grid is smooth and has continuous first and second derivatives (and possibly also continuous higher-order derivatives), grid lines are appropriately stretched or clustered, and grid lines are orthogonal or nearly so over most of the grid domain. The source terms in the grid-generating PDEs (hereafter called "defining" PDEs) make it possible for the grid to satisfy requirements for clustering and orthogonality properties in the vicinity of specific surfaces in three dimensions or in the vicinity of specific lines in two dimensions. The grid parameters in question are decay parameters that appear in the source terms of the inhomogeneous defining PDEs. The decay parameters are characteristic lengths in exponential- decay factors that express how the influences of the boundaries decrease with distance from the boundaries. These terms govern the rates at which distance between adjacent grid lines change with distance from nearby boundaries. Heretofore, users have arbitrarily specified decay parameters. However, the characteristic lengths are coupled with the strengths of the source terms, such that arbitrary specification could lead to conflicts among parameter values. Moreover, the manual insertion of decay parameters is cumbersome for static grids and infeasible for dynamically changing grids. In the present method, manual insertion and user specification of decay parameters are neither required nor allowed. Instead, the decay parameters are determined automatically as part of the solution of the defining PDEs. Depending on the shape of the boundary segments and the physical nature of the problem to be solved on the grid, the solution of the defining PDEs may provide for rates of decay to vary along and among the boundary segments and may lend itself to interpretation in terms of one or more physical quantities associated with the problem.

Kaul, Upender K.↗

NASA NTRS2007

Enhanced Elliptic Grid Generation

Method and system for generating an elliptic grid in generalized coordinates in two or three dimensions, where one or more decay parameters near a boundary segment of a grid are determined as part of the grid solution, rather than being prescribed initially by a user. The decay parameters may vary with one or more generalized coordinates and determine the rate(s) at which separation distances between adjacent grid lines change as one moves toward or away from a grid boundary segment.

Kaul, Upender K.↗

NASA NTRS2005

FIDDLE: A Computer Code for Finite Difference Development of Linear Elasticity in Generalized Curvilinear Coordinates

A three-dimensional numerical solver based on finite-difference solution of three-dimensional elastodynamic equations in generalized curvilinear coordinates has been developed and used to generate data such as radial and tangential stresses over various gear component geometries under rotation. The geometries considered are an annulus, a thin annular disk, and a thin solid disk. The solution is based on first principles and does not involve lumped parameter or distributed parameter systems approach. The elastodynamic equations in the velocity-stress formulation that are considered here have been used in the solution of problems of geophysics where non-rotating Cartesian grids are considered. For arbitrary geometries, these equations along with the appropriate boundary conditions have been cast in generalized curvilinear coordinates in the present study.

Kaul, Upender K.↗

NASA NTRS2002

Elliptic Grid Generation of Spiral-Bevel Pinion Gear Typical of OH-58 Helicopter Transmission

This paper discusses the source term treatment in the numerical solution of elliptic partial differential equations for an interior grid generation problem in generalized curvilinear coordinates. The geometry considered is that of a planar cross-section of a generic spiral-bevel gear tooth typical of a pinion in the OH-58 helicopter transmission. The source terms used are appropriate for an interior grid domain where all the boundaries are prescribed via a combination of Dirichlet and Neumann boundary conditions. New constraints based on the Green's Theorem are derived which uniquely determine the coefficients in the source terms. These constraints are designed for boundary clustered grids where gradients in physical quantities need to be resolved adequately. However, it is seen that the present formulation works satisfactorily for mild clustering also. Thus, a fully automated elliptic grid generation technique is made possible where there is no need for a parametric study of these parameters since the new relations fix these free parameters uniquely.

Kaul, Upender K.↗

NASA NTRS2002

New Boundary Constraints for Elliptic Systems used in Grid Generation Problems

This paper discusses new boundary constraints for elliptic partial differential equations as used in grid generation problems in generalized curvilinear coordinate systems. These constraints, based on the principle of local conservation of thermal energy in the vicinity of the boundaries, are derived using the Green's Theorem. They uniquely determine the so called decay parameters in the source terms of these elliptic systems. These constraints' are designed for boundary clustered grids where large gradients in physical quantities need to be resolved adequately. It is observed that the present formulation also works satisfactorily for mild clustering. Therefore, a closure for the decay parameter specification for elliptic grid generation problems has been provided resulting in a fully automated elliptic grid generation technique. Thus, there is no need for a parametric study of these decay parameters since the new constraints fix them uniquely. It is also shown that for Neumann type boundary conditions, these boundary constraints uniquely determine the solution to the internal elliptic problem thus eliminating the non-uniqueness of the solution of an internal Neumann boundary value grid generation problem.

Kaul, Upender K.↗

NASA NTRS1991

Computated Turbulent Flow In A Turnaround Duct

Kinetic energy of turbulence near wall predicted accurately. Report describes numerical simulation of two-dimensional turbulent flow in duct having 180 degree bend. Study contributes to fund of information needed to design advanced engines containing turbulent flows.

Kaul, Upender K.↗

NASA NTRS1991

Computed Hypersonic Flow About A Sharp Cone

Report describes mostly computational study of hypersonic flow about sharp cone. Part of continuing effort to understand and predict complicated aerodynamic and heat-transfer phenomena in flows over forebodies of such hypersonic vehicles as Space Shuttle and proposed National Aerospace Plane. Focuses on transition from laminar to turbulent flow, with particular reference to heat-transfer predictions.

Kaul, Upender K.↗

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 NTRS1989

Turbulent flow in a 180 deg bend: Modeling and computations

A low Reynolds number k-epsilon turbulence model was presented which yields accurate predictions of the kinetic energy near the wall. The model is validated with the experimental channel flow data of Kreplin and Eckelmann. The predictions are also compared with earlier results from direct simulation of turbulent channel flow. The model is especially useful for internal flows where the inflow boundary condition of epsilon is not easily prescribed. The model partly derives from some observations based on earlier direct simulation results of near-wall turbulence. The low Reynolds number turbulence model together with an existing curvature correction appropriate to spinning cylinder flows was used to simulate the flow in a U-bend with the same radius of curvature as the Space Shuttle Main Engine (SSME) Turn-Around Duct (TAD). The present computations indicate a space varying curvature correction parameter as opposed to a constant parameter as used in the spinning cylinder flows. Comparison with limited available experimental data is made. The comparison is favorable, but detailed experimental data is needed to further improve the curvature model.

Kaul, Upender K.↗

NASA NTRS1989

Prediction of transitional (laminar-turbulent) hypersonic flows using the parabolized Navier-Stokes equations

Hypersonic transitional flow predictions have been made using the parabolized Navier-Stokes equations with an algebraic transition/turbulence model by appropriately modulating the turbulent viscosity with the available intermittency functions for incompressible and compressible flows. A comparison between the predictions with and without a low Reynolds number correction has also been made. The predictions are compared with the available experimental data and with the theory over a range of Mach number. A simple Re(theta)/M(delta) criterion is shown to satisfactorily predict the meridional variation of the onset location of transition on a cone at a small angle of attack, whereas none of the correlations discussed can do that. Various available correlations are discussed vis-a-vis the predictions as to the locations of the onset and the end of transition.

Kaul, Upender K.↗

NASA NTRS1988

Do large structures control their own growth in a mixing layer? - An assessment

Two different two-dimensional free shear layers, the T-layer developing in time from an initial tangential velocity discontinuity separating the two half-spaces, and the S-layer which develops downstream of the origin where two uniform streams of unequal velocity are brought into tangential contact, are compared. Calculations are performed in order to determine to what extent the perturbations induced upstream by large concentrations of vorticity found downstream hasten or retard the subharmonic instability that leads to the formations of these large structures. The results show that the elliptic influence, or the feedback, in a mixing layer is relatively small for small velocity ratios.

Kaul, Upender K.↗