Search NASA⌕ Search

SEARCH · Search NASA

Results for “Higher order calculations”

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 73 records · Page 4

Effects of Bifurcations on Aft-Fan Engine Nacelle Noise

Aft-fan engine nacelle noise is a significant factor in the increasingly important issue of aircraft community noise. The ability to predict such noise within complex duct geometries is a valuable tool in studying possible noise attenuation methods. A recent example of code development for such predictions is the ducted fan noise propagation and radiation code CDUCT-LaRC. This work focuses on predicting the effects of geometry changes (i.e. bifurcations, pylons) on aft fan noise propagation. Beginning with simplified geometries, calculations show that bifurcations lead to scattering of acoustic energy into higher order modes. In addition, when circumferential mode number and the number of bifurcations are properly commensurate, bifurcations increase the relative importance of the plane wave mode near the exhaust plane of the bypass duct. This is particularly evident when the bypass duct surfaces include acoustic treatment. Calculations involving more complex geometries further illustrate that bifurcations and pylons clearly affect modal content, in both propagation and radiation calculations. Additionally, results show that consideration of acoustic radiation results may provide further insight into acoustic treatment effectiveness for situations in which modal decomposition may not be straightforward. The ability of CDUCT-LaRC to handle complex (non-axisymmetric) multi-block geometries, as well as axially and circumferentially segmented liners, allows investigation into the effects of geometric elements (bifurcations, pylons).

Nark, Douglas M.↗

Improved solution for potential flow about arbitrary axisymmetric bodies by the use of a higher-order surface source method. Part 1. Theory and results

The surface-source method of calculating potential flow is improved by refining the underlying numerical analysis. The analysis uses parabolic elements and linearly-varying source density. The result is a large increase in computing speed and accuracy. The theory is described, and the effectiveness of the modification is illustrated.

Hess, J. L.↗

Search for extreme ultraviolet radiation from selected stars

Results are reported for a search for EUV radiation between 140 and 430 A from selected objects which was conducted with a rocket-borne EUV filter photometer. The candidate objects included the bright planetary nebula NGC 2392 and the stars gamma Gem, alpha CMa, omicron-prime CMa, and gamma-2 Vel. Upper limits (2-sigma level) to the flux near 200 A are determined for the five objects, and estimates of the expected 200-A helium-continuum flux are derived for two of the stars on the basis of the hot-corona model. It is found that the calculated upper limits are at least two orders of magnitude higher than plausible flux estimates based on the hot-corona model. However, it is noted that uncertainties in the flux estimates are great enough to make the upper limits useful.

Riegler, G. R.↗

Higher-order derivative correlations and the alignment of small-scale structures in isotropic numerical turbulence

In a three dimensional simulation higher order derivative correlations, including skewness and flatness factors, are calculated for velocity and passive scalar fields and are compared with structures in the flow. The equations are forced to maintain steady state turbulence and collect statistics. It is found that the scalar derivative flatness increases much faster with Reynolds number than the velocity derivative flatness, and the velocity and mixed derivative skewness do not increase with Reynolds number. Separate exponents are found for the various fourth order velocity derivative correlations, with the vorticity flatness exponent the largest. Three dimensional graphics show strong alignment between the vorticity, rate of strain, and scalar-gradient fields. The vorticity is concentrated in tubes with the scalar gradient and the largest principal rate of strain aligned perpendicular to the tubes. Velocity spectra, in Kolmogorov variables, collapse to a single curve and a short minus 5/3 spectral regime is observed.

Kerr, R. A.↗

Using a microgravity environment to probe wave turbulence

The experimental key to observing stochasticity or turbulence in a distribution of interacting propagating waves is the achievement of high amplitude and the use of a medium with a large coefficient of nonlinearity. The research indicates that capillary waves are the best means of observing this phenomenon; however, gravitational modifications of the capillary wave dispersion law greatly reduce the large coefficient of nonlinearity. Thus, a search for wave turbulence in a large drop of fluid that is positioned in a microgravity experiment was conducted. Capillary waves that run around the surface of the drop are excited, and their power spectrum and higher order correlations are analyzed for wave turbulence. The theoretical calculations indicate that modulations of the power spectrum should propagate as second sound waves. These issues have consequences for signal processing and plasma confinement.

Putterman, S.↗

A symmetric modal formulation of fluid-structure interaction, including a static approximation to higher order fluid modes

A method of fluid-structure coupling which provides symmetrical matrix equations of standard form solved by existing finite element computer programs is presented. The method postulates that the uncoupled vibration modes of the fluid or the structure be calculated before the coupled analysis. A numerical solution of vibration modes in an axisymmetric container demonstrated that a static approximation to higher order fluid modes can improve the accuracy of dynamic response computations using modal methods.

Macneal, R. H.↗

Application of higher-order numerical methods to the boundary-layer equations

