SEARCH · Search NASA
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.
Rudy and Bushnell's (1972) mixing length theory for the prediction of mean velocity fields in turbulent shear layers is extended to include supersonic free turbulent mixing. The normal momentum equation is coupled with the conventional equations of motion and is solved iteratively. Results show that transverse static pressure variation has very little effect on the mean flow variables in compressible free turbulent shear layers.
The service life of aircraft structural components undergoing random stress cycling was analyzed by the application of fracture mechanics. The initial crack sizes at the critical stress points for the fatigue-crack growth analysis were established through proof load tests. The fatigue-crack growth rates for random stress cycles were calculated using the half-cycle method. A new equation was developed for calculating the number of remaining flights for the structural components. The number of remaining flights predicted by the new equation is much lower than that predicted by the conventional equation.
The service life of aircraft structural components undergoing random stress cycling was analyzed by the application of fracture mechanics. The initial crack sizes at the critical stress points for the fatigue crack growth analysis were established through proof load tests. The fatigue crack growth rates for random stress cycles were calculated using the half-cycle method. A new equation was developed for calculating the number of remaining flights for the structural components. The number of remaining flights predicted by the new equation is much lower than that predicted by the conventional equation. This report describes the application of fracture mechanics and the half-cycle method to calculate the number of remaining flights for aircraft structural components.
A general equation governing aerodynamic sound generation in the presence of solid boundaries is derived. It is shown that all the theories in the literature appear as special cases of this general equation. Derived special equations for propeller and fan noise are likewise shown to be more general than the conventional equations in that they make allowance for variation in retarded time over the blade surfaces.
Measurements of the drag and of the nacelle internal pressures on a wing and nacelle that housed a horizontally opposed piston engine were made in the 40- by 80-Foot Wind Tunnel at Ames Research Center. These tests are follow-ons to earlier tests made with the same wing and nacelle but in which the engine was replaced with an electric motor and an adjustable orifice plate. In the initial tests the orifice plate was used to control the rate of cooling-air flow through the nacelle and thereby to simulate a range of gasoline engine types. Good agreement was found between the results of those tests and of the test reported here. Also, the upper and lower plenum pressure and cooling-air flow rate were found to be related by conventional equations used to represent the flow through orifices. Tests were run with three cooling air inlet sizes over a free-stream velocity range from 50 to 150 knots, an angle of attack range from 0 deg to 10 deg, and a cowl-flap deflection range from 0 deg to 30 deg. The data were analyzed by computing a flow coefficient similar to that used in the analysis of orifices. It was found that all of the flow coefficient values fell within a band that varied linearly with inlet area. The linear mean line through this band provides an estimate of the relationship between cooling-air flow rate and upper plenum pressure over a wide range of test conditions.
Conventionally kinematical constraints in multibody systems are treated similar to geometrical constraints and are modeled by constraint reaction forces which are perpendicular to constraint surfaces. However, in reality, one may want to achieve the desired kinematical conditions by control forces having different directions in relation to the constraint surfaces. The conventional equations of motion for multibody systems subject to kinematical constraints are generalized by introducing general direction control forces. Conditions for the selections of the control force directions are also discussed. A redundant robotic system subject to prescribed end-effector motion is analyzed to illustrate the methods proposed.
The predictive capability of computational fluid dynamics (CFD) codes for turbulent flow through curved ducts is of significant importance to the design and performance analysis of modern rocket engine flowpaths. Code calibration and validation studies for this class of flow are desireable to estimate the performance margin and operating range of components designed using Navier-Stokes methods. Parametric experimental studies such as that of Weske (NACA ARR W-39) provided a wealth of performance data for the design of single- and compound elbow configurations with various cross-sections, curvature and aspect ratios at varying Reynolds numbers. In that work, the majority of data is presented in the form of loss coefficients, characterizing pressure losses due to duct curvature, and including losses due to wall friction. Using measured friction coefficients, losses of equivalent straight lengths of duct are subtracted, resulting in performance curves useful for design computations. These data are currently used in a CFD-based parametric study covering a broad range of operating conditions. Of particular interest for the accuracy of CFD predictions are the effects on pressure loss due to inlet boundary layer thickness (dependent on upstream development length), and the wall treatment for the turbulence equations (conventional wall functions vs. wall integration using a two-layer model). The experimental data are reassessed in the form of an error analysis, and are compared with CFD predictions for 18 computational cases. Grid-independence, grid spacing, and convergence requirements of the cases are discussed. Conclusions regarding the relative importance of the parametric variables will be presented.
Bending actuators of a proposed type would partly resemble ordinary bending actuators, but would include simple additional components that would render them capable of exerting large forces at small displacements. Like an ordinary bending actuator, an actuator according to the proposal would include a thin rectangular strip that would comprise two bonded layers (possibly made of electroactive polymers with surface electrodes) and would be clamped at one end in the manner of a cantilever beam. Unlike an ordinary bending actuator, the proposed device would include a rigid flat backplate that would support part of the bending strip against backward displacement; because of this feature, the proposed device is called a backed bending actuator. When an ordinary bending actuator is inactive, the strip typically lies flat, the tip displacement is zero, and the force exerted by the tip is zero. During activation, the tip exerts a transverse force and undergoes a bending displacement that results from the expansion or contraction of one or more of the bonded layers. The tip force of an ordinary bending actuator is inversely proportional to its length; hence, a long actuator tends to be weak. The figure depicts an ordinary bending actuator and the corresponding backed bending actuator. The bending, the tip displacement (d(sub t)), and the tip force (F) exerted by the ordinary bending actuator are well approximated by the conventional equations for the loading and deflection of a cantilever beam subject to a bending moment which, in this case, is applied by the differential expansion or contraction of the bonded layers. The bending, displacement, and tip force of the backed bending actuator are calculated similarly, except that it is necessary to account for the fact that the force F(sub b) that resists the displacement of the tip could be sufficient to push part of the strip against the backplate; in such a condition, the cantilever beam would be effectively shortened (length L*) and thereby stiffened and, hence, made capable of exerting a greater tip force for a given degree of differential expansion or contraction of the bonded layers. Taking all of these effects into account, the cantilever-beam equations show that F(sub b) would be approximately inversely proportional to d(sup 1/2) for d less than a calculable amount, denoted the transition displacement (dt). For d less than d(sub t), part of the strip would be pressed against the backplate. Therefore, the force F(sub b) would be very large for d at or near zero and would decrease as d increases toward d(sub t). At d greater than d(sub t), none of the strip would be pressed against the backplate and F(sub b) would equal the tip force F of the corresponding ordinary bending actuator. The advantage of the proposal is that a backed bending actuator could be made long to obtain large displacement when it encountered little resistance but it could also exert a large zero-displacement force, so that it could more easily start the movement of a large mass, throw a mechanical switch, or release a stuck mechanism.
The n-dimensional scattering problem is solved by means of a k-space formulation of the field equations, thereby replacing the conventional integral equation formulation by a set of two algebraic equations in two unknowns in two spaces (the constitutive equation being an algebraic equation in x-space). These equations are solved by an iterative method with the aid of the fast Fourier transform (FFT) algorithm connecting the two spaces, requiring very simple initial approximations. Since algebraic and FFT equations are used, the number of arithmetic multiple-add operations and storage allocations required for a numerical solution are reduced from the order of N sq (for solving the matrix equations resulting from the conventional integral equations) to the order of N(log base 2 of N) and N, respectively (where N is the number of data points required for the specification of the problem). The advantage gained in speed and storage is thus of the order of N/log base 2 of N and N, respectively. This method is thus considerably more efficient than the conventional matrix method, and permits exact numerical solutions for much larger problems. Arguments are presented toward the view that the field equations are more fundamental in k-space. The details and some numerical results of the application of this method to the three-dimensional electromagnetic scattering problems are presented as an example.
Numerical predictions, flow visualization experiments and time-mean velocity measurements were obtained for six basic nonreacting flowfields (with inlet swirl vane angles of 0 (swirler removed), 45 and 70 degrees and sidewall expansion angles of 90 and 45 degrees) in an idealized axisymmetric combustor geometry. A flowfield prediction computer program was developed which solves appropriate finite difference equations including a conventional two equation k-epsilon eddy viscosity turbulence model. The wall functions employed were derived from previous swirling flow measurements, and the stairstep approximation was employed to represent the sloping wall at the inlet to the test chamber. Recirculation region boundaries have been sketched from the entire flow visualization photograph collection. Tufts, smoke, and neutrally buoyant helium filled soap bubbles were employed as flow tracers. A five hole pitot probe was utilized to measure the axial, radial, and swirl time mean velocity components.
A mathematical representation has been developed for the electromagnetic force field, fluid flow field, and solute concentration field of levitation-melted metal specimens. The governing equations consist of the conventional transport equations combined with the appropriate expressions for the electromagnetic force field. The predictions obtained by solving the governing equations numerically on a digital computer are in good agreement with lifting force and average temperature measurements reported in the literature.
We implemented a two- and three-dimensional finite difference/front tracking technique to solve liquid-vapor phase change problems. The mathematical and the numerical features of the method were explained in great detail in our previous reports, Briefly, we used a single formula representation which incorporated jump conditions into the governing equations. The interfacial terms were distributed as singular terms using delta functions so that the governing equations would be the same as conventional conservation equations away from the interface and in the vicinity of the interface they would provide correct jump conditions. We used a fixed staggered grid to discretize these equations and an unstructured grid to explicitly track the front. While in two dimensions the front was simply a connection of small line segments, in three dimensions it was represented by a connection of small triangular elements. The equations were written in conservative forms and during the course of computations we used regriding to control the size of the elements of the unstructured grid. Moreover, we implemented a coalescence in two dimensions which allowed the merging of different fronts or two segments of the same front when they were sufficiently close. We used our code to study thermocapillary migration of bubbles, burst of bubbles at a free surface, buoyancy-driven interactions of bubbles, evaporation of drops, rapid evaporation of an interface, planar solidification of an undercooled melt, dendritic solidification, and a host of other problems cited in the reference.
An equation that relates the pressure flow characteristics of hydraulic spool valves was developed. The dependent variable is valve output pressure, and the independent variables are spool position and flow. This causal form of equation is preferred in applications that simulate the effects of hydraulic line dynamics. Results from this equation are compared with those from the conventional valve equation, whose dependent variable is flow. A computer program of the valve equations includes spool stops, leakage spool clearances, and dead-zone characteristics of overlap spools.
The region of validity of conventional rate equation formulations used to describe the excitation of flame gases for spectroscopic diagnostic purposes is examined by direct comparison with the density matrix formulation. The density matrix equations are formulated for an ensemble of atoms or molecules in the gas state. These equations are reduced to the steady-state limit for a two-level system under the effect of single-mode and multimode optical fields. The steady-state equations are shown to be valid only when the laser pulse rises slowly with respect to the collisional dephasing rate for single-mode excitation, and for certain special cases under multimode excitation.
Nonlinear strain displacement relations for three-dimensional elasticity are determined in orthogonal curvilinear coordinates. To develop a two-dimensional theory, the displacements are expressed by trigonometric series representation through-the-thickness. The nonlinear strain-displacement relations are expanded into series which contain all first and second degree terms. In the series for the displacements only the first few terms are retained. Insertion of the expansions into the three-dimensional virtual work expression leads to nonlinear equations of equilibrium for laminated and thick plates and shells that include the effects of transverse shearing. Equations of equilibrium and buckling equations are derived for flat plates and cylindrical shells. The shell equations reduce to conventional transverse shearing shell equations when the effects of the trigonometric terms are omitted and to classical shell equations when the trigonometric terms are omitted and the shell is assumed to be thin.
The paper describes the application of techniques based on computational fluid dynamics to the simulation of jet flowfields. A solution code for the Reynolds-averaged Navier-Stokes equations is supplemented by conventional two-equation turbulence models based on the Boussinesq approximation. The axisymmetric SCIPVIS code is enhanced with the PARCH and CRAFT codes to examine plug-jet flowfields and imperfectly expanded axisymmetric free round jets. The sensitivity of shock/boundary layer interactions is observed in simulations of the plug case, and the adaptive gridding in the disk region and turbulence levels generated at the triple point are identified as areas in the Mach case that require improvement. Jet-wave structure in the region beyond the first several shock cells can be predicted, and turbulence modeling can be undertaken with respect to improving compressibility, length scale, vorticity, and energy budget. The mean flow structure of imperfectly expanded jets can be studied to develop related noise suppression concepts for the High-Speed Civilian Transport (HSCT).
Photoionization plays a central role in the development of streamer discharges and other non-equilibrium plasma phenomena. It creates seed electrons, which are essential for positive streamer propagation, allowing the ionization front to move forward. Because of this, accurate modeling of photoionization is very important for predicting streamer behavior and plasma evolution. The photoionization process in air (N 2 – O 2 mixture) is often described by the Zheleznyak model (1982). This model is usually solved through Helmholtz-type equations that approximate the Zheleznyak photoionization model (Zheleznyak et al. 1982) as Partial Differential Equations (PDEs). Conventional numerical methods, such as the Finite Difference Method (FDM) or Finite Volume Method (FVM), are widely used to solve these equations. Although they are prevalent, the computational cost of these methods due to their need for matrix operations and iterative solver is demanding. To address this challenge, this work develops a spectral solver based on the Fast Fourier Transform (FFT) combined with Discrete Cosine Transform (DCT) and Discrete Sine Transform (DST) to calculate the photoionization rate efficiently in an axisymmetric cylindrical domain. This method naturally satisfies the boundary conditions used in the model and converts the PDE into algebraic ones in spectral space. Thus, avoids the need for iterative matrix solvers. When compared with FDM results, it is demonstrated that the new solver not only maintains accuracy, but also reduces the computational cost, showing a performance increase of approximately 100 compared to FDM over a wide range of problem sizes. The method is parallelized using Message Passing Interface (MPI) and has been integrated into a fluid plasma model for streamer simulation. Here, this FFT-based approach provides a fast and reliable alternative for calculating photoionization in fluid models, helping large-scale plasma simulations run faster and efficiently, and allows higher-resolution simulation without extra computational cost.
This paper is one of a series of papers describing the development of a new numerical method for the Navier-Stokes equations. Unlike conventional numerical methods, the current method concentrates on the discrete simulation of both the integral and differential forms of the Navier-Stokes equations. Conservation of mass, momentum, and energy in space-time is explicitly provided for through a rigorous enforcement of both the integral and differential forms of the governing conservation laws. Using local polynomial expansions to represent the discrete primitive variables on each cell, fluxes at cell interfaces are evaluated and balanced using exact functional expressions. No interpolation or flux limiters are required. Because of the generality of the current method, it applies equally to the steady and unsteady Navier-Stokes equations. In this paper, we generalize and extend the authors' 2-D, steady state implicit scheme. A general closure methodology is presented so that all terms up through a given order in the local expansions may be retained. The scheme is also extended to nonorthogonal Cartesian grids. Numerous flow fields are computed and results are compared with known solutions. The high accuracy of the scheme is demonstrated through its ability to accurately resolve developing boundary layers on coarse grids. Finally, we discuss applications of the current method to the unsteady Navier-Stokes equations.