Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite difference methods”

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 1,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Vibrations of cantilevered shallow cylindrical shells of rectangular planform

A cantilevered, shallow shell of circular cylindrical curvature and rectangular planform exhibits free vibration behavior which differs considerably from that of a cantilevered beam or of a flat plate. Some numerical results can be found for the problem in the previously published literature, mainly obtained by using various finite element methods. The present paper is the first definitive study of the problem, presenting accurate non-dimensional frequency parameters for wide ranges of aspect ratio, shallowness ratio and thickness ratio. The analysis is based upon shallow shell theory. Numerical results are obtained by using the Ritz method, with algebraic polynomial trial functions for the displacements. Convergence is investigated, with attention being given both to the number of terms taken for each co-ordinate direction and for each of the three components of displacement. Accuracy of the results is also established by comparison with finite element results for shallow shells and with other accurate flat plate solutions.

Leissa, A. W.↗
NASA NTRS1985

An evaluation of fracture analysis methods

The results of an experimental and predictive round robin on the applications of fracture analysis methods are presented. The objective of the round robin was to verify whether fracture analysis methods currently in use can or cannot predict failure loads on complex structural components containing cracks. Fracture results from tests on a number of compact specimens were used to make the predictions. The accuracy of the prediction methods was evaluated in terms of the variation in the ratio of predicted to experimental failure loads, and the predictions methods are ranked in order of minimum standard error. A range of applicability of the different methods was also considered in assessing their usefulness. For 7075-T651 aluminum alloy, the best methods were: the effective K sub R curve; the critical crack-tip opening displacement (CTOD) criterion using a finite element analysis; and the K sub R curve with the Dugdale model. For the 2024-T351 aluminum alloy, the best methods included: the two-parameter fracture criterion (TPFC); the CTOD parameter using finite element analysis; the K-curve with the Dugdale model; the deformation plasticity failure assessment diagram (DPFAD); and the effective K sub R curve with a limit load condition. For 304 stainless steel, the best methods were the limit load analysis; the CTOD criterion using finite-element analysis TPFC and DPFAD. Some sample experimental results are given in an appendix.

Newman, J. C., Jr.↗
NASA NTRS1987

Fracture of composites using singular finite elements

Six noded quarter point natural isoparametric triangular elements are employed to obtain displacement and stress distributions around the crack tip in unidirectional graphite epoxy. The performance of the quarter point element for predicting crack growth direction in unidirectional composite materials is studied for different loading conditions and fiber orientations. The finite element results are compared to analytical solutions of anisotropic elasticity where available. Notched Iosipescu shear specimens are analyzed to study fracture behavior under shear loading. Two methods for determining stress intensity factors are presented.

Choksi, G. N.↗
NASA NTRS1988

An elastic-plastic finite-element analysis of the J-resistance curve using a CTOD criterion

The fracture process in compact and bend specimens was simulated using a two-dimensional finite-element analysis of the J-resistance (J-R) curve and a CTOD criterion. The J-R curves were calculated from the numerical results for each specimen type using several different methods. In general, the J-R curves obtained for the bend specimens were found to be higher than those for the compact specimens, especially beyond maximum load. However, below the maximum load, the modified deformation theory of plasticity and the contour-integral J(r) method resulted in very similar J-R curves for both specimen types.

Newman, J. C., Jr.↗
NASA NTRS1989

Low-frequency sound propagation modeling over a locally-reacting boundary using the parabolic approximation

There is substantial interest in the analytical and numerical modeling of low-frequency, long-range atmospheric acoustic propagation. Ray-based models, because of frequency limitations, do not always give an adequate prediction of quantities such as sound pressure or intensity levels. However, the parabolic approximation method, widely used in ocean acoustics, and often more accurate than ray models for lower frequencies of interest, can be applied to acoustic propagation in the atmosphere. Modifications of an existing implicit finite-difference implementation for computing solutions to the parabolic approximation are discussed. A locally-reacting boundary is used together with a one-parameter impedance model. Intensity calculations are performed for a number of flow resistivity values in both quiescent and windy atmospheres. Variations in the value of this parameter are shown to have substantial effects on the spatial variation of the acoustic signal.

