Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coordinate transformations”

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 631 records · Page 35

Development of a numerical procedure to map a general 3-d body onto a near-circle

Conformal mapping is a classical technique utilized for solving problems in aerodynamics and hydrodynamics. Conformal mapping is utilized in the construction of grids around airfoils, engine inlets and other aircraft configurations. These shapes are transformed onto a near-circle image for which the equations of fluid motion are discretized on the mapped plane and solved numerically by utilizing the appropriate techniques. In comparison to other grid-generation techniques such as algerbraic or differential type, conformal mapping offers an analytical and accurate form even if the grid deformation is large. One of the most appealing features is that the grid can be constrained to remain orthogonal to the body after the transformation. Hence, the grid is suitable for analyzing the supersonic flow past a blunt object. The associated shock as a coordinate surface adjusts its position in the course of computation until convergence is reached. The present work applied conformal mapping to 3-D bodies with no axis of symmetry such as the Aerobraking Flight Experiment (AFE) vehicle, transforming the AFE shape onto a near-circle image. A numerical procedure and code are used to generate grids around the AFE body.

Hommel, M. J.↗

A well posed boundary value problem in transonic gas dynamics

A boundary value problem for the Tricomi equation was studied in connection with transonic gas dynamics. The transformed equation delta u plus 1/3Y u sub Y equals 0 in canonical coordinates was considered in the complex domain of two independent complex variables. A boundary value problem was then set by prescribing the real part of the solution on the boundary of the real unit circle. The Dirichlet problem in the upper unit semicircle with vanishing values of the solution at Y = 0 was solved explicitly in terms of the hypergeometric function for the more general Euler-Poisson-Darboux equation. An explicit representation of the solution was also given for a mixed Dirichlet and Neumann problem for the same equation and domain.

Sanz, J. M.↗

Geometric methods in computational fluid dynamics

General methods for the construction of geometric computational fluid dynamic algorithms are presented which simulate a variety of flow fields in various nontrivial regions. Included are: basic developments with tensors; various forms for the equations of motion; generalized numerical methods and boundary conditions; and methods for mesh generation to meet the strong geometric constraints of turbomachines. Coordinate generation is shown generally to yield mesh descriptions from one or more transformations that are smoothly joined together to form a composite mesh.

Eiseman, P. R.↗

Further development of XTRAN3S computer program

This report describes modifications and enhancements to XTRAN3S, a computer program for aerodynamic, static aeroelastic, and dynamic aeroelastic analysis of three-dimensional wings in the transonic speed regime. Modifications to the program include incorporation of a viscous boundary layer, modified coefficient and generalized force integration, direct input of airfoil surface coordinates and slopes, variable dimensions of computational arrays and a modified grid mapping transformation to improve computational stability for highly swept and tapered planforms. Results obtained with the modified program are included. Modifications discussed but not incorporated include a state matrix integration method, enhanced vectorization, and use of a cartesian physical mesh.

Borland, C. J.↗

Manipulator interactive design with interconnected flexible elements

This paper describes the development of an analysis tool for the interactive design of control systems for manipulators and similar electro-mechanical systems amenable to representation as structures in a topological chain. The chain consists of a series of elastic bodies subject to small deformations and arbitrary displacements. The bodies are connected by hinges which permit kinematic constraints, control, or relative motion with six degrees of freedom. The equations of motion for the chain configuration are derived via Kane's method, extended for application to interconnected flexible bodies with time-varying boundary conditions. A corresponding set of modal coordinates has been selected. The motion equations are imbedded within a simulation that transforms the vector-dyadic equations into scalar form for numerical integration. The simulation also includes a linear, time-invariant controler specified in transfer function format and a set of sensors and actuators that interface between the structure and controller. The simulation is driven by an interactive set-up program resulting in an easy-to-use analysis tool.

Singh, R. P.↗

On a new method for calculating the potential flow past a body of revolution

A new method is presented for obtaining the velocity potential of the flow about a body of revolution moving uniformly in the direction of its axis of symmetry in a fluid otherwise at rest. This method is based essentially on the fact that the form of the differential equation for the velocity potential is invariant with regard to conformal transformation of the meridian plane. By means of the conformal transformation of the meridian profile into a circle a system of orthogonal curvilinear coordinates is obtained, the main feature of which is that one of the coordinate lines is the meridian profile itself. The use of this type of coordinate system yields a simple expression of the boundary condition at the surface of the solid and leads to a rational process of iteration for the solution of the differential equation for the velocity potential. It is shown that the velocity potential for an arbitrary body of revolution may be expressed in terms of universal functions which, although not normal, are obtainable by means of simple quadratures. The general results are applied to a body of revolution obtained by revolving a symmetrical Joukowski profile about its axis of symmetry. A numerical example further serves to illustrate the theory.

