SEARCH · Engineering Papers

Results for “Formulation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Incompressible flow computations based on the vorticity-stream function and velocity-pressure formulations

Finite element procedures and computations based on the velocity-pressure and vorticity-stream function formulations of incompressible flows are presented. Two new multistep velocity-pressure formulations are proposed and compared with the vorticity-stream function and one-step formulations. The example problems chosen are the standing vortex problem and flow past a circular cylinder. Benchmark quality computations are performed for the cylinder problem. The numerical results indicate that the vorticity-stream function formulation and one of the two new multistep formulations involve much less numerical dissipation than the one-step formulation.

Tezduyar, T. E.↗

NASA NTRS1991

A novel approach in formulation of special transition elements: Mesh interface elements

The objective of this research program is in the development of more accurate and efficient methods for solution of singular problems encountered in various branches of mechanics. The research program can be categorized under three levels. The first two levels involve the formulation of a new class of elements called 'mesh interface elements' (MIE) to connect meshes of traditional elements either in three dimensions or in three and two dimensions. The finite element formulations are based on boolean sum and blending operators. MEI are being formulated and tested in this research to account for the steep gradients encountered in aircraft and space structure applications. At present, the heat transfer and structural analysis problems are being formulated from uncoupled theory point of view. The status report: (1) summarizes formulation for heat transfer and structural analysis; (2) explains formulation of MEI; (3) examines computational efficiency; and (4) shows verification examples.

Sarigul, Nesrin↗

NASA NTRS1995

Completed Beltrami-Michell formulation for analyzing mixed boundary value problems in elasticity

In elasticity, the method of forces, wherein stress parameters are considered as the primary unknowns, is known as the Beltrami-Michell formulation (BMF). The existing BMF can only solve stress boundary value problems; it cannot handle the more prevalent displacement of mixed boundary value problems of elasticity. Therefore, this formulation, which has restricted application, could not become a true alternative to the Navier's displacement method, which can solve all three types of boundary value problems. The restrictions in the BMF have been alleviated by augmenting the classical formulation with a novel set of conditions identified as the boundary compatibility conditions. This new method, which completes the classical force formulation, has been termed the completed Beltrami-Michell formulation (CBMF). The CBMF can solve general elasticity problems with stress, displacement, and mixed boundary conditions in terms of stresses as the primary unknowns. The CBMF is derived from the stationary condition of the variational functional of the integrated force method. In the CBMF, stresses for kinematically stable structures can be obtained without any reference to the displacements either in the field or on the boundary. This paper presents the CBMF and its derivation from the variational functional of the integrated force method. Several examples are presented to demonstrate the applicability of the completed formulation for analyzing mixed boundary value problems under thermomechanical loads. Selected example problems include a cylindrical shell wherein membrane and bending responses are coupled, and a composite circular plate.

Patnaik, Surya N.↗

NASA NTRS2001

Stress Formulation in Three-Dimensional Elasticity

The theory of elasticity evolved over centuries through the contributions of eminent scientists like Cauchy, Navier, Hooke Saint Venant, and others. It was deemed complete when Saint Venant provided the strain formulation in 1860. However, unlike Cauchy, who addressed equilibrium in the field and on the boundary, the strain formulation was confined only to the field. Saint Venant overlooked the compatibility on the boundary. Because of this deficiency, a direct stress formulation could not be developed. Stress with traditional methods must be recovered by backcalculation: differentiating either the displacement or the stress function. We have addressed the compatibility on the boundary. Augmentation of these conditions has completed the stress formulation in elasticity, opening up a way for a direct determination of stress without the intermediate step of calculating the displacement or the stress function. This Completed Beltrami-Michell Formulation (CBMF) can be specialized to derive the traditional methods, but the reverse is not possible. Elasticity solutions must be verified for the compliance of the new equation because the boundary compatibility conditions expressed in terms of displacement are not trivially satisfied. This paper presents the variational derivation of the stress formulation, illustrates the method, examines attributes and benefits, and outlines the future course of research.

