Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “delta method”

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 595 records · Page 33

Computation of wing-vortex interaction in transonic flow using implicit finite difference algorithm

An implicit delta form finite difference algorithm for Euler equations in conservation law form was used in preliminary calculations of three dimensional wing vortex interaction. Both steady and unsteady transonic flow wing vortex interactions are computed. The computations themselves are meant to guide upcoming wind tunnel experiments of the same flow field. Various modifications to the numerical method that are intended to improve computational efficiency are also described and tested in both two and three dimensions. Combination of these methods can reduce the overall computational time by a factor of 4.

Srinivasan, G.↗

Crack tip field and fatigue crack growth in general yielding and low cycle fatigue

Fatigue life consists of crack nucleation and crack propagation periods. Fatigue crack nucleation period is shorter relative to the propagation period at higher stresses. Crack nucleation period of low cycle fatigue might even be shortened by material and fabrication defects and by environmental attack. In these cases, fatigue life is largely crack propagation period. The characteristic crack tip field was studied by the finite element method, and the crack tip field is related to the far field parameters: the deformation work density, and the product of applied stress and applied strain. The cyclic carck growth rates in specimens in general yielding as measured by Solomon are analyzed in terms of J-integral. A generalized crack behavior in terms of delta is developed. The relations between J and the far field parameters and the relation for the general cyclic crack growth behavior are used to analyze fatigue lives of specimens under general-yielding cyclic-load. Fatigue life is related to the applied stress and strain ranges, the deformation work density, crack nucleus size, fracture toughness, fatigue crack growth threshold, Young's modulus, and the cyclic yield stress and strain. The fatigue lives of two aluminum alloys correlate well with the deformation work density as depicted by the derived theory. The general relation is reduced to Coffin-Manson low cycle fatigue law in the high strain region.

Minzhong, Z.↗

Solutions for hypersonic viscous flow over delta wings

A conical Navier-Stokes solver is developed to study hypersonic flow over delta wings. The differencing scheme used is van Leer flux vector splitting with MUSCL differencing. A flux limiter which involves a second difference in pressure was found to provide sharp, nonoscillatory shocks in windward and leeward flow regions. The governing equations are cast in finite-volume form and are solved using an explicit multistage method. A von Neumann stability analysis of a one-dimensional flux split model is used to predict the stability characteristics of the conical solver. The stability behavior of the scheme is found to strongly depend on the stage coefficients of the multistage scheme and on the order of accuracy of the upwind discretization.

Ruffin, Stephen M.↗

Dynamical approach study of spurious steady-state numerical solutions of nonlinear differential equations. Part 1: The ODE connection and its implications for algorithm development in computational fluid dynamics

Spurious stable as well as unstable steady state numerical solutions, spurious asymptotic numerical solutions of higher period, and even stable chaotic behavior can occur when finite difference methods are used to solve nonlinear differential equations (DE) numerically. The occurrence of spurious asymptotes is independent of whether the DE possesses a unique steady state or has additional periodic solutions and/or exhibits chaotic phenomena. The form of the nonlinear DEs and the type of numerical schemes are the determining factor. In addition, the occurrence of spurious steady states is not restricted to the time steps that are beyond the linearized stability limit of the scheme. In many instances, it can occur below the linearized stability limit. Therefore, it is essential for practitioners in computational sciences to be knowledgeable about the dynamical behavior of finite difference methods for nonlinear scalar DEs before the actual application of these methods to practical computations. It is also important to change the traditional way of thinking and practices when dealing with genuinely nonlinear problems. In the past, spurious asymptotes were observed in numerical computations but tended to be ignored because they all were assumed to lie beyond the linearized stability limits of the time step parameter delta t. As can be seen from the study, bifurcations to and from spurious asymptotic solutions and transitions to computational instability not only are highly scheme dependent and problem dependent, but also initial data and boundary condition dependent, and not limited to time steps that are beyond the linearized stability limit.

Yee, H. C.↗

Dynamical approach study of spurious steady-state numerical solutions of nonlinear differential equations. I - The dynamics of time discretization and its implications for algorithm development in computational fluid dynamics

