Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXACT SOLUTIONS”

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 91 records · Page 5

An Exact Dual Adjoint Solution Method for Turbulent Flows on Unstructured Grids

An algorithm for solving the discrete adjoint system based on an unstructured-grid discretization of the Navier-Stokes equations is presented. The method is constructed such that an adjoint solution exactly dual to a direct differentiation approach is recovered at each time step, yielding a convergence rate which is asymptotically equivalent to that of the primal system. The new approach is implemented within a three-dimensional unstructured-grid framework and results are presented for inviscid, laminar, and turbulent flows. Improvements to the baseline solution algorithm, such as line-implicit relaxation and a tight coupling of the turbulence model, are also presented. By storing nearest-neighbor terms in the residual computation, the dual scheme is computationally efficient, while requiring twice the memory of the flow solution. The scheme is expected to have a broad impact on computational problems related to design optimization as well as error estimation and grid adaptation efforts.

Nielsen, Eric J.↗

Feedline dynamic effects on shuttle POGO stability

The transmission parameters for the dynamic characteristics of a feedline were approximated using both power and product series expansions. The feedline transfer functions of a shuttle orbiter feedline configuration were obtained using power and product series approximations of 60th, 120th, 180th, and 240th, order. Bode plots using the above polynomial approximations were obtained and the results compared with the exact solution. The exact solution to the feedline transfer function was obtained by using the transcendental terms appearing in the transmission parameters. The results show that the shuttle orbiter feedline may be modeled adequately by using polynomial approximations for the transcendental functions appearing in the transmission parameters. The power series approach was shown to be preferable to the product series method.

Dimaggio, O. D.↗

Exact analytical solution to a transient conjugate heat-transfer problem

An exact analytical solution is found for laminar, constant-property, slug flow over a thin plate which is also convectively cooled from below. The solution is found by means of two successive Laplace transformations when a transient in the plate and the fluid is initiated by a step change in the fluid inlet temperature. The exact solution yields the transient fluid temperature, surface heat flux, and surface temperature distributions. The results of the exact transient solution for the surface heat flux are compared to the quasi-steady values, and a criterion for the validity of the quasi-steady results is found. Also the effect of the plate coupling parameter on the surface heat flux are investigated.

Sucec, J.↗

An improved solution of the simplified two-phase plume

An exact solution of the two-phase, gas-particle, solid propellant rocket exhaust plume flowfield is presented. The good agreement between the exact solution and the experimental data proves that the accuracy of the present technique is acceptable for, at least, preliminary design calculations. A good comparison between the plume impingement effects predicted by the exact solution and those by the best computer code available has also been found for various propellant formulations and rocket sizes. This shows that the present exact solution does describe the basic characteristics of the two-phase plume. In case of a parametric analysis for obtaining an optimum design, the present exact solution will be a more convenient tool than any one of the numerical techniques especially when many parameters are involved in the optimization study. Several typical results are shown herein and their physical significance is discussed.

Wang, S. Y.↗

Reversed flow above a plate with suction

A minor extension is developed for Kovasznay's (1948) solution, and an exact solution is presented for Navier-Stokes equations that may represent the reversed flow about a flat plate with suction. This exact solution is the result of a search for a solution that may serve as an exact basic flow for a model stability analysis of nonparallel flows involving flow reversal.

Lin, S. P.↗

Hierarchic models for laminated plates

Structural plates and shells are three-dimensional bodies, one dimension of which happens to be much smaller than the other two. Thus, the quality of a plate or shell model must be judged on the basis of how well its exact solution approximates the corresponding three-dimensional problem. Of course, the exact solution depends not only on the choice of the model but also on the topology, material properties, loading and constraints. The desired degree of approximation depends on the analyst's goals in performing the analysis. For these reasons models have to be chosen adaptively. Hierarchic sequences of models make adaptive selection of the model which is best suited for the purposes of a particular analysis possible. The principles governing the formulation of hierarchic models for laminated plates are presented. The essential features of the hierarchic models described models are: (1) the exact solutions corresponding to the hierarchic sequence of models converge to the exact solution of the corresponding problem of elasticity for a fixed laminate thickness; and (2) the exact solution of each model converges to the same limit as the exact solution of the corresponding problem of elasticity with respect to the laminate thickness approaching zero. The formulation is based on one parameter (beta) which characterizes the hierarchic sequence of models, and a set of constants whose influence was assessed by a numerical sensitivity study. The recommended selection of these constants results in the number of fields increasing by three for each increment in the power of beta. Numerical examples analyzed with the proposed sequence of models are included and good correlation with the reference solutions was found. Results were obtained for laminated strips (plates in cylindrical bending) and for square and rectangular plates with uniform loading and with homogeneous boundary conditions. Cross-ply and angle-ply laminates were evaluated and the results compared with those of MSC/PROBE. Hierarchic models make the computation of any engineering data possible to an arbitrary level of precision within the framework of the theory of elasticity.

