Search NASA⌕ Search

SEARCH · Search NASA

Results for “generation equation”

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 109 records · Page 6

A multistage mesh generator for solving the average-passage equation system

One means of numerically simulating the 3-D flow field within a multistage turbomachine is through the solution of the average-passage equation system. One requirement of a current algorithm used to solve this system of equations has been the ability to generate multiple blade row meshes which satisfy specific geometrical constraints. In addition to meeting this criterion, one desires a mesh generation code which requires minimal user input, utilizes variable mesh control parameters, generates diagnostics helpful to the user, and possesses the capability to handle widely varying geometries. A mesh generation code with these features was written and has been used in solving the inviscid form of the average-passing equation system for both ducted and unducted multiple blade row geometries. This paper serves as a user reference guide, with a description of the mesh generation algorithm, a sample input file, and examples of typical meshes generated.

Mulac, Richard A.↗

Grid generation in three dimensions by Poisson equations with control of cell size and skewness at boundary surfaces

An algorithm for generating computational grids about arbitrary three-dimensional bodies is developed. The elliptic partial differential equation (PDE) approach developed by Steger and Sorenson and used in the NASA computer program GRAPE is extended from two to three dimensions. Forcing functions which are found automatically by the algorithm give the user the ability to control mesh cell size and skewness at boundary surfaces. This algorithm, as is typical of PDE grid generators, gives smooth grid lines and spacing in the interior of the grid. The method is applied to a rectilinear wind-tunnel case and to two body shapes in spherical coordinates.

Sorenson, R. L.↗

A viscous flow analysis of the tip vortex generation process

A viscous primary-secondary flow analysis designed for numerical solution by forward marching integration is applied to the tip vortex generation problem. Equations governing streamwise momentum, streamwise vorticity, and a gas law are solved in conjunction with equations governing scalar and vector surface potentials for the secondary velocities. The numerical method used combines a consistently split linearized block implicit (LBI) scheme for parabolic equations and a scalar iterative ADI scheme for elliptic equations. Computed results are presented for a wing of rectangular planform immersed in a high Reynolds number stream at 6 degrees incidence. The physical mechanism for generation of the tip vortex, as revealed by the analysis and computations, is discussed. Detailed flow results for the case of 6 degree incidence are included.

Shamroth, S. J.↗

Approach guidance logic for a tilt-rotor aircraft

The distinctive feature of a tilt-rotor aircraft is that the pilot can change the rotor mast angles to go from a helicopter configuration for take-off and landing to an airplane configuration for high cruise speeds and vice-versa. An approach path for such an aircraft is proposed and the logic required to fly along this path in the presence of wind is determined. The main contribution of this work is an efficient and, to my knowledge, new method for generating the nominal state and control histories taking into account an estimate of the mean wind velocity and direction. The method requires the solution of algebraic (mostly linear) equations to generate a 'universal nominal', and feedforward and feedback gains. Then, in flight the additional state and control corrections due to deviation in descent rate, deceleration, and flight in a steady wind are obtained by multiplying simple precalculated functions of time by descent rate, deceleration or sine and cosine components of the mean wind vector. Simulations of approach flights for different wind conditions, assuming perfect state information in the feedback signal, indicated satisfactory performance.

Beser, J.↗

Multibody dynamics model building using graphical interfaces

In recent years, the extremely laborious task of manually deriving equations of motion for the simulation of multibody spacecraft dynamics has largely been eliminated. Instead, the dynamicist now works with commonly available general purpose dynamics simulation programs which generate the equations of motion either explicitly or implicitly via computer codes. The user interface to these programs has predominantly been via input data files, each with its own required format and peculiarities, causing errors and frustrations during program setup. Recent progress in a more natural method of data input for dynamics programs: the graphical interface, is described.

Macala, Glenn A.↗

Excitation of turbulence by density waves

