Search NASA⌕ Search

SEARCH · Search NASA

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 235 records · Page 13

Evaluation of candidate working fluid formulations for the electrothermal - chemical wind tunnel

Various candidate chemical formulations are evaluated as a precursor for the working fluid to be used in the electrothermal hypersonic test facility which was under study at the NASA LaRC Hypersonic Propulsion Branch, and the formulations which would most closely satisfy the goals set for the test facility are identified. Out of the four tasks specified in the original proposal, the first two, literature survey and collection of kinetic data, are almost completed. The third task, work on a mathematical model of the ET wind tunnel operation, was started and concentrated on the expansion in the nozzle with finite rate kinetics.

Akyurtlu, Jale F.↗

Double simple-harmonic-oscillator formulation of the thermal equilibrium of a fluid interacting with a coherent source of phonons

A formulation is given for a collection of phonons (sound) in a fluid at a non-zero temperature which uses the simple harmonic oscillator twice; one to give a stochastic thermal 'noise' process and the other which generates a coherent Glauber state of phonons. Simple thermodynamic observables are calculated and the acoustic two point function, 'contrast' is presented. The role of 'coherence' in an equilibrium system is clarified by these results and the simple harmonic oscillator is a key structure in both the formulation and the calculations.

Defacio, B.↗

Analytical formulation of optimum rotor interdisciplinary design with a three-dimensional wake

An analytical formulation of optimum rotor interdisciplinary design is presented. A finite-state aeroelastic rotor model, coupling generalized dynamic wake with blade finite elements, is applied to perform the optimum rotor blade design for improved aerodynamic performance and vehicle vibration, while a feasible direction nonlinear optimizer, CONMIN, provides the optimization algorithm. The approach features a systematic rotor aeroelastic model which offers an efficient analytical tool, and retains necessary aerodynamic and blade dynamic building blocks for a sufficient rotor dynamic response analysis. The formulation is well suited for an efficient design sensitivity computation without resorting to finite difference, and thus provides a practical design tool. The results show improved rotor aerodynamic performance and reduced hub vibratory loads for the optimized blade as compared to the advanced rotor of reference design.

He, Chengjian↗

Improved solar models constructed with a formulation of convection for stellar structure and evolution calculations without the mixing-length theory approximations

We have updated a previous attempt to incorporate within a solar model a treatment of convection based upon numerical simulations of convection rather than mixing-length theory (MLT). We have modified our formulation of convection for a better treatment of the kinetic energy flux. Our solar model has been updated to include a complete range of OPAL opacities, the Debye-Hueckel correction to the equation of state, helium diffusion due to gravitational settling, and atmospheres by Kurucz. We construct a series of models using both MLT and our revised formulation of convection and the compared results to measurements of the solar radius, the solar luminosity, and the depth of the solar convection zone as inferred from helioseismology. We find X(solar) = 0.702 +/- 0.005, Y(solar) = 0.278 +/- 0.005, and Z(solar) = 0.0193 +/- 0.0005.

Lydon, Thomas J.↗

Stiffness, thermal expansion, and thermal bending formulation of stiffened, fiber-reinforced composite panels

A method is presented for formulating stiffness terms and thermal coefficients of stiffened, fiber-reinforced composite panels. The method is robust enough to handle panels with general cross sectional shapes, including those which are unsymmetric and/or unbalanced. Nonlinear, temperature and load dependent constitutive material data of each laminate are used to 'build-up' the stiffened panel membrane, bending, and membrane-bending coupling stiffness terms and thermal coefficients. New thermal coefficients are introduced to quantify panel response from through-the-thickness temperature gradients. A technique of implementing this capability with a single plane of shell finite elements using the MSC/NASTRAN analysis program (FEA) is revealed that provides accurate solutions of entire airframes or engines with coarsely meshed models. An example of a composite, hat-stiffened panel is included to demonstrate errors that occur when an unsymmetric panel is symmetrically formulated as traditionally done. The erroneous results and the correct ones produced from this method are compared to analysis from discretely meshed three-dimensional FEA.

Collier, Craig S.↗

A comparison of two formulations for high-order accurate essentially non-oscillatory schemes

The finite-volume and finite-difference implementations of high-order accurate essentially nonoscillatory shock-capturing schemes are discussed and compared. Results obtained with fourth-order accurate algorithms based on both formulations are examined for accuracy, sensitivity to grid irregularities, resolution of waves that are oblique to the mesh, and computational efficiency. Some algorithm modifications that may be required for a given application are suggested. Conclusions that pertain to the relative merits of both formulations are drawn, and some circumstances for which each might be useful are noted.

Casper, J.↗

