Search NASA⌕ Search

SEARCH · Search NASA

Results for “generalized finite differences”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1978

An evaluation of the sandwich beam in four-point bending as a compressive test method for composites

The experimental phase of the study included compressive tests on HTS/PMR-15 graphite/polyimide, 2024-T3 aluminum alloy, and 5052 aluminum honeycomb at room temperature, and tensile tests on graphite/polyimide at room temperature, -157 C, and 316 C. Elastic properties and strength data are presented for three laminates. The room temperature elastic properties were generally found to differ in tension and compression with Young's modulus values differing by as much as twenty-six percent. The effect of temperature on modulus and strength was shown to be laminate dependent. A three-dimensional finite element analysis predicted an essentially uniform, uniaxial compressive stress state in the top flange test section of the sandwich beam. In conclusion, the sandwich beam can be used to obtain accurate, reliable Young's modulus and Poisson's ratio data for advanced composites; however, the ultimate compressive stress for some laminates may be influenced by the specimen geometry.

Shuart, M. J.↗
NASA NTRS1979

A program to compute three-dimensional subsonic unsteady aerodynamic characteristics using the doublet lattice method, L216 (DUBFLEX). Volume 2: Supplemental system design and maintenance document

The information necessary for execution of the digital computer program L216 on the CDC 6600 is described. L216 characteristics are based on the doublet lattice method. Arbitrary aerodynamic configurations may be represented with combinations of nonplanar lifting surfaces composed of finite constant pressure panel elements, and axially summetric slender bodies composed of constant pressure line elements. Program input consists of configuration geometry, aerodynamic parameters, and modal data; output includes element geometry, pressure difference distributions, integrated aerodynamic coefficients, stability derivatives, generalized aerodynamic forces, and aerodynamic influence coefficient matrices. Optionally, modal data may be input on magnetic field (tape or disk), and certain geometric and aerodynamic output may be saved for subsequent use.

Harrison, B. A.↗
NASA NTRS1983

The efficient solution of transonic wing flow fields

An evaluation of the transonic-wing-analysis computer code TWING is presented. TWING utilizes a fully implicit, approximate-factorization iteration scheme to solve the full-potential equation in conservative form. A numerical elliptic-solver grid-generation scheme is used to generate the required finite-difference mesh. Several wing configurations have been analyzed, and comparisons of computed results have been made with available experimental data. Results indicate that the code is robust, accurate (when significant viscous effects are not present), and efficient. TWING generally produces solutions an order of magnitude faster than other conservative, full-potential codes using successive-line overrelaxation. The present method is applicable to a wide range of isolated wing configurations, including high-aspect-ratio transport wings and low-aspect-ratio, high-sweep, fighter configurations.

Holst, T. L.↗
NASA NTRS1991

Sensitivity analysis of a wing aeroelastic response

A variation of Sobieski's Global Sensitivity Equations (GSE) approach is implemented to obtain the sensitivity of the static aeroelastic response of a three-dimensional wing model. The formulation is quite general and accepts any aerodynamics and structural analysis capability. An interface code is written to convert one analysis's output to the other's input, and visa versa. Local sensitivity derivatives are calculated by either analytic methods or finite difference techniques. A program to combine the local sensitivities, such as the sensitivity of the stiffness matrix or the aerodynamic kernel matrix, into global sensitivity derivatives is developed. The aerodynamic analysis package FAST, using a lifting surface theory, and a structural package, ELAPS, implementing Giles' equivalent plate model are used.

Kapania, Rakesh K.↗
NASA NTRS2005

Emissions Prediction and Measurement for Liquid-Fueled TVC Combustor with and without Water Injection

An investigation is performed to evaluate the performance of a computational fluid dynamics (CFD) tool for the prediction of the reacting flow in a liquid-fueled combustor that uses water injection for control of pollutant emissions. The experiment consists of a multisector, liquid-fueled combustor rig operated at different inlet pressures and temperatures, and over a range of fuel/air and water/fuel ratios. Fuel can be injected directly into the main combustion airstream and into the cavities. Test rig performance is characterized by combustor exit quantities such as temperature and emissions measurements using rakes and overall pressure drop from upstream plenum to combustor exit. Visualization of the flame is performed using gray scale and color still photographs and high-frame-rate videos. CFD simulations are performed utilizing a methodology that includes computer-aided design (CAD) solid modeling of the geometry, parallel processing over networked computers, and graphical and quantitative post-processing. Physical models include liquid fuel droplet dynamics and evaporation, with combustion modeled using a hybrid finite-rate chemistry model developed for Jet-A fuel. CFD and experimental results are compared for cases with cavity-only fueling, while numerical studies of cavity and main fueling was also performed. Predicted and measured trends in combustor exit temperature, CO and NOx are in general agreement at the different water/fuel loading rates, although quantitative differences exist between the predictions and measurements.

