Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analytical solution”

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 307 records · Page 17

A defect stream function formulation for compressible turbulent boundary layers

Progress to date on the development of a method for turbulent, wall-bounded flow which uses the defect stream function formulation in the outer layer and an analytic law of the wall and wake formulation in the inner region is reviewed. This two-formulation approach avoids the need to computationally resolve the high-gradient inner layer. One of the most appealing recent developments is the transformation of the compressible governing equation for the defect stream function into a linear, second-order differential equation which has analytic solutions for many problems of practical interest. Numerical and analytic results for incompressible and compressible flows are shown to be in excellent agreement with experimental results. In this paper the two-formulation approach is applied to primitive-variable computations. Excellent comparisons with experiment are presented for two compressible flat plate flows.

Barnwell, Richard W.↗

Electron transport and energy degradation in the ionosphere: Evaluation of the numerical solution, comparison with laboratory experiments and auroral observations

Auroral electron transport calculations are a critical part of auroral models. We evaluate a numerical solution to the transport and energy degradation problem. The numerical solution is verified by reproducing simplified problems to which analytic solutions exist, internal self-consistency tests, comparison with laboratory experiments of electron beams penetrating a collision chamber, and by comparison with auroral observations, particularly the emission ratio of the N2 second positive to N2(+) first negative emissions. Our numerical solutions agree with range measurements in collision chambers. The calculated N(2)2P to N2(+)1N emission ratio is independent of the spectral characteristics of the incident electrons, and agrees with the value observed in aurora. Using different sets of energy loss cross sections and different functions to describe the energy distribution of secondary electrons that emerge from ionization collisions, we discuss the uncertainties of the solutions to the electron transport equation resulting from the uncertainties of these input parameters.

Lummerzheim, D.↗

Slow Crack Growth Analysis of Brittle Materials with Finite Thickness Subjected to Constant Stress-Rate Flexural Loading

A two-dimensional, numerical analysis of slow crack growth (SCG) was performed for brittle materials with finite thickness subjected to constant stress-rate ("dynamic fatigue") loading in flexure. The numerical solution showed that the conventional, simple, one-dimensional analytical solution can be used with a maximum error of about 5% in determining the SCG parameters of a brittle material with the conditions of a normalized thickness (a ratio of specimen thickness to initial crack size) T > 3.3 and of a SCG parameter n > 10. The change in crack shape from semicircular to elliptical configurations was significant particularly at both low stress rate and low T, attributed to predominant difference in stress intensity factor along the crack front. The numerical solution of SCG parameters was supported within the experimental range by the data obtained from constant stress-rate flexural testing for soda-lime glass microslides at ambient temperature.

Chio, S. R.↗

Observations on the Computation of Eigenvalue and Eigenvector Jacobians

Many scientific and engineering problems benefit from analytic expressions for eigenvalue and eigenvector derivatives with respect to the elements of the parent matrix. While there exists extensive literature on the calculation of these derivatives, which take the form of Jacobian matrices, there are a variety of deficiencies that have yet to be addressed — including the need for both left and right eigenvectors, limitations on the matrix structure, and issues with complex eigenvalues and eigenvectors. This work addresses these deficiencies by proposing a new analytic solution for the eigenvalue and eigenvector derivatives. The resulting analytic Jacobian matrices are numerically efficient to compute and are valid for the general complex case. It is further shown that this new general result collapses to previously known relations for the special cases of real symmetric matrices and real diagonal matrices. Finally, the new Jacobian expressions are validated using forward finite differencing and performance is compared with another technique.

Jacobian↗

The general solution to the classical problem of the finite Euler-Bernoulli beam

An analytical solution is obtained for the problem of free and forced vibrations of a finite Euler-Bernoulli beam with arbitrary (partially fixed) boundary conditions. The effects of linear viscous damping, Winkler foundation, constant axial tension, a concentrated mass, and an arbitrary forcing function are included in the analysis. No restriction is placed on the values of the parameters involved, and the solution presented here contains all cited previous solutions as special cases.

Hussaini, M. Y.↗

