Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite difference 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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

Assessment of existing analytic methods for prediction of high angle-of-attack loads on delta wings at supersonic speeds

An assessment of the applicability of four loading prediction methods to high angle-of-attack conditions for simplified wing-body configurations is provided. The methods are: The tangent wedge approximation, the linear theory methods of Middleton and Woodward, and a shock-fitting finite-difference technique. Estimates obtained by these methods were compared with experimental pressure data on delta wings to examine the effects of Mach number, camber, sweep angle, and angle of attack. Results indicate that all of the methods provided reasonable estimates at moderate angles of attack. At these moderate angles of attack, the methods of Middleton and Woodward provided good estimates at Mach numbers higher than those usually associated with linear theory. Only the finite-difference method provided reasonable load estimates at high angles of attack.

Landrum, E. J.↗
NASA NTRS1992

Two-dimensional hybrid continuum/particle approach for rarefied flows

A hybrid numerical technique previously developed for one-dimensional rarefied gas flows is generalized to two dimensions. The method is based on the fact that the flowfield that develops near a body in a rarefied gas typically contains local regions of continuum, transitional and free-molecular flow. By utilizing the solution technique most appropriate for each region and coupling the techniques in an interface region where both are applicable, more computationally efficient solutions can be obtained, or equivalently, more complex flowfields can be analyzed. The present method combines finite difference solution of the Navier-Stokes equations in the continuum regions, with Direct Simulation Monte Carlo in the more rarefied regions. The two schemes are coupled interactively via a general conservative flux boundary condition. The method is tested by application to the model problem of pressure-driven rarefied flow through a slit. Results show the hybrid scheme offers a speedup of a factor of nearly two for the nominal conditions considered, due to the decrease in the size of the Direct Simulation domain.

Wadsworth, Dean C.↗
NASA NTRS1990

Numerical Methods For Chemically Reacting Flows

Issues related to numerical stability, accuracy, and resolution discussed. Technical memorandum presents issues in numerical solution of hyperbolic conservation laws containing "stiff" (relatively large and rapidly changing) source terms. Such equations often used to represent chemically reacting flows. Usually solved by finite-difference numerical methods. Source terms generally necessitate use of small time and/or space steps to obtain sufficient resolution, especially at discontinuities, where incorrect mathematical modeling results in unphysical solutions.

Leveque, R. J.↗
NASA NTRS1977

Analysis of simple 2-D and 3-D metal structures subjected to fragment impact

Theoretical methods were developed for predicting the large-deflection elastic-plastic transient structural responses of metal containment or deflector (C/D) structures to cope with rotor burst fragment impact attack. For two-dimensional C/D structures both, finite element and finite difference analysis methods were employed to analyze structural response produced by either prescribed transient loads or fragment impact. For the latter category, two time-wise step-by-step analysis procedures were devised to predict the structural responses resulting from a succession of fragment impacts: the collision force method (CFM) which utilizes an approximate prediction of the force applied to the attacked structure during fragment impact, and the collision imparted velocity method (CIVM) in which the impact-induced velocity increment acquired by a region of the impacted structure near the impact point is computed. The merits and limitations of these approaches are discussed. For the analysis of 3-d responses of C/D structures, only the CIVM approach was investigated.

Witmer, E. A.↗
NASA NTRS1991

Aerodynamic design optimization using sensitivity analysis and computational fluid dynamics

A new and efficient method is presented for aerodynamic design optimization, which is based on a computational fluid dynamics (CFD)-sensitivity analysis algorithm. The method is applied to design a scramjet-afterbody configuration for an optimized axial thrust. The Euler equations are solved for the inviscid analysis of the flow, which in turn provides the objective function and the constraints. The CFD analysis is then coupled with the optimization procedure that uses a constrained minimization method. The sensitivity coefficients, i.e. gradients of the objective function and the constraints, needed for the optimization are obtained using a quasi-analytical method rather than the traditional brute force method of finite difference approximations. During the one-dimensional search of the optimization procedure, an approximate flow analysis (predicted flow) based on a first-order Taylor series expansion is used to reduce the computational cost. Finally, the sensitivity of the optimum objective function to various design parameters, which are kept constant during the optimization, is computed to predict new optimum solutions. The flow analysis of the demonstrative example are compared with the experimental data. It is shown that the method is more efficient than the traditional methods.

Baysal, Oktay↗
NASA NTRS1979

Evaluation of finite element formulations for transient conduction forced-convection analysis

Conventional versus upwind convective finite elements, and lumped versus consistent formulations for practical conduction/forced convection analysis are evaluated on the basis of numerical studies, with finite element and finite difference lumped-parameter temperatures compared to closed-form analytical solutions for convection problems. Attention is given to two practical combined conduction and forced convection applications, stressing that the finite element method, showing superior accuracy, is competitive with the finite difference lumped-parameter method. Also considered are the computational time savings offered by the zero capacitance nodes procedure and comparative finite element and finite difference lumped-parameter computer times. The present study has reference to the design of actively cooled engine and airframe structures for hypersonic flight.

Thornton, E. A.↗
NASA NTRS1996

On the Sensitivity of Piezoceramics and Piezopolymers in Structural Integrity Monitoring of Large Trusses

