Engineering Papers⌕ Search

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.

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

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS1991

Parallel computing using a Lagrangian formulation

A new Lagrangian formulation of the Euler equation is adopted for the calculation of 2-D supersonic steady flow. The Lagrangian formulation represents the inherent parallelism of the flow field better than the common Eulerian formulation and offers a competitive alternative on parallel computers. The implementation of the Lagrangian formulation on the Thinking Machines Corporation CM-2 Computer is described. The program uses a finite volume, first-order Godunov scheme and exhibits high accuracy in dealing with multidimensional discontinuities (slip-line and shock). By using this formulation, a better than six times speed-up was achieved on a 8192-processor CM-2 over a single processor of a CRAY-2.

Liou, May-Fun↗
NASA NTRS1992

Parallel computing using a Lagrangian formulation

This paper adopts a new Lagrangian formulation of the Euler equation for the calculation of two dimensional supersonic steady flow. The Lagrangian formulation represents the inherent parallelism of the flow field better than the common Eulerian formulation and offers a competitive alternative on parallel computers. The implementation of the Lagrangian formulation on the Thinking Machines Corporation CM-2 Computer is described. The program uses a finite volume, first-order Godunov scheme and exhibits high accuracy in dealing with multidimensional discontinuities (slip-line and shock). By using this formulation, we have achieved better than six times speed-up on a 8192-processor CM-2 over a single processor of a CRAY-2.

Liou, May-Fun↗
NASA NTRS1995

Non-parabolic hydrodynamic formulations for the simulation of inhomogeneous semiconductor devices

Hydrodynamic models are becoming prevalent design tools for small scale devices and other devices in which high energy effects can dominate transport. Most current hydrodynamic models use a parabolic band approximation to obtain fairly simple conservation equations. Interest in accounting for band structure effects in hydrodynamic device simulation has begun to grow since parabolic models can not fully describe the transport in state of the art devices due to the distribution populating non-parabolic states within the band. This paper presents two different non-parabolic formulations of the hydrodynamic model suitable for the simulation of inhomogeneous semiconductor devices. The first formulation uses the Kane dispersion relationship (hk)(exp 2)/2m = W(1 + alpha(W)). The second formulation makes use of a power law ((hk)(exp 2)/2m = xW(sup y)) for the dispersion relation. Hydrodynamic models which use the first formulation rely on the binomial expansion to obtain moment equations with closed form coefficients. This limits the energy range over which the model is valid. The power law formulation readily produces closed form coefficients similar to those obtained using the parabolic band approximation. However, the fitting parameters (x,y) are only valid over a limited energy range. The physical significance of the band non-parabolicity is discussed as well as the advantages/disadvantages and approximations of the two non-parabolic models. A companion paper describes device simulations based on the three dispersion relationships: parabolic, Kane dispersion, and power low dispersion.

Smith, Arlynn W.↗
NASA NTRS1996

A priori testing of subgrid-scale models for the velocity-pressure and vorticity-velocity formulations

Subgrid-scale models for Large Eddy Simulation (LES) in both the velocity-pressure and the vorticity-velocity formulations were evaluated and compared in a priori tests using spectral Direct Numerical Simulation (DNS) databases of isotropic turbulence: 128(exp 3) DNS of forced turbulence (Re(sub(lambda))=95.8) filtered, using the sharp cutoff filter, to both 32(exp 3) and 16(exp 3) synthetic LES fields; 512(exp 3) DNS of decaying turbulence (Re(sub(Lambda))=63.5) filtered to both 64(exp 3) and 32(exp 3) LES fields. Gaussian and top-hat filters were also used with the 128(exp 3) database. Different LES models were evaluated for each formulation: eddy-viscosity models, hyper eddy-viscosity models, mixed models, and scale-similarity models. Correlations between exact versus modeled subgrid-scale quantities were measured at three levels: tensor (traceless), vector (solenoidal 'force'), and scalar (dissipation) levels, and for both cases of uniform and variable coefficient(s). Different choices for the 1/T scaling appearing in the eddy-viscosity were also evaluated. It was found that the models for the vorticity-velocity formulation produce higher correlations with the filtered DNS data than their counterpart in the velocity-pressure formulation. It was also found that the hyper eddy-viscosity model performs better than the eddy viscosity model, in both formulations.

