SEARCH · Engineering Papers

Results for “Moment equations”

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 937 records · Page 52

All sources NASA NTRS DOE OSTI

NASA NTRS1973

An economical method of analyzing transient motion of gas-lubricated rotor-bearing systems.

A method of economically evaluating the hydrodynamic forces generated in a gas-lubricated tilting-pad bearing is presented. The numerical method consists of solving the case of the infinite width bearing and then converting this solution to the case of the finite bearing by accounting for end leakage. The approximate method is compared to the finite-difference solution of Reynolds equation and yields acceptable accuracy while running about one-hundred times faster. A mathematical model of a gas-lubricated tilting-pad vertical rotor systems is developed. The model is capable of analyzing a two-bearing-rotor system in which the rotor center of mass is not at midspan by accounting for gyroscopic moments. The numerical results from the model are compared to actual test data as well as analytical results of other investigators.

Falkenhagen, G. L.↗

NASA NTRS1975

Aeroelastic stability of trimmed helicopter blades in forward flight

Equations for moderately large amplitude coupled flap-lag motion of a torsionally rigid hingeless elastic helicopter blade in forward flight are derived. Quasi-steady aerodynamic loads are considered and the effects of reversed flow are included. By using Galerkin's method the spatial dependence of the problem is eliminated and the equations are linearized about a time dependent equilibrium position determined from the trimmed equilibrium position of the rotor in forward flight. In the first trim procedure the rotor is maintained at a fixed value of thrust coefficient with forward flight and horizontal and vertical force equilibrium is satisfied in addition to maintaining zero pitch and roll moments. The second trim procedure maintains only zero pitch and roll moment simulating conditions under which a rotor would be tested in the wind tunnel.

Friedmann, P.↗

NASA NTRS1990

Calculation of turbulent flow in complex geometries with a second-moment closure model

A full Reynolds stress turbulence model is applied to calculate the flow in various channels with complex cross-section. The model is shown to be more robust and to produce more consistent results than algebraic stress models. The relationship between various versions of the latter is provided. It is shown that differential or algebraic stress models derived from a direct application of a popular linear approximation for the pressure-strain terms in the Reynolds stress equations will always underpredict the secondary motion, and hence the accompanying distortion of other flow properties. A simple approximation which leads to better predictions is presented and related to previous practice.

Demuren, A. O.↗

NASA NTRS1992

Feautrier solution of the electron transport equation

An efficient, accurate, and stable solution to the electron transport problem is presented. The two-stream approximation is discussed in detail, but the solution is developed for the general multistream case. Applied to the photoelectron problem, it is shown that accurate escape fluxes are obtained using only 10 altitudes distributed over the region 100 - 1000 km. The local equilibrium, two-stream, and multistream solutions to the photoelectron problem are compared. In addition, a recursion relation for the moments of the elastic angular scattering phase function is presented.

Link, R.↗

NASA NTRS1959

Modified Matrix Method for Calculating Steady-State Span Loading on Flexible Wings in Subsonic Flight

A method is presented for shortening the computations required to determine the steady-state span loading on flexible wings in subsonic flight. The method makes use of tables of downwash factors to find the necessary aerodynamic-influence coefficients for the application of lifting-line theory. Explicit matrix equations of equilibrium are converted into a matrix power series with a finite number of terms by utilizing certain characteristic properties of matrices. The number of terms in the series is determined by a trial-and-error process dependent upon the required accuracy of the solution. Spanwise distributions of angle of attack, airload, shear, bending moment, and pitching moment are readily obtained as functions of qm(sub R) where q denotes the dynamic pressure and mR denotes the lift-curve slope of a rigid wing. This method is intended primarily to make it practical to solve steady-state aeroelastic problems on the ordinary manually operated desk calculators, but the method is also readily adaptable to automatic computing equipment.

Gainer, Patrick A.↗

NASA NTRS1990

The semidiscrete Galerkin finite element modeling of compressible viscous flow past an airfoil

The primary project was the numerical simulation, by a finite element/finite difference method, of the viscous flow about an airfoil. The secondary project involved the numerical simulation of the three-dimensional separated and vortex-dominated flow about a hemispherically capped cylinder in the transonic regime. Preliminary calculations were started for the hemisphere-cylinder at 0 and 5 degree angle of attack. The solution of the flow field about airfoils and wings is required to determine the important parameters of lift, moment, and drag. Viscous effects must be accounted for if the drag is to be accurately calculated. At present there are basically two approaches to the numerical simulation of the flow field, the use of fully viscous models and the inviscid/viscous models. The fully viscous models require the solution of an approximation of the Navier-Stokes equations and therefore should simulate most of the physical mechanisms. A fast, accurate, and computationally efficient inviscid flow solver was recently developed by Hartwich. It is thought that Hartwich's program coupled to a fast, accurate, and computationally efficient boundary layer code, will make an excellent tool for airfoil design. The purpose of the primary project was to develop a compressible boundary layer code using the semidiscrete Galerkin finite element method. The numerical scheme employed used the combination of a Dorodnitsyn formulation of the boundary layer equations, with a finite difference/finite element procedure (semidiscrete Galerkin method), in the solution of the compressible two-dimensional boundary layer equations. A laminar compressible boundary layer code was developed and tested for a NACA 0012 airfoil at a Mach number of 0.5, a Reynolds number of 5000, and zero angle of attack. At present the boundary layer program solves up to, but not beyond, separation.