A nonlinear system describes the microdynamical state of turbulence that is excited by density waves. It consists of an equation of propagation and a master equation. A group-scaling generates the scaled equations of many interacting groups of distribution functions. The two leading groups govern the transport processes of evolution and eddy diffusivity. The remaining sub-groups represent the relaxation for the approach of diffusivity to equilibrium. In strong turbulence, the sub-groups disperse themselves and the ensemble acts like a medium that offers an effective damping to close the hierarchy. The kinetic equation of turbulence is derived. It calculates the eddy viscosity and identifies the effective damping of the assumed medium self-consistently. It formulates the coupling mechanism for the intensification of the turbulent energy at the expense of the wave energy, and the transfer mechanism for the cascade. The spectra of velocity and density fluctuations find the power law k sup-2 and k sup-4, respectively.

Tichen, C. M.↗

Three-dimensional elliptic grid generation for an F-16

A case history depicting the effort to generate a computational grid for the simulation of transonic flow about an F-16 aircraft at realistic flight conditions is presented. The flow solver for which this grid is designed is a zonal one, using the Reynolds averaged Navier-Stokes equations near the surface of the aircraft, and the Euler equations in regions removed from the aircraft. A body conforming global grid, suitable for the Euler equation, is first generated using 3-D Poisson equations having inhomogeneous terms modeled after the 2-D GRAPE code. Regions of the global grid are then designated for zonal refinement as appropriate to accurately model the flow physics. Grid spacing suitable for solution of the Navier-Stokes equations is generated in the refinement zones by simple subdivision of the given coarse grid intervals. That grid generation project is described, with particular emphasis on the global coarse grid.

Reese L. Sorenson↗

Evaluation of the Accuracy of the Load Prediction Equations of Low-Load Balance Calibration Data

Accuracy and reliability of the load prediction equations of a low-load calibration data set of a force balance were investigated. First, independent load prediction equations were generated from the data of a full-load and low-load machine calibration. Then, the low- load data set was processed as a set of precision check loads for the load prediction equations that were obtained from the full-load and low-load calibrations. Finally, the load prediction equations of the two calibrations were applied to four manual check load data sets that were recorded between 2009 and 2022. The residuals of the predicted calibration and check loads were compared in all cases. No systematic improvement of the load prediction accuracy was observed when the load prediction equations of the low-load calibration were applied to the check load data sets. Therefore, it is recommended to apply a full-load calibration to a six-component balance even if the balance is not used across its entire design envelope. This approach has the advantage that the resulting load prediction equations are less likely to be applied outside of the calibration load ranges. In addition, it is expected that the numerical estimates of the gage sensitivities are more reliable assuming that the uncertainties of the applied loads and measured outputs of the calibration data are more or less constant across the load range of the balance.

strain-gage balance↗

Evaluation of the Accuracy of the Load Prediction Equations of Low-Load Balance Calibration Data

Accuracy and reliability of the load prediction equations of a low-load calibration data set of a force balance were investigated. First, independent load prediction equations were generated from the data of a full-load and low-load machine calibration. Then, the low-load data set was processed as a set of precision check loads for the load prediction equations that were obtained from the full-load and low-load calibrations. Finally, the load prediction equations of the two calibrations were applied to four manual check load data sets that were recorded between 2009 and 2022. The residuals of the predicted calibration and check loads were compared in all cases. No systematic improvement of the load prediction accuracy was observed when the load prediction equations of the low-load calibration were applied to the check load data sets. Therefore, it is recommended to apply a full-load calibration to a six-component balance even if the balance is not used across its entire design envelope. This approach has the advantage that the resulting load prediction equations are less likely to be applied outside of the calibration load ranges. In addition, it is expected that the numerical estimates of the bridge sensitivities are more reliable assuming that the uncertainties of the applied loads and measured outputs of the calibration data are more or less constant across the load range of the balance.

strain-gage balance↗

A computer program to generate two-dimensional grids about airfoils and other shapes by the use of Poisson's equation

A method for generating two dimensional finite difference grids about airfoils and other shapes by the use of the Poisson differential equation is developed. The inhomogeneous terms are automatically chosen such that two important effects are imposed on the grid at both the inner and outer boundaries. The first effect is control of the spacing between mesh points along mesh lines intersecting the boundaries. The second effect is control of the angles with which mesh lines intersect the boundaries. A FORTRAN computer program has been written to use this method. A description of the program, a discussion of the control parameters, and a set of sample cases are included.