Actis, Ricardo Luis↗

A new class of exact, nonlinear solutions to the Grad-Shafranov equation

We have constructed a new class of exact, nonlinear solutions to the Grad-Shafranov equation, representing force-free magnetic fields with translational symmetry. These exact solutions are pertinent to the study of magnetic structures in the solar corona that are subjected to photospheric shearing motions.

Roumeliotis, George↗

The Pekeris problem, and the interior of the pulsating sphere, first success of exact ray theory, and what rays?

This paper presents what is probably the first application of exact ray theory to provide nontrivial exact solutions of the Helmholtz equation. First, the Pekeris problem is considered. This is followed by consideration of the problem of the field maintained within a uniformly pulsating sphere. In the Pekeris problem, an iterative solution provides the exact solution on the second iteration. In the second problem, it is shown that the exact ray theory eliminates the focus of the geometric theory.

Mcaninch, Gerry L.↗

Theoretical and lidar studies of the density response of the mesospheric sodium layer to gravity wave perturbations

The density response of atmospheric layers to gravity waves is developed in two forms, an exact solution and a perturbation series solution. The degree of nonlinearity in the layer density response is described by the series solution whereas the exact solution gives insight into the nature of the responses. Density perturbation in an atmospheric layer are shown to be substantially greater than the atmospheric density perturbation associated with the propagation of a gravity wave. Because of the density gradients present in atmospheric layers, interesting effects were observed such as a phase reversal in the linear layer response which occurs near the layer peak. Once the layer response is understood, the sodium layer can be used as a tracer of atmospheric wave motions. A two dimensional digital signal processing technique was developed. Both spatial and temporal filtering are utilized to enhance the resolution by decreasing shot noise by more han 10 dB. Many of the features associated with a layer density response to gravity waves were observed in high resolution density profiles of the mesospheric sodium layer. These include nonlinearities as well as the phase reversal in the linear layer response.

Shelton, J. D.↗

Verification of RadCAD: Specular Capabilities

As part of the RadCAD's development process, it is necessary to compare RadCAD's results with other radiation tools and exact solutions when and where possible. Form factor algorithms have been previously verified with exact solutions. This paper will consider RadCAD's specular capabilities. First, radiation exchange factors will be compared against exact solutions and results from TRASYS for various geometries. Critical dimensions and optical properties are changed for each geometry. Second, a specular adjunct plate system will be used to verify absorbed heat fluxes. This particular geometric problem has had some attention in the literature. Previous authors have used this problem to validate software results with exact analytical solution. This paper will compare absorbed heat rates against the exact solution and other published results from other thermal radiation tools. The agreement between RadCAD and the exact solutions is good. The maximum error for both specular and diffuse exchange factors for both geometries and all optical properties was 3%. The absorbed fluxes differed by a maximum of 4% for the adjunct plate problem.

Lucas, Sam↗

Hierarchic models for laminated composites

The principles governing the formulation of hierarchic models for laminated composites are discussed. The essential features of the hierarchic models described herein are: (1) the exact solutions corresponding to the hierarchic sequence of models converge to the exact solution of the corresponding problem of elasticity for a fixed laminate thickness; and (2) the exact solution of each model converges to the same limit as the exact solution of the corresponding problem of elasticity with respect to the laminate thickness approaching zero. Hierarchic models make the computation of any engineering data possible to an arbitrary level of precision within the framework of the theory of elasticity.

Babuska, Ivo↗

Numerical Uncertainty Analysis for Computational Fluid Dynamics using Student T Distribution -- Application of CFD Uncertainty Analysis Compared to Exact Analytical Solution

Computational Fluid Dynamics (CFD) is the standard numerical tool used by Fluid Dynamists to estimate solutions to many problems in academia, government, and industry. CFD is known to have errors and uncertainties and there is no universally adopted method to estimate such quantities. This paper describes an approach to estimate CFD uncertainties strictly numerically using inputs and the Student-T distribution. The approach is compared to an exact analytical solution of fully developed, laminar flow between infinite, stationary plates. It is shown that treating all CFD input parameters as oscillatory uncertainty terms coupled with the Student-T distribution can encompass the exact solution.

Groves, Curtis E.↗

Exact vacuum solution to conformal Weyl gravity and galactic rotation curves