Meade, Andrew J.↗

NASA NTRS1993

Conical Euler analysis and active roll suppression for unsteady vortical flows about rolling delta wings

A conical Euler code was developed to study unsteady vortex-dominated flows about rolling, highly swept delta wings undergoing either forced motions or free-to-roll motions that include active roll suppression. The flow solver of the code involves a multistage, Runge-Kutta time-stepping scheme that uses a cell-centered, finite-volume, spatial discretization of the Euler equations on an unstructured grid of triangles. The code allows for the additional analysis of the free to-roll case by simultaneously integrating in time the rigid-body equation of motion with the governing flow equations. Results are presented for a delta wing with a 75 deg swept, sharp leading edge at a free-stream Mach number of 1.2 and at 10 deg, 20 deg, and 30 deg angle of attack alpha. At the lower angles of attack (10 and 20 deg), forced-harmonic analyses indicate that the rolling-moment coefficients provide a positive damping, which is verified by free-to-roll calculations. In contrast, at the higher angle of attack (30 deg), a forced-harmonic analysis indicates that the rolling-moment coefficient provides negative damping at the small roll amplitudes. A free-to-roll calculation for this case produces an initially divergent response, but as the amplitude of motion grows with time, the response transitions to a wing-rock type of limit cycle oscillation, which is characteristic of highly swept delta wings. This limit cycle oscillation may be actively suppressed through the use of a rate-feedback control law and antisymmetrically deflected leading-edge flaps. Descriptions of the conical Euler flow solver and the free-to roll analysis are included in this report. Results are presented that demonstrate how the systematic analysis of the forced response of the delta wing can be used to predict the stable, neutrally stable, and unstable free response of the delta wing. These results also give insight into the flow physics associated with unsteady vortical flows about delta wings undergoing forced motions and free-to-roll motions, including the active suppression of the wing-rock type phenomenon. The conical Euler methodology developed is directly extend able to three-dimensional calculations.

Lee-Rausch, Elizabeth M.↗

NASA NTRS2008

Deflection-Compensating Beam for use inside a Cylinder

A design concept for a beam for a specific application permits variations and options for satisfying competing requirements to minimize certain deflections under load and to minimize the weight of the beam. In the specific application, the beam is required to serve as a motion-controlled structure for supporting a mirror for optical testing in the lower third portion of a horizontal, cylindrical vacuum chamber. The cylindrical shape of the chamber is fortuitous in that it can be (and is) utilized as an essential element of the deflection-minimizing design concept. The beam is, more precisely, a table-like structure comprising a nominally flat, horizontal portion with vertical legs at its ends. The weights of the beam and whatever components it supports are reacted by the contact forces between the lower ends of the legs and the inner cylindrical chamber wall. Whereas the bending moments arising from the weights contribute to a beam deflection that is concave with its lowest point at midlength, the bending moments generated by the contact forces acting on the legs contribute to a beam deflection that is convex with its highest point at midlength. In addition, the bending of the legs in response to the weights causes the lower ends of the legs to slide downward on the cylindrical wall. By taking the standard beam-deflection equations, combining them with the geometric relationships among the legs and the horizontal portion of the beam, and treating the sliding as a component of deflection, it is possible to write an equation for the net vertical deflection as a function of the load and of position along the beam. A summary of major conclusions drawn from the equation characterization is included.

Goodman, Dwight↗

NASA NTRS2015

An Analytical Solution for Yaw Maneuver Optimization on the International Space Station and Other Orbiting Space Vehicles