Patnaik, Surya N.↗

NASA NTRS1999

Verification and Analysis of Formulation 4 of Langley for the Study of Noise From High Speed Surfaces

There are several approaches to the prediction of the noise from sources on high speed surfaces. Two of these are the Kirchhoff and the Ffowcs williams-Hawkings methods. It can be shown that both of these methods depend on the solution of the wave equation with mathematically similar inhomogeneous source terms. Two subsonic solutions known as Formulation 1 and 1A of Langley are simple and efficient for noise prediction. The supersonic solution known as Formulation 3 is very complicated and difficult to code. Because of the complexity of the result, the computation time is longer than the subsonic formulas. Furthermore, it is difficult to assess the accuracy of noise prediction. We have been searching for a new and simpler supersonic formulation without these shortcomings. In the last AIAA Aeroacoustics Conference in Toulouse, Farassat, Dunn and Brentner presented a paper in which such a result was presented and called Formulation 4 of Langley. In this paper we will present two analytic tests of the validity this Formulation: 1) the noise from dipole distribution on the unit circle whose strength varies radially with the square of the distance from the center and 2) the noise from dipole distribution on the unit sphere whose strength varies with the cosine of the angle from the polar axis. We will discuss the question of singularities of Formulation 4.

Farassat, F.↗

NASA NTRS2004

Thermoelastic Formulation of Stiffened, Unsymmetric Composite Panels for Finite Element Analysis of High Speed Aircraft

An emerging technology need for capturing 3-D panel thermoelastic response with 2-D planar finite element models (FEMs) is aided with an equivalent plate stiffness and thermal coefficient formulation. The formulation is general and applies to all panel concepts. Included with the formulation is the ability to provide membrane-bending coupling of unsymmetric sections and calculation of all thermal expansion and bending responses from in-plane and through-the-thickness temperature gradients. Thermal residual strains for both the laminates and plies are included. The general formulation is defined and then applied to a hat-shaped, corrugated stiffened panel. Additional formulations are presented where required to include all of the hat's unique characteristics. Each formulation is validated independently with 3-D FEA.

Collier, Craig S.↗

NASA NTRS2007

Decision-Tree Formulation With Order-1 Lateral Execution

A compact symbolic formulation enables mapping of an arbitrarily complex decision tree of a certain type into a highly computationally efficient multidimensional software object. The type of decision trees to which this formulation applies is that known in the art as the Boolean class of balanced decision trees. Parallel lateral slices of an object created by means of this formulation can be executed in constant time considerably less time than would otherwise be required. Decision trees of various forms are incorporated into almost all large software systems. A decision tree is a way of hierarchically solving a problem, proceeding through a set of true/false responses to a conclusion. By definition, a decision tree has a tree-like structure, wherein each internal node denotes a test on an attribute, each branch from an internal node represents an outcome of a test, and leaf nodes represent classes or class distributions that, in turn represent possible conclusions. The drawback of decision trees is that execution of them can be computationally expensive (and, hence, time-consuming) because each non-leaf node must be examined to determine whether to progress deeper into a tree structure or to examine an alternative. The present formulation was conceived as an efficient means of representing a decision tree and executing it in as little time as possible. The formulation involves the use of a set of symbolic algorithms to transform a decision tree into a multi-dimensional object, the rank of which equals the number of lateral non-leaf nodes. The tree can then be executed in constant time by means of an order-one table lookup. The sequence of operations performed by the algorithms is summarized as follows: 1. Determination of whether the tree under consideration can be encoded by means of this formulation. 2. Extraction of decision variables. 3. Symbolic optimization of the decision tree to minimize its form. 4. Expansion and transformation of all nested conjunctive-disjunctive paths to a flattened conjunctive form composed only of equality checks when possible. If each reduced conjunctive form contains only equality checks and all of these forms use the same variables, then the decision tree can be reduced to an order-one operation through a table lookup. The speedup to order one is accomplished by distributing each decision variable over a surface of a multidimensional object by mapping the equality constant to an index