A fourth-order method is presented for calculating numerical solutions to parabolic, partial differential equations in two variables or ordinary differential equations. The method is the natural extension of the second-order Keller Box Scheme to fourth order and is demonstrated with application to the incompressible, laminar and turbulent boundary-layer equations for both attached and separated flows. The efficiency of the present method is compared with other higher-order methods; namely, the Keller Box Scheme with Richardson extrapolation, the method of deferred corrections, the three-point spline methods, and a modified finite-element method. For equivalent accuracy, numerical results show the present method to be more efficient than the other higher-order methods for both laminar and turbulent flows.

Wornom, S. F.↗

Apparatus, Method and Program Storage Device for Determining High-Energy Neutron/Ion Transport to a Target of Interest

An apparatus, method and program storage device for determining high-energy neutron/ion transport to a target of interest. Boundaries are defined for calculation of a high-energy neutron/ion transport to a target of interest; the high-energy neutron/ion transport to the target of interest is calculated using numerical procedures selected to reduce local truncation error by including higher order terms and to allow absolute control of propagated error by ensuring truncation error is third order in step size, and using scaling procedures for flux coupling terms modified to improve computed results by adding a scaling factor to terms describing production of j-particles from collisions of k-particles; and the calculated high-energy neutron/ion transport is provided to modeling modules to control an effective radiation dose at the target of interest.

Wilson, John W.↗

Critical study of higher order numerical methods for solving the boundary-layer equations

A fourth order box method is presented for calculating numerical solutions to parabolic, partial differential equations in two variables or ordinary differential equations. The method, which is the natural extension of the second order box scheme to fourth order, was demonstrated with application to the incompressible, laminar and turbulent, boundary layer equations. The efficiency of the present method is compared with two point and three point higher order methods, namely, the Keller box scheme with Richardson extrapolation, the method of deferred corrections, a three point spline method, and a modified finite element method. For equivalent accuracy, numerical results show the present method to be more efficient than higher order methods for both laminar and turbulent flows.

Wornom, S. F.↗

Investigations of medium wavelength magnetic anomalies in the eastern Pacific using MAGSAT data

The computer program package for dealing with MAGSAT data was refined and a program listing and deck of cards was forwarded for use by interested investigators. Instead of doing equivalent source calculations, it is proposed to use spherical harmonics to model the observed data since spherical harmonic functions are orthogonal and addition of higher degree and order harmonics can be done without recalculation of the harmonics of lower degree and order already calculated. These functions are much more easy to manipulate in that the height of the term of the observation can be taken care of with great facility. Downward calculation (after suitable filtering) can also be performed very simply.

Harrison, C. G. A.↗

A comparison between onestep and other multiconic trajectory propagation methods

The paper presents an optimized version of the one-step multiconic method of spacecraft trajectory propagation in the presence of two attracting bodies. This method is compared to eleven other methods ranging from simple conic propagation to full integration of the equations of motion by applying the methods to a wide variety of test cases taken from planned space missions. This optimize one-step method is shown to be an order of magnitude better than other single-step methods for interplanetary and most lunar trajectories. Multistep multiconics methods and integration of the state offer another order of magnitude improvement, but at a proportionally higher cost and calculation time.

Kiedron, Krystyna↗

Dynamic Analysis and Control of Lightweight Manipulators with Flexible Parallel Link Mechanisms

The objective is the theoretical analysis and the experimental verification of dynamics and control of a two link flexible manipulator with a flexible parallel link mechanism. Nonlinear equations of motion of the lightweight manipulator are derived by the Lagrangian method in symbolic form to better understand the structure of the dynamic model. The resulting equation of motion have a structure which is useful to reduce the number of terms calculated, to check correctness, or to extend the model to higher order. A manipulator with a flexible parallel link mechanism is a constrained dynamic system whose equations are sensitive to numerical integration error. This constrained system is solved using singular value decomposition of the constraint Jacobian matrix. Elastic motion is expressed by the assumed mode method. Mode shape functions of each link are chosen using the load interfaced component mode synthesis. The discrepancies between the analytical model and the experiment are explained using a simplified and a detailed finite element model.

Lee, Jeh Won↗

Toward a formal verification of a floating-point coprocessor and its composition with a central processing unit

Discussed here is work to formally specify and verify a floating point coprocessor based on the MC68881. The HOL verification system developed at Cambridge University was used. The coprocessor consists of two independent units: the bus interface unit used to communicate with the cpu and the arithmetic processing unit used to perform the actual calculation. Reasoning about the interaction and synchronization among processes using higher order logic is demonstrated.

Pan, Jing↗

Parametric studies of stitching effectiveness for preventing substructure disbond