The Effect of Rotating a Faraday Disc Perpendicular to an Applied Magnetic Field Theory and Experiment

A magnetohydrodynamic model that examines the effect of rotating an electrically conducting cylinder with a uniform external magnetic field applied orthogonal to its axis is presented. Noting a simple geometry, it can be classified as a fundamental dynamo problem. For the case of an infinitely long cylinder, an analytical solution is obtained and analyzed in detail. A semi-analytical model was developed that considers a finite cylinder. Experimental data from a spinning brass wheel in the presence of Earth's magnetic field were compared to the proposed theory and found to fit well.

Mazuruk, Konstantin↗

Combined Numerical/Analytical Perturbation Solutions of the Navier-Stokes Equations for Aerodynamic Ejector/Mixer Nozzle Flows

In spite of rapid advances in both scalar and parallel computational tools, the large number of variables involved in both design and inverse problems make the use of sophisticated fluid flow models impractical, With this restriction, it is concluded that an important family of methods for mathematical/computational development are reduced or approximate fluid flow models. In this study a combined perturbation/numerical modeling methodology is developed which provides a rigorously derived family of solutions. The mathematical model is computationally more efficient than classical boundary layer but provides important two-dimensional information not available using quasi-1-d approaches. An additional strength of the current methodology is its ability to locally predict static pressure fields in a manner analogous to more sophisticated parabolized Navier Stokes (PNS) formulations. To resolve singular behavior, the model utilizes classical analytical solution techniques. Hence, analytical methods have been combined with efficient numerical methods to yield an efficient hybrid fluid flow model. In particular, the main objective of this research has been to develop a system of analytical and numerical ejector/mixer nozzle models, which require minimal empirical input. A computer code, DREA Differential Reduced Ejector/mixer Analysis has been developed with the ability to run sufficiently fast so that it may be used either as a subroutine or called by an design optimization routine. Models are of direct use to the High Speed Civil Transport Program (a joint government/industry project seeking to develop an economically.viable U.S. commercial supersonic transport vehicle) and are currently being adopted by both NASA and industry. Experimental validation of these models is provided by comparison to results obtained from open literature and Limited Exclusive Right Distribution (LERD) sources, as well as dedicated experiments performed at Texas A&M. These experiments have been performed using a hydraulic/gas flow analog. Results of comparisons of DREA computations with experimental data, which include entrainment, thrust, and local profile information, are overall good. Computational time studies indicate that DREA provides considerably more information at a lower computational cost than contemporary ejector nozzle design models. Finally. physical limitations of the method, deviations from experimental data, potential improvements and alternative formulations are described. This report represents closure to the NASA Graduate Researchers Program. Versions of the DREA code and a user's guide may be obtained from the NASA Lewis Research Center.

DeChant, Lawrence Justin↗

Asymptotic solution for solar electric low thrust orbit raising with eclipse penalty

This paper describes an analytic solution for the solar electric low thrust orbit raising problem derived by application of the two-variable asymptotic expansion method. Included in the formulation are the effects of (1) high initial eccentricity, (2) large change of inclination, (3) variable mass, and (4) earth shadow penalty including the rotation of the shadow axis relative to the osculating orbit. An advantage of this mathematical strategy is that the effects of various steering policies can be evaluated readily from a physical point of view.

Flandro, G. A.↗

Simplified solution of diffraction from a Lyot system

This paper presents a derivation of a simplified analytical solution of diffraction from a Lyot (1939) system designed for observation of the solar corona outside of eclipses. Applying the theorem of Papoulis (1986) to simplify the calculations, a simplified solution is derived which is found to be in reasonable agreement with the exact solution. The simplified solution suffices for the preliminary evaluation of the amount of diffraction reduction needed to meet certain system requirements and also serves as a guideline for further apodization.

Wang, Yaujen↗

On the statistical uncertainties associated with line profile fitting