James, Mark↗

NASA NTRS2014

On Formulations of Discontinuous Galerkin and Related Methods for Conservation Laws

A formulation for the discontinuous Galerkin (DG) method that leads to solutions using the differential form of the equation (as opposed to the standard integral form) is presented. The formulation includes (a) a derivative calculation that involves only data within each cell with no data interaction among cells, and (b) for each cell, corrections to this derivative that deal with the jumps in fluxes at the cell boundaries and allow data across cells to interact. The derivative with no interaction is obtained by a projection, but for nodal-type methods, evaluating this derivative by interpolation at the nodal points is more economical. The corrections are derived using the approximate (Dirac) delta functions. The formulation results in a family of schemes: different approximate delta functions give rise to different methods. It is shown that the current formulation is essentially equivalent to the flux reconstruction (FR) formulation. Due to the use of approximate delta functions, an energy stability proof simpler than that of Vincent, Castonguay, and Jameson (2011) for a family of schemes is derived. Accuracy and stability of resulting schemes are discussed via Fourier analyses. Similar to FR, the current formulation provides a unifying framework for high-order methods by recovering the DG, spectral difference (SD), and spectral volume (SV) schemes. It also yields stable, accurate, and economical methods.

Huynh, H. T.↗

NASA NTRS2017

A New Formulation of Time Domain Boundary Integral Equation for Acoustic Wave Scattering in the Presence of a Uniform Mean Flow

It has been well-known that under the assumption of a constant uniform mean flow, the acoustic wave propagation equation can be formulated as a boundary integral equation, in both the time domain and the frequency domain. Compared with solving partial differential equations, numerical methods based on the boundary integral equation have the advantage of a reduced spatial dimension and, hence, requiring only a surface mesh. However, the constant uniform mean flow assumption, while convenient for formulating the integral equation, does not satisfy the solid wall boundary condition wherever the body surface is not aligned with the uniform mean flow. In this paper, we argue that the proper boundary condition for the acoustic wave should not have its normal velocity be zero everywhere on the solid surfaces, as has been applied in the literature. A careful study of the acoustic energy conservation equation is presented that shows such a boundary condition in fact leads to erroneous source or sink points on solid surfaces not aligned with the mean flow. A new solid wall boundary condition is proposed that conserves the acoustic energy and a new time domain boundary integral equation is derived. In addition to conserving the acoustic energy, another significant advantage of the new equation is that it is considerably simpler than previous formulations. In particular, tangential derivatives of the solution on the solid surfaces are no longer needed in the new formulation, which greatly simplifies numerical implementation. Furthermore, stabilization of the new integral equation by Burton-Miller type reformulation is presented. The stability of the new formulation is studied theoretically as well as numerically by an eigenvalue analysis. Numerical solutions are also presented that demonstrate the stability of the new formulation.

Hu, Fang↗

NASA NTRS2023

A General Procedure to Formulate 3D Elements for Finite Element Applications

This paper presents a general procedure to formulate and implement 3D elements of arbitrary order in meshes with multiple element types. This procedure includes obtaining shape functions and integration quadrature and establishing an approach for checking the generated element’s compatibility with adjacent elements’ surfaces. This procedure was implemented in Matlab, using its symbolic and graphics toolbox, and complied as a GUI interface named ShapeGen3D to provide finite element users with a tool to tailor elements according to their analysis needs. ShapeGen3D also outputs files with the element formulation needed to enable users to implement the generated elements in other programming languages or through user elements in commercial finite element software. Currently, finite element (FE) users are limited to employing element formulation available in the literature, commercial software, or existing element libraries. Thus, the developed procedure implemented in ShapeGen3D offers FEM users the possibility to employ elements beyond those readily available. The procedure was tested by generating the formulation for a brick element, a brick transition element, and higher-order hexahedron and tetrahedron elements that can be used in a spectral finite element analysis. The formulation obtained for the 20-node element was in perfect agreement with the formulation available in the literature. In addition, the results showed that the interpolation condition was met for all the generated elements, which provides confidence in the implementation of the process. Researchers and educators can use this procedure to efficiently develop and illustrate three-dimensional elements.

