Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical techniques”

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 109 records · Page 6

Undergraduate research report

Studies are reported of the probability distribution of yield strength for mild steel, and the application of numerical techniques to solve mechanical vibration problems. Tensile experiments were conducted for 50 cold rolled and 38 hot rolled mild steel rods. Chi-square, Kolmogro-Smirnov and goodness-of-fit methods are used to test the normality or lognormality of the experimental results. As an example of the application of numerical techniques, a nonlinear dash-pot-speing system subjected to seismic excitation is analyzed by using the fourth-order Runge-Kutta integrated method.

Yeh, H. Y.↗

Numerical marching techniques for fluid flows with heat transfer

The finite difference formulation and method of solution is presented for a wide variety of fluid flow problems with associated heat transfer. Only a few direct results from these formulations are given as examples, since the book is intended primarily to serve a discussion of the techniques and as a starting point for further investigations; however, the formulations are sufficiently complete that a workable computer program may be written from them. In the appendixes a number of topics are discussed which are of interest with respect to the finite difference equations presented. These include a very rapid method for solving certain sets of linear algebraic equations, a discussion of numerical stability, the inherent error in flow rate for confined flow problems, and a method for obtaining high accuracy with a relatively small number of mesh points.

Hornbeck, R. W.↗

Modeling of the float zone process

The interaction of heat, mass, and momentum transport in the floating zone method for growing single crystals from the melt is examined. Methods for detailed numerical simulation of the transport phenomena in a floating zone are developed. Results of the calculations are combined with experiments to determine the effects of solidification induced, surface tension driven, and buoyancy driven convection in establishing dopant redistribution in the melt and the roles of heat transfer in crystal and melt and melt/solid interface shape in determining crystal quality. State of the art finite element techniques were developed for calculating the influence of natural convection in the melt on the shape of a melt/crystal interface and dopant segregation in the crystal. These techniques are demonstrated for solidification by the Bridgman technique. Numerical techniques are developed that calculate the shapes of both the melt/solid and melt/gas interfaces simultaneously with the thermal fields in melt and solid. Models for the fluid flows due to the rotation of the feed and crystal rods are completed and the effects of these flows on dopant segregation are studied, especially in the case of zones longer than can be achieved on Earth.

Brown, R. A.↗

Analytical approximation of a distorted reflector surface defined by a discrete set of points

Reflector antennas on Earth orbiting spacecrafts generally cannot be described analytically. The reflector surface is subjected to a large temperature fluctuation and gradients, and is thus warped from its true geometrical shape. Aside from distortion by thermal stresses, reflector surfaces are often purposely shaped to minimize phase aberrations and scanning losses. To analyze distorted reflector antennas defined by discrete surface points, a numerical technique must be applied to compute an interpolatory surface passing through a grid of discrete points. In this paper, the distorted reflector surface points are approximated by two analytical components: an undistorted surface component and a surface error component. The undistorted surface component is a best fit paraboloid polynomial for the given set of points and the surface error component is a Fourier series expansion of the deviation of the actual surface points, from the best fit paraboloid. By applying the numerical technique to approximate the surface normals of the distorted reflector surface, the induced surface current can be obtained using physical optics technique. These surface currents are integrated to find the far field radiation pattern.

Acosta, Roberto J.↗

CAVE: A computer code for two-dimensional transient heating analysis of conceptual thermal protection systems for hypersonic vehicles

A digital computer code CAVE (Conduction Analysis Via Eigenvalues), which finds application in the analysis of two dimensional transient heating of hypersonic vehicles is described. The CAVE is written in FORTRAN 4 and is operational on both IBM 360-67 and CDC 6600 computers. The method of solution is a hybrid analytical numerical technique that is inherently stable permitting large time steps even with the best of conductors having the finest of mesh size. The aerodynamic heating boundary conditions are calculated by the code based on the input flight trajectory or can optionally be calculated external to the code and then entered as input data. The code computes the network conduction and convection links, as well as capacitance values, given basic geometrical and mesh sizes, for four generations (leading edges, cooled panels, X-24C structure and slabs). Input and output formats are presented and explained. Sample problems are included. A brief summary of the hybrid analytical-numerical technique, which utilizes eigenvalues (thermal frequencies) and eigenvectors (thermal mode vectors) is given along with aerodynamic heating equations that have been incorporated in the code and flow charts.

Rathjen, K. A.↗

Recommendations for improving rigor and reproducibility in site specific characterization