Winckelmans, G. S.↗
NASA NTRS1996

Non-Parabolic Hydrodynamic Formulations for the Simulation of Inhomogeneous Semiconductor Devices

Hydrodynamic models are becoming prevalent design tools for small scale devices and other devices in which high energy effects can dominate transport. Most current hydrodynamic models use a parabolic band approximation to obtain fairly simple conservation equations. Interest in accounting for band structure effects in hydrodynamic device simulation has begun to grow since parabolic models cannot fully describe the transport in state of the art devices due to the distribution populating non-parabolic states within the band. This paper presents two different non-parabolic formulations or the hydrodynamic model suitable for the simulation of inhomogeneous semiconductor devices. The first formulation uses the Kane dispersion relationship ((hk)(exp 2)/2m = W(1 + alphaW). The second formulation makes use of a power law ((hk)(exp 2)/2m = xW(exp y)) for the dispersion relation. Hydrodynamic models which use the first formulation rely on the binomial expansion to obtain moment equations with closed form coefficients. This limits the energy range over which the model is valid. The power law formulation readily produces closed form coefficients similar to those obtained using the parabolic band approximation. However, the fitting parameters (x,y) are only valid over a limited energy range. The physical significance of the band non-parabolicity is discussed as well as the advantages/disadvantages and approximations of the two non-parabolic models. A companion paper describes device simulations based on the three dispersion relationships; parabolic, Kane dispersion and power law dispersion.

Smith, A. W.↗
NASA NTRS2000

Higher Order Modeling In the BEM/FEM Hybrid Formulation

Hybrid formulations using low order curl-conforming bases to represent the total electric field within a finite element region and low order divergence-conforming bases to represent equivalent electric and magnetic currents on the boundary are well known. However, higher-order divergence and curl-conforming bases have been shown to provide significant benefits in convergence rates and accuracy when employed in strictly integral equation and strictly finite element formulations. In this paper, a hybrid electric field formulation employing higher order bases is presented. The paper addresses benefits and issues associated with using higher order divergence-and curl-conforming bases in the hybrid finite element/boundary element electric field formulation. The method of singularity subtraction may be used to compute the self terms of the boundary integral when the bases are of low order. But this method becomes laborious and requires great care when the divergence conforming bases are of higher order. In order to handle these singularities simply and accurately, a generalized Gaussian quadrature method is employed in which the expansion functions account for the singularity. In preliminary tests of the higher order hybrid formulation, the equivalent electric current induced by scattering of a plane wave from a square dielectric cylinder is examined. Accurate results are obtained using only a two-triangle mesh when the current basis is of order 4 or 5. Additional results are presented comparing the error obtained using higher order bases to that obtained using lower order bases when the number of unknowns is approximately equal. Also, convergence rates obtained with higher order bases are compared to those obtained with lower order bases for selected sample problems.

Fink, Patrick W.↗
NASA NTRS2002

Slow Crack Growth of Brittle Materials With Exponential Crack-Velocity Formulation: Analysis - Part 1

Extensive slow-crack-growth (SCG) analysis was made using a primary exponential crack-velocity formulation under three widely used load configurations: constant stress rate, constant stress, and cyclic stress. Although the use of the exponential formulation in determining SCG parameters of a material requires somewhat inconvenient numerical procedures, the resulting solutions presented gave almost the same degree of simplicity in both data analysis and experiments as did the power-law formulation. However, the fact that the inert strength of a material should be known in advance to determine the corresponding SCG parameters was a major drawback of the exponential formulation as compared with the power-law formulation.

Choi, Sung R.↗
NASA NTRS2002

Slow Crack Growth of Brittle Materials With Exponential Crack-Velocity Formulation: Constant Stress and Cyclic Stress Experiments - Part 3

The previously determined life prediction analysis based on an exponential crack-velocity formulation was examined using a variety of experimental data on advanced structural ceramics tested under constant stress and cyclic stress loading at ambient and elevated temperatures. The data fit to the relation between the time to failure and applied stress (or maximum applied stress in cyclic loading) was very reasonable for most of the materials studied. It was also found that life prediction for cyclic stress loading from data of constant stress loading in the exponential formulation was in good agreement with the experimental data, resulting in a similar degree of accuracy as compared with the power-law formulation. The major limitation in the exponential crack-velocity formulation, however, was that the inert strength of a material must be known a priori to evaluate the important slow-crack-growth (SCG) parameter n, a significant drawback as compared with the conventional power-law crack-velocity formulation.

Choi, Sung R.↗
NASA NTRS2005

Completed Beltrami-Michell Formulation in Polar Coordinates

A set of conditions had not been formulated on the boundary of an elastic continuum since the time of Saint-Venant. This limitation prevented the formulation of a direct stress calculation method in elasticity for a continuum with a displacement boundary condition. The missed condition, referred to as the boundary compatibility condition, is now formulated in polar coordinates. The augmentation of the new condition completes the Beltrami-Michell formulation in polar coordinates. The completed formulation that includes equilibrium equations and a compatibility condition in the field as well as the traction and boundary compatibility condition is derived from the stationary condition of the variational functional of the integrated force method. The new method is illustrated by solving an example of a mixed boundary value problem for mechanical as well as thermal loads.

Patnaik, Surya N.↗
NASA NTRS2012

A New and General Formulation of the Parametric HFGMC Micromechanical Method for Three-Dimensional Multi-Phase Composites

The recent two-dimensional (2-D) parametric formulation of the high fidelity generalized method of cells (HFGMC) reported by the authors is generalized for the micromechanical analysis of three-dimensional (3-D) multiphase composites with periodic microstructure. Arbitrary hexahedral subcell geometry is developed to discretize a triply periodic repeating unit-cell (RUC). Linear parametric-geometric mapping is employed to transform the arbitrary hexahedral subcell shapes from the physical space to an auxiliary orthogonal shape, where a complete quadratic displacement expansion is performed. Previously in the 2-D case, additional three equations are needed in the form of average moments of equilibrium as a result of the inclusion of the bilinear terms. However, the present 3-D parametric HFGMC formulation eliminates the need for such additional equations. This is achieved by expressing the coefficients of the full quadratic polynomial expansion of the subcell in terms of the side or face average-displacement vectors. The 2-D parametric and orthogonal HFGMC are special cases of the present 3-D formulation. The continuity of displacements and tractions, as well as the equilibrium equations, are imposed in the average (integral) sense as in the original HFGMC formulation. Each of the six sides (faces) of a subcell has an independent average displacement micro-variable vector which forms an energy-conjugate pair with the transformed average-traction vector. This allows generating symmetric stiffness matrices along with internal resisting vectors for the subcells which enhances the computational efficiency. The established new parametric 3-D HFGMC equations are formulated and solution implementations are addressed. Several applications for triply periodic 3-D composites are presented to demonstrate the general capability and varsity of the present parametric HFGMC method for refined micromechanical analysis by generating the spatial distributions of local stress fields. These applications include triply periodic composites with inclusions in the form of a cavity, spherical inclusion, ellipsoidal inclusion, discontinuous aligned short fiber. A 3-D repeating unit-cell for foam material composite is simulated.

Haj-Ali, Rami↗
NASA NTRS2015

Compact Assumption Applied to the Monopole Term of Farassat's Formulations

Farassat's formulations provide an acoustic prediction at an observer location provided a source surface, including motion and flow conditions. This paper presents compact forms for the monopole term of several of Farassat's formulations. When the physical surface is elongated, such as the case of a high aspect ratio rotorcraft blade, compact forms can be derived which are shown to be a function of the blade cross sectional area by reducing the computation from a surface integral to a line integral. The compact forms of all formulations are applied to two example cases: a short span wing with constant airfoil cross section moving at three forward flight Mach numbers and a rotor at two advance ratios. Acoustic pressure time histories and power spectral densities of monopole noise predicted from the compact forms of all the formulations at several observer positions are shown to compare very closely to the predictions from their non-compact counterparts. A study on the influence of rotorcraft blade shape on the high frequency portion of the power spectral density shows that there is a direct correlation between the aspect ratio of the airfoil and the error incurred by using the compact form. Finally, a prediction of pressure gradient from the non-compact and compact forms of the thickness term of Formulation G1A shows that using the compact forms results in a 99.6% improvement in computation time, which will be critical when noise is incorporated into a design environment.

Lopes, Leonard V.↗