Brankovic, A.↗
NASA NTRS1995

A kinematically driven anisotropic viscoelastic constitutive model applied to tires

Aircraft tires are composite structures manufactured with viscoelastic materials such as carbon black filled rubber and nylon cords. When loaded they experience large deflections and moderately large strains. Detailed structural models of tires require the use of either nonlinear shell or nonlinear three dimensional solid finite elements. Computational predictions of the dynamic response of tires must consider the composite viscoelastic material behavior in a realistic fashion. We describe a modification to a nonlinear anisotropic shell finite element so it can be used to model viscoelastic stresses during general deformations. The model is developed by introducing internal variables of the type used to model elastic strain energy. The internal variables are strains, curvatures, and transverse shear angles which are in a one-to-one correspondence with the generalized coordinates used to model the elastic strain energy for nonlinear response. A difference-relaxation equation is used to relate changes in the observable strain field to changes in the internal strain field. The internal stress state is introduced into the equilibrium equations by converting it to nodal loads associated with the element's displacement degrees of freedom. In this form the tangent matrix in the Newton-Raphson solution algorithm is not modified from its form for the nonlinear statics problem. Only the gradient vector is modified and the modification is not computationally costly. The existing finite element model for the Space Shuttle nose gear tire is used to provide examples of the algorithm. In the first example, the tire's rim is displaced at a constant rate up to a fixed value. In the second example, the tire's rim is enforced to follow a saw tooth load and unload curve to generate hysteresis loops.

Johnson, Arthur R.↗
NASA NTRS1973

Study of hypervelocity meteoroid impact on orbital space stations

Structural damage resulting in hypervelocity impact of a meteorite on a spacecraft is discussed. Of particular interest is the backside spallation caused by such a collision. To treat this phenomenon two numerical schemes were developed in the course of this study to compute the elastic-plastic flow fracture of a solid. The numerical schemes are a five-point finite difference scheme and a four-node finite element scheme. The four-node finite element scheme proved to be less sensitive to the type of boundary conditions and loadings. Although further development work is needed to improve the program versatility (generalization of the network topology, secondary storage for large systems, improving of the coding to reduce the run time, etc.), the basic framework is provided for a utilitarian computer program which may be used in a wide variety of situations. Analytic results showing the program output are given for several test cases.

Leimbach, K. R.↗
NASA NTRS1981

A numerical simulation of hypersonic viscous flow over arbitrary geometries at angle of attack

An implicit conservative, noniterative, finite-difference algorithm that predicts the supersonic, laminar or turbulent viscous flow about arbitrary geometries at large angles of attack is presented. The three-dimensional parabolized form of the thin-layer Navier-Stokes equations are written in generalized coordinates. These equations are solved using the delta form of the Beam-Warming implicit algorithm. Flow field simulations have been obtained for a blunt biconic with windward and leeward cuts and an X-24C lifting body for both laminar and turbulent flow at various Mach numbers and angles of attack. When compared with experiment or with previous theories, these computational predictions show good agreement.

Chaussee, D. S.↗
NASA NTRS1982

Postbuckling behavior of selected flat stiffened graphite-epoxy panels loaded in compression

Results of an experimental study of the postbuckling behavior of selected flat stiffened graphite-epoxy panels loaded in compression are presented. The postbuckling response and failure characteristics of undamaged panels and panels damaged by low-speed impact are described. Each panel had four equally-spaced I-spaced stiffeners and 16- or 24-ply quasi-isotropic skins. Panels with three different stiffener spacings were tested. Some undamaged specimens supported as much as three times their initial buckling load before failing. Failure of all panels initiated in a skin-stiffener interface region. Analytical results obtained from a nonlinear general shell finite element analysis computer code correlate well with typical postbuckling test results up to failure. The analytical modeling detail necessary to predict accurately the response of a panel is described. Test results show that low-speed impact damage can reduce the postbuckling strength of a stiffened panel and that the skin-stiffener interface region is more sensitive to impact damage than the skin midway between stiffeners.