This paper suggests a new method for optimizing yaw maneuvers on the International Space Station (ISS). Yaw rotations are the most common large maneuvers on the ISS often used for docking and undocking operations, as well as for other activities. When maneuver optimization is used, large maneuvers, which were performed on thrusters, could be performed either using control moment gyroscopes (CMG), or with significantly reduced thruster firings. Maneuver optimization helps to save expensive propellant and reduce structural loads ‐ an important factor for the ISS service life. In addition, optimized maneuvers reduce contamination of the critical elements of the vehicle structure, such as solar arrays. This paper presents an analytical solution for optimizing yaw attitude maneuvers. Equations describing pitch and roll motion needed to counteract the major torques during a yaw maneuver are obtained. A yaw rate profile is proposed. Also the paper describes the physical basis of the suggested optimization approach. In the obtained optimized case, the torques are significantly reduced. This torque reduction was compared to the existing optimization method which utilizes the computational solution. It was shown that the attitude profiles and the torque reduction have a good match for these two methods of optimization. The simulations using the ISS flight software showed similar propellant consumption for both methods. The analytical solution proposed in this paper has major benefits with respect to computational approach. In contrast to the current computational solution, which only can be calculated on the ground, the analytical solution does not require extensive computational resources, and can be implemented in the onboard software, thus, making the maneuver execution automatic. The automatic maneuver significantly simplifies the operations and, if necessary, allows to perform a maneuver without communication with the ground. It also reduces the probability of command errors. The suggested analytical solution provides a new method of maneuver optimization which is less complicated, automatic and more universal. A maneuver optimization approach, presented in this paper, can be used not only for the ISS, but for other orbiting space vehicles.

Dobrinskaya, Tatiana↗

DOE OSTI2022

Hybrid simulation of energetic particles interacting with magnetohydrodynamics using a slow manifold algorithm and GPU acceleration

The hybrid method combining particle-in-cell and magnetohydrodynamics can be used to study the interaction between energetic particles and global plasma modes. In this paper we introduce the M3D-C1-K code, which is developed based on the M3D-C1 finite element code solving the magnetohydrodynamics equations, with a newly developed kinetic module simulating energetic particles. The particle pushing is done using a new algorithm by applying the Boris pusher to the classical Pauli particles to simulate the slow-manifold of particle orbits, with long-term accuracy and fidelity. The particle pushing can be accelerated using GPUs with a significant speedup. The moments of the particles are calculated using the δƒ method, and are coupled into the magnetohydrodynamics simulation through pressure or current coupling schemes. Several linear simulations of magnetohydrodynamics modes driven by energetic particles have been conducted using M3D-C1-K with the δƒ method, including fishbone, toroidal Alfvén eigenmodes and reversed shear Alfvén eigenmodes. Good agreement with previous results from other eigenvalue, kinetic and hybrid codes have been achieved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1973

High-energy proton radiation belt.

The experiments and theories to explain the high-energy protons trapped in the earth's radiation belt are reviewed. The theory of cosmic ray albedo neutron decay injection of protons into the radiation belt is discussed. Radial diffusion and change in the earth's dipole moment are considered along with losses of protons by ionization and nuclear collision. It is found that the measured albedo neutron escape current is sufficient to supply trapped protons above 30 MeV. The theoretical calculations of the trapped protons are in agreement with the measurements for L less than or equal to 1.7 both on and off the equator. For L greater than or equal to 1.7, additional trapped proton differential energy measurements should be made before the theory can be adequately tested. It appears that an additional loss mechanism such as pitch angle scattering may be required.

White, R. S.↗

NASA NTRS1987

Calculation of viscous effects on transonic flow for oscillating airfoils and comparisons with experiment

A method is described for calculating unsteady transonic flow with viscous interaction by coupling a steady integral boundary-layer code with an unsteady, transonic, inviscid small-disturbance computer code in a quasi-steady fashion. Explicit coupling of the equations together with viscous -inviscid iterations at each time step yield converged solutions with computer times about double those required to obtain inviscid solutions. The accuracy and range of applicability of the method are investigated by applying it to four AGARD standard airfoils. The first-harmonic components of both the unsteady pressure distributions and the lift and moment coefficients have been calculated. Comparisons with inviscid calcualtions and experimental data are presented. The results demonstrate that accurate solutions for transonic flows with viscous effects can be obtained for flows involving moderate-strength shock waves.

Howlett, James T.↗

NASA NTRS2008

An Alternate Set of Basis Functions for the Electromagnetic Solution of Arbitrarily-Shaped, Three-Dimensional, Closed, Conducting Bodies Using Method of Moments

In this work, we present an alternate set of basis functions, each defined over a pair of planar triangular patches, for the method of moments solution of electromagnetic scattering and radiation problems associated with arbitrarily-shaped, closed, conducting surfaces. The present basis functions are point-wise orthogonal to the pulse basis functions previously defined. The prime motivation to develop the present set of basis functions is to utilize them for the electromagnetic solution of dielectric bodies using a surface integral equation formulation which involves both electric and magnetic cur- rents. However, in the present work, only the conducting body solution is presented and compared with other data.

Mackenzie, Anne I.↗

NASA NTRS1981

Interaction of upstream flow distortions with high Mach number cascades