Sorenson, R. L.↗

Zonal-Mean Temperature Variations Inferred from SABER Measurements on TIMED Compared with UARS Observations

In the Numerical Spectral Model (NSM, Mayr et al., 2003), small-scale gravity waves propagating in the north/south direction can generate zonal mean (m = 0) meridional wind oscillations with periods between 2 and 4 months. These oscillations tend to be confined to low latitudes and have been interpreted to be the meridional counterpart of the wave-driven Quasi Biennial Oscillation in the zonal circulation. Wave driven meridional winds across the equator should generate, due to dynamical heating and cooling, temperature oscillations with opposite phase in the two hemispheres. We have analyzed SABER temperature measurements in the altitude range between 55 and 95 km to investigate the existence such variations. Because there are also strong tidal signatures (up to approximately 20 K) in the data, our algorithm estimates both mean values and tides together from the data. Based on SABER temperature data, the intra-annual variations with periods between 2 and 4 months can have amplitudes up to 5 K or more, depending on the altitude. Their amplitudes are in qualitative agreement with those inferred Erom UARS data (from different years). The SABER temperature variations also reveal pronounced hemispherical asymmetries, which are qualitatively consistent with wave driven meridional wind oscillations across the equator. Oscillations with similar periods have been seen in the meridional winds based on UARS data (Huang and Reber, 2003).

Huang, Frank T.↗

Approximate-factorization schemes for solving the transonic full-potential equation

The present paper provides a general discussion of approximate-factorization techniques applied to the transonic full-potential equation. Giving particular attention to the AF2 approximate-factorization scheme. This scheme was first introduced by Ballhaus and Steger (1975) for solving the low-frequency (unsteady), transonic small-disturbance equation. The full-potential equation algorithm is examined, taking into account the governing equations, grid generation, the artificial density scheme (spatial differencing), the alternating direction implicit scheme, the AF2 iteration scheme, temporal damping, and boundary conditions. Computed results are also presented. It is shown that fast, fully-implicit algorithms of the approximate-factorization variety are both efficient and reliable for solving the conservative full-potential equation.

Holst, T. L.↗

Trees, bialgebras and intrinsic numerical algorithms

Preliminary work about intrinsic numerical integrators evolving on groups is described. Fix a finite dimensional Lie group G; let g denote its Lie algebra, and let Y(sub 1),...,Y(sub N) denote a basis of g. A class of numerical algorithms is presented that approximate solutions to differential equations evolving on G of the form: dot-x(t) = F(x(t)), x(0) = p is an element of G. The algorithms depend upon constants c(sub i) and c(sub ij), for i = 1,...,k and j is less than i. The algorithms have the property that if the algorithm starts on the group, then it remains on the group. In addition, they also have the property that if G is the abelian group R(N), then the algorithm becomes the classical Runge-Kutta algorithm. The Cayley algebra generated by labeled, ordered trees is used to generate the equations that the coefficients c(sub i) and c(sub ij) must satisfy in order for the algorithm to yield an rth order numerical integrator and to analyze the resulting algorithms.

Crouch, Peter↗

Intra-seasonal Oscillations Inferred from SABER (TIMED) and MLS (UARS) Temperature Measurements

