Search NASA⌕ Search

SEARCH · Search NASA

Results for “Semi-analytical 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2026

Developing a Pyrolysis Gas Thermal Blocking Model for Reentry Demise

In NASA’s Object Reentry Survival Analysis Tool (ORSAT), aerodynamic drag and aerothermal heating coefficients are computed for each of the free-molecular, continuum, and transitional flow regimes using analytical and semi-analytical methods. These heating coefficients were derived for typical metallic materials that melt and do not have a strong gas-phase contribution to the flow in the boundary layer. Modern satellites typically feature fiber-reinforced polymer (FRP) components, such as solar array booms, facesheets of sandwich panels, or overwraps for composite-overwrapped pressure vessels (COPV). These FRP materials do not behave the same as metals in the reentry environment, but instead will pyrolyze and develop significant volumes of gas into the boundary layer. Accurately predicting the reentry demise of FRP components is critical to assessing the reentry casualty risk for modern spacecraft. Research in recent years has shown that this demisability can depend heavily on how the expulsion of gaseous pyrolysis products through the outer surface of the material affects the heat flux at the surface. The ODPO has been developing a reduced-order model of the effect of pyrolysis gas blowing on the heat flux based on correlations between a blowing factor and a non-dimensional heat flux to be incorporated in the upcoming version 7.3 of the Object Reentry Survivability Analysis Tool (ORSAT). This presentation discusses the progress of this development project and the challenges remaining for generalizing the model across families of FRP materials.

Benton Greene↗
NASA NTRS1989

Error analysis for semi-analytic displacement derivatives with respect to shape and sizing variables

Sensitivity analysis is fundamental to the solution of structural optimization problems. Consequently, much research has focused on the efficient computation of static displacement derivatives. As originally developed, these methods relied on analytical representations for the derivatives of the structural stiffness matrix (K) with respect to the design variables (b sub i). To extend these methods for use with complex finite element formulations and facilitate their implementation into structural optimization programs using the general finite element method analysis codes, the semi-analytic method was developed. In this method the matrix the derivative of K/the derivative b sub i is approximated by finite difference. Although it is well known that the accuracy of the semi-analytic method is dependent on the finite difference parameter, recent work has suggested that more fundamental inaccuracies exist in the method when used for shape optimization. Another study has argued qualitatively that these errors are related to nonuniform errors in the stiffness matrix derivatives. The accuracy of the semi-analytic method is investigated. A general framework was developed for the error analysis and then it is shown analytically that the errors in the method are entirely accounted for by errors in delta K/delta b sub i. Furthermore, it is demonstrated that acceptable accuracy in the derivatives can be obtained through careful selection of the finite difference parameter.

Fenyes, Peter A.↗
NASA NTRS1992

An analytical evaluation of mascon effects in a semi-analytic theory

A semianalytic method to predict the orbit of a spacecraft orbiting a planet and perturbed by oblateness and a gravity anomaly (mascon) is presented. The Hamiltonian is first expressed in closed form in terms of the nonsingular canonical set of Poincare elements where the mascon is represented by a point-mass. The planet global potential is limited to the second and third harmonics. The mascon potential in harmonics is limited to second degree and first order. The original Lie-Hori method is then applied to transform the canonical equations of motion from the osculating element space to the mean element space. Numerical integration of the new equations is used.

Salama, Ahmed H.↗
NASA NTRS1990

Acoustic radiation from lifting airfoils in compressible subsonic flow

The far field acoustic radiation from a lifting airfoil in a three-dimensional gust is studied. The acoustic pressure is calculated using the Kirchhoff method, instead of using the classical acoustic analogy approach due to Lighthill. The pressure on the Kirchhoff surface is calculated using an existing numerical solution of the unsteady flow field. The far field acoustic pressure is calculated in terms of these values using Kirchhoff's formula. The method is validated against existing semi-analytical results for a flat plate. The method is then used to study the problem of an airfoil in a harmonic three-dimensional gust, for a wide range of Mach numbers. The effect of variation of the airfoil thickness and angle of attack on the acoustic far field is studied. The changes in the mechanism of sound generation and propagation due to the presence of steady loading and nonuniform mean flow are also studied.

Atassi, Hafiz M.↗
NASA NTRS1990

Acoustic radiation from lifting airfoils in compressible subsonic flow

