SEARCH · Engineering Papers

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

SourceResearch connectionFrom yearThrough year

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

At least 667 records · Page 37

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Three dimensional flow simulation with application to aeroelastic analysis

The three-dimensional flowfield about realistic launch vehicle configurations is simulated using the Reynolds-averaged Navier-Stokes equations. Turbulent mixing is accounted for by means of the two-layer Baldwin and Lomax (1978) algebraic eddy viscosity model. The Beam and Warming (1976) implicit approximate factorization algorithm is used for the solution of the finite difference equations. Applications include the study of the flowfield about a hemisphere-cylinder configuration both in the subsonic and supersonic flight regimes, and about two hammerhead payload configurations at transonic speeds. A method is also described which permits incorporating this flow solver into a complete algorithm to perform time-domain aeroelastic stability analyses. The vehicle is modeled as a free-free beam and modal superposition techniques are used for the structural-dynamic formulation. Aeroelastic analyses were performed for the hammerhead configurations.

Azevedo, Joao Luiz F.↗

NASA NTRS2013

Shape and Stress Sensing of Multilayered Composite and Sandwich Structures Using an Inverse Finite Element Method

The marked increase in the use of composite and sandwich material systems in aerospace, civil, and marine structures leads to the need for integrated Structural Health Management systems. A key capability to enable such systems is the real-time reconstruction of structural deformations, stresses, and failure criteria that are inferred from in-situ, discrete-location strain measurements. This technology is commonly referred to as shape- and stress-sensing. Presented herein is a computationally efficient shape- and stress-sensing methodology that is ideally suited for applications to laminated composite and sandwich structures. The new approach employs the inverse Finite Element Method (iFEM) as a general framework and the Refined Zigzag Theory (RZT) as the underlying plate theory. A three-node inverse plate finite element is formulated. The element formulation enables robust and efficient modeling of plate structures instrumented with strain sensors that have arbitrary positions. The methodology leads to a set of linear algebraic equations that are solved efficiently for the unknown nodal displacements. These displacements are then used at the finite element level to compute full-field strains, stresses, and failure criteria that are in turn used to assess structural integrity. Numerical results for multilayered, highly heterogeneous laminates demonstrate the unique capability of this new formulation for shape- and stress-sensing.

Cerracchio, Priscilla↗

NASA NTRS1991

An exact algebraic solution of the infimum in H-infinity optimization with output feedback

This paper presents a simple and noniterative procedure for the computation of the exact value of the infimum in the standard H-infinity-optimal control with output feedback. The problem formulation is general and does not place any restrictions on the direct feedthrough terms between the control input and the controlled output variables, and between the disturbance input and the measurement output variables. The method is applicable to systems that satisfy (1) the transfer function from the control input to the controlled output is right-invertible and has no invariant zeros on the j(w) axis and, (2) the transfer function from the disturbance to the measurement output is left-invertible and has no invariant zeros on the j(w) axis. A set of necessary and sufficient conditions for the solvability of H-infinity-almost disturbance decoupling problem via measurement feedback with internal stability is also given.

Chen, Ben M.↗

NASA NTRS1984

An implict LU scheme for the Euler equations applied to arbitrary cascades

An implicit scheme for solving the Euler equations is derived and demonstrated. The alternating-direction implicit (ADI) technique is modified, using two implicit-operator factors corresponding to lower-block-diagonal (L) or upper-block-diagonal (U) algebraic systems which can be easily inverted. The resulting LU scheme is implemented in finite-volume mode and applied to 2D subsonic and transonic cascade flows with differing degrees of geometric complexity. The results are presented graphically and found to be in good agreement with those of other numerical and analytical approaches. The LU method is also 2.0-3.4 times faster than ADI, suggesting its value in calculating 3D problems.

Buratynski, E. K.↗

NASA NTRS1992

Maxwell's second- and third-order equations of transfer for non-Maxwellian gases