Symmetric Galerkin boundary formulations employing curved elements

Accounts of the symmetric Galerkin approach to boundary element analysis (BEA) have recently been published. This paper attempts to add to the understanding of this method by addressing a series of fundamental issues associated with its potential computational efficiency. A new symmetric Galerkin theoretical formulation for both the (harmonic) heat conduction and the (biharmonic) elasticity problem that employs regularized singular and hypersingular boundary integral equations (BIEs) is presented. The novel use of regularized BIEs in the Galerkin context is shown to allow straightforward incorporation of curved, isoparametric elements. A symmetric reusable intrinsic sample point (RISP) numerical integration algorithm is shown to produce a Galerkin (i.e., double) integration strategy that is competitive with its counterpart (i.e., singular) integration procedure in the collocation BEA approach when the time saved in the symmetric equation solution phase is also taken into account. This new formulation is shown to be capable of employing hypersingular BIEs while obviating the requirement of C 1 continuity, a fact that allows the employment of the popular continuous element technology. The behavior of the symmetric Galerkin BEA method with regard to both direct and iterative equation solution operations is also addressed. A series of example problems are presented to quantify the performance of this symmetric approach, relative to the more conventional unsymmetric BEA, in terms of both accuracy and efficiency. It is concluded that appropriate implementations of the symmetric Galerkin approach to BEA indeed have the potential to be competitive with, if not superior to, collocation-based BEA, for large-scale problems.

Kane, J. H.↗

A comparison of two formulations for high-order accurate essentially non-oscillatory schemes

The finite-volume and finite-difference implementations of high-order accurate essentially non-oscillatory shock-capturing schemes are discussed and compared. Results obtained with fourth-order accurate algorithms based on both formulations are examined for accuracy, sensitivity to grid irregularities, resolution of waves that are oblique to the mesh, and computational efficiency. Some algorithm modifications that may be required for a given application are suggested. Conclusions that pertain to the relative merits of both formulations are drawn, and some circumstances for which each might be useful are noted.

Casper, Jay↗

Cryogenic Gellant and Fuel Formulation for Metallized Gelled Propellants: Hydrocarbons and Hydrogen with Aluminum

An experimental program to determine the viability of nanoparticulate gellant materials for gelled hydrocarbons and gelled liquid hydrogen was conducted. The gellants included alkoxides (BTMSE and BTMSH) and silica-based materials. Hexane, ethane, propane and hydrogen were gelled with the newly-formulated materials and their rheological properties were determined: shear stress versus shear rate and their attendant viscosities. Metallized hexane with aluminum particles was also rheologically characterized. The propellant and gellant formulations were selected for the very high surface area and relatively-high energy content of the gellants. These new gellants can therefore improve rocket engine specific impulse over that obtained with traditional cryogenic-fuel gellant materials silicon dioxide, frozen methane, or frozen ethane particles. Significant reductions in the total mass of the gellant were enabled in the fuels. In gelled liquid hydrogen, the total mass of gellant was reduced from 10-40 wt percent of frozen hydrocarbon particles to less that 8 wt percent with the alkoxide.

Wong, Wing↗

A survey of the core-congruential formulation for geometrically nonlinear TL finite elements

This article presents a survey of the core-congruential formulation (CCF) for geometrically nonlinear mechanical finite elements based on the total Lagrangian (TL) kinematic description. Although the key ideas behind the CCF can be traced back to Rajasekaran and Murray in 1973, it has not subsequently received serious attention. The CCF is distinguished by a two-phase development of the finite element stiffness equations. The initial phase developed equations for individual particles. These equations are expressed in terms of displacement gradients as degrees of freedom. The second phase involves congruential-type transformations that eventually binds the element particles of an individual element in terms of its node-displacement degrees of freedom. Two versions of the CCF, labeled direct and generalized, are distinguished. The direct CCF (DCCF) is first described in general form and then applied to the derivation of geometrically nonlinear bar, and plane stress elements using the Green-Lagrange strain measure. The more complex generalized CCF (GCCF) is described and applied to the derivation of 2D and 3D Timoshenko beam elements. Several advantages of the CCF, notably the physically clean separation of material and geometric stiffnesses, and its independence with respect to the ultimate choice of shape functions and element degrees of freedom, are noted. Application examples involving very large motions solved with the 3D beam element display the range of applicability of this formulation, which transcends the kinematic limitations commonly attributed to the TL description.

Felippa, Carlos A.↗

A three-dimensional nonlinear Timoshenko beam based on the core-congruential formulation