Spurious stable as well as unstable steady state numerical solutions, spurious asymptotic numerical solutions of higher period, and even stable chaotic behavior can occur when finite difference methods are used to solve nonlinear differential equations (DE) numerically. The occurrence of spurious asymptotes is independent of whether the DE possesses a unique steady state or has additional periodic solutions and/or exhibits chaotic phenomena. The form of the nonlinear DEs and the type of numerical schemes are the determining factor. In addition, the occurrence of spurious steady states is not restricted to the time steps that are beyond the linearized stability limit of the scheme. In many instances, it can occur below the linearized stability limit. Therefore, it is essential for practitioners in computational sciences to be knowledgeable about the dynamical behavior of finite difference methods for nonlinear scalar DEs before the actual application of these methods to practical computations. It is also important to change the traditional way of thinking and practices when dealing with genuinely nonlinear problems. In the past, spurious asymptotes were observed in numerical computations but tended to be ignored because they all were assumed to lie beyond the linearized stability limits of the time step parameter delta t. As can be seen from the study, bifurcations to and from spurious asymptotic solutions and transitions to computational instability not only are highly scheme dependent and problem dependent, but also initial data and boundary condition dependent, and not limited to time steps that are beyond the linearized stability limit.

Yee, H. C.↗

Surface Pressure Distribution at Hypersonic Speeds for Blunt Delta Wings at Angle of Attack

Surface pressures were measured over a blunt 60 deg delta wing with extended trailing edge at a Mach number of 5.7, a free-stream Reynolds number of 20,000 per inch, and angles of attack from -10 to +10 deg. Aft of four leading-edge thicknesses the pressure distributions evidenced no appreciable three-dimensional effects and were predicted qualitatively by a method described herein for calculation of pressure distribution in two-dimensional flow. Results of tests performed elsewhere on blunt triangular wings were found to substantiate the near two-dimensionality of the flow and were used to extend the range of applicability of the method of surface pressure predictions to Mach numbers of 11.5 in air and 13.3 in helium.

Creager, Marcus O.↗

A first-principles derivation of Doppler noise expected from solar wind density fluctuations

The level of Doppler noise (DN) expected from solar wind (SW) density fluctuations (DF) is derived beginning with the expression for refractive index variations. The calculation takes account of up and downlink paths and of the method actually used to produce the DN values. The usual assumptions that the DF are frozen in, that the large-scale radial variation can be separated from the DF, that the DF power spectrum is a power law with outer scale k sub 0 and that the DF are homogeneous on scales less than 2c delta t, delta t = sample time, are made. The result agrees quite well with previous observations of DN. Corrections for the finite number of points used in the actual algorithm are discussed.

Callahan, P. S.↗

Three-dimensional time-marching aeroelastic analyses using an unstructured-grid Euler method

The modifications to a 3D implicit, upwind, unstructured-grid Euler code for aeroelastic analysis of complete aircraft configurations described involve the addition of structural equations of motion for their simultaneous time integration with the governing flow equations. A detailed description of the time-marching aeroelastic procedure and comparisons with experimental data are presented, to provide an assessment of capabilities. Flutter results are shown for an isolated 45-deg sweptback wing and an SST with a fuselage, clipped delta wing, and two identical rearward-mounted nacelles. Comparisons between computed and experimental flutter characteristics show good agreement.

Rausch, Russ D.↗

Longitudinal Stability Characteristics of a 1/40-Scale Model of a Proposed Configuration of the XF-91 Airplane Measured by the Wing-Flow Method

This report presents the results of an investigation in the transonic speed range of the longitudinal stability characteristics of a proposed configuration for the Republic XF-91 airplane. The tests covered a Mach number range of 0.55 to 1.05 and a Reynolds number range from 400,000 to 1,375,000. Lift, pitching-moment, and rolling-moment characteristics of the half model and the hinge moments on the all-moving tail were measured. The downwash factor delta x epsilon / delta x alpha at the tail was determined from the pitching-moment data. A calculation of the elevator deflection and stick force required for trim was also made. It was found that the variation of force and moment coefficients was linear through the test angle-of-attack range of -1 deg to 8 deg at any Mach number; that the stability increased markedly at Mach numbers above 0.85; that the effectiveness of the tail in producing pitching moments decreased about one-third with increasing Mach numbers and that the value of the downwash factor, delta x epsilon / delta x alpha, at the tail decreased from about 0.35 at a Mach number of 0.85 to about zero at a Mach number near 0.95 and became slightly negative at higher Mach numbers. The calculated values of stick force per g and elevator deflection per g, assuming no aerodynamic balance, increased rapidly above a Mach number of 0.85.

