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,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Dynamic delamination fracture toughness of a graphite/epoxy laminate under impact

Dynamic delamination fracture toughness in a (90/0)5s T300/934 graphite/epoxy laminate was investigated using impact loading. Delamination cracks of three different sizes were embedded at the midplane of the composite specimen. The threshold impact velocity that causes propagation of the delamination crack was used in the dynamic analysis with the finite element method. From the finite element solution, the time-history of the strain energy release rate was calculated. The critical strain energy release rate was taken to equal the maximum value of the response history.

Sun, C. T.↗
NASA NTRS1991

High-order ENO schemes applied to two- and three-dimensional compressible flow

High order essentially non-oscillatory (ENO) finite difference schemes are applied to the 2-D and 3-D compressible Euler and Navier-Stokes equations. Practical issues, such as vectorization, efficiency of coding, cost comparison with other numerical methods, and accuracy degeneracy effects, are discussed. Numerical examples are provided which are representative of computational problems of current interest in transition and turbulence physics. These require both nonoscillatory shock capturing and high resolution for detailed structures in the smooth regions and demonstrate the advantage of ENO schemes.

Shu, Chi-Wang↗
NASA NTRS1991

Thermal postbuckling of thin-walled composite stiffeners

A study is made of the thermal postbuckling response of composite stiffeners subjected to prescribed edge displacement and a temperature rise. The flanges and web of the stiffeners are modeled by using two-dimensional plate finite elements. A mixed formulation is used with the fundamental unknowns consisting of the generalized displacements and the stress resultants of the plate. A reduction method is used in conjunction with mixed finite element models for determining the postbuckling response of the stiffeners. Sensitivity derivatives are evaluated and used to study the effects of variations in the different lamination and material parameters of the stiffeners on their postbuckling response characteristics. Numerical studies are presented for anisotropic stiffeners with Zee and channel sections.

Noor, Ahmed K.↗
NASA NTRS1996

A Full Multi-Grid Method for the Solution of the Cell Vertex Finite Volume Cauchy-Riemann Equations

The system of inhomogeneous Cauchy-Riemann equations defined on a square domain and subject to Dirichlet boundary conditions is considered. This problem is discretised by using the cell vertex finite volume method on quadrilateral meshes. The resulting algebraic problem is overdetermined and the solution is defined in a least squares sense. By this approach a consistent algebraic problem is obtained which differs from the original one by O(h(exp 2)) perturbations of the right-hand side. A suitable cell-based convergent smoothing iteration is presented which is naturally linked to the least squares formulation. Hence, a standard multi-grid algorithm is reported which combines the given smoother and cell-based transfer operators. Some remarkable reduction properties of these operators are shown. A full multi-grid method is discussed which solves the discrete problem to the level of truncation error by employing one multi-grid cycle at each current level of discretisation. Experiments and applications of the full multi-grid scheme are presented.

Borzi, A.↗
NASA NTRS2011

Numerical Solution of the Radiative Transfer Equation: X-Ray Spectral Formation from Cylindrical Accretion onto a Magnetized Neutron Star

Predicting the emerging X-ray spectra in several astrophysical objects is of great importance, in particular when the observational data are compared with theoretical models. This requires developing numerical routines for the solution of the radiative transfer equation according to the expected physical conditions of the systems under study. Aims. We have developed an algorithm solving the radiative transfer equation in the Fokker-Planck approximation when both thermal and bulk Comptonization take place. The algorithm is essentially a relaxation method, where stable solutions are obtained when the system has reached its steady-state equilibrium. Methods. We obtained the solution of the radiative transfer equation in the two-dimensional domain defined by the photon energy E and optical depth of the system pi using finite-differences for the partial derivatives, and imposing specific boundary conditions for the solutions. We treated the case of cylindrical accretion onto a magnetized neutron star. Results. We considered a blackbody seed spectrum of photons with exponential distribution across the accretion column and for an accretion where the velocity reaches its maximum at the stellar surface and at the top of the accretion column, respectively. In both cases higher values of the electron temperature and of the optical depth pi produce flatter and harder spectra. Other parameters contributing to the spectral formation are the steepness of the vertical velocity profile, the albedo at the star surface, and the radius of the accretion column. The latter parameter modifies the emerging spectra in a specular way for the two assumed accretion profiles. Conclusions. The algorithm has been implemented in the XPEC package for X-ray fitting and is specifically dedicated to the physical framework of accretion at the polar cap of a neutron star with a high magnetic field (approx > 10(exp 12) G). This latter case is expected to be of typical accreting systems such as X-ray pulsars and supergiant fast X ray transients.

Fairnelli, R.↗
NASA NTRSDate not supplied

Nonlinear Homogenization of Finitely Deformed Viscoelastic-Viscoplastic Composites Using Mechanics of Structure Genome

The objective of this paper is to develop a micromechanics approach to homogenizing finitely deformed viscoelastic-viscoplastic composites using the mechanics of structure genome. The incremental constitutive relation for glassy polymers, formulated in the spatial configuration, is implemented in the present approach.This involves (1) pulling-back the constitutive model to the material configuration and (2)choosing the deformation gradient tensor and the first Piola–Kirchhoff stress tensor as the strain and the stress measures during homogenization, respectively. An Euler–Newton predictor–corrector method is developed for homogenization. Each step involves formulating a variational statement using the mechanics of structure genome, discretizing the statement in a finite-dimensional space, and solving the problem using an Euler/multilevel Newton method. The present approach is demonstrated by homogenizing fiber- and particle-reinforced composites undergoing uniaxial, biaxial, or shear deformation, at different stain rates.

Multi-scale modeling, High Strain Composites, Visc↗
NASA NTRS1978

Numerical method for boundary layers with blowing - The exponential box scheme

The paper describes a new numerical scheme based on exponential difference operator concepts combined with Keller's (1968) box scheme approach to produce a stable second-order accurate finite-difference scheme for convection-diffusion problems arising in boundary layer flows in the presence of massive injection through a porous surface. The technique is demonstrated by application to the self-similar boundary layer equations with massive blowing at the surface.

El-Mistikawy, T. M.↗
NASA NTRS2001

Computational Aeroacoustics by the Space-time CE/SE Method