Kaplan, Carl↗

Two-Dimensional Inlet Simulation Using a Diagonal Implicit Algorithm

A modification of an implicit approximate-factorization finite-difference algorithm applied to the two-dimensional Euler and Navier-Stokes equations in general curvilinear coordinates is presented for supersonic freestream flow about and through inlets. The modification transforms the coupled system of equations Into an uncoupled diagonal form which requires less computation work. For steady-state applications the resulting diagonal algorithm retains the stability and accuracy characteristics of the original algorithm. Solutions are given for inviscid and laminar flow about a two-dimensional wedge inlet configuration. Comparisons are made between computed results and exact theory.

Chaussee, D.S.↗

Role of Turbulent Prandtl Number on Heat Flux at Hypersonic Mach Number

Present simulation of turbulent flows involving shock wave/boundary layer interaction invariably overestimates heat flux by almost a factor of two. One possible reason for such a performance is a result of the fact that the turbulence models employed make use of Morkovin's hypothesis. This hypothesis is valid for non-hypersonic Mach numbers and moderate rates of heat transfer. At hypersonic Mach numbers, high rates of heat transfer exist in regions where shock wave/boundary layer interactions are important. As a result, one should not expect traditional turbulence models to yield accurate results. The goal of this investigation is to explore the role of a variable Prandtl number formulation in predicting heat flux in flows dominated by strong shock wave/boundary layer interactions. The intended applications involve external flows in the absence of combustion such as those encountered in supersonic inlets. This can be achieved by adding equations for the temperature variance and its dissipation rate. Such equations can be derived from the exact Navier-Stokes equations. Traditionally, modeled equations are based on the low speed energy equation where the pressure gradient term and the term responsible for energy dissipation are ignored. It is clear that such assumptions are not valid for hypersonic flows. The approach used here is based on the procedure used in deriving the k-zeta model, in which the exact equations that governed k, the variance of velocity, and zeta, the variance of vorticity, were derived and modeled. For the variable turbulent Prandtl number, the exact equations that govern the temperature variance and its dissipation rate are derived and modeled term by term. The resulting set of equations are free of damping and wall functions and are coordinate-system independent. Moreover, modeled correlations are tensorially consistent and invariant under Galilean transformation. The final set of equations will be given in the paper.

Xiao, X.↗

Three-Dimensional Piecewise-Continuous Class-Shape Transformation of Wings

Class-Shape Transformation (CST) is a popular method for creating analytical representations of the surface coordinates of various components of aerospace vehicles. A wide variety of two- and three-dimensional shapes can be represented analytically using only a modest number of parameters, and the surface representation is smooth and continuous to as fine a degree as desired. This paper expands upon the original two-dimensional representation of airfoils to develop a generalized three-dimensional CST parametrization scheme that is suitable for a wider range of aircraft wings than previous formulations, including wings with significant non-planar shapes such as blended winglets and box wings. The method uses individual functions for the spanwise variation of airfoil shape, chord, thickness, twist, and reference axis coordinates to build up the complete wing shape. An alternative formulation parameterizes the slopes of the reference axis coordinates in order to relate the spanwise variation to the tangents of the sweep and dihedral angles. Also discussed are methods for fitting existing wing surface coordinates, including the use of piecewise equations to handle discontinuities, and mathematical formulations of geometric continuity constraints. A subsonic transport wing model is used as an example problem to illustrate the application of the methodology and to quantify the effects of piecewise representation and curvature constraints.

Olson, Erik D.↗

Simulation Development and Analysis of Crew Vehicle Ascent Abort