Heterogeneous catalysis is driven by the interaction of reactant molecules and the catalyst surface. The locus of this interaction as well as the surrounding ensemble of atoms is referred to as the catalyst active site. Active site characterization attempts to distinguish active catalytic sites from inactive surface sites, to elucidate the structural and chemical nature of active sites, and to quantify active site concentration. Numerous techniques have been demonstrated to provide compositional and structural information about the active sites within a catalyst. However, each technique has its own limitations and experimental pitfalls that can lead to data misinterpretation or irreproducible results. Further, this work aims to provide an overview of the types of data that can be collected, to outline common experimental challenges and how to avoid them, and to assemble relevant references for the most used active site characterization techniques. More broadly, we aim to outline best practices for researchers to collect, interpret, and report active site characterization data in a way that provides the most benefit to the broader catalysis community. Increasing the rigor and reproducibility of active site characterization offers a strategy to better link properties with catalytic performance and to enable the community to develop consensus concerning these relationships.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Calculation of free fall trajectories based on numerical optimization techniques

The development of a means of computing free-fall (nonthrusting) trajectories from one specified point in the solar system to another specified point in the solar system in a given amount of time was studied. The problem is that of solving a two-point boundary value problem for which the initial slope is unknown. Two standard methods of attack exist for solving two-point boundary value problems. The first method is known as the initial value or shooting method. The second method of attack for two-point boundary value problems is to approximate the nonlinear differential equations by an appropriate linearized set. Parts of both boundary value problem solution techniques described above are used. A complete velocity history is guessed such that the corresponding position history satisfies the given boundary conditions at the appropriate times. An iterative procedure is then followed until the last guessed velocity history and the velocity history obtained from integrating the acceleration history agree to some specified tolerance everywhere along the trajectory.

Source record↗

A mesh gradient technique for numerical optimization

A class of successive-improvement optimization methods in which directions of descent are defined in the state space along each trial trajectory are considered. The given problem is first decomposed into two discrete levels by imposing mesh points. Level 1 consists of running optimal subarcs between each successive pair of mesh points. For normal systems, these optimal two-point boundary value problems can be solved by following a routine prescription if the mesh spacing is sufficiently close. A spacing criterion is given. Under appropriate conditions, the criterion value depends only on the coordinates of the mesh points, and its gradient with respect to those coordinates may be defined by interpreting the adjoint variables as partial derivatives of the criterion value function. In level 2, the gradient data is used to generate improvement steps or search directions in the state space which satisfy the boundary values and constraints of the given problem.

Willis, E. A., Jr.↗

Comparison of numerical integration techniques for orbital applications

The present work gives a brief comparison of the performance of programs for integrating differential equations for orbital applications. The evaluation criteria and the method of testing are described, and the results of the test problem set are included. Integration methods that were chosen for comparison include high-order Runge-Kutta methods; a rational extrapolation method (Bulirsch and Stoer, 1966); a variable step, variable order, multistep method (Krogh, 1969); classical multistep methods of Adams and Cowell; and modified multistep methods. The high-order Runge-Kutta methods used in the comparison include RKF 7(8) and RKF 8(9) (Fehlberg, 1968), and RKS 8-10 (Shanks, 1966).

Moore, H.↗

Light beams and radiation patterns in random media

The purpose of this paper is to apply a numerical technique in solving problems involving light beams in a random medium. The technique starts by generating numerically sample media with prescribed statistical properties. Rays are then traced in these sample media, and the ray statistics are compiled. These statistics are utilized to find the mean square displacement and distribution of a beam. Problems on beam broadening and distortion of radiation patterns are considered. Fluctuations in wave amplitude and phase are also investigated. When possible, the numerical results are compared with the analytical results and experimental results.

Youakim, M. Y.↗

Numerical optimization techniques for bound circulation distribution for minimum induced drag of nonplanar wings: Basic formulations

A theoretical method was developed for determining the optimum span load distribution for minimum induced drag for subsonic nonplanar configurations. The undistorted wing wake is assumed to have piecewise linear variation of shed vortex sheet strength, resulting in a quadratic variation of bound circulation and span load. The optimum loading is obtained either through a direct technique, whereby derivatives of the drag expression are calculated analytically in terms of the unknown wake vortex sheet strengths. Both techniques agree well with each other and with available exact solutions for minimum induced drag.

Kuhlman, J.↗

Numerical optimization techniques for bound circulation distribution for minimum induced drag of Nonplanar wings: Computer program documentation

A two dimensional advanced panel far-field potential flow model of the undistorted, interacting wakes of multiple lifting surfaces was developed which allows the determination of the spanwise bound circulation distribution required for minimum induced drag. This model was implemented in a FORTRAN computer program, the use of which is documented in this report. The nonplanar wakes are broken up into variable sized, flat panels, as chosen by the user. The wake vortex sheet strength is assumed to vary linearly over each of these panels, resulting in a quadratic variation of bound circulation. Panels are infinite in the streamwise direction. The theory is briefly summarized herein; sample results are given for multiple, nonplanar, lifting surfaces, and the use of the computer program is detailed in the appendixes.

Kuhlman, J. M.↗

Numerical evaluation of the performance of active noise control systems