A three-dimensional, geometrically nonlinear two-node Timoshenkoo beam element based on the total Larangrian description is derived. The element behavior is assumed to be linear elastic, but no restrictions are placed on magnitude of finite rotations. The resulting element has twelve degrees of freedom: six translational components and six rotational-vector components. The formulation uses the Green-Lagrange strains and second Piola-Kirchhoff stresses as energy-conjugate variables and accounts for the bending-stretching and bending-torsional coupling effects without special provisions. The core-congruential formulation (CCF) is used to derived the discrete equations in a staged manner. Core equations involving the internal force vector and tangent stiffness matrix are developed at the particle level. A sequence of matrix transformations carries these equations to beam cross-sections and finally to the element nodal degrees of freedom. The choice of finite rotation measure is made in the next-to-last transformation stage, and the choice of over-the-element interpolation in the last one. The tangent stiffness matrix is found to retain symmetry if the rotational vector is chosen to measure finite rotations. An extensive set of numerical examples is presented to test and validate the present element.

Crivelli, Luis A.↗

Coupling equivalent plate and finite element formulations in multiple-method structural analyses

A coupled multiple-method analysis procedure for use late in conceptual design or early in preliminary design of aircraft structures is described. Using this method, aircraft wing structures are represented with equivalent plate models, and structural details such as engine/pylon structure, landing gear, or a 'stick' model of a fuselage are represented with beam finite element models. These two analysis methods are implemented in an integrated multiple-method formulation that involves the assembly and solution of a combined set of linear equations. The corresponding solution vector contains coefficients of the polynomials that describe the deflection of the wing and also the components of translations and rotations at the joints of the beam members. Two alternative approaches for coupling the methods are investigated; one using transition finite elements and the other using Lagrange multipliers. The coupled formulation is applied to the static analysis and vibration analysis of a conceptual design model of a fighter aircraft. The results from the coupled method are compared with corresponding results from an analysis in which the entire model is composed of finite elements.

VIBRATION ANALYSIS↗

Hybrid fuel formulation and technology development

The objective was to develop an improved hybrid fuel with higher regression rate, a regression rate expression exponent close to 0.5, lower cost, and higher density. The approach was to formulate candidate fuels based on promising concepts, perform thermomechanical analyses to select the most promising candidates, develop laboratory processes to fabricate fuel grains as needed, fabricate fuel grains and test in a small lab-scale motor, select the best candidate, and then scale up and validate performance in a 2500 lbf scale, 11-inch diameter motor. The characteristics of a high performance fuel have been verified in 11-inch motor testing. The advanced fuel exhibits a 15% increase in density over an all hydrocarbon formulation accompanied by a 50% increase in regression rate, which when multiplied by the increase in density yields a 70% increase in fuel mass flow rate; has a significantly lower oxidizer-to-fuel (O/F) ratio requirement at 1.5; has a significantly decreased axial regression rate variation making for more uniform propellant flow throughout motor operation; is very clean burning; extinguishes cleanly and quickly; and burns with a high combustion efficiency.

Dean, D. L.↗

Constrained Multipoint Aerodynamic Shape Optimization Using an Adjoint Formulation and Parallel Computers

An aerodynamic shape optimization method that treats the design of complex aircraft configurations subject to high fidelity computational fluid dynamics (CFD), geometric constraints and multiple design points is described. The design process will be greatly accelerated through the use of both control theory and distributed memory computer architectures. Control theory is employed to derive the adjoint differential equations whose solution allows for the evaluation of design gradient information at a fraction of the computational cost required by previous design methods. The resulting problem is implemented on parallel distributed memory architectures using a domain decomposition approach, an optimized communication schedule, and the MPI (Message Passing Interface) standard for portability and efficiency. The final result achieves very rapid aerodynamic design based on a higher order CFD method. In order to facilitate the integration of these high fidelity CFD approaches into future multi-disciplinary optimization (NW) applications, new methods must be developed which are capable of simultaneously addressing complex geometries, multiple objective functions, and geometric design constraints. In our earlier studies, we coupled the adjoint based design formulations with unconstrained optimization algorithms and showed that the approach was effective for the aerodynamic design of airfoils, wings, wing-bodies, and complex aircraft configurations. In many of the results presented in these earlier works, geometric constraints were satisfied either by a projection into feasible space or by posing the design space parameterization such that it automatically satisfied constraints. Furthermore, with the exception of reference 9 where the second author initially explored the use of multipoint design in conjunction with adjoint formulations, our earlier works have focused on single point design efforts. Here we demonstrate that the same methodology may be extended to treat complete configuration designs subject to multiple design points and geometric constraints. Examples are presented for both transonic and supersonic configurations ranging from wing alone designs to complex configuration designs involving wing, fuselage, nacelles and pylons.

