Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 1,171 records · Page 65

Effect of airfoil (trailing-edge) thickness on the numerical solution of panel methods based on the Dirichlet boundary condition

The practical limit of airfoil thickness ratio for which acceptable engineering results are obtainable with the Dirichlet boundary-condition-based numerical methods is investigated. This is done by studying the effect of thickness on the calculated pressure distribution near the trailing edge and by comparing the aerodynamic coefficients with available exact solutions. The first objective of this study, owing to the wide use of such computational methods, is to demonstrate the numerical symptoms that occur when the body or wing thickness approaches zero and to increase the awareness of potential users of these methods. Additionally, an effort is made to obtain the practical limits of the trailing-edge thickness where such problems will appear in the flow solution, and to propose some possible cures for very thin airfoils or those with cusped trailing edges.

Yon, Steven↗

Spectral solution of inviscid supersonic flows over wedges and axisymmetric cones

A shock-fitted multidomain spectral collocation method is used to solve both steady and unsteady inviscid supersonic flows over bodies. New aspects of the method include two subdomain interface types and a zonal solution procedure to get efficient convergence to steady-state in supersonic regions. Steady-state examples include flow over a sharp cone, a hyperbolic cone and a hyperbolic wedge. For the cone, the exact flow solution is used to show that the method is spectrally accurate. As an example of an unsteady flow, a calculation of the interaction of a free-stream hot ring with the flow over a sharp cone is presented.

Kopriva, David A.↗

On development of a finite dynamic element and solution of associated eigenproblem by a block Lanczos procedure

The paper first presents the details of the development of a new six-noded plane triangular finite dynamic element. A block Lanczos algorithm is developed next for the accurate and efficient solution of the quadratic matrix eigenvalue problem associated with the finite dynamic element formulation. The resulting computer program fully exploits matrix sparsity inherent in such a discretization and proves to be most efficient for the extraction of the usually required first few roots and vectors, including repeated ones. Most importantly, the present eigenproblem solution is shown to be comparable to that of the corresponding finite element analysis, thereby rendering the associated dynamic element method rather attractive owing to superior convergence characteristics of such elements, presented herein.

Gupta, K. K.↗

Modifications to Encke's method for long arc orbit determination solutions

An expanded model is developed that permits the extrapolation of Encke's method for the determination of orbits with long arcs, and the model is used to determine a solution for the Lageos trajectory. Encke's method is reviewed emphasizing the nature of the growth of the Encke ratio and the reliability of extrapolated reference trajectories. The reference-orbit formulation is improved by including parameters that accommodate drag and large-amplitude perodic variations in the orbital elements. The proposed Long Arc Model is expected to provide a maximum Encke ratio that is an order of magnitude more reliable than that given by the secularly precessing ellipse. The computational cost of using the Long Arc Model is shown to compare favorably with that of the true force model, and the long arc solutions are useful for current orbit-determination needs.

Lundberg, J. B.↗

Schroedinger's radial equation - Solution by extrapolation

A high-accuracy numerical method for the solution of a 1D Schroedinger equation that is suitable for a diatomic molecule, obtained by combining a finite-difference method with iterative extrapolation to the limit, is presently shown to have several advantages over more conventional methods. Initial guesses for the term values are obviated, and implementation of the algorithm is straightforward. The method is both less sensitive to round-off error, and faster than conventional methods for equivalent accuracy. These advantages are illustrated through the solution of Schroedinger's equation for a Morse potential function suited for HCl and a numerically derived Rydberg-Klein-Rees potential function for the X 1Sigma(+) state of CO.

Goorvitch, D.↗

A solution-adaptive mesh algorithm for dynamic/static refinement of two and three dimensional grids