Condensed algebraic forms for Maxwell's second- and third-order equations of transfer are developed for the case of molecules described by either elastic hard spheres, inverse-power potentials, or by Bird's variable hard-sphere model. These hardly reduced, yet exact, equations provide a new point of origin, when using the moment method, in seeking approximate solutions in the kinetic theory of gases for molecular models that are physically more realistic than that provided by the Maxwell model. An important by-product of the analysis when using these second- and third-order relations is that a clear mathematical connection develops between Bird's variable hard-sphere model and that for the inverse-power potential.

Baganoff, D.↗

NASA NTRS1997

Multilevel methods for elliptic problems on unstructured grids

Multilevel methods on unstructured grids for elliptic problems are reviewed. The advantages of these techniques are the flexible approximation of the boundaries of complicated physical domains and the ability to adapt the grid to the resolution of fine scaled structures. Multilevel methods, which include multigrid methods and domain decomposition methods, depend on the correct splitting of appropriate finite element spaces. The standard splittings used in the structured grid case cannot be directly extended to unstructured grids due to their requirement for a hierarchical grid structure. Issues related to the application of multilevel methods to unstructured grids are discussed, including how the coarse spaces and transfer operators are defined and how different types of boundary conditions are treated. An obvious way to generate a coarse mesh is to regrid the physical domain several times. Several alternatives are proposed and discussed: node nested coarse spaces, agglomerated coarse spaces and algebraically generated coarse spaces.

Chan, Tony F.↗

NASA NTRS2003

Summary of EASM Turbulence Models in CFL3D With Validation Test Cases

This paper summarizes the Explicit Algebraic Stress Model in k-omega form (EASM-ko) and in k-epsilon form (EASM-ke) in the Reynolds-averaged Navier-Stokes code CFL3D. These models have been actively used over the last several years in CFL3D, and have undergone some minor modifications during that time. Details of the equations and method for coding the latest versions of the models are given, and numerous validation cases are presented. This paper serves as a validation archive for these models.

Rumsey, Christopher L.↗

NASA NTRS2014

Experimental Aeroheating Study of Mid-L/D Entry Vehicle Geometries: NASA LaRC 20-Inch Mach 6 Air Tunnel Test 6966

Aeroheating data on mid lift-to-drag ratio entry vehicle configurations has been obtained through hypersonic wind tunnel testing. Vehicles of this class have been proposed for high-mass Mars missions, such as sample return and crewed exploration, for which the conventional sphere-cone entry vehicle geometries of previous Mars missions are insufficient. Several configurations were investigated, including elliptically-blunted cylinders with both circular and elliptical cross sections, biconic geometries based on launch vehicle dual-use shrouds, and parametrically-optimized analytic geometries. Testing was conducted at Mach 6 over a range of Reynolds numbers sufficient to generate laminar, transitional, and turbulent flow. Global aeroheating data were obtained using phosphor thermography. Both stream-wise and cross-flow transition occured on different configurations. Comparisons were made with laminar and turbulent computational predictions generated with an algebraic turbulence model. Predictions were generally in good agreement in regions of laminar or fully-turbulent flow; however for transitional cases, the lack of a transition onset prediction capability produced less accurate comparisons. The data obtained in this study are intended to be used for prelimary mission design studies and the development and validation of computational methods.

Hollis, Brian R.↗

NASA NTRS1977

Perturbation of a discontinuous transonic flow

The main difficulty in perturbing a discontinuous transonic flow is in the representation of the shift in the location of the discontinuity (shock wave). Herein presented is a method of overcoming this difficulty by using a distorted airfoil as the initial case rather than the real physical airfoil; the distortion is chosen such that the shock location is unchanged by the perturbation. The distorted airfoil is obtained by the use of a strained coordinate system. A direct consequence of the theory is the derivation of an algebraic similarity relation between related airfoils with shock waves at differing locations. Results for simple examples are shown.

Nixon, D.↗

NASA NTRS1985

Computational aspects of unsteady flows