NASA's Commercial Crew Program is an integral step in its journey to Mars as it would expedite development of space technologies and open up partnership with U.S. commercial companies. NASA reviews and independent assessment of Commercial Crew Program is fundamental to its success, and being able to model a commercial crew vehicle in a simulation rather than conduct a live test would be a safer, faster, and less expensive way to assess and certify the capabilities of the vehicle. To this end, my project was to determine the feasibility of using a simulation tool named SOMBAT version 2.0 to model a multiple parachute system for Commercial Crew Program simulation. The main tasks assigned to me were to debug and test the main parachute system model, (capable of simulating one to four main parachute bodies), and to utilize a graphical program to animate the simulation results. To begin tackling the first task, I learned how to use SOMBAT by familiarizing myself with its mechanics and by understanding the methods used to tweak its various parameters and outputs. I then used this new knowledge to set up, run, and analyze many different situations within SOMBAT in order to explore the limitations of the parachute model. Some examples of parameters that I varied include the initial velocity and orientation of the falling capsule, the number of main parachutes, and the location where the parachutes were attached to the capsule. Each parameter changed would give a different output, and in some cases, would expose a bug or limitation in the model. A major bug that I discovered was the inability of the model to handle any number of parachutes other than three. I spent quite some time trying to debug the code logically, but was unable to figure it out until my mentor taught me that digital simulation limitations can occur when some approximations are mistakenly assumed for certain in a physical system. This led me to the realization that unlike in all of the programming classes I have taken thus far that focus on pure logic, simulation code focuses on mimicking the physical world with some approximation and can have inaccuracies or numerical instabilities. Learning from my mistake, I adopted new methods to analyze these different simulations. One method the student used was to numerically plot various physical parameters using MATLAB to confirm the mechanical behavior of the system in addition to comparing the data to the output from a separate simulation tool called FAST. By having full control over what was being outputted from the simulation, I could choose which parameters to change and to plot as well as how to plot them, allowing for an in depth analysis of the data. Another method of analysis was to convert the output data into a graphical animation. Unlike the numerical plots, where all of the physical components were displayed separately, this graphical display allows for a combined look at the simulation output that makes it much easier for one to see the physical behavior of the model. The process for converting SOMBAT output for EDGE graphical display had to be developed. With some guidance from other EDGE users, I developed a process and created a script that would easily allow one to display simulations graphically. Another limitation with the SOMBAT model was the inability for the capsule to have the main parachutes instantly deployed with a large angle between the air speed vector and the chutes drag vector. To explore this problem, I had to learn about different coordinate frames used in Guidance, Navigation & Control (J2000, ECEF, ENU, etc.) to describe the motion of a vehicle and about Euler angles (e.g. Roll, Pitch, Yaw) to describe the orientation of the vehicle. With a thorough explanation from my mentor about the description of each coordinate frame, as well as how to use a directional cosine matrix to transform one frame to another, I investigated the problem by simulating different capsule orientations. In the end, I was able to show that this limitation could be avoided if the capsule is initially oriented antiparallel to its velocity vector.

Wong, Chi S.↗

Laser Tracker Utilization Methodology in Measuring Truth Trajectories for INS Testing on 6 Degree of Freedom Table at the Marshall Space Flight Center's Contact Dynamics Simulation Laboratory with Lessons Learned

When performing Inertial Navigation System (INS) testing at the Marshall Space Flight Center's (MSFC) Contact Dynamics Simulation Laboratory (CDSL) early in 2017, a Leica Geosystems AT901 Laser Tracker system (LLT) measured the twist & sway trajectories as generated by the 6 Degree Of Freedom (6DOF) Table in the CDSL. These LLT measured trajectories were used in the INS software model validation effort. Several challenges were identified and overcome during the preparation for the INS testing, as well as numerous lessons learned. These challenges included determining the position and attitude of the LLT with respect to an INS-shared coordinate frame using surveyed monument locations in the CDSL and the accompanying mathematical transformation, accurately measuring the spatial relationship between the INS and a 6DOF tracking probe due to lack of INS visibility from the LLT location, obtaining the data from the LLT during a test, determining how to process the results for comparison with INS data in time and frequency domains, and using a sensitivity analysis of the results to verify the quality of the results. While many of these challenges were identified and overcome before or during testing, a significant lesson on test set-up was not learned until later in the data analysis process. It was found that a combination of trajectory-dependent gimbal locking and environmental noise introduced non-negligible noise in the angular measurements of the LLT that spanned the evaluated frequency spectrum. The lessons learned in this experiment may be useful for others performing INS testing in similar testing facilities.

Leggett, Jared O.↗

Geometrical approach to neural net control of movements and posture