In the zonal mean meridional winds of the upper mesosphere, intra-seasonal oscillations with periods between 1 and 4 months have been inferred from UARS measurements and independently predicted with the Numerical Spectral Model WSM). The wind oscillations tend to be confined to low latitudes and appear to be driven, at least in part, by small-scale gravity waves propagating in the meridional direction. Winds across the equator should generate, due to dynamical heating and cooling, temperature oscillations with opposite phase in the two hemispheres. Investigating this phenomenon, we have analyzed SABER temperatures from TIMED in the altitude range between 55 and 95 km to delineate with an empirical model, the year-long variability of the migrating tides and zonal mean components. The inferred seasonal variations of the diurnal tide, characterized by amplitude maxima near equinox, are in substantial agreement with UARS observations and results from the NSM. For the zonal mean, the dominant seasonal variations in the SABER temperatures, with annual (12 months) and semiannual (6 months) periodicities, agree well with those derived from UARS measurements. The intra-seasonal variations with periods between 2 and 4 months have amplitudes close to 2 K, almost half as large as those for the dominant seasonal variations. Their amplitudes are in qualitative agreement with the corresponding values inferred from UARS during different years. The SABER and UARS temperature variations reveal pronounced hemispherical asymmetries, consistent with meridional wind oscillations across the equator. The phase of the semi-annual temperature oscillations from the NSM agrees with the observations from UARS and SABER. But the amplitudes are systematically smaller, which may indicate that planetary waves are more important than is allowed for in the model. For the shorter-period intra-seasonal variations, which can be generated by gravity wave drag, the model results are generally in better agreement with the observations.

Huang, F. T.↗

On numerical nonlinear analysis of highly flexible spinning cantilevers

The general nonlinear discretized equations of motion of spinning elastic solids and structures are derived as a set of nonlinear ordinary differential equations for the case when the strain-displacement and velocity-displacement relations are nonlinear up to the second order. It is shown that the cost of generation of such equations is proportional to the fourth power of the number of degrees of freedom. A computer program is written to automatically generate the equations for the case of spinning cantilevers with initial imperfections. The types and the number of the coordinate functions used in the trial solution are parameters of the program.

Utku, S.↗

Grid generation for general 2-D regions using hyperbolic equations

A method for applying a hyperbolic grid generation scheme to the construction of meshes in general 2-D regions has been developed. This approach, which follows the theory developed by Steger and Chaussee (1980) and the algorithm outlined by Kinsey and Barth (1984), is based on improving local grid control. This is accomplished by adding an angle control source term to the equations and using a new algorithm for computing the volume source term. These modifications lead to superior methods for fixing the 'local' problems of hyperbolic grid generation, namely, propagation of initial discontinuities and formation of grid shocks (crossing grid lines). More importantly, a method for solving the global problem of constraining the grid with more than one boundary (internal grid generation) has been developed. These algorithms have been implemented in an interactive grid generation program and the results for several geometries are presented and discussed.

Cordova, Jeffrey Q.↗

Evaluation of Delamination Growth Characterization Methods Under Mode I Fatigue Loading

Reliable delamination characterization data for laminated composites are needed for input to analytical models of structures to predict delamination. The double-cantilevered beam (DCB) specimen is used with laminated composites to measure fracture toughness, G(sub Ic), delamination onset strain energy release rate, and growth rate data under cyclic loading. In the current study, DCB specimens of IM7/8552 graphite/epoxy supplied by two different manufacturers were tested in static and fatigue to compare the measured characterization data from the two sources, and to evaluate a proposed ASTM standard for generating Paris Law equations. Static results were used to generate compliance calibration constants for the fatigue data, and a delamination resistance curve, G(sub IR), which was used to determine the effects of fiber-bridging on delamination growth. Specimens were tested in fatigue at a cyclic G(sub Imax) level equal to 50, 40 or 30% of G(sub Ic), to determine a delamination onset curve and delamination growth rate. The delamination onset curve equations had similar exponents and the same trends. Delamination growth rate was calculated by fitting a Paris Law to the da/dN versus G(sub Imax) data. Both a 2-point and a 7-point data reduction method were used and the Paris Law equations were compared. To determine the effects of fiber-bridging, growth rate results were normalized by the delamination resistance curve for each material and compared to the non-normalized results. Paris Law exponents were found to decrease by 31% to 37% due to normalizing the growth data. Normalizing the data also greatly reduced the amount of scatter between the different specimens. Visual data records from the fatigue testing were used to calculate individual compliance calibration constants from the fatigue data for some of the specimens. The resulting da/dN versus G(sub Imax) plots showed much improved repeatability between specimens. Gretchen

Murri, Gretchen B.↗