An adaptive grid algorithm has been developed in two and three dimensions that can be used dynamically with a solver or as part of a grid refinement process. The algorithm employs a transformation from the Cartesian coordinate system to a general coordinate space, which is defined as a parallelepiped in three dimensions. A weighting function, independent for each coordinate direction, is developed that will provide the desired refinement criteria in regions of high solution gradient. The adaptation is performed in the general coordinate space and the new grid locations are returned to the Cartesian space via a simple, one-step inverse mapping. The algorithm for relocation of the mesh points in the parametric space is based on the center of mass for distributed weights. Dynamic solution-adaptive results are presented for laminar flows in two and three dimensions.

Benson, Rusty A.↗

Solution of nonlinear flow equations for complex aerodynamic shapes

Solution-adaptive CFD codes based on unstructured methods for 3-D complex geometries in subsonic to supersonic regimes were investigated, and the computed solution data were analyzed in conjunction with experimental data obtained from wind tunnel measurements in order to assess and validate the predictability of the code. Specifically, the FELISA code was assessed and improved in cooperation with NASA Langley and Imperial College, Swansea, U.K.

Djomehri, M. Jahed↗

Analytical solutions to constrained hypersonic flight trajectories

The flight trajectory of aerospace vehicles subject to a class of path constraints is considered. The constrained dynamics is shown to be a natural two-time-scale system. Asymptotic analytical solutions are obtained. Problems of trajectory optimization and guidance can be dramatically simplified with these solutions. Applications in trajectory design for an aerospace plane strongly support the theoretical development.

Lu, Ping↗

Methods for predicting properties and tailoring salt solutions for industrial processes

An algorithm developed at Oak Ridge National Laboratory accurately and quickly predicts thermodynamic properties of concentrated aqueous salt solutions. This algorithm is much simpler and much faster than other modeling schemes and is unique because it can predict solution behavior at very high concentrations and under varying conditions. Typical industrial applications of this algorithm would be in manufacture of inorganic chemicals by crystallization, thermal storage, refrigeration and cooling, extraction of metals, emissions controls, etc.

Ally, Moonis R.↗

ERS-1 scatterometer calibration and validation activities at ECMWF. B: From radar backscatter characteristics to wind vector solutions

Calibration and validation activities for the ERS-1 scatterometer were carried out at ECMWF (European Center for Medium range Weather Forecast) complementary to the 'Haltenbanken' field campaign off the coast of Norway. At a Numerical Weather Prediction (NWP) center a wealth of verifying data is available both in time and space. This data is used to redefine the wind retrieval procedure given the instrumental characteristics. It was found that a maximum likelihood estimation procedure to obtain the coefficients of a reformulated sigma deg to wind relationship should use radar measurements in logarithmic rather than physical space, and use winds as the wind components rather than wind speed and direction. Doing this, a much more accurate transfer function than the one currently operated by ESA was derived. Sigma deg measurement space shows no signature of a separation in an upwind solution cone and a downwind solution cone. As such signature was anticipated in ESA's wind direction ambiguity removal algorithm, reconsideration of the procedure is necessary. Despite the fact that revisions have to be made in the process of wind retrieval; a grid potential is shown for scatterometry in meteorology and climatology.

Stoffelen, AD↗

A brief description of a new numerical framework for solving conservation laws: The method of space-time conservation element and solution element

A new numerical method for solving conservation laws is being developed. It differs substantially from the well established methods, i.e., finite difference, finite volume, finite element, and spectral methods, in both concept and methodology. It is much simpler than a typical high resolution method. No flux limiter or any technique related to characteristics is involved. No artificial viscosity or smoothing is introduced, and no moving mesh is used. Yet this method is capable of generating highly accurate shock tube solutions. The slight numerical overshoot and/or oscillations generated can be removed if a simple averaging formula initially used is replaced by a weighted formula. This modification has little effect on other parts of the solution. Because of its simplicity, generalization of this new method for multi-dimensional problems is straightforward.

Chang, Sin-Chung↗

On the accurate long-time solution of the wave equation in exterior domains: Asymptotic expansions and corrected boundary conditions