This paper presents a generalized numerical technique for evaluating the optimal performance of active noise controllers. In this technique, the indirect BEM numerical procedures are used to derive the active noise controllers for optimal control of enclosed harmonic sound fields where the strength of the noise sources or the description of the enclosure boundary may not be known. The performance prediction for a single-input single-output system is presented, together with the analysis of the stability and observability of an active noise-control system employing detectors. The numerical procedures presented can be used for the design of both the physical configuration and the electronic components of the optimal active noise controller.

Mollo, C. G.↗

Accuracy Assessment of Response Surface Approximations for Supersonic Turbine Design

There is a growing trend to employ CFD tools to supply the necessary information for design optimization of fluid dynamics components/systems. Such results are prone to uncertainties due to reasons including discretization. errors, incomplete convergence of computational procedures, and errors associated with physical models such as turbulence closures. Based on this type of information, gradient-based optimization algorithms often suffer from the noisy calculations, which can seriously compromise the outcome. Similar problems arise from the experimental measurements. Global optimization techniques, such as those based on the response surface (RS) concept are becoming popular in part because they can overcome some of these barriers. However, there are also fundamental issues related to such global optimization technique such as RS. For example, in high dimensional design spaces, typically only a small number of function evaluations are available due to computational and experimental costs. On the other hand, complex features of the design variables do not allow one to model the global characteristics of the design space with simple quadratic polynomials. Consequently a main challenge is to reduce the size of the region where we fit the RS, or make it more accurate in the regions where the optimum is likely to reside. Response Surface techniques using either polynomials or and Neural Network (NN) methods offer designers alternatives to conduct design optimization. The RS technique employs statistical and numerical techniques to establish the relationship between design variables and objective/constraint functions, typically using polynomials. In this study, we aim at addressing issues related to the following questions: (1) How to identify outliers associated with a given RS representation and improve the RS model via appropriate treatments? (2) How to focus on selected design data so that RS can give better performance in regions critical to design optimization? (3) How to combine NN and polynomial techniques for improving the accuracy of the RS model?

Papila, Nilay↗

Solar Activity Forecasting for use in Orbit Prediction

Orbital prediction for satellites in low Earth orbit (LEO) or low planetary orbit depends strongly on exospheric densities. Solar activity forecasting is important in orbital prediction, as the solar UV and EUV inflate the upper atmospheric layers of the Earth and planets, forming the exosphere in which satellites orbit. Geomagnetic effects also relate to solar activity. Because of the complex and ephemeral nature of solar activity, with different cycles varying in strength by more than 100%, many different forecasting techniques have been utilized. The methods range from purely numerical techniques (essentially curve fitting) to numerous oddball schemes, as well as a small subset, called 'Precursor techniques.' The situation can be puzzling, owing to the numerous methodologies involved, somewhat akin to the numerous ether theories near the turn of the last century. Nevertheless, the Precursor techniques alone have a physical basis, namely dynamo theory, which provides a physical explanation for why this subset seems to work. I discuss this solar cycle's predictions, as well as the Sun's observed activity. I also discuss the SODA (Solar Dynamo Amplitude) index, which provides the user with the ability to track the Sun's hidden, interior dynamo magnetic fields. As a result, one may then update solar activity predictions continuously, by monitoring the solar magnetic fields as they change throughout the solar cycle. This paper ends by providing a glimpse into what the next solar cycle (#24) portends.

Schatten, Kenneth↗

Anomalous Results from PO Applied to Reflector Antennas: The importance of Near Field Computations

There is no doubt among the designers of reflector antennas that the physical optics (PO) analysis technique is the most popular numerical technique. Powerful computer codes are available for the analysis of single or multi reflector antenna systems. Additionally, ever increasingly demand on the antenna performance necessitates the computation of antenna far field patterns under various situations. For example, in using multi reflector antennas such as, Gregorian or Cassegrain, it may become necessary to determine the total fields including the feed radiation pattern, subreflector scattered pattern and the main reflector scattered pattern. In these situations, the common practice is to sum up various scattered fields and the incident field contributions to obtain the desired total field, it is the purpose of this paper to demonstrate that the typical approach based on the far field pattern of the feed would result into erroneous result and special care must be exercised to obtain the correct result. This will be demonstrated through a detailed investigation of a representative test case.

Rahmat-Samii, Yahya↗

Estimating QSVT angles for matrix inversion with large condition numbers

Quantum Singular Value Transformation (QSVT) is a state-of-the-art, near-optimal quantum algorithm that can be used for matrix inversion. The QSVT circuit is parameterized by a sequence of angles that must be pre-calculated classically, with the number of angles increasing as the matrix condition number grows. Computing QSVT angles for ill-conditioned problems is a numerically challenging task. Here, we propose a numerical technique for estimating QSVT angles for large condition numbers. This technique allows one to avoid expensive numerical computations of QSVT angles and to emulate QSVT circuits for solving ill-conditioned problems.

97 MATHEMATICS AND COMPUTING↗