The paper considers the statistical uncertainties in the best-fit parameters associated with the fitting of observational data by model line profiles of specific form. The parameters involved are amplitude, center, and width. A general procedure is described for any profile type of parametric form. The case of Gaussian profiles is treated in detail, and numerical solutions are presented for equally spaced data of equal weight. Approximate analytical solutions are derived that give good general agreement with the numerical solutions over the whole range of interest. Approximate solutions are also obtained for the case of Lorentzian profiles. It is shown that the Lorentzian results can be used in conjunction with the Gaussian results to estimate the corresponding parameter uncertainties in the case of the general Voigt profile.

Landman, D. A.↗

Numerical solutions for heat flow in adhesive lap joints

The present formulation for the modeling of heat transfer in thin, adhesively bonded lap joints precludes difficulties associated with large aspect ratio grids required by standard FEM formulations. This quasi-static formulation also reduces the problem dimensionality (by one), thereby minimizing computational requirements. The solutions obtained are found to be in good agreement with both analytical solutions and solutions from standard FEM programs. The approach is noted to yield a more accurate representation of heat-flux changes between layers due to a disbond.

Howell, P. A.↗

Analysis of gas absorption to a thin liquid film in the presence of a zero-order chemical reaction

The paper presents a detailed theoretical analysis of the process of gas absorption to a thin liquid film adjacent to a horizontal rotating disk. The film is formed by the impingement of a controlled liquid jet at the center of the disk and subsequent radial spreading of liquid along the disk. The chemical reaction between the gas and the liquid film can be expressed as a zero-order homogeneous reaction. The process was modeled by establishing equations for the conservation of mass, momentum, and species concentration and solving them analytically. A scaling analysis was used to determine dominant transport processes. Appropriate boundary conditions were used to solve these equations to develop expressions for the local concentration of gas across the thickness of the film and distributions of film height, bulk concentration, and Sherwood number along the radius of the disk. The partial differential equation for species concentration was solved using the separation of variables technique along with the Duhamel's theorem and the final analytical solution was expressed using confluent hypergeometric functions. Tables for eigenvalues and eigenfunctions are presented for a number of reaction rate constants. A parametric study was performed using Reynolds number, Ekman number, and dimensionless reaction rate as parameters. At all radial locations, Sherwood number increased with Reynolds number (flow rate) as well as Ekman number (rate of rotation). The enhancement of mass transfer due to chemical reaction was found to be small when compared to the case of no reaction (pure absorption), but the enhancement factor was very significant when compared to pure absorption in a stagnant liquid film. The zero-order reaction processes considered in the present investigation included the absorption of oxygen in aqueous alkaline solutions of sodiumdithionite and rhodium complex catalyzed carbonylation of methanol. Present analytical results were compared to previous theoretical results for limiting conditions, and were found to have very good agreement.

Rajagopalan, S.↗

Continuum calculations of continental deformation in transcurrent environments

A thin viscous sheet approximation is used to investigate continental deformation near a strike-slip boundary. The vertically averaged velocity field is calculated for a medium characterized by a power law rheology with stress exponent n. Driving stresses include those applied along boundaries of the sheet and those arising from buoyancy forces related to lateral differences in crustal thickness. Exact and approximate analytic solutions for a region with a sinusoidal strike-slip boundary condition are compared with solutions for more geologically relevant boundary conditions obtained using a finite element technique. The across-strike length scale of the deformation is approximately 1/4pi x sq rt n times the dominant wavelength of the imposed strike-slip boundary condition for both the analytic and the numerical solutions; this result is consistent with length scales observed in continental regions of large-scale transcurrent faulting. An approximate, linear relationship between displacement and rotation is found that depends only on the deformation length scale and the rheology. Calculated displacements, finite rotations, and distribution of crustal thicknesses are consistent with those observed in the region of the Pacific-North America plate boundary in California.

Sonder, L. J.↗

Optimal rendezvous in the neighborhood of a circular orbit

The minimum velocity change rendezvous solutions, when the motion may be linearized about a circular orbit, fall into two separate regions; the phase-for-free region and the general region. Phase-for-free solutions are derived from the optimum transfer solutions, require the same velocity change expenditure, but may not be unique. Analytic solutions are presented in two of the three subregions. An algorithm is presented for determining the unique solutions in the general region. Various sources of initial conditions are discussed and three examples presented.

Jones, J. B.↗