Crane, Harold L.↗

Two Carrier Analysis of Persistent Photoconductivity in Modulation-Doped Structures

A simultaneous fit of Hall and conductivity data gives quantitative results on the carrier concentration and mobility in both the quantum well and the parallel conduction channel. In this study this method was applied to reveal several new findings on the effect of persistent photoconductivity (PPC) on free-carrier concentrations and mobilities. The increase in the two-dimensional electron-gas (2DEG) concentration is significantly smaller than the apparent one derived from single carrier analysis of the Hall coefficient. In the two types of structures investigated, delta doped and continuously doped barrier, the apparent concentration almost doubles following illumination, while analysis reveals an increase of about 20% in the 2DEG. The effect of PPC on mobility depends on the structure. For the sample with a continuously doped barrier the mobility in the quantum well more than doubles. This increase is attributed to the effective screening of the ionized donors by the large electron concentration in the barrier. In the delta doped barrier sample the mobility is reduced by almost a factor of 2. This decrease is probably caused by strong coupling between the two wells, as is demonstrated by self-consistent analysis.

Schacham, S. E.↗

Flight test evaluation of a method to determine the level flight performance of a propeller-driven aircraft

A flight test method has been developed for determining the level flight drag and propulsive efficiency of propeller-driven aircraft. The overall drag of the aircraft is expressed in terms of the measured increment of power required to overcome a corresponding known increment of drag, which is generated by a towed drogue. The simplest form of the governing equations, D = delta D SHP/delta SHP (where D is drag and SHP is the horsepower delivered to the propeller) is such that all of the parameters on the right side of the equation can be measured in flight. An evaluation of the governing equations has been performed using data generated by flight test of a Beechcraft T-34B.

Bridges, P. D.↗

Improved vortex methods for three-dimensional flows

Robust numerical methods are developed for three-dimensional incompressible vortical flows, using Lagrangian vortex elements. A successful scheme must be able to handle regions of intense vortex stretching and vortex reconnection with reasonable accuracy (without diverging). Here, consideration is given to vortex particles, also commonly called vortons or vortex sticks. The following issues are discussed: (1) use of delta-function elements and weak solutions of the vorticity equation; (2) use of smoothed elements and the choice of the smoothing function; (3) representation of viscous effects and the redistribution of element strength; and (4) conservation laws (are they satisfied?). The various proposed schemes have been tested on flows involving a strong interaction between two vortex rings.

Winckelmans, G.↗

Subsonic aerodynamic and flutter characteristics of several wings calculated by the SOUSSA P1.1 panel method

Several applications of the steady, oscillatory, and unsteady subsonic and supersonic aerodynamics (SOUSSA) computer program to wings with steady and oscillatory motion, including flutter, are discussed. The program employs a generalized Green's function to the full, time-dependent potential-flow equation to obtain an integral equation for the velocity potential at any point in a flow, even points on a body or whole bodies in a flow. Aerodynamic calculations are provided for two rectangular wings, a clipped-tip delta wing, and two swept wings with and without a fuselage. The number and distribution of the finite element panels are varied in order to demonstrate the convergence of the results. The results are shown to be close to those of lifting-surface theory, and further applications with bodies having deformities, arbitrary shapes, motions, and deformations are indicated.

Yates, E. C., Jr.↗

Application of unstructured grid methods to steady and unsteady aerodynamic problems

The purpose is to describe the development of unstructured grid methods which have several advantages when compared to methods which make use of structured grids. Unstructured grids, for example, easily allow the treatment of complex geometries, allow for general mesh movement for realistic motions and structural deformations of complete aircraft configurations which is important for aeroelastic analysis, and enable adaptive mesh refinement to more accurately resolve the physics of the flow. Steady Euler calculations for a supersonic fighter configuration to demonstrate the complex geometry capability; unsteady Euler calculations for the supersonic fighter undergoing harmonic oscillations in a complete-vehicle bending mode to demonstrate the general mesh movement capability; and vortex-dominated conical-flow calculations for highly-swept delta wings to demonstrate the adaptive mesh refinement capability are discussed. The basic solution algorithm is a multi-stage Runge-Kutta time-stepping scheme with a finite-volume spatial discretization based on an unstructured grid of triangles in 2D or tetrahedra in 3D. The moving mesh capability is a general procedure which models each edge of each triangle (2D) or tetrahedra (3D) with a spring. The resulting static equilibrium equations which result from a summation of forces are then used to move the mesh to allow it to continuously conform to the instantaneous position or shape of the aircraft. The adaptive mesh refinement procedure enriches the unstructured mesh locally to more accurately resolve the vortical flow features. These capabilities are described in detail along with representative results which demonstrate several advantages of unstructured grid methods. The applicability of the unstructured grid methodology to steady and unsteady aerodynamic problems and directions for future work are discussed.