The far field acoustic radiation from a lifting airfoil in a three-dimensional gust is studied. The acoustic pressure is calculated using the Kirchhoff method, instead of using the classical acoustic analogy approach due to Lighthill. The pressure on the Kirchhoff surface is calculated using an existing numerical solution of the unsteady flow field. The far field acoustic pressure is calculated in terms of these values using Kirchhoff's formula. The method is validated against existing semi-analytical results for a flat plate. The method is then used to study the problem of an airfoil in a harmonic three-dimensional gust, for a wide range of Mach numbers. The effect of variation of the airfoil thickness and angle of attack on the acoustic far field is studied. The changes in the mechanism of sound generation and propagation due to the presence of steady loading and non-uniform mean flow are also studied.

Atassi, Hafiz M.↗
NASA NTRS1979

Determination of the ground albedo and the index of absorption of atmospheric particulates by remote sensing. I - Theory

A statistical technique is developed for inferring the optimum values of the ground albedo and the effective imaginary term of the complex refractive index of atmospheric particulates. The procedure compares measurements of the ratio of the hemispheric diffuse to directly transmitted solar flux density at the earth's surface with radiative transfer computations of the same as suggested by Herman et al. (1975). A detailed study is presented which shows the extent to which the ratio of diffuse to direct solar radiation is sensitive to many of the radiative transfer parameters. Results indicate that the optical depth and size distribution of atmospheric aerosol particles are the two parameters which uniquely specify the radiation field to the point where ground albedo and index of absorption can be inferred. Varying the real part of the complex refractive index of atmospheric particulates as well as their vertical distribution is found to have a negligible effect on the diffuse-direct ratio. The statistical procedure utilizes a semi-analytic gradient search method from least-squares theory and includes a detailed error analysis.

King, M. D.↗
NASA NTRS1984

Quantification of Saint-Venant's Principle for a General Prismatic Member

A semi-analytic finite element method for finding the decay rates of edge effects in anisotropic layered plates was extended to the treatment of prismatic members of an arbitrary cross section composed of any number of anisotropic elements. The displacement fields corresponding to the decay rates with positive real part were determined for a homogenous graphite/epoxy plate two units thick. The behavior of a rectangular beam is also considered.

Source record↗
NASA NTRS1985

Quantification of Saint-Venant's principle for a general prismatic member

A method to quantify Saint-Venant's principle for general prismatic members is described. Finite elements are used to model the member cross-section to allow the analysis of arbitrary geometries, material properties, and nonhomogeneities. The variation in the structural response along the member's length, away from the loaded end, is modelled with an exponential decay which has been established in the literature to be the correct form. Many analytic results exist in the literature for special cases. Because of its ability to handle arbitrary cross sectional geometries, this semi-analytic finite element method can provide a firm footing for the utilization of Saint-Venant's principle for whole new classes of problems. Results for several typical problems are given.

Goetschel, D. B.↗
NASA NTRS1988

Accuracy problems associated with semi-analytical derivatives of static response

The semianalytical method is widely used for calculating derivatives of static response with respect to design variables for structures modeled by finite elements. This paper shows that the method can have serious accuracy problems for shape design variables in structures modeled by beam elements. The problem was first discovered in the analysis of a complex car model which is described. Next, it is shown that the same accuracy problem occurs for the simplest beam structures. Finally, the problem is shown to be associated with a high ratio of rigid body rotation to elastic deformation, and a test is developed to diagnose the severity of the problem in a given structure. Several examples are presented to show the validity of this test.

Barthelemy, Bruno↗
NASA NTRS1991

A method for the geometrically nonlinear analysis of compressively loaded prismatic composite structures

A method was developed for the geometrically nonlinear analysis of the static response of thin-walled stiffened composite structures loaded in uniaxial or biaxial compression. The method is applicable to arbitrary prismatic configurations composed of linked plate strips, such as stiffened panels and thin-walled columns. The longitudinal ends of the structure are assumed to be simply supported, and geometric shape imperfections can be modeled. The method can predict the nonlinear phenomena of postbuckling strength and imperfection sensitivity which are exhibited by some buckling-dominated structures. The method is computer-based and is semi-analytic in nature, making it computationally economical in comparison to finite element methods. The method uses a perturbation approach based on the use of a series of buckling mode shapes to represent displacement contributions associated with nonlinear response. Displacement contributions which are of second order in the model amplitudes are incorported in addition to the buckling mode shapes. The principle of virtual work is applied using a finite basis of buckling modes, and terms through the third order in the model amplitudes are retained. A set of cubic nonlinear algebraic equations are obtained, from which approximate equilibrium solutions are determined. Buckling mode shapes for the general class of structure are obtained using the VIPASA analysis code within the PASCO stiffened-panel design code. Thus, subject to some additional restrictions in loading and plate anisotropy, structures which can be modeled with respect to buckling behavior by VIPASA can be analyzed with respect to nonlinear response using the new method. Results obtained using the method are compared with both experimental and analytical results in the literature. The configurations investigated include several different unstiffened and blade-stiffening panel configurations, featuring both homogeneous, isotropic materials, and laminated composite material.

