Character of the Equations of Elasto-Plastic Flow in Three Independent Variables
Elastoplastic flow equations character in three independent variables, noting solution methods for initial and boundary value problems
Swedlow, J. L.↗
SEARCH · Search NASA
Results for “FLOW 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.
Interplanetary gas dynamics.
Observations of the interplanetary medium are considered together with the structure of the solar atmosphere and the kinetic properties of the interplanetary plasma. Methods of gas kinetics and continuum flows are examined, taking into account microscopic and macroscopic representations, the kinetic flow equations, the formal solution of the Vlasov equation, and the continuum flow equations. The collective particle behavior of that interplanetary gas is discussed along with the hydrodynamic coronal expansion, a wave-pump problem, the free expansion phenomenon, and a generalized free expansion problem.
Liu, V. C.↗
An investigation on the effect of second-order additional thickness distributions to the upper surface of an NACA 64 sub 1-212 airfoil
An investigation was conducted on a CDC 7600 digital computer to determine the effects of additional thickness distributions to the upper surface of an NACA 64 sub 1 - 212 airfoil. Additional thickness distributions employed were in the form of two second-order polynomial arcs which have a specified thickness at a given chordwise location. The forward arc disappears at the airfoil leading edge, the aft arc disappears at the airfoil trailing edge. At the juncture of the two arcs, x = x, continuity of slope is maintained. The effect of varying the maximum additional thickness and its chordwise location on airfoil lift coefficient, pitching moment, and pressure distribution was investigated. Results were obtained at a Mach number of 0.2 with an angle-of-attack of 6 degrees on the basic NACA 64 sub 1 - 212 airfoil, and all calculations employ the full potential flow equations for two dimensional flow. The relaxation method of Jameson was employed for solution of the potential flow equations.
Hague, D. S.↗
An investigation on the effect of second-order additional thickness distributions to the upper surface of an NACA 64-206 airfoil
An investigation was conducted on a CDC 7600 digital computer to determine the effects of additional thickness distributions to the upper surface of an NACA 64-206 airfoil. Additional thickness distributions employed were in the form of two second-order polynomial arcs which have a specified thickness at a given chordwise location. The forward arc disappears at the airfoil leading edge, the aft arc disappears at the airfoil trailing edge. At the juncture of the two arcs, x = x, continuity of slope is maintained. The effect of varying the maximum additional thickness and its chordwise location on airfoil lift coefficient, pitching moment, and pressure distribution was investigated. Results were obtained at a Mach number of 0.2 with an angle-of-attack of 6 degrees on the basic NACA 64-206 airfoil, and all calculations employ the full potential flow equations for two dimensional flow. The relaxation method of Jameson was employed for solution of the potential flow equations.
Merz, A. W.↗
Theoretical effect of modifications to the upper surface of two NACA airfoils using smooth polynomial additional thickness distributions which emphasize leading edge profile and which vary quadratically at the trailing edge
An investigation was conducted on a CDC 7600 digital computer to determine the effects of additional thickness distributions to the upper surface of the NACA 64-206 and 64 sub 1 - 212 airfoils. The additional thickness distribution had the form of a continuous mathematical function which disappears at both the leading edge and the trailing edge. The function behaves as a polynomial of order epsilon sub 1 at the leading edge, and a polynomial of order epsilon sub 2 at the trailing edge. Epsilon sub 2 is a constant and epsilon sub 1 is varied over a range of practical interest. The magnitude of the additional thickness, y, is a second input parameter, and the effect of varying epsilon sub 1 and y on the aerodynamic performance of the airfoil was investigated. Results were obtained at a Mach number of 0.2 with an angle-of-attack of 6 degrees on the basic airfoils, and all calculations employ the full potential flow equations for two dimensional flow. The relaxation method of Jameson was employed for solution of the potential flow equations.
Merz, A. W.↗
Theoretical effect of modifications to the upper surface of two NACA airfoils using smooth polynomial additional thickness distributions which emphasize leading edge profile and which vary linearly at the trailing edge
An investigation was conducted on a CDC 7600 digital computer to determine the effects of additional thickness distributions to the upper surface of airfoils. The additional thickness distribution had the form of a continuous mathematical function which disappears at both the leading edge and the trailing edge. Results were obtained at a Mach number of 0.2 with an angle of attack of 6 deg. All calculations employed the full potential flow equations for two dimensional flow. The relaxation method of Jameson was used for solution of the potential flow equations. It is shown that increasing the thickness and variations in shape increases the lift and the adverse pitching moment coefficients.
Hague, D. S.↗
A two-point difference scheme for computing steady-state solutions to the conservative one-dimensional Euler equations
An implicit finite-difference method is presented for obtaining steady-state solutions to the time-dependent, conservative Euler equations for flows containing shocks. The method uses a two-point central-difference scheme for the flux derivatives with dissipation added at supersonic points via the retarded density concept. Application of the method to 1-dimensional nozzle flow equations for various combinations of subsonic and supersonic boundary conditions show the method to be very efficient. Residuals are typically reduced to machine zero in approximately 35 time steps for 50 mesh points. For 1-dimensional Euler calculations, it is shown that the scheme offers two advantages over the more widely-used three-point schemes. The first is in regard to application of boundary conditions, and the second relates to the fact that the two-point algorithm is well-conditioned for large time steps.
Wornom, S. F.↗
Some Experiences with High Order Discretization Schemes for Compressible Fluid Flow
It is widely recognized that the structure and complexity of compressible fluid flow places high demands on numerical discretization techniques for the fluid flow equations. Fluid flows arising in external aerodynamics often contain both flow field discontinuities and fluid boundary-layers. Both must be accurately resolved to provide useful information to aerodynamic design and analysis engineers. These accuracy requirements motivated the present author to examine a class of finite-volume techniques on arbitrary triangulated domains based on linear or quadratic reconstruction of integral-averaged data followed by upwind flux function evaluation and small time evolution. More recently, we have considered some new upwind techniques which yield compact discretizations while maintaining higher order accuracy. In the mini-symposium talk we will discuss both of these techniques as well as demonstrate the relative merits of each method by computing a number of aerodynamic flows containing shock waves and boundary-layers.
Barth, Timothy J.↗
Simplified Modeling of Oxidation of Hydrocarbons
A method of simplified computational modeling of oxidation of hydrocarbons is undergoing development. This is one of several developments needed to enable accurate computational simulation of turbulent, chemically reacting flows. At present, accurate computational simulation of such flows is difficult or impossible in most cases because (1) the numbers of grid points needed for adequate spatial resolution of turbulent flows in realistically complex geometries are beyond the capabilities of typical supercomputers now in use and (2) the combustion of typical hydrocarbons proceeds through decomposition into hundreds of molecular species interacting through thousands of reactions. Hence, the combination of detailed reaction- rate models with the fundamental flow equations yields flow models that are computationally prohibitive. Hence, further, a reduction of at least an order of magnitude in the dimension of reaction kinetics is one of the prerequisites for feasibility of computational simulation of turbulent, chemically reacting flows. In the present method of simplified modeling, all molecular species involved in the oxidation of hydrocarbons are classified as either light or heavy; heavy molecules are those having 3 or more carbon atoms. The light molecules are not subject to meaningful decomposition, and the heavy molecules are considered to decompose into only 13 specified constituent radicals, a few of which are listed in the table. One constructs a reduced-order model, suitable for use in estimating the release of heat and the evolution of temperature in combustion, from a base comprising the 13 constituent radicals plus a total of 26 other species that include the light molecules and related light free radicals. Then rather than following all possible species through their reaction coordinates, one follows only the reduced set of reaction coordinates of the base. The behavior of the base was examined in test computational simulations of the combustion of heptane in a stirred reactor at various initial pressures ranging from 0.1 to 6 MPa. Most of the simulations were performed for stoichiometric mixtures; some were performed for fuel/oxygen mole ratios of 1/2 and 2.
Bellan, Josette↗
Deepak Condenser Model (DeCoM)
Development of the DeCoM comes from the requirement of analyzing the performance of a condenser. A component of a loop heat pipe (LHP), the condenser, is interfaced with the radiator in order to reject heat. DeCoM simulates the condenser, with certain input parameters. Systems Improved Numerical Differencing Analyzer (SINDA), a thermal analysis software, calculates the adjoining component temperatures, based on the DeCoM parameters and interface temperatures to the radiator. Application of DeCoM is (at the time of this reporting) restricted to small-scale analysis, without the need for in-depth LHP component integrations. To efficiently develop a model to simulate the LHP condenser, DeCoM was developed to meet this purpose with least complexity. DeCoM is a single-condenser, single-pass simulator for analyzing its behavior. The analysis is done based on the interactions between condenser fluid, the wall, and the interface between the wall and the radiator. DeCoM is based on conservation of energy, two-phase equations, and flow equations. For two-phase, the Lockhart- Martinelli correlation has been used in order to calculate the convection value between fluid and wall. Software such as SINDA (for thermal analysis analysis) and Thermal Desktop (for modeling) are required. DeCoM also includes the ability to implement a condenser into a thermal model with the capability of understanding the code process and being edited to user-specific needs. DeCoM requires no license, and is an open-source code. Advantages to DeCoM include time dependency, reliability, and the ability for the user to view the code process and edit to their needs.
Patel, Deepak↗
The structure of protostellar accretion disks and the origin of bipolar flows
Equations are obtained which govern the disk-wind structure and identify the physical parameters relevant to circumstellar disks. The system of equations is analyzed in the thin-disk approximation, and it is shown that the system can be consistently reduced to a set of ordinary differential equations in z. Representative solutions are presented, and it is shown that the apparent paradox discussed by Shu (1991) is resolved when the finite thickness of the disk is taken into account. Implications of the results for the origin of bipolar flows in young stellar objects and possible application to active galactic nuclei are discussed.
Wardle, Mark↗
Finite Span Wings in Compressible Flow
Equations are developed using the source distribution method for the velocity potential function and pressure on thin wings in steady and unsteady motion. Closed form solutions are given for harmonically oscillating wings of general plan form including the effect of the wing wake. Some useful examples are presented in an appendix for arrow, semielliptical, and hexagonal plan form wings.
WINGS, ELLIPTICAL↗
Chart procedures for design velocity distri- bution
Charts for solution of radial distribution flow equations used in axial flow compressor design
RADIAL DISTRIBUTION↗
Equations of motion and two-equation turbulence model for plane or axisymmetric turbulent flows in body-oriented orthogonal curvilinear coordinates and mass-averaged dependent variables
The full Navier-Stokes time-dependent, compressible, turbulent, mean-flow equations in mass-averaged variables for plane or axisymmetric flow are presented. The equations are derived in a body-oriented, orthogonal, curvilinear coordinate system. Turbulence is modelled by a system of two equations for mass-averaged turbulent kinetic energy and dissipation rate proposed. These equations are rederived and some new features are discussed. A system of second order boundary layer equations is then derived which includes the effects of longitudinal curvature and the normal pressure gradient. The Wilcox and Chambers approach is used in considering effects of streamline curvature on turbulence phenomena in turbulent boundary layer type flows. Their two-equation turbulence model with curvature terms are rederived for the cases considered in the present report. The derived system equations serves as a basis for an investigation of problems where streamline curvature is of the order of the characteristic length in the longitudinal direction.
Sislian, J. P.↗
Integral equations for flows in wind tunnels
This paper surveys recent work on the use of integral equations for the calculation of wind tunnel interference. Due to the large number of possible physical situations, the discussion is limited to two-dimensional subsonic and transonic flows. In the subsonic case, the governing boundary value problems are shown to reduce to a class of Cauchy singular equations generalizing the classical airfoil equation. The theory and numerical solution are developed in some detail. For transonic flows nonlinear singular equations result, and a brief discussion of the work of Kraft and Kraft and Lo on their numerical solution is given. Some typical numerical results are presented and directions for future research are indicated.
Fromme, J. A.↗
Comparisons of Flutter Analyses for an Experimental Fan
Two propulsion aeroelasticity codes were used to model the aeroelastic characteristics of an experimental forward-swept fan that encountered flutter during wind tunnel testing. Both of these three-dimensional codes model the unsteady flowfield due to blade vibrations using the Navier-Stokes equations. In the first approach, the unsteady flow equations are solved using an implicit time-marching approach. In the second approach, the unsteady flow equations are converted to a harmonic balance form and solved using a pseudo-time marching method. This paper describes the flutter calculations and compares the results to experimental measurements.
Bakhle, Milind A.↗
Harmonic Balance Computations of Fan Aeroelastic Stability
A harmonic balance (HB) aeroelastic analysis, which has been recently developed, was used to determine the aeroelastic stability (flutter) characteristics of an experimental fan. To assess the numerical accuracy of this HB aeroelastic analysis, a time-domain aeroelastic analysis was also used to determine the aeroelastic stability characteristics of the same fan. Both of these three-dimensional analysis codes model the unsteady flowfield due to blade vibrations using the Reynolds-averaged Navier-Stokes (RANS) equations. In the HB analysis, the unsteady flow equations are converted to a HB form and solved using a pseudo-time marching method. In the time-domain analysis, the unsteady flow equations are solved using an implicit time-marching approach. Steady and unsteady computations for two vibration modes were carried out at two rotational speeds: 100 percent (design) and 70 percent (part-speed). The steady and unsteady results obtained from the two analysis methods compare well, thus verifying the recently developed HB aeroelastic analysis. Based on the results, the experimental fan was found to have no aeroelastic instability (flutter) at the conditions examined in this study.
Bakhle, Milind A.↗
Linear/nonlinear behavior in unsteady transonic aerodynamics
The accurate calculation of the aerodynamic forces in unsteady transonic flow requires the solution of the nonlinear flow equations. The aeroelastician, on the other hand, seeks to treat his problems (flutter, for example) by means of linear equations whenever possible. He may do this, even when the underlying flow is nonlinear, if the perturbation forces are linear over some (perhaps small) range of unsteady amplitude of motion. This paper assesses the range of parameters over which linear behavior occurs. In particular calculations are made for an NACA 64A006 airfoil oscillating in pitch over a range of amplitudes, frequencies, and Mach numbers. The primary aerodynamic method used is the well known LTRAN2 code of Ballhaus and Goorjian that provides a finite-difference solution to the low frequency, small disturbance, two-dimensional potential flow equation. Comparisons are made with linear subsonic theory, local linearization, and, for steady flow, with the full potential equation code of Bauer, Garabedian, and Korn.
Dowell, E. H.↗