Robertson, J. S.↗
NASA NTRS1988

Three dimensional flow simulation with application to aeroelastic analysis

The three-dimensional flowfield about realistic launch vehicle configurations is simulated using the Reynolds-averaged Navier-Stokes equations. Turbulent mixing is accounted for by means of the two-layer Baldwin and Lomax (1978) algebraic eddy viscosity model. The Beam and Warming (1976) implicit approximate factorization algorithm is used for the solution of the finite difference equations. Applications include the study of the flowfield about a hemisphere-cylinder configuration both in the subsonic and supersonic flight regimes, and about two hammerhead payload configurations at transonic speeds. A method is also described which permits incorporating this flow solver into a complete algorithm to perform time-domain aeroelastic stability analyses. The vehicle is modeled as a free-free beam and modal superposition techniques are used for the structural-dynamic formulation. Aeroelastic analyses were performed for the hammerhead configurations.

Azevedo, Joao Luiz F.↗
NASA NTRS1990

Recent advances, trends and new perspectives via enthalpy-based finite element formulations for applications to solidification problems

The present paper describes recent advances and trends in finite element developments and applications for solidification problems. In particular, in comparison to traditional methods of approach, new enthalpy-based architectures based on a generalized trapezoidal family of representations are presented which provide different perspectives, physical interpretation and solution architectures for effective numerical simulation of phase change processes encountered in solidification problems. Various numerical test models are presented and the results support the proposition for employing such formulations for general phase change applications.

Tamma, Kumar K.↗
NASA NTRS1991

Lamb wave scattering by a surface-breaking crack in a plate

An NDE method based on finite-element representation and modal expansion has been developed for solving the scattering of Lamb waves in an elastic plate waveguide. This method is very powerful for handling discontinuities of arbitrary shape, weldments of different orientations, canted cracks, etc. The advantage of the method is that it can be used to study the scattering of Lamb waves in anisotropic elastic plates and in multilayered plates as well.

Datta, S. K.↗
NASA NTRS1994

Spectral solution of the viscous blunt body problem. 2: Multidomain approximation

We present steady solutions of high speed viscous flows over blunt bodies using a multidomain Chebyshev spectral collocation method. The region with the shock layer is divided into subdomains so that internal layers can be well-resolved. In the interiors of the subdomains, the solution is approximated by Chebyshev collocation. At interfaces between subdomains, the advective terms are upwinded and the viscous terms are treated by a penalty method. The method is applied to five flows, the Mach number range 5-25 and Reynolds number range 2,000 - 83,000, based on nose radius. Results are compared to experimental data and to a finite difference result.

Kopriva, David A.↗
NASA NTRS2001

Effect of Under-Resolved Grids on High Order Methods

There has been much discussion on verification and validation processes for establishing the credibility of CFD simulations. Since the early 1990s, many of the aeronautical and mechanical engineering related reference journals mandated that any accepted articles in numerical simulations (without known solutions to compared with) need to perform a minimum of one level of grid refinement and time step reduction. Due to the difficulty in analysis, the effect of under-resolved grids and the nonlinear behavior of available spatial discretizations, are scarcely discussed in the literature. Here, an under-resolved numerical simulation is one where the grid spacing being used is too coarse to resolve the smallest physically relevant scales of the chosen continuum governing equations that are of interest to the numerical modeler. With the advent of new developments in fourth-order or higher spatial schemes, it has become common to regard high order schemes as more accurate, reliable and require less grid points. The danger comes when one tries to perform computations with the coarsest grid possible while still hoping to maintain numerical results sufficiently accurate for complex flows, and especially, data-limited problems. On one hand, high order methods when applies to highly coupled multidimensional complex nonlinear problems might have different stability, convergence and reliability behavior than their well studied low order counterparts, especially for nonlinear schemes such as TVD, MUSCL with limiters, ENO, WENO and discrete Galerkin. On the other hand, high order methods involve more operation counts and systematic grid convergence study can be time consuming and prohibitively expansive. At the same time it is difficult to fully understand or categorize the different nonlinear behavior of finite discretizations, especially at the limits of under-resolution when different types of bifurcation phenomena might occur, depending on the combination of grid spacings, time steps, initial conditions and numerical treatments of boundary conditions.