Starnes, J. H., Jr.↗
NASA NTRS1988

Fatigue damage modeling for coated single crystal superalloys

A high temperature, low-cycle fatigue life prediction method for coated single crystal nickel-base superalloys is being developed. The method is being developed for use in predicting crack initiation life of coated single crystal turbine airfoils. Although the models are being developed using coated single crystal PWA 1480, they should be readily adaptable to other coated nickel-base single crystal materials. The coatings choosen for this effort were of two generic types: a low pressure plasma sprayed NiCoCrAlY overlay, designated PWA 286, and an aluminide diffusion, designated PWA 273. In order to predict the useful crack initiation life of airfoils, the constitutive and failure behavior of the coating/substrate combination must be taken into account. Coatings alter the airfoil surface microstructure and are a primary source from which cracks originate. The adopted life prediction approach addresses this complexity by separating the coating and single crystal crack initiation regimes. This provides a flexible means for using different life model formulations for the coating and single crystal materials. At the completion of this program, all constitutive and life model formulations will be available in equation form and as software. The software will use the MARC general purpose finite element code to drive the constitutive models and calculate life parameters.

Nissley, David M.↗
DOE OSTI2017

A general CFD framework for fault-resilient simulations based on multi-resolution information fusion

We develop a general CFD framework for multi-resolution simulations to target multiscale problems but also resilience in exascale simulations, where faulty processors may lead to gappy, in space-time, simulated fields. We combine approximation theory and domain decomposition together with statistical learning techniques, e.g. coKriging, to estimate boundary conditions and minimize communications by performing independent parallel runs. To demonstrate this new simulation approach, we consider two benchmark problems. First, we solve the heat equation (a) on a small number of spatial “patches” distributed across the domain, simulated by finite differences at fine resolution and (b) on the entire domain simulated at very low resolution, thus fusing multi-resolution models to obtain the final answer. Second, we simulate the flow in a lid-driven cavity in an analogous fashion, by fusing finite difference solutions obtained with fine and low resolution assuming gappy data sets. We investigate the influence of various parameters for this framework, including the correlation kernel, the size of a buffer employed in estimating boundary conditions, the coarseness of the resolution of auxiliary data, and the communication frequency across different patches in fusing the information at different resolution levels. In addition to its robustness and resilience, the new framework can be employed to generalize previous multiscale approaches involving heterogeneous discretizations or even fundamentally different flow descriptions, e.g. in continuum-atomistic simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Generalized and High-Order Multiscale Finite-Element Methods for Seismic Wave Propagation [Slides]

We have developed a series of multiscale finite-element methods for simulating seismic wave propagation in complex media. The most important feature of multiscale finite-element methods is using medium- and mesh-dependent multiscale basis functions. For different wave equations, or for wave equation modeling in different domains, we design different local spectral/linear problems to compute multiscale basis functions. Multiscale basis functions can accurately capture the fine-scale medium information for coarse-scale modeling problems. Numerical examples have demonstrated the efficiency and accuracy of our multiscale finite-element method framework for wave equation modeling. Future works aim at further improving the efficiency using more levels of information encapsulation, and address more realistic physics problems.

58 GEOSCIENCES↗
DOE OSTI2022

Machine-Learning of Nonlocal Kernels for Anomalous Subsurface Transport from Breakthrough Curves

Anomalous behavior is ubiquitous in subsurface solute transport due to the presence of high degrees of heterogeneity at different scales in the media. Although fractional models have been extensively used to describe the anomalous transport in various subsurface applications, their application is hindered by computational challenges. Simpler nonlocal models characterized by integrable kernels and finite interaction length represent a computationally feasible alternative to fractional models; yet, the informed choice of their kernel functions still remains an open problem. We propose a general data-driven framework for the discovery of optimal kernels on the basis of very small and sparse data sets in the context of anomalous subsurface transport. Using spatially sparse breakthrough curves recovered from fine-scale particle-density simulations, we learn the best coarse-scale nonlocal model using a nonlocal operator regression technique. Predictions of the breakthrough curves obtained using the optimal nonlocal model show good agreement with fine-scale simulation results even at locations and time intervals different from the ones used to train the kernel, confirming the excellent generalization properties of the proposed algorithm. A comparison with trained classical models and with black-box deep neural networks confirms the superiority of the predictive capability of the proposed model.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1994