In recent years, a new numerical methodology for conservation laws-the Space-Time Conservation Element and Solution Element Method (CE/SE), was developed by Dr. Chang of NASA Glenn Research Center and collaborators. In nature, the new method may be categorized as a finite volume method, where the conservation element (CE) is equivalent to a finite control volume (or cell) and the solution element (SE) can be understood as the cell interface. However, due to its rigorous treatment of the fluxes and geometry, it is different from the existing schemes. The CE/SE scheme features: (1) space and time treated on the same footing, the integral equations of conservation laws are solve( for with second order accuracy, (2) high resolution, low dispersion and low dissipation, (3) novel, truly multi-dimensional, simple but effective non-reflecting boundary condition, (4) effortless implementation of computation, no numerical fix or parameter choice is needed, an( (5) robust enough to cover a wide spectrum of compressible flow: from weak linear acoustic waves to strong, discontinuous waves (shocks) appropriate for linear and nonlinear aeroacoustics. Currently, the CE/SE scheme has been developed to such a stage that a 3-13 unstructured CE/SE Navier-Stokes solver is already available. However, in the present paper, as a general introduction to the CE/SE method, only the 2-D unstructured Euler CE/SE solver is chosen as a prototype and is sketched in Section 2. Then applications of the CE/SE scheme to linear, nonlinear aeroacoustics and airframe noise are depicted in Sections 3, 4, and 5 respectively to demonstrate its robustness and capability.

Loh, Ching Y.↗
NASA NTRS1984

Efficient linear and nonlinear heat conduction with a quadrilateral element

A method is presented for performing efficient and stable finite element calculations of heat conduction with quadrilaterals using one-point quadrature. The stability in space is obtained by using a stabilization matrix which is orthogonal to all linear fields and its magnitude is determined by a stabilization parameter. It is shown that the accuracy is almost independent of the value of the stabilization parameter over a wide range of values; in fact, the values 3, 2 and 1 for the normalized stabilization parameter lead to the 5-point finite difference, 9-point finite difference and fully integrated finite element operators, respectively, for rectangular meshes; numerical experiments reported here show that the three have identical rates of convergence in the L2 norm. Eigenvalues of the element matrices, which are needed for stability limits, are also given. Numerical applications are used to show that the method yields accurate solutions with large increases in efficiency, particularly in nonlinear problems.

Liu, W. K.↗
NASA NTRS1977

A numerical study of the temperature field in a cooled radial turbine rotor

The three dimensional temperature distribution in the cooled rotor of a radial inflow turbine is determined numerically using the finite element method. Through this approach, the complicated geometries of the hot rotor and coolant passage surfaces are handled easily, and the temperatures are determined without loss of accuracy at these convective boundaries. Different cooling techniques with given coolant to primary flow ratios are investigated, and the corresponding rotor temperature fields are presented for comparison.

Hamed, A.↗
NASA NTRS1977

Over-relaxation applied to the MacCormack finite-difference scheme

An over-relaxation procedure is applied to the MacCormack finite-difference scheme in order to reduce the computation time required to obtain a steady-state solution. The implementation of this acceleration procedure to an existing computer program using the regular MacCormack method is extremely simple and does not require additional storage. The over-relaxation procedure does not alter the steady-state solution, which is second-order accurate. The method is first applied to Burgers' equation. A stability condition and an expression for the increase in the rate of convergence are derived. The method is then applied to the calculation of the hypersonic viscous flow over a flat plate, using the complete Navier-Stokes equations, and the inviscid flow over a wedge. Reductions in computing time by factors of 3 and 1.5, respectively, are obtained by over-relaxation.

Desideri, J.-A.↗
NASA NTRS1983

Simplified method for nonlinear structural analysis

A simplified inelastic analysis computer program was developed for predicting the stress-strain history of a thermomechanically cycled structure from an elastic solution. The program uses an iterative and incremental procedure to estimate the plastic strains from the material stress-strain properties and a simulated plasticity hardening model. The simplified method was exercised on a number of problems involving uniaxial and multiaxial loading, isothermal and nonisothermal conditions, and different materials and plasticity models. Good agreement was found between these analytical results and nonlinear finite element solutions for these problems. The simplified analysis program used less than 1 percent of the CPU time required for a nonlinear finite element analysis.

Kaufman, A.↗
NASA NTRS1984

Error reduction program: A progress report

Five finite differences schemes were evaluated for minimum numerical diffusion in an effort to identify and incorporate the best error reduction scheme into a 3D combustor performance code. Based on this evaluated, two finite volume method schemes were selected for further study. Both the quadratic upstream differencing scheme (QUDS) and the bounded skew upstream differencing scheme two (BSUDS2) were coded into a two dimensional computer code and their accuracy and stability determined by running several test cases. It was found that BSUDS2 was more stable than QUDS. It was also found that the accuracy of both schemes is dependent on the angle that the streamline make with the mesh with QUDS being more accurate at smaller angles and BSUDS2 more accurate at larger angles. The BSUDS2 scheme was selected for extension into three dimensions.

Syed, S. A.↗
NASA NTRS1987

Viscous shock-layer analysis of hypersonic flows over long slender bodies

Numerical solutions from the time-steady viscous shock-layer equations are presented for the hypersonic laminar and turbulent flow of a perfect gas over long slender bodies. These results are obtained from a spatial-marching implicit finite-difference technique. Detailed comparisons have been made with other predictions and experimental data to assess the accuracy of the present numerical technique, especially for slender-body flows. Results from the present method show that coupling the normal momentum and continuity equations and the use of the Vigneron pressure condition in the subsonic nose region give quite accurate and stable results. These results compare (with some exceptions) favorably with those obtained using PNS and other VSL methods.

Gupta, R. N.↗
NASA NTRS1988

Improved method for stress and compatibility analysis of multicomponent rotating systems

An improved method of analyzing multicomponent rotating assemblies for the determination of operating stresses and component compatibility has been developed. In this method, a single finite element model is developed which contains all of the separate components in the rotating assembly. This is made possible by using gap elements to simulate the contact surfaces between components. The MARC finite element computer program is then used to perform the analysis. This improved method is less time consuming and more reliable than the conventional method of analyzing such systems. Results are presented for two different stress-compatibility analyses of a six-component axial flow compressor rotor. The results for the previously used flexibility analysis method are compared with those for the improved analysis method. The stresses predicted by each method compare quite well with each other. The predictions of the component compatibility, as well as the magnitude of the forces at the contact surfaces, also compare well for these two analysis procedures. It is therefore recommended that the improved analysis method be used to determine the stress-compatibility characterstics of multicomponent rotating systems.

Carek, Gerald A.↗
NASA NTRS1988

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.↗
NASA NTRS1989

Computation of hypersonic flow fields

The present discussion of the theory, solution methods, and development status of chemically reactive flow CFD codes pertaining to the flow field around a hypersonic vehicle encompasses the formulation of multicomponent viscous equations, grid-generation techniques, and finite-difference algorithms. The Navier-Stokes equations presented focus on such particularities of high speed vehicles as their aerodynamic configuration, lee and wake flow, bow shock, and chemistry and low density effects. Issues for which further development is desirable encompass boundary treatments, grid quality, shock capturing, and the decoupling of chemistry from flow variables.

Li, Chien-Peng↗
NASA NTRS1991

Computer modeling of inversion layer MOS solar cells and arrays

A two dimensional numerical model of the inversion layer metal insulator semiconductor (IL/MIS) solar cell is proposed by using the finite element method. The two-dimensional current flow in the device is taken into account in this model. The electrostatic potential distribution, the electron concentration distribution, and the hole concentration distribution for different terminal voltages are simulated. The results of simple calculation are presented. The existing problems for this model are addressed. Future work is proposed. The MIS structures are studied and some of the results are reported.

Ho, Fat Duen↗