higher-order elements↗

NASA NTRS1980

MCC level C formulation requirements. Entry guidance and entry autopilot, optional TAEM targeting

The level C software formulations requirements for the entry guidance and the simplified autopilot to be used by the Mission Control Center (MCC) entry processor are presented. The modifications required to functionally simulate optional Terminal Area Energy Management (TAEM) targeting capability (OTT) are incorporated. Implementation of this logic in the MCC must be coordinated with flight software OTT implementation and MCC TAEM guidance OTT. The entry guidance logic is based on the orbiter avionics entry guidance software. Descriptions of the entry guidance detailed formulation requirements, the detailed autopilot formulation requirements, and the targeting routine are given. Also included are a definition of coordinate systems, a list of parameter definitions for the software formulations, and a set of formulation flow charts.

Harpold, J. C.↗

NASA NTRS1980

Formulation and implementation of a practical algorithm for parameter estimation with process and measurement noise

A new formulation is proposed for the problem of parameter estimation of dynamic systems with both process and measurement noise. The formulation gives estimates that are maximum likelihood asymptotically in time. The means used to overcome the difficulties encountered by previous formulations are discussed. It is then shown how the proposed formulation can be efficiently implemented in a computer program. A computer program using the proposed formulation is available in a form suitable for routine application. Examples with simulated and real data are given to illustrate that the program works well.

Maine, R. E.↗

NASA NTRS1981

Steady-state boundary lubrication with formulated C-ethers to 260 C

Steady state wear and friction studies were made at boundary lubrication conditions in a pin on disk (pure iron on rotating CVM M 50 steel) sliding friction apparatus with five C ether formulated fluids (modified polyphenyl ether containing phosphrous ester, organic acid, and other additives). Conditions included 20, 150, and 260 C disk temperatures, dry air test atmosphere, 1 kilogram load, 50 rpm disk speed, and test times to 130 minutes. Results were compared with those obtained with a formulated MIL L 27502 candidate ester and the C ether base fluid. Three of the C ether formulations gave better lubrication than both reference fluids under most conditions. The other two C ether formulations yielded higher wear rates and friction coefficients than the C ether base fluid for most of the temperature range. Only one C ether formulation showed consistently higher steady state wear rates than the ester.

Loomis, W. R.↗

NASA NTRS1983

Polyimides formulated from a partially fluorinated diamine for aerospace tribological applications

Preliminary tribological studies on polyimides formulated from the diamine 2,2-bis 4-(4-aminophenoxy)phenyl hexafluorapane (4-BDAF) indicate that polyimides formulated from this diamine have excellent potential for high temperature tribological applications. The dianhydrides used to make the polyimides were pyromellitic (PMDA) and benzophenonetetracarboxylic acid (BTDA). Friction and wear studies at 25 and 200 C indicate that polyimides formulated using 50 mole percent of the PMDA dianhydride and 50 mole percent of the BTDA dianhydride perform better than polyimides formulated solely with the BTDA dianhydride. Graphite fiber reinforced polyimide composites were formulated with the polyimide made from the BTDA dianhydride, both graphitic and non-graphitic fibers were evaluated. Graphitic fibers produced better tribological results, since thin, flowing, "layer-like' transfer films were produced which did not build-up with long sliding durations. Non-graphitic fibers did not produce this type of transfer.

Fusaro, R. L.↗

NASA NTRS1984

Polyimides formulated from a partially fluorinated diamine for aerospace tribological applications