The calculation of unsteady flows and the development of numerical methods for solving unsteady boundary layer equations and their application to the flows around important configurations such as oscillating airfoils are presented. A brief review of recent work is provided with emphasis on the need for numerical methods which can overcome possible problems associated with flow reversal and separation. The zig-zag and characteristic box schemes are described in this context, and when embodied in a method which permits interaction between solutions of inviscid and viscous equations, the characteristic box scheme is shown to avoid the singularity associated with boundary layer equations and prescribed pressure gradient. Calculations were performed for a cylinder started impulsively from rest and oscillating airfoils. The results are presented and discussed. It is conlcuded that turbulence models based on an algebraic specification of eddy viscosity can be adequate, that location of translation is important to the calculation of the location of flow separation and, therefore, to the overall lift of an oscillating airfoil.

Cebeci, T.↗

NASA NTRS1988

A numerical method for computing unsteady 2-D boundary layer flows

A numerical method for computing unsteady two-dimensional boundary layers in incompressible laminar and turbulent flows is described and applied to a single airfoil changing its incidence angle in time. The solution procedure adopts a first order panel method with a simple wake model to solve for the inviscid part of the flow, and an implicit finite difference method for the viscous part of the flow. Both procedures integrate in time in a step-by-step fashion, in the course of which each step involves the solution of the elliptic Laplace equation and the solution of the parabolic boundary layer equations. The Reynolds shear stress term of the boundary layer equations is modeled by an algebraic eddy viscosity closure. The location of transition is predicted by an empirical data correlation originating from Michel. Since transition and turbulence modeling are key factors in the prediction of viscous flows, their accuracy will be of dominant influence to the overall results.

Krainer, Andreas↗

NASA NTRS1994

Application of a Third Order Upwind Scheme to Viscous Flow over Clean and Iced Wings

A 3-D compressible Navier-Stokes solver has been developed and applied to 3-D viscous flow over clean and iced wings. This method uses a third order accurate finite volume scheme with flux difference splitting to model the inviscid fluxes, and second order accurate symmetric differences to model the viscous terms. The effects of turbulence are modeled using a Kappa-epsilon model. In the vicinity of the sold walls the kappa and epsilon values are modeled using Gorski's algebraic model. Sampling results are presented for surface pressure distributions, for untapered swept clean and iced wings made of NACA 0012 airfoil sections. The leading edge of these sections is modified using a simulated ice shape. Comparisons with experimental data are given.

Bangalore, A.↗

NASA NTRS1996

Research in Computational Aeroscience Applications Implemented on Advanced Parallel Computing Systems

Improving the numerical linear algebra routines for use in new Navier-Stokes codes, specifically Tim Barth's unstructured grid code, with spin-offs to TRANAIR is reported. A fast distance calculation routine for Navier-Stokes codes using the new one-equation turbulence models is written. The primary focus of this work was devoted to improving matrix-iterative methods. New algorithms have been developed which activate the full potential of classical Cray-class computers as well as distributed-memory parallel computers.

Wigton, Larry↗

NASA NTRS2017

Aerodynamic Design of the Hybrid Wing Body Propulsion-Airframe Integration

A hybrid wingbody (HWB) concept is being considered by NASA as a potential subsonic transport aircraft that meets aerodynamic, fuel, emission, and noise goals in the time frame of the 2030s. While the concept promises advantages over conventional wing-and-tube aircraft, it poses unknowns and risks, thus requiring in-depth and broad assessments. Specifically, the configuration entails a tight integration of the airframe and propulsion geometries; the aerodynamic impact has to be carefully evaluated. With the propulsion nacelle installed on the (upper) body, the lift and drag are affected by the mutual interference effects between the airframe and nacelle. The static margin for longitudinal stability is also adversely changed. We develop a design approach in which the integrated geometry of airframe (HWB) and propulsion is accounted for simultaneously in a simple algebraic manner, via parameterization of the planform and airfoils at control sections of the wingbody. In this paper, we present the design of a 300-passenger transport that employs distributed electric fans for propulsion. The trim for stability is achieved through the use of the wingtip twist angle. The geometric shape variables are determined through the adjoint optimization method by minimizing the drag while subject to lift, pitch moment, and geometry constraints. The design results clearly show the influence on the aerodynamic characteristics of the installed nacelle and trimming for stability. A drag minimization with the trim constraint yields a reduction of 10 counts in the drag coefficient.