Batina, John T.↗

Flutter analysis of low aspect ratio wings

Several very low aspect ratio flat plate wing configurations are analyzed for their aerodynamic instability (flutter) characteristics. All of the wings investigated are delta planforms with clipped tips, made of aluminum alloy plate and cantilevered from the supporting vehicle body. Results of both subsonic and supersonic NASTRAN aeroelastic analyses as well as those from another version of the program implementing the supersonic linearized aerodynamic theory are presented. Results are selectively compared with the experimental data; however, supersonic predictions of the Mach Box method in NASTRAN are found to be erratic and erroneous, requiring the use of a separate program.

Parnell, L. A.↗

Observation of radiation induced changes in stress and electrical properties in MOS devices

Strain measurements using X-ray diffraction were performed on irradiated commercial and radiation-hardened metal gate CMOS devices in addition to polysilicon gate NMOS devices. I-V curves were taken and V(ot) and V(it) were separated using the subthreshold slope method for all devices. A correlation has been shown to exist between physical strain relaxation and the electrical properties as a function of radiation dose and recovery. Data shown suggest that the physical response (strain relaxation) in the silicon at the oxide interface is a measure of the type of damage induced and the recovery mechanism. Postradiation measurements of Delta V(it) and Delta V(ot) taken immediately after irradiation support the conclusions of V. Zekeriya and T.-P. Ma (1983) and K. Kasama et al. (1986, 1987); compressive stress at the silicon/SiO2 interface does reduce radiation damage in the device.

Shaw, D. C.↗

Flight Hardware Development and Research at MSFC for Optimizing Success on the International Space Station

To optimize biological crystallization success in microgravity in-house personnel at the MSFC are working on the development of innovative flight hardware such as Delta-L and the Iterative Biological Crystallization (IBC) apparatus as well as troubleshooting the performance of existing hardware. Delta-L will provide a diagnostic hardware to examine the relationship between crystal growth characteristics and crystal quality improvement in microgravity. IBC is a new hardware being designed to allow iteration of crystal growth experiments in microgravity using innovative lab on a chip technology. While being built to obtain scientific data of benefit to the scientific community, the design methods involved in the development of these hardware have directly benefited other groups within NASA and keep NASA at the forefront of innovation.

Source record↗

Crack Instability Predictions Using a Multi-Term Approach

Present crack instability analysis for fracture critical flight hardware is normally performed using a single parameter, K(sub C), fracture toughness value obtained from standard ASTM 2D geometry test specimens made from the appropriate material. These specimens do not sufficiently match the boundary conditions and the elastic‐plastic constraint characteristics of the hardware component, and also, the crack instability of most commonly used aircraft and aerospace structural materials have some amount of stable crack growth before fracture which makes the normal use of a K(sub C) single parameter toughness value highly approximate. In the past, extensive studies have been conducted to improve the single parameter (K or J controlled) approaches by introducing parameters accounting for the geometry or in‐plane constraint effects. Using 'J‐integral' and 'A' parameter as a measure of constraint is one of the most accurate elastic‐plastic crack solutions currently available. In this work the feasibility of the J‐A approach for prediction of the crack instability was investigated first by ignoring the effects of stable crack growth i.e. using a critical J and A and second by considering the effects of stable crack growth using the corrected J‐delta a using the 'A' parameter. A broad range of initial crack lengths and a wide range of specimen geometries including C(T), M(T), ESE(T), SE(T), Double Edge Crack (DEC), Three‐Hole‐Tension (THT) and NC (crack from a notch) manufactured from Al7075 were studied. Improvements in crack instability predictions were observed compared to the other methods available in the literature.

Zanganeh, Mohammad↗