Stoll, Frederick↗
NASA NTRS1990

Computational aspects of sensitivity calculations in linear transient structural analysis

A study was performed focusing on the calculation of sensitivities of displacements, velocities, accelerations, and stresses in linear, structural, transient response problems. One significant goal of the study was to develop and evaluate sensitivity calculation techniques suitable for large-order finite element analyses. Accordingly, approximation vectors such as vibration mode shapes are used to reduce the dimensionality of the finite element model. Much of the research focused on the accuracy of both response quantities and sensitivities as a function of number of vectors used. Two types of sensitivity calculation techniques were developed and evaluated. The first type of technique is an overall finite difference method where the analysis is repeated for perturbed designs. The second type of technique is termed semi-analytical because it involves direct, analytical differentiation of the equations of motion with finite difference approximation of the coefficient matrices. To be computationally practical in large-order problems, the overall finite difference methods must use the approximation vectors from the original design in the analyses of the perturbed models. In several cases this fixed mode approach resulted in very poor approximations of the stress sensitivities. Almost all of the original modes were required for an accurate sensitivity and for small numbers of modes, the accuracy was extremely poor. To overcome this poor accuracy, two semi-analytical techniques were developed. The first technique accounts for the change in eigenvectors through approximate eigenvector derivatives. The second technique applies the mode acceleration method of transient analysis to the sensitivity calculations. Both result in accurate values of the stress sensitivities with a small number of modes and much lower computational costs than if the vibration modes were recalculated and then used in an overall finite difference method.

Greene, William H.↗
NASA NTRS1980

An Economical Semi-Analytical Orbit Theory for Retarded Satellite Motion About an Oblate Planet

Brouwer and Brouwer-Lyddanes' use of the Von Zeipel-Delaunay method is employed to develop an efficient analytical orbit theory suitable for microcomputers. A succinctly simple pseudo-phenomenologically conceptualized algorithm is introduced which accurately and economically synthesizes modeling of drag effects. The method epitomizes and manifests effortless efficient computer mechanization. Simulated trajectory data is employed to illustrate the theory's ability to accurately accommodate oblateness and drag effects for microcomputer ground based or onboard predicted orbital representation. Real tracking data is used to demonstrate that the theory's orbit determination and orbit prediction capabilities are favorably adaptable to and are comparable with results obtained utilizing complex definitive Cowell method solutions on satellites experiencing significant drag effects.

Gordon, R. A.↗
NASA NTRS1992

Sensitivity analysis of dynamic aeroelastic responses

This paper summarizes ongoing research on the sensitivity analysis of dynamic aeroelastic response of wings. Two approaches are being used to express the unsteady aerodynamic loads: (1) the frequency-domain approach, and (2) the state-space approach. The frequency-domain approach is demonstrated on a three-dimensional box wing and the state-space domain approach is demonstrated on a simple two-dimensional sectional model. Three different methods are used to find the sensitivities: (1) a purely finite difference approach, (2) a semi-analytical approach in which an analytical expression is used for calculating the sensitivity of an eigenvalue of the complex valued aeroelastic matrix, however, the derivatives of the components of the matrix are obtained using finite difference, and (3) a semi-analytic approach that differs from (2) in the sense that the sensitivity of the aerodynamic matrix is now obtained analytically. A good agreement is seen between the three sets of results. For the two-dimensional sectional model the results for the sensitivities of the flutter speed with respect to the various parameters are obtained.

Kapania, Rakesh K.↗
NASA NTRS1989

An algebraic homotopy method for generating quasi-three-dimensional grids for high-speed configurations

A fast and versatile procedure for algebraically generating boundary conforming computational grids for use with finite-volume Euler flow solvers is presented. A semi-analytic homotopic procedure is used to generate the grids. Grids generated in two-dimensional planes are stacked to produce quasi-three-dimensional grid systems. The body surface and outer boundary are described in terms of surface parameters. An interpolation scheme is used to blend between the body surface and the outer boundary in order to determine the field points. The method, albeit developed for analytically generated body geometries is equally applicable to other classes of geometries. The method can be used for both internal and external flow configurations, the only constraint being that the body geometries be specified in two-dimensional cross-sections stationed along the longitudinal axis of the configuration. Techniques for controlling various grid parameters, e.g., clustering and orthogonality are described. Techniques for treating problems arising in algebraic grid generation for geometries with sharp corners are addressed. A set of representative grid systems generated by this method is included. Results of flow computations using these grids are presented for validation of the effectiveness of the method.