Yee, H. C.↗
NASA NTRS2008

Factors Influencing Progressive Failure Analysis Predictions for Laminated Composite Structure

Progressive failure material modeling methods used for structural analysis including failure initiation and material degradation are presented. Different failure initiation criteria and material degradation models are described that define progressive failure formulations. These progressive failure formulations are implemented in a user-defined material model for use with a nonlinear finite element analysis tool. The failure initiation criteria include the maximum stress criteria, maximum strain criteria, the Tsai-Wu failure polynomial, and the Hashin criteria. The material degradation model is based on the ply-discounting approach where the local material constitutive coefficients are degraded. Applications and extensions of the progressive failure analysis material model address two-dimensional plate and shell finite elements and three-dimensional solid finite elements. Implementation details are described in the present paper. Parametric studies for laminated composite structures are discussed to illustrate the features of the progressive failure modeling methods that have been implemented and to demonstrate their influence on progressive failure analysis predictions.

Knight, Norman F., Jr.↗
NASA NTRSDate not supplied

Developing a Spectral Correlation based Method for mitigating Angular Mismatch Effects in Intercalibration Using a Benchmark Hyperspectral Sensor

The CLARREO Pathfinder (CPF) mission will implement a state-of-the-art intercalibration method for transferring CPF’s in-orbit Système Internationale (SI)-traceable reference to the shortwave channel of the Clouds and the Earth’s Radiance Energy System (CERES) and the reflective solar bands of the Visible Infrared Imaging Radiometer Suite (VIIRS) aboard the NOAA-20 satellite platform with a targeted intercalibration methodology uncertainty of 0.3% (k=1). In order to achieve such a high intercalibration accuracy, the CPF intercalibration measurements to be scheduled need to closely match those from CERES and VIIRS in time, space, angles, and wavelength. Despite’s CPF’s two-axis pointing capability to match its boresight line of sight to that of a target sensor, there will be finite residual differences in angular samplings of the two instruments. The potential angular mismatch between CPF and the target sensors can introduce systematic errors in the intercalibration results, and therefore needs to be corrected using a rigorously designed algorithm. The CPF intercalibration team has been developing a correction method for mitigating the impact of these angular anisotropic effects in the CPF-VIIRS and CPF-CERES intercalibration. The method explores the spectral correlation relationship between the reflected solar radiances from the same surface target that would be measured by CPF at two adjacent angles. Our studies have shown that the spectral information based on CPF measurements can be used to accurately predict the spectral radiance or reflectance difference due to a given mismatch in the viewing and solar geometry. The hyper-spectral information can also be used to provide scene stratification without using any auxiliary data, which is critical to reduce the angular correction uncertainty. The angular correction relationship can be well established using simulated CPF-like top-of-atmosphere spectral radiances observed at all sorts of viewing geometry and solar angles and for different scenes. Intensive radiative transfer simulations have been completed using a state-of-art radiative transfer model developed by the CPF team members. The implementation of the algorithm on high-fidelity event simulation data and the characterization for the angular adjustment uncertainty will be presented.

Wan Wu↗
NASA NTRS1981

Scattering of Pruppacher-Pitter raindrops at 30 GHz