The complete, exact exterior solution for a static, spherically symmetric source in locally conformal invariant Weyl gravity is presented. The solution includes the familiar exterior Schwarzschild solution as a special case and contains an extra gravitational potential term which grows linearly with distance. The obtained solution provides a potential explanation for observed galactic rotation curves without the need for dark matter. The solution also has some interesting implications for cosmology.

Mannheim, Philip D.↗

Second Computational Aeroacoustics (CAA) Workshop on Benchmark Problems

The proceedings of the Second Computational Aeroacoustics (CAA) Workshop on Benchmark Problems held at Florida State University are the subject of this report. For this workshop, problems arising in typical industrial applications of CAA were chosen. Comparisons between numerical solutions and exact solutions are presented where possible.

Tam, C. K. W.↗

Conjugate quasilinear Dirichlet and Neumann problems and a posteriori error bounds

Quasilinear Dirichlet and Neumann problems on a rectangle D with boundary D prime are considered. Using these concepts, conjugate problems, that is, a pair of one Dirichlet and one Neumann problem, the minima of the energies of which add to zero, are introduced. From the concept of conjugate problems, two-sided bounds for the energy of the exact solution of any given Dirichlet or Neumann problem are constructed. These two-sided bounds for the energy at the exact solution are in turn used to obtain a posteriori error bounds for the norm of the difference of the approximate and exact solutions of the problem. These bounds do not involve the unknown exact solution and are easily constructed numerically.

Lavery, J. E.↗

Nonlinear structure formation in flat cosmological models

This paper describes the formation of nonlinear structure in flat (zero curvature) Friedmann cosmological models. We consider models with two components: the usual nonrelativistic component that evolves under gravity and eventually forms the large-scale structure of the universe, and a uniform dark matter component that does not clump under gravity, and whose energy density varies with the scale factor a(t) like a(t)(sup -n), where n is a free parameter. Each model is characterized by two parameters: the exponent n and the present density parameter Omega(sub 0) of the nonrelativistic component. The linear perturbation equations are derived and solved for these models, for the three different cases n = 3, n is greater than 3, and n is less than 3. The case n = 3 is relevant to model with massive neutrinos. The presence of the uniform component strongly reduces the growth of the perturbation compared with the Einstein-de Sitter model. We show that the Meszaros effect (suppression of growth at high redshift) holds not only for n = 4, radiation-dominated models, but for all models with n is greater than 3. This essentially rules out any such model. For the case n is less than 3, we numerically integrate the perturbation equations from the big bang to the present, for 620 different models with various values of Omega(sub 0) and n. Using these solutions, we show that the function f(Omega(sub 0), n) = (a/delta(sub +))d(delta)(sub +)/da, which enters in the relationship between the present density contrast delta(sub 0) and peculiar velocity field u(sub 0) is essentially independent of n. We derive approximate solutions for the second-order perturbation equations. These second-order solutions are tested against the exact solutions and the Zel'dovich approximation for spherically symmetric perturbations in the marginally nonlinear regime (the absolute value of delta is less than or approximately 1). The second-order and Zel'dovich solutions have comparable accuracy, significantly higher than the accuracy of the linear solutions. We then investigate the dependence of the delta(sub 0) - u(sub 0) relationship upon the value of n in the nonlinear regime using the second-order solutions for marginally nonlinear, general perturbations, and the exact solutions for strongly nonlinear, spherically symmetric perturbations. In both cases, we find that the delta(sub 0) - u(sub 0) relationship remains independent of n. We speculate that this result extends to strongly nonlinear, general perturbations as well. This eliminates any hope to determine the presence of the uniform component or the value of n using dynamical methods. Finally, we compute the nonlinear evolution of the skewness of the distribution of values of delta, assuming Gaussian initial conditions. We find that the skewness is not only independent of n, but also of Omega(sub 0). Thus the skewness cannot be used to discriminate among various models with Gaussian initial conditions. However, it can be used for testing the Gaussianity of the initial conditions themselves. We conclude that the uniform component leaves no observable signature in the present large-scale structure of the universe. To determine its presence and nature, we must investigate the relationship between the past and present universe, using redshift-dependent tests.

Martel, Hugo↗

Direct analytical solutions to non-uniform beam problems

The direct analytical solution to the vibration of non-uniform beams with and without discontinuities and with various boundary conditions is presented. Results are compared to results from the exact solution for certain cases where the exact solution has been obtained. It is shown that the direct solution converges to the exact solution, in fact, with 'indefinite accuracy' just as Hamilton stated that it would.

Bailey, C. D.↗