In one approach to modeling brain function, sensorimotor integration is described as geometrical mapping among coordinates of non-orthogonal frames that are intrinsic to the system; in such a case sensors represent (covariant) afferents and motor effectors represent (contravariant) motor efferents. The neuronal networks that perform such a function are viewed as general tensor transformations among different expressions and metric tensors determining the geometry of neural functional spaces. Although the non-orthogonality of a coordinate system does not impose a specific geometry on the space, this "Tensor Network Theory of brain function" allows for the possibility that the geometry is non-Euclidean. It is suggested that investigation of the non-Euclidean nature of the geometry is the key to understanding brain function and to interpreting neuronal network function. This paper outlines three contemporary applications of such a theoretical modeling approach. The first is the analysis and interpretation of multi-electrode recordings. The internal geometries of neural networks controlling external behavior of the skeletomuscle system is experimentally determinable using such multi-unit recordings. The second application of this geometrical approach to brain theory is modeling the control of posture and movement. A preliminary simulation study has been conducted with the aim of understanding the control of balance in a standing human. The model appears to unify postural control strategies that have previously been considered to be independent of each other. Third, this paper emphasizes the importance of the geometrical approach for the design and fabrication of neurocomputers that could be used in functional neuromuscular stimulation (FNS) for replacing lost motor control.

Review↗

Three-dimensional, shock-on-shock interaction problem

The unsteady, three-dimensional flowfield resulting from the interaction of a plane shock with a cone-shaped vehicle traveling supersonically is determined, using a second-order, shock-capturing, finite-difference approach. The time-dependent, inviscid gasdynamic equations are transformed to include the self-similar property of the flow, to align various coordinate surfaces with known shock waves, and to cluster points in the vicinity of the intersection of the transmitted incident shock and the surface of the vehicle. The governing partial differential equations in conservation-law form are then solved iteratively using MacCormack's (1969) algorithm.

Kutler, P.↗

Multispectral histogram normalization contrast enhancement

A multispectral histogram normalization or decorrelation enhancement which achieves effective color composites by removing interband correlation is described. The enhancement procedure employs either linear or nonlinear transformations to equalize principal component variances. An additional rotation to any set of orthogonal coordinates is thus possible, while full histogram utilization is maintained by avoiding the reintroduction of correlation. For the three-dimensional case, the enhancement procedure may be implemented with a lookup table. An application of the enhancement to Landsat multispectral scanning imagery is presented.

Soha, J. M.↗

A diagonal form of an implicit approximate-factorization algorithm with application to a two dimensional inlet

A modification of an implicit approximate-factorization finite-difference algorithm applied to the two dimensional Euler and Navier-Stokes equations in general curvilinear coordinates is presented for supersonic free stream flow about and through inlets. The modification transforms the coupled system of equations into an uncoupled diagonal form which requires less computation work. For steady-state applications the resulting diagonal algorithm retains the stability and accuracy characteristics of the original algorithm. Solutions are given for inviscid and laminar flow about a two dimensional wedge inlet configuration. Comparisons are made between computed results and exact theory.

Chaussee, D. S.↗

A unified lunar control network: April 1991

This program was designed to combine and transform various control networks of the Moon into a common center-of-mass coordinate system. The first phase, dealing with the near side, was completed and published. This report contains coordinates of 1166 points on the near side of the Moon.

Davies, Merton E.↗

Insights into Tropospheric Ozone from the INTEX Ozonesonde Network Study (IONS)

Ozone profile data from soundings integrate models, aircraft and other ground-based measurements for better interpretation of atmospheric chemistry and dynamics. A well-designed network of ozonesonde stations, with consistent sampling, can answer questions not possible with short campaigns or current satellite technology. The SHADOZ (Southern Hemisphere Additional Ozonesondes) project, for example, has led to these findings about tropical ozone: definition of the zonal tropospheric wave-one pattern in equatorial ozone, characterization of the "Atlantic ozone paradox" and establishment of a link between tropical Atlantic and Indian Ocean pollution. Building on the SHADOZ concept, a short-term ozone network was formed in July-August 2004 to coordinate ozonesonde launches during the ICARTT/INTEX/NEAQS (International Consortium on Atmospheric Research on Transport and Transformation)/Intercontinental Transport Experiment/New England Air Quality Study. In IONS (INTEX Ozonesonde Network Study), more than 250 soundings, with daily frequency at half the sites, were launched from eleven North American stations and an oceanographic ship in the Gulf of Maine. Although the goal was to examine pollution influences under stable high-pressure systems and transport associated with "warm conveyor belt" flows, the INTEX study region was dominated by a series of weak frontal system that mixed aged pollution with stratospheric ozone in the middle troposphere. Deconvoluting ozone sources provides new insights into ozone in the transition between mid-latitude and polar air.

Thompson, Anne M.↗