Reuther, James↗

A Mixed Multi-Field Finite Element Formulation for Thermopiezoelectric Composite Shells

Analytical formulations are presented which account for the coupled mechanical, electrical, and thermal response of piezoelectric composite shell structures. A new mixed multi-field laminate theory is developed which combines "single layer" assumptions for the displacements along with layerwise fields for the electric potential and temperature. This laminate theory is formulated using curvilinear coordinates and is based on the principles of linear thermopiezoelectricity. The mechanics have the inherent capability to explicitly model both the active and sensory responses of piezoelectric composite shells in thermal environment. Finite element equations are derived and implemented for an eight-noded shell element. Numerical studies are conducted to investigate both the sensory and active responses of piezoelectric composite shell structures subjected to thermal loads. Results for a cantilevered plate with an attached piezoelectric layer are com- pared with corresponding results from a commercial finite element code and a previously developed program. Additional studies are conducted on a cylindrical shell with an attached piezoelectric layer to demonstrate capabilities to achieve thermal shape control on curved piezoelectric structures.

Lee, Ho-Jun↗

Numerical Modeling of Crystal of ZnSe by Physical Vapor Transport - Towards a more Comprehensive Formulations

Crystal growth from the vapor phase has various advantages over melt growth. The main advantage is from a lower processing temperature which makes the process more amenable in instances where the melting temperature of the crystal is high. Other benefits stem from the inherent purification mechanism in the process due to differences in the vapor pressures of the native elements and impurities, and the enhanced interfacial morphological stability during the growth process. Further, the implementation of PVT growth in closed ampoules affords experimental simplicity with minimal needs for complex process control which makes it an ideal candidate for space investigations in systems where gravity tends to have undesirable effects on the growth process. Bulk growth of wide band gap II-VI semiconductors by physical vapor transport has been developed and refined over the past several years at NASA MSFC. Results from a modeling study of PVT crystal growth of ZnSe are reported in this paper. The PVT process is numerically investigated using both two-dimensional and fully three-dimensional formulation of the governing equations and associated boundary conditions. Both the incompressible Boussinesq approximation and the compressible model are tested to determine the influence of gravity on the process and to discern the differences between the two approaches. The influence of a residual gas is included in the models. The results show that both the incompressible and compressible approximations provide comparable results and the presence of a residual gas tends to measurably reduce the mass flux in the system. Detailed flow, thermal and concentration profiles will be provided in the final manuscript along with computed heat and mass transfer rates. Comparisons with the 1-D model will also be provided. The effect of gravity on the process from numerical computations shows subtle effects although experimental evidence from vertically and horizontally grown samples show dramatic evidence of gravitational effects. The shortcomings of the problem formulation will be discussed and a framework will be provided leading up towards a more comprehensive model of PVT systems.

Ramachandran, N.↗

Fuel Plume Image Mixing Analysis Formulation With Proper Treatment of Non-Constant Velocity Flowfields

A previously developed technique allows an estimate of integral mixing to be obtained from an image of laser scattered light from particle seeded fuel in the hypervelocity flow through a scramjet combustor. This previous mixing analysis formulation contains an assumption of a constant velocity flowfield across the plane of the fuel plume image. For high-speed scramjet combustors, the velocity flowfield is quite uniform and an assumption of constant velocity works well. Applying this same mixing analysis technique to fuel plume images obtained from a mid-speed scramjet combustor makes it desirable to remove the constant velocity assumption. This is due to the non-uniform velocity flowfields present in mid-speed scramjet combustors. A new formulation of the mixing analysis methodology is developed and presented so that the technique can be applied to a mid-speed scramjet combustor without the need to assume a constant velocity flowfield.

Mekkes, Gregory L.↗

Review of the Toxicity of the Metallic Compounds Proposed for Use in the New Green Missile Formulation

In an effort to modernize and minimize hazards posed by the toxic components of missile propellant, the USACHPPM has been tasked to provide a comparison of the toxicity of compounds currently in use as missile propellants and the suite of compounds proposed to replace them. This report deals with the portion of this work concerning the toxicity of the organometallic compounds used in these formulations. Toxicity assessments of the organic compounds used in these formulations are published elsewhere. In general, toxicity data were available for all the metal compounds of concern or for closely related compounds that can serve as surrogates for the assessment of toxicity. We have high confidence in the reliability of these comparisons. This report is organized by element to provide the reader with an in-depth assessment with a minimum of redundancy. The narrative will first describe general concepts about the toxicity of each metal and then provide a summary of the toxicological information available for the specific compound.

Major, Michael A.↗