An analytical assessment has been made of the reliability of using integrated microactuators and sensors in the form of piezoceramics and piezopolymers as joint integrity monitors in trussed systems. The concept is first implemented for a simple structure which consists of two truss members with a 45 deg lift angle joined at the apex. A piezoceramic patch (or piezopolymer film) bonded on the surface of one of the members at a location near the joint is used as a collocated actuator/sensor. The overall structural dynamic response under an excitation was modeled by finite element method. Different degrees of nodal constraints at the joints representing various degrees of joint integrity are employed. The resulting dynamic response showed distinct responses for varying joint stiffnesses. Parallel experimental work on a truss model using a multichannel data acquisition system and a digital signal analyzer confirms the results from analysis. We further studied the sensitivity of the micro-sensors to the behavior of joints of large arch truss structure. Results obtained for large trusses with many degrees of freedom indicate optimum locations of sensors for which the dynamic response signatures are distinct and distinguishable for relatively small changes in joint integrity and/or structural geometry. Computations based on finite element modeling show that locating the single actuator/sensor at the joint corresponding to the first loss of static stability appear optimal. Hence, static stability analysis of complex trusses can give us a good indication of the optimum placement of sensors for maximum response. This observation is important if few distributed sensors and actuators are available for placement in constructed facilities made from large trusses with many degrees of freedom. As an extension of this work a dynamic response signature identification technique to monitor in-service degradation of joints is under development for application to the monitoring of the integrity of adhesive joints in composite structures.

Abatan, A. O.↗
NASA NTRSDate not supplied

Thermal Modelling of the in-Situ Consolidation of Automated Fiber Placement of Thermoplastic Composites

NASA is developing carbon fiber reinforced thermoplastic composites processing methods under the Hi-rate Composites Aircraft Manufacturing (HiCAM) project. The in-situ consolidation automated fiber placement (AFP) of thermoplastics (ICAT) process has high potential to increase manufacturing throughput and lower costs as it combines the repeatability and fast lay down rates of AFP with the thermoformability of thermoplastic materials to achieve part consolidation out of the autoclave. Physics-based process models are employed in the ICAT process development to understand, fundamentally, the thermal response of the carbon fiber (CF)/polyaryletherketone (PAEK) material during the rapid heating and cooling associated with laser-assisted AFP. A one-dimensional (1-D) closed-form model was developed capable of predicting the temperature profile through-the-thickness of the tape material. In addition, a two-dimensional (2-D) model was developed to predict the heat transfer through-the-thickness and in the laydown x-direction of the AFP head motion. The solution of the Fourier heat-transfer equations in the 2-D model is approximated using the explicit finite difference numerical method. The temperature profile during placement of PAEK slit-tape materials at various laser target temperatures and placement speeds were measured during ICAT process development trials at Electroimpact®, Inc. †† and the resulting experimental data is compared with the predictions of the 1-D and 2-D models.

Thermoplastics↗
NASA NTRSDate not supplied

Thermal Modelling of the In-Situ Consolidation of Automated Fiber Placement Of Thermoplastic Composites

NASA is developing carbon fiber reinforced thermoplastic composites processing methods under the Hi-rate Composites Aircraft Manufacturing (HiCAM) project. The in-situ consolidation automated fiber placement (AFP) of thermoplastics (ICAT) process has high potential to increase manufacturing throughput and lower costs as it combines the repeatability and fast lay down rates of AFP with the thermoformability of thermoplastic materials to achieve part consolidation out of the autoclave. Physics-based process models are employed in the ICAT process development to understand, fundamentally, the thermal response of the carbon fiber (CF)/polyaryletherketone (PAEK) material during the rapid heating and cooling associated with laser-assisted AFP. A one-dimensional (1-D) closed-form model was developed capable of predicting the temperature profile through-the-thickness of the tape material. In addition, a two-dimensional (2-D) model was developed to predict the heat transfer through-the-thickness and in the laydown x-direction of the AFP head motion. The solution of the Fourier heat-transfer equations in the 2-D model is approximated using the explicit finite difference numerical method. The temperature profile during placement of PAEK slit-tape materials at various laser target temperatures and placement speeds were measured during ICAT process development trials at Electroimpact®, Inc. †† and the resulting experimental data is compared with the predictions of the 1-D and 2-D models.

Christopher J. Stelter↗
NASA NTRS1990

The semidiscrete Galerkin finite element modeling of compressible viscous flow past an airfoil

The primary project was the numerical simulation, by a finite element/finite difference method, of the viscous flow about an airfoil. The secondary project involved the numerical simulation of the three-dimensional separated and vortex-dominated flow about a hemispherically capped cylinder in the transonic regime. Preliminary calculations were started for the hemisphere-cylinder at 0 and 5 degree angle of attack. The solution of the flow field about airfoils and wings is required to determine the important parameters of lift, moment, and drag. Viscous effects must be accounted for if the drag is to be accurately calculated. At present there are basically two approaches to the numerical simulation of the flow field, the use of fully viscous models and the inviscid/viscous models. The fully viscous models require the solution of an approximation of the Navier-Stokes equations and therefore should simulate most of the physical mechanisms. A fast, accurate, and computationally efficient inviscid flow solver was recently developed by Hartwich. It is thought that Hartwich's program coupled to a fast, accurate, and computationally efficient boundary layer code, will make an excellent tool for airfoil design. The purpose of the primary project was to develop a compressible boundary layer code using the semidiscrete Galerkin finite element method. The numerical scheme employed used the combination of a Dorodnitsyn formulation of the boundary layer equations, with a finite difference/finite element procedure (semidiscrete Galerkin method), in the solution of the compressible two-dimensional boundary layer equations. A laminar compressible boundary layer code was developed and tested for a NACA 0012 airfoil at a Mach number of 0.5, a Reynolds number of 5000, and zero angle of attack. At present the boundary layer program solves up to, but not beyond, separation.

Meade, Andrew J.↗