A methodology is desired that will allow a designer to select appropriate amounts of through-thickness reinforcement needed to meet design requirements. The goal is to use a relatively simple analysis to minimize the amount of testing that needs to be performed, and to make test results from simple configurations applicable to more general structures. Using this methodology one should be able to optimize the selection of stitching materials, the weight of the yarn, and the stitching density. The analysis approach is to treat substructure disbond as a crack propagation problem. In this approach, the stitches have little influence until a delamination begins to grow. Once the delamination reaches, or extends beyond a stitch, the stitch serves to reduce the strain-energy-release-rate (G) at the crack tip for a given applied load. The reduced G can then be compared to the unstitched materials toughness to predict the load required to further extend the crack. The current model treats the stitch as a simple spring which responds to displacements in the vertical (through-thickness) direction. In concept, this approach is similar to that proposed by other authors. Test results indicate that the model should be refined to include the shearing stiffness of the stitch. The strain-energy-release-rate calculations are performed using a code which uses interconnected higher-order plates to model built-up composite cross-sections. When plates are stacked vertically, the interfacial tractions between the plates can be computed. The plate differential equations are solved in closed-form. The code, called SUBLAM, was developed as part of this section in one dimension. Because of this limitation, rows of stitches are treated as a two-dimensional sheet. The spring stiffness of a row of stitches can be estimated from the stitch material, weight, and density. As a practical and conservative approach, we can assume that the stitch is bonded until a crack passes the stitch location. After the crack passes, it is fully bonded. A series of tests were performed to exercise this methodology and incorporated an attached flange such that the sudden change in thickness initiated a delamination. The analysis was used to estimate the materials' critical G from that of the unstitched specimens. With this data, a prediction was made for the load required to delaminate the stitched specimens. Using the methodology, design charts have been created for simplified geometries. These charts give stitch force, along with G(sub 1) and G(sub 2) as as function of the stitch spring stiffness. Using the charts, it should be possible to determine the stitch spring stiffness and strength required to reduce the G to a desired level. From these parameters, the actual stitching material, weight, and density can be computed.

Flanagan, Gerry↗

Higher order moments of the matter distribution in scale-free cosmological simulations with large dynamic range

We calculate reduced moments (xi bar)(sub q) of the matter density fluctuations, up to order q = 5, from counts in cells produced by particle-mesh numerical simulations with scale-free Gaussian initial conditions. We use power-law spectra P(k) proportional to k(exp n) with indices n = -3, -2, -1, 0, 1. Due to the supposed absence of characteristic times or scales in our models, all quantities are expected to depend on a single scaling variable. For each model, the moments at all times can be expressed in terms of the variance (xi bar)(sub 2), alone. We look for agreement with the hierarchical scaling ansatz, according to which ((xi bar)(sub q)) proportional to ((xi bar)(sub 2))(exp (q - 1)). For n less than or equal to -2 models, we find strong deviations from the hierarchy, which are mostly due to the presence of boundary problems in the simulations. A small, residual signal of deviation from the hierarchical scaling is however also found in n greater than or equal to -1 models. The wide range of spectra considered and the large dynamic range, with careful checks of scaling and shot-noise effects, allows us to reliably detect evolution away from the perturbation theory result.

Lucchin, Francesco↗

A critical study of higher-order numerical methods for solving the boundary-layer equations

A fourth-order box method is presented for calculating numerical solutions to parabolic, partial differential equations in two variables or ordinary differential equations. The method is the natural extension of the second-order Keller Box Scheme to fourth order and is demonstrated with application to the incompressible, laminar and turbulent boundary-layer equations. The efficiency of the present method is compared with other two-point and three-point higher-order methods; namely, the Keller Box Scheme with Richardson extrapolation, the method of deferred corrections, and the three-point spline methods. For equivalent accuracy, numerical results show the present method to be more efficient than the other higher-order methods for both laminar and turbulent flows.

Wornom, S. F.↗

Postbuckling response of long thick isotropic plates loaded in compression including higher order transverse shearing effects

Buckling and postbuckling results for aluminum plates loaded in compression are presented. The buckling results were plotted to show the effects of thickness on the stress coefficient. Buckling results are given for various length-to-width ratios. Postbuckling results for plates with transverse shearing flexibility are compared to results from classical theory for various width-to-thickness ratios. The plates are considered to be long with side edges simply supported, with edges free of stress and the plates are subjected to longitudinal compressive displacement. Characteristic curves indicating the average longitudinal direct stress resultant as a function of the applied displacements are calculated based on four different theories: Classical von Karman, first-order shear deformation, higher-order shear deformation, and three-dimensional flexibility. Present results indicate that the three-dimensional flexibility theory gives the lowest and therefore, most accurate results. The higher-order shear deformation theory has fewer unknowns than the three-dimensional flexibility but is not as accurate. The figures presented show that small differences occur in the maximum stress resultants and the transverse displacements calculated when the effects of transverse shear are included.

Stein, Manuel↗

Corrections to vibrational transition probabilities calculated from a three-dimensional model.

Corrections to the collision-induced vibration transition probability calculated by Hansen and Pearson from a three-dimensional semiclassical model are examined. These corrections come from the retention of higher order terms in the expansion of the interaction potential and the use of the actual value of the deflection angle in the calculation of the transition probability. It is found that the contribution to the transition cross section from previously neglected potential terms can be significant for short range potentials and for the large relative collision velocities encountered at high temperatures. The correction to the transition cross section obtained from the use of actual deflection angles will not be appreciable unless the change in the rotational quantum number is large.

Stallcop, J. R.↗