hybrid wing↗

NASA NTRSDate not supplied

Expedited Model Building Method for Rocket Nozzle Side Loading

Statistical analysis of test data is currently the preferred method of creating models for uncharacterized side loading in rocket nozzles. However, statistical analysis requires multiple data sets to be reliable; getting these sets requires building and testing multiple rocket motors, which is not feasible given the incredible time and cost required to do so. By using a combination of geometry, simplified finite element analysis (FEA) analogs, and algebra we can create a limited predictive model from only one test data set. Once created, the model requires only the input from one strain gauge located at each yaw and pitch attachment bracket for the actuation system. This process is not computationally intensive, therefore saving an abundance of time and budget over using full 3D FEA models to develop a similar predictive model. To develop the model, conversion factors between the strain at the gauge locations and corresponding actuator load is derived from the data. These conversion factors are used to create a “conversion ellipse” that bounds the constitutive properties of the nozzle as one lumped factor. Using a 2D cross-sectional analog of the nozzle exit plane, we create a function that allows us to determine the azimuthal location of the unexpected load based on the ratio between yaw and pitch strain readings. Combining this with the “conversion ellipse” yields a function that given the strain reading outputs an induced actuator load, regardless of the azimuthal location of the unexpected load. As previously mentioned, this predictive model is limited in its application as there is inaccuracy introduced by the use of data by only one test as well as simplifying assumptions. This approach makes a useful tool to quickly determine whether uncharacterized actuator loading events can damage actuation systems or induce loading that exceeds design limits.

Seetha A Kolli↗

NASA NTRS1986

An approach to simultaneous system design. I - Semialgebraic geometric methods

This paper introduces semialgebraic parameterization as an approach to analyze simultaneous stabilization and pole placement problems. Rational families of plants of a given McMillan degree, that are simultaneously stabilizable by a fixed family of compensators, are parameterized. For a discrete family of plants, the parameterization problem reduces to the simultaneous stabilization or the pole placement problem of a r-tuple of multi input multi output plants by a nonswitching compensator. It is shown that by removing a semialgebraic subset of a proper algebraic set, the 'space of plants' can be decomposed into components that are either simultaneously stabilizable or simultaneously unstabilizable. Under special cases, explicit parameterization of the semialgebraic set is obtained. Finally a necessary condition for the simultaneous stabilization of single input or single output plants is obtained.

Ghosh, B. K.↗

NASA NTRS1980

Finite difference grid generation by multivariate blending function interpolation

The General Interpolants Method (GIM) code which solves the multidimensional Navier-Stokes equations for arbitrary geometric domains is described. The geometry module in the GIM code generates two and three dimensional grids over specified flow regimes, establishes boundary condition information and computes finite difference analogs for use in the GIM code numerical solution module. The technique can be classified as an algebraic equation approach. The geometry package uses multivariate blending function interpolation of vector-values functions which define the shapes of the edges and surfaces bounding the flow domain. By employing blending functions which conform to the cardinality conditions the flow domain may be mapped onto a unit square (2-D) or unit cube (3-D), thus producing an intrinsic coordinate system for the region of interest. The intrinsic coordinate system facilitates grid spacing control to allow for optimum distribution of nodes in the flow domain.

Anderson, P. G.↗

NASA NTRS1982

Three dimensional flow field inside compressor rotor, including blade boundary layers

The Reynolds stress equation, pressure strain correlation, and dissipative terms and diffusion are discussed in relation to turbulence modelling using the Reynolds stress model. Algebraic modeling of Reynolds stresses and calculation of the boundary layer over an axial cylinder are examined with regards to the kinetic energy model for turbulence modelling. The numerical analysis of blade and hub wall boundary layers, and an experimental study of rotor blade boundary layer in an axial flow compressor rotor are discussed. The Patankar-Spalding numerical method for two dimensional boundary layers is included.

Galmes, J. M.↗