Optimum design of modern ground-satellite communication systems requires the knowledge of rain-induced differential attenuation, differential phase shift, and cross polarization factors. Different available analytical techniques for raindrop scattering problems were assessed. These include: (1) geometrical theory of diffraction; (2) method of moment; (3) perturbation method; (4) point matching methods; (5) extended boundary condition method; and (6) global-local finite element method. The advantages and disadvantages of each are listed. The extended boundary condition method, which was determined to yield the most scattering results, is summarized. The scattered fields for Pruppacher-Pitter raindrops with sizes ranging from 0.5 mm to 3.5 mm at 20 C and at 30 GHz for several incidence angles are tabulated.

Yeh, C.↗
NASA NTRS1981

Determination of two-stroke engine exhaust noise by the method of characteristics

A computational technique was developed for the method of characteristics solution of a one-dimensional flow in a duct as applied to the wave action in an engine exhaust system. By using the method, it was possible to compute the unsteady flow in both straight pipe and tuned expansion chamber exhaust systems as matched to the flow from the cylinder of a small two-stroke engine. The radiated exhaust noise was then determined by assuming monopole radiation from the tailpipe outlet. Very good agreement with experiment on an operation engine was achieved in the calculation of both the third octave radiated noise and the associated pressure cycles at several locations in the different exhaust systems. Of particular interest is the significance of nonlinear behavior which results in wave steepening and shock wave formation. The method computes the precise paths on the x-t plane of a finite number of C(sub +), C(sub -) and P characteristics, thereby obtaining high accuracy in determining the tailpipe outlet velocity and the radiated noise.

Jones, A. D.↗
NASA NTRS1985

Finite element method for nonlinear forced vibrations of circular plates

Geometric nonlinearities for large amplitude free and forced vibrations of circular plates are investigated. Inplane displacement and inertia are included in the formulation. Finite element method is used. Harmonic force matrix for nonlinear forced vibration analysis is introduced and derived. Various out-of-plane and inplane boundary conditions are considered. The relations of amplitude - frequency ratio for different boundary conditions and various loads conditions are presented.

Decha-Umphai, K.↗
NASA NTRS1985

Numerical Techniques in Acoustics

This is the compilation of abstracts of the Numerical Techniques in Acoustics Forum held at the ASME's Winter Annual Meeting. This forum was for informal presentation and information exchange of ongoing acoustic work in finite elements, finite difference, boundary elements and other numerical approaches. As part of this forum, it was intended to allow the participants time to raise questions on unresolved problems and to generate discussions on possible approaches and methods of solution.

Baumeister, K. J.↗
NASA NTRS1989

Involute composite design evaluation using global design sensitivity derivatives

An optimization capability for involute structures has been developed. Its key feature is the use of global material geometry variables which are so chosen that all combinations of design variables within a set of lower and upper bounds correspond to manufacturable designs. A further advantage of global variables is that their number does not increase with increasing mesh density. The accuracy of the sensitivity derivatives has been verified both through finite difference tests and through the successful use of the derivatives by an optimizer. The state of the art in composite design today is still marked by point design algorithms linked together using ad hoc methods not directly related to a manufacturing procedure. The global design sensitivity approach presented here for involutes can be applied to filament wound shells and other composite constructions using material form features peculiar to each construction. The present involute optimization technology is being applied to the Space Shuttle SRM nozzle boot ring redesigns by PDA Engineering.

Hart, J. K.↗
NASA NTRS1989

Three-dimensional viscous rotor flow calculations using boundary-layer equations

A three-dimensional viscous-inviscid interaction analysis has been developed to predict the performance of rotors in hover and forward flight at subsonic and transonic tip speeds. The analysis solves the full-potential and boundary-layer equations by finite-difference numerical procedures. Calculations were made for several different model rotor configurations in hover and forward flight at subsonic and transonic tip speeds. The results were compared with predictions from a two-dimensional integral method and with experimental data. The comparisons show good agreement between test data and predictions.

Chen, Ching S.↗