Preliminary pribological studies on polyimides formulated from the diamine 2,2-bis 4-(4-aminophenoxy)phenyl hexafluoropane (4-BDAF) indicate that polyimides formulated from this diamine have excellent potential for high temperature tribological applications. The dianhydrides used to make the polyimides were pyromellitic (PMDA) and benzophenonetetracarboxylic acid (BTDA). Friction and wear studies at 25 and 200 C indicate that polyimides formulated using 50 mole percent of the PMDA dianhydride and 50 mole percent of the BTDA dianhydride perform better than polyimides formulated solely with the BTDA dianhydride. Graphite fiber reinforced polyimide composites were formulated with the polyimide made from the BTDA dianhydride and both graphitic and non-graphitic fibers were evaluated. Graphitic fibers produced better tribological results, since thin, flowing, 'layer-like' transfer films were produced which did not build up with long sliding durations. Non-graphitic fibers did not produce this type of transfer. Previously announced in STAR as N83-22423

Fusaro, R. L.↗

NASA NTRS1987

Comparison of two transonic noise prediction formulations using the aircraft noise prediction program

This paper addresses recently completed work on using Farassat's Formulation 3 noise prediction code with the Aircraft Noise Prediction Program (ANOPP). Software was written to link aerodynamic loading generated by the Propeller Loading (PLD) module within ANOPP with formulation 3. Included are results of comparisons between Formulation 3 with ANOPP's existing noise prediction modules, Subsonic Propeller Noise (SPN) and Transonic Propeller Noise (TPN). Four case studies are investigated. Results of the comparison studies show excellent agreement for the subsonic cases. Differences found in the comparisons made under transonic conditions are strictly numerical and can be explained by the way in which the time derivative is calculated in Formulation 3. Also included is a section on how to execute Formulation 3 with ANOPP.

Spence, Peter L.↗

NASA NTRS1986

An analysis of finite-difference and finite-volume formulations of conservation laws

Finite-difference and finite-volume formulations are analyzed in order to clear up the confusion concerning their application to the numerical solution of conservation laws. A new coordinate-free formulation of systems of conservation laws is developed, which clearly distinguishes the role of physical vectors from that of algebraic vectors which characterize the system. The analysis considers general types of equations--potential, Euler, and Navier-Stokes. Three-dimensional unsteady flows with time-varying grids are described using a single, consistent nomeclature for both formulations. Grid motion due to a non-inertial reference frame as well as flow adaptation is covered. In comparing the two formulations, it is found useful to distinguish between differences in numerical methods and differences in grid definition. The former plays a role for non-Cartesian grids, and results in only cosmetic differences in the manner in which geometric terms are handled. The differences in grid definition for the two formulations is found to be more important, since it affects the manner in which boundary conditions, zonal procedures, and grid singularities are handled at computational boundaries. The proper interpretation of strong and weak conservation-law forms for quasi-one-dimensional and axisymmetric flows is brought out.

Vinokur, Marcel↗

NASA NTRS1989

An analysis of finite-difference and finite-volume formulations of conservation laws

Finite-difference and finite-volume formulations are analyzed in order to clear up the confusion concerning their application to the numerical solution of conservation laws. A new coordinate-free formulation of systems of conservation laws is developed, which clearly distinguishes the role of physical vectors from that of algebraic vectors which characterize the system. The analysis considers general types of equations: potential, Euler, and Navier-Stokes. Three-dimensional unsteady flows with time-varying grids are described using a single, consistent nomenclature for both formulations. Grid motion due to a non-inertial reference frame as well as flow adaptation is covered. In comparing the two formulations, it is found useful to distinguish between differences in numerical methods and differences in grid definition. The former plays a role for non-Cartesian grids, and results in only cosmetic differences in the manner in which geometric terms are handled. The differences in grid definition for the two formulations is found to be more important, since it affects the manner in which boundary conditions, zonal procedures, and grid singularities are handled at computational boundaries. The proper interpretation of strong and weak conservation-law forms for quasi-one-dimensional and axisymmetric flows is brought out.

Vinokur, Marcel↗