Moitra, Anutosh↗
NASA NTRS1990

A hybrid perturbation-Galerkin technique for partial differential equations

A two-step hybrid perturbation-Galerkin technique for improving the usefulness of perturbation solutions to partial differential equations which contain a parameter is presented and discussed. In the first step of the method, the leading terms in the asymptotic expansion(s) of the solution about one or more values of the perturbation parameter are obtained using standard perturbation methods. In the second step, the perturbation functions obtained in the first step are used as trial functions in a Bubnov-Galerkin approximation. This semi-analytical, semi-numerical hybrid technique appears to overcome some of the drawbacks of the perturbation and Galerkin methods when they are applied by themselves, while combining some of the good features of each. The technique is illustrated first by a simple example. It is then applied to the problem of determining the flow of a slightly compressible fluid past a circular cylinder and to the problem of determining the shape of a free surface due to a sink above the surface. Solutions obtained by the hybrid method are compared with other approximate solutions, and its possible application to certain problems associated with domain decomposition is discussed.

Geer, James F.↗
NASA NTRS1992

Elastic Analysis of Three-Dimensional Solids with Fiber Inclusions by Bem

A very efficient method of analysis for the elastic behavior of three-dimensional solids with a large number of embedded fiber inclusions has been developed. For specific applications to the analysis of cylindrical shaped fiber inclusions significant gain in efficiency can be achieved by defining these inclusions by a system of curvilinear line elements with a prescribed diameter and by assuming a variation in the traction and displacement field in the circumferential directions in terms of a trigonometric shape function together with a linear or quadratic variation along the longitudinal direction. The resulting integrals are then treated semi-analytically. A number of examples of the analysis which has been developed for composite elements are described.

Banerjee, P. K.↗
NASA NTRS1968

Internal low Reynolds number turbulent heat transfer

The results of a semi-analytical and experimental study of adiabatic tube flow and an analytical and experimental investigation of the thermal entrance region for gas flowing through electrically heated circular tubes are presented. Emphasis is placed on the low Reynolds number turbulent flow regime--defined as fully turbulent flow at bulk Reynolds numbers from 3,000 to about 15,000. Adiabatic air velocity and friction data and localized heat transfer measurements for air and helium, at low heating rates, are presented for this range. The adiabatic data were obtained in a 1.61 inch 1D tube for flow at bulk Reynolds numbers from 3,000 to 15,000. A continuous, Reynolds number-dependent, profile is developed from the data by using a modification of Reichardts' wall and middle law eddy diffusivity expressions. The velocity profile satisfies continuity. It is valid for all Reynolds numbers in excess of 3,000 for which the flow is fully turbulent and the Blasius friction factor expression is valid. The thermal entrance problem for a fully developed velocity profile is solved analytically by the method of Sparrow, Hallman, and Siegel. The solution is based on the profile developed from the velocity study. Tabular values of the eigenvalues and normalized Nusselt numbers for gases are presented for it range of Reynolds numbers from 3,000 to 50,000. The axial variation of Nusselt number is found to be correlated by [Nu/Nu_(∞)] = 1 + 0.8 (1+ 70,000 Re^(-3/2)) ((x/D)^(-1)) to within ± 5 per cent for x/D ≥2. The fully developed value agrees with the Dittus - Boelter correlation, Nu_(∞) = 0.021 Re^(0.8) Pr^(0.4) For the eigenvalues, λ^(2)_(n), and the associated constants, A_(n), correlations of the form λ^(2)_(n) =A_(1,n) Re^(-b_(1,n) + C_(1,n) Re^(-d_(1,n)) A_(n) = -A_(2,n) Re^(-b_(2,n) + C_(2,n) Re^(-d_(2,n)) are obtained. The coefficients and powers are presented in tabular form. Heat transfer data are presented, primarily for helium, for the low Reynolds number turbulent range. A one-quarter inch, resistively heated, vertical, circular tube was used for the study. The data cover an axial range from 1.2 to 96 diameters; wall-to-bulk temperature ratios vary from 1 to 1.4. In the low Reynolds number turbulent regime, these data clearly support the present analytical solution rather than the prediction obtained by applying the eddy diffusivity distribution used by Sparrow, Hallman, and Siegel.

Harold C. Reynolds, Jr.↗