Features of the interaction of flow distortions, such as gusts and wakes with blade rows of advance type fans and compressors having high tip Mach numbers are modeled. A typical disturbance was assumed to have harmonic time dependence and was described, at a far upstream location, in three orthogonal spatial coordinates by a double Fourier series. It was convected at supersonic relative to a linear cascade described as an unrolled annulus. Conditions were selected so that the component of this velocity parallel to the axis of the turbomachine was subsonic, permitting interaction between blades through the upstream as well as downstream flow media. A strong, nearly normal shock was considered in the blade passages which was allowed curvature and displacement. The flows before and after the shock were linearized relative to uniform mean velocities in their respective regions. Solution of the descriptive equations was by adaption of the Wiener-Hopf technique, enabling a determination of distortion patterns through and downstream of the cascade as well as pressure distributions on the blade and surfaces. Details of interaction of the disturbance with the in-passage shock were discussed. Infuences of amplitude, wave length, and phase of the disturbance on lifts and moments of cascade configurations are presented. Numerical results are clarified by reference to an especially orderly pattern of upstream vertical motion in relation to the cascade parameters.

Englert, G. W.↗

NASA NTRS1982

Mars' atmospheric behavior from Viking infra-red thermal mapper measurements

Data from the 15 micron band of CO2 readings with the two Viking IR thermal mappers are discussed. Contrasts were observed to be strong between clear and dusty conditions, with a latitudinal gradient and a diurnal amplitude variation in the winter southern hemisphere. Consistency was found in zonal mean temperatures in the absence of dust, with a peak temperature of 180 K at the poles and a diurnal amplitude of 15 K at the equator. Large temperature increases occur in dusty conditions, with global dust storms being present in the northern, but not southern, hemisphere. Estimations of the surface and atmospheric temperatures are calculated in order to derive optical depths from the IR measurements of atmospheric opacity. The optical depth around the whole planet is found to be relatively uniform at any given moment. Finally, the diurnal behavior of the brightness temperature is outlined for 1.4 Mars years.

Martin, T. Z.↗

NASA NTRS2016

Application of a Full Reynolds Stress Model to High Lift Flows

A recently developed second-moment Reynolds stress model was applied to two challenging high-lift flows: (1) transonic flow over the ONERA M6 wing, and (2) subsonic flow over the DLR-F11 wing-body configuration from the second AIAA High Lift Prediction Workshop. In this study, the Reynolds stress model results were contrasted with those obtained from one- and two{equation turbulence models, and were found to be competitive in terms of the prediction of shock location and separation. For an ONERA M6 case, results from multiple codes, grids, and models were compared, with the Reynolds stress model tending to yield a slightly smaller shock-induced separation bubble near the wing tip than the simpler models, but all models were fairly close to the limited experimental surface pressure data. For a series of high-lift DLR{F11 cases, the range of results was more limited, but there was indication that the Reynolds stress model yielded less-separated results than the one-equation model near maximum lift. These less-separated results were similar to results from the one-equation model with a quadratic constitutive relation. Additional computations need to be performed before a more definitive assessment of the Reynolds stress model can be made.

Lee-Rausch, E. M.↗

NASA NTRSDate not supplied

FUN3D Mesh Adaptation Simulations of 4th AIAA High-Lift Prediction Workshop Case 1a

This paper describes in detail the FUN3D Mesh adaptation simulation results for Case 1a of the 4th AIAA High-Lift Prediction Workshop (HLPW-4). The Stabilized Finite Element library with FUN3D is used along with the Negative Spalart-Allmaras One-Equation Model to simulate fully turbulent free-air flow around the High Lift version of the Common Research Model at three flap settings. Meshes are adapted using multi-scale and goal-oriented metrics. The results are presented and the discussion focuses on the predicted differences between the three flap settings. Coefficients of lift, drag, and pitching moment are predicted within 7% of values measured in the QinetiQ 5m Wind Tunnel at each flap setting. Predicted surface pressures and skin friction distributions show inconsistent discrepancies with measurements across the flap settings. Consequently, flap increments for these figures of merit are not well predicted. These predictions are consistent with submissions to the HLPW-4 from other RANS-SA methods. Visualizations of Liutex vortex cores colored by angular velocity are presented as possible new quantitative and qualitative perspective for solution evaluation.

High-Lift Prediction Workshop↗

NASA NTRS1990

Numerical solutions of the linearized Euler equations for unsteady vortical flows around lifting airfoils

A linearized unsteady aerodynamic analysis is presented for unsteady, subsonic vortical flows around lifting airfoils. The analysis fully accounts for the distortion effects of the nonuniform mean flow on the imposed vortical disturbances. A frequency domain numerical scheme which implements this linearized approach is described, and numerical results are presented for a large variety of flow configurations. The results demonstrate the effects of airfoil thickness, angle of attack, camber, and Mach number on the unsteady lift and moment of airfoils subjected to periodic vortical gusts. The results show that mean flow distortion can have a very strong effect on the airfoil unsteady response, and that the effect depends strongly upon the reduced frequency, Mach number, and gust wave numbers.

Scott, James R.↗