We consider the solution of scattering problems for the wave equation using approximate boundary conditions at artificial boundaries. These conditions are explicitly viewed as approximations to an exact boundary condition satisfied by the solution on the unbounded domain. We study the short and long term behavior of the error. It is provided that, in two space dimensions, no local in time, constant coefficient boundary operator can lead to accurate results uniformly in time for the class of problems we consider. A variable coefficient operator is developed which attains better accuracy (uniformly in time) than is possible with constant coefficient approximations. The theory is illustrated by numerical examples. We also analyze the proposed boundary conditions using energy methods, leading to asymptotically correct error bounds.

Hagstrom, Thomas↗

Implications of new GALLEX results for the Mikheyev-Smirnov-Wolfenstein solution of the solar neutrino problem

We compare the implications for Be-7 and pp neutrinos of the two Mikheyev-Smirnov-Wolfenstein fits to the new GALLEX solar neutrino measurements. Small-mixing-angle solutions tend to suppress the former as electron neutrinos, but not the latter, and large-angle solutions tend to reduce both by about a factor of two. The consequences for BOREXINO and similar solar neutrino-electron scattering experiments are discussed.

Gelb, James M.↗

Finite element solution for the structural behavior of a scientific balloon

A finite element solution for the structural behavior of a scientific balloon has been obtained using a non-linear finite element code. The pneumatic skin is modelled by shell elements that are given a small artificial bending stiffness to overcome numerical problems yet the membrane solution remains unaffected. Validation of the analysis approach is provided through strain measurements on a small scale balloon that exhibits all essential features of a full scale balloon.

Schur, W. W.↗

Efficient yet accurate solution of the linear transport equation in the presence of internal sources - The exponential-linear-in-depth approximation

The present solutions to the linear transport equation pertain to monoenergetic particles' interaction with a multiple scattering/absorbing layered medium with a general anisotropic internal source term. Attention is given to a novel exponential-linear approximation to the internal source, as a function of scattering depth, which furnishes an at-once efficient and accurate solution to the linear transport equation through its reduction of the spatial mesh size. The great superiority of the proposed method is demonstrated by the numerical results obtained in the illustrative cases of (1) an embedded thermal source and (2) a rapidly varying beam pseudosource.

Kylling, Arve↗

On convexity of H-infinity Riccati solutions

The authors revealed several important eigen properties of the stabilizing solutions of the two H-infinity Riccati equations and their product. Among them, the most prominent one is that the spectral radius of the product of these two Riccati solutions is a continuous, nonincreasing, convex function of gamma in the domain of interest. Based on these properties, quadratically convergent algorithms are developed to compute the optimal H-infinity norm. Two examples are used to illustrate the algorithms.

Li, X. P.↗

Unstructured grid solutions to a wing/pylon/store configuration using VGRID3D/USM3D

The purpose of this paper is to validate an inviscid flow solution package based on a new unstructured grid methodology using experimental data on a wing/pylon/store configuration. The solution package consists of an advancing front unstructured grid generator, VGRID3D, and an efficient Euler equation solver, USM3D. Comparisons of computed data versus experimental data are made for two free-stream Mach numbers at five store locations relative to the wing. Both rigid body aerodynamics and mutual interference effects are explored. A very good agreement is observed between computed and wind tunnel data.

Parikh, Paresh↗

Finite difference solution for transient radiative cooling of a conducting semitransparent square region

Transient solutions were obtained for a square region of heat conducting semitransparent material cooling by thermal radiation. The region is in a vacuum environment, so energy is dissipated only by radiation from within the medium leaving through its boundaries. The effect of heat conduction during the transient is to partially equalize the internal temperature distribution. As the optical thickness of the region is increased, the temperature gradients increase near the boundaries and corners, unless heat conduction is large. The solution procedure must provide accurate temperature distributions in these regions to prevent error in the calculated radiation losses. Two-dimensional numerical Gaussian integration is used to obtain the local radiative source term. A finite difference procedure with variable space and time increments is used to solve the transient energy equation. Variable spacing was used to concentrate grid points in regions with large temperature gradients.

Siegel, R.↗