Interlaminar stress analysis of dropped-ply laminated plates and shells by a mixed method

A mixed method of approximation based on Reissner's variational principle is developed for the linear analysis of interlaminar stresses in laminated composites, with special interest in laminates that contain terminated internal plies (dropped-ply laminates). Two models are derived, one for problems of generalized plane deformation and the other for the axisymmetric response of shells of revolution. A layerwise approach is taken in which the stress field is assumed with an explicit dependence on the thickness coordinate in each layer. The dependence of the stress field on the thickness coordinate is determined such that the three-dimensional equilibrium equations are satisfied by the approximation. The solution domain is reduced to one dimension by integration through the thickness. Continuity of tractions and displacements between layers is imposed. The governing two-point boundary value problem is composed of a system of both differential and algebraic equations (DAE's) and their associated boundary conditions. Careful evaluation of the system of DAE's was required to arrive at a form that allowed application of a one-step finite difference approximation. A two-stage Gauss implicit Runge-Kutta finite difference scheme was used for the solution because of its relatively high degree of accuracy. Patch tests of the two models revealed problems with solution accuracy for the axisymmetric model of a cylindrical shell loaded by internal pressure. Parametric studies of dropped-ply laminate characteristics and their influence on the interlaminar stresses were performed using the generalized plane deformation model. Eccentricity of the middle surface of the laminate through the ply drop-off was found to have a minimal effect on the interlaminar stresses under longitudinal compression, transverse tension, and in-plane shear. A second study found the stiffness change across the ply termination to have a much greater influence on the interlaminar stresses.

Harrison, Peter N.↗
DOE OSTI2016

Smoothed Boundary Method for simulating bulk and grain boundary transport in complex polycrystalline microstructures

Grain boundaries have a major impact on material properties, but explicit consideration of the complex geometries of grain structures in simulations poses a challenge. In this paper, we present a general method for incorporating the effect of grain boundaries based on the Smoothed Boundary Method (SBM). By using multiple domain parameters to define the domains of different grains, this method circumvents time-consuming mesh generation steps that are associated with finite element calculations involving complex microstructures. To validate the approach, we evaluate the accuracy of the SBM against the sharp interface method. The capabilities of this approach were demonstrated through simulations of surface and grain boundary diffusion, as well as those of electrochemical impedance spectroscopy. We conclude this method is applicable to many material systems in which grain boundaries play a crucial role.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Surface plasmon transmission through discontinuous conducting surfaces: Plasmon amplitude modulation by grazing scattered fields

We have studied numerically the diffraction of a surface plasmon polariton (SPP) when it encounters a wide multi-wavelength slit in conducting films. As a jump process a SPP is excited beyond the slit by wave scattering at the second slit edge. The exciting radiation is produced when the incident SPP collapses at the first slit edge. We have found that the transmitted SPP supports inherent and unavoidable interference with grazing scattered radiation; the spatial modulation extends to the fields in the diffraction region where a series of low intensity spots arises. We demonstrate that the SPP generated on the second slab depends on the frequency but not on the wave vector of the collapsed SPP; a SPP is transmitted even when the two metals forming the slit are different. The numerical results were obtained using the Finite Difference Time Domain (FDTD) method with a grid size λ/100.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1972

Flow field analysis of aircraft configurations using a numerical solution to the three-dimensional unified supersonic/hypersonic small disturbance equations, part 1

The unified small disturbance equations are numerically solved using the well-known Lax-Wendroff finite difference technique. The method allows complete determination of the inviscid flow field and surface properties as long as the flow remains supersonic. Shock waves and other discontinuities are accounted for implicity in the numerical method. This technique was programed for general application to the three-dimensional case. The validity of the method is demonstrated by calculations on cones, axisymmetric bodies, lifting bodies, delta wings, and a conical wing/body combination. Part 1 contains the discussion of problem development and results of the study. Part 2 contains flow charts, subroutine descriptions, and a listing of the computer program.

Gunness, R. C., Jr.↗