Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical mathematics”

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 667 records · Page 37

An Automated Method to Compute Orbital Re-entry Trajectories with Heating Constraints

Determining how to properly manipulate the controls of a re-entering re-usable launch vehicle (RLV) so that it is able to safely return to Earth and land involves the solution of a two-point boundary value problem (TPBVP). This problem, which can be quite difficult, is traditionally solved on the ground prior to flight. If necessary, a nearly unlimited amount of time is available to find the 'best' solution using a variety of trajectory design and optimization tools. The role of entry guidance during flight is to follow the pre- determined reference solution while correcting for any errors encountered along the way. This guidance method is both highly reliable and very efficient in terms of onboard computer resources. There is a growing interest in a style of entry guidance that places the responsibility of solving the TPBVP in the actual entry guidance flight software. Here there is very limited computer time. The powerful, but finicky, mathematical tools used by trajectory designers on the ground cannot in general be converted to do the job. Non-convergence or slow convergence can result in disaster. The challenges of designing such an algorithm are numerous and difficult. Yet the payoff (in the form of decreased operational costs and increased safety) can be substantiaL This paper presents an algorithm that incorporates features of both types of guidance strategies. It takes an initial RLV orbital re-entry state and finds a trajectory that will safely transport the vehicle to Earth. During actual flight, the computed trajectory is used as the reference to be flown by a more traditional guidance method.

Zimmerman, Curtis↗

A Dynamic Approach to Monitoring Particle Fallout in a Cleanroom Environment

This slide presentation discusses a mathematical model to monitor particle fallout in a cleanroom. "Cleanliness levels" do not lead to increases with regards to cleanroom type or time because the levels are not linear. Activity level, impacts the cleanroom class. The numerical method presented leads to a simple Class-hour formulation, that allows for dynamic monitoring of the particle using a standard air particle counter.

Perry, Radford L., III↗

The Sizing and Optimization Language, (SOL): Computer language for design problems

The Sizing and Optimization Language, (SOL), a new high level, special purpose computer language was developed to expedite application of numerical optimization to design problems and to make the process less error prone. SOL utilizes the ADS optimization software and provides a clear, concise syntax for describing an optimization problem, the OPTIMIZE description, which closely parallels the mathematical description of the problem. SOL offers language statements which can be used to model a design mathematically, with subroutines or code logic, and with existing FORTRAN routines. In addition, SOL provides error checking and clear output of the optimization results. Because of these language features, SOL is best suited to model and optimize a design concept when the model consits of mathematical expressions written in SOL. For such cases, SOL's unique syntax and error checking can be fully utilized. SOL is presently available for DEC VAX/VMS systems. A SOL package is available which includes the SOL compiler, runtime library routines, and a SOL reference manual.

Lucas, Stephen H.↗

Evaluation of massless-spring modeling of suspension-line elasticity during the parachute unfurling process

A general theory on mathematical modeling of elastic parachute suspension lines during the unfurling process was developed. Massless-spring modeling of suspension-line elasticity was evaluated in detail. For this simple model, equations which govern the motion were developed and numerically integrated. The results were compared with flight test data. In most regions, agreement was satisfactory. However, poor agreement was obtained during periods of rapid fluctuations in line tension.

Poole, L. R.↗

Similarity laws of lunar and terrestrial volcanic flows

A mathematical model of a one dimensional, steady duct flow of a mixture of a gas and small solid particles (rock) was analyzed and applied to the lunar and the terrestrial volcanic flows under geometrically and dynamically similar conditions. Numerical results for the equilibrium two phase flows of lunar and terrestrial volcanoes under similar conditions are presented. The study indicates that: (1) the lunar crater is much larger than the corresponding terrestrial crater; (2) the exit velocity from the lunar volcanic flow may be higher than the lunar escape velocity but the exit velocity of terrestrial volcanic flow is much less than that of the lunar case; and (3) the thermal effects on the lunar volcanic flow are much larger than those of the terrestrial case.

Pai, S. I.↗

Control of the errors of discretization and idealization in finite element analysis

Understanding of the basic principles which control errors of discretization in finite element analysis has increased very substantially since 1980. The main milestones were: (1) development of the theoretical basis of p-extensions (1981); (2) understanding of the proper interplay between mesh design and assignment of polynomial degree to elements. Practical realization of exponential convergence rates, independently of the smoothness of the exact solution (1984); and (3) industrial experience with the new finite element technology known as the p- or hp-version of the finite element method: General Dynamics reported thirty- to forty-fold savings in terms of human time and large savings in computer time (1986). Lockheed reported favorably on their evaluation of error estimation and quality control capabilities of the p-version in industrial settings (1987). The gains in our understanding of how to control the errors of discretization represent only half of the control necessary to ensure that a numerical model is in fact an accurate representation of the corresponding physical system. Control of the errors of idealization is equally important. A brief overview of the main ideas of how to ensure the quality and reliability of mathematical models of structural systems is presented.

Szabo, Barna A.↗

Vibrational relaxation in hypersonic flow fields

Mathematical formulations of vibrational relaxation are derived from first principles for application to fluid dynamic computations of hypersonic flow fields. Relaxation within and immediately behind shock waves is shown to be substantially faster than that described in current numerical codes. The result should be a significant reduction in nonequilibrium radiation overshoot in shock layers and in radiative heating of hypersonic vehicles; these results are precisely the trends needed to bring theoretical predictions more in line with flight data. Errors in existing formulations are identified and qualitative comparisons are made.

Meador, Willard E.↗

Experimental Verification of Computational Models for Laminated Composites

The objective of the research reported herein is to develop a progressive damage methodology capable of predicting the residual strength of continuous fiber-reinforced, laminated, polymer matrix composites with through-penetration damage. The fracture behavior of center-notch tension panels with thin crack-like slits was studied. Since fibers are the major load-carrying constituent in polymer matrix composites, predicting the residual strength of a laminate requires a criterion for fiber fracture. The effects on fiber strain due to other damage mechanisms such as matrix cracking and delaminations must also be modeled. Therefore, the research herein examines the damage mechanisms involved in translaminate fracture and identifies the toughening mechanisms responsible for damage growth resistance in brittle epoxy matrix systems. The mechanics of matrix cracking and fiber fracture are discussed as is the mathematical framework for the progressive damage model developed by the authors. The progressive damage analysis algorithms have been implemented into a general purpose finite element code developed by NASA, the Computational Structural Mechanics Testbed (COMET). Damage growth is numerically simulated and the analytical residual strength predictions are compared to experimental results for a variety of notched panel configurations and materials systems.

Harris, Charles E.↗

Limits to randomness in paleobiologic models: the case of Phanerozoic species diversity

The question of how random, or unconstrained, paleobiologic models should be is examined with a case study: Signor's (1982, 1985) inverse calculation of levels of marine species diversity through the Phanerozoic. His calculation involved an ingenious model that estimated species numbers and species abundances in the world oceans of the past by correcting known numbers of fossil species for variations in sedimentary rocks available for sampling and in effort paleontologists might devote to sampling. The model proves robust to changes in possible shapes of species-abundance distributions, but it is sensitive to alterations in the assumption that paleontologists collect fossils at random. If it is assumed that ease of collecting varies with age of sediment (with the Cenozoic offering easy sampling) or that paleontologists tend to seek out rarer fossils, results of the inverse calculation change. In particular, the magnitude of the calculated Cenozoic diversity increase always declines from the factor of about seven as originally reported to something considerably smaller. This leaves open the problem of the magnitude of Cenozoic increase in marine species diversity, awaiting better empirical data and, perhaps, more exacting models, random or otherwise.

Non-NASA Center↗

Comparison of EKF and UKF for Robotic Missions to Mars

The batched Kalman lter works well for a wide variety of orbit determination scenarios but requires linearized approximations to the underlying dynamical system. Furthermore, implementation depends on the computation of partial derivatives for all measurements and parameters used in the lter, a non-trivial task even when analytical derivations are possible. To address these limitations, many extensions and modi - cations to the linear Kalman lter have been proposed over the years since Kalman's initial publication. Our investigation examines two popular variants of the Kalman lter, namely the Extended Kalman Filter (EKF) and Unscented Kalman Filter (UKF), within the context of spacecraft missions to Mars. The EKF relies on the same mathematical basis as the classic version of the lter, however it sequentially updates the linearization as the ltering process is conducted. In contrast, the UKF is built upon the unscented transform which samples the nearby solution space and numerically propagates perturbed trajectories as part of the update / estimation process. We assess the performance of the ltering algorithms under a variety of mission scenarios with the goal of establishing decision criteria for which implementation to use under di ering circumstances.

Ely, Todd↗

On propagation of long waves in curved ducts.

A two-dimensional detailed study of the behavior of long waves in curved ducts and in junctions between straight and curved ducts will be given. The mathematical treatment of the problem utilizes the method of separation variables. Solutions and expressions for principal mode of the wave are obtained by using the linearized equation of motion solved for its characteristic values. The unavoidable approximations in the numerical solutions of the cylindrical functions are due to use of series expansion of Bessel functions and from restrictions necessary to solve infinite matrices.

Rostafinski, W.↗

Band models and correlations for infrared radiation

Absorption of infrared radiation by various line and band models are briefly reviewed. Narrow band model relations for absorptance are used to develop 'exact' formulations for total absorption by four wide band models. Application of a wide band model to a particular gas largely depends upon the spectroscopic characteristic of the absorbing-emitting molecule. Seven continuous correlations for the absorption of a wide band model are presented and each one of these is compared with the exact (numerical) solutions of the wide band models. Comparison of these results indicate the validity of a correlation for a particular radiative transfer application. In radiative transfer analyses, use of continuous correlations for total band absorptance provides flexibilities in various mathematical operations.

Tiwari, S. N.↗

Similar explosive eruptions of lunar and terrestrial volcanoes

A mathematical model of one-dimensional, steady duct flow of a mixture of gas and small solid particles (rock), has been applied to the earth and the moon under geometrically and dynamically similar conditions to simulate explosive eruption conditions. Numerical results for equilibrium two-phase flows for lunar and terrestrial explosive eruptions under similar conditions indicate that: (1) a lunar explosive vent is much larger than the corresponding terrestrial vent, (2) the exit velocity from the lunar explosive flow may be higher than the lunar escape velocity, but the exit velocity of terrestrial explosive flow is much less than that of the lunar case (this result supports the hypothesis that Australian tektites came from the moon as a stream of a mixture of rock and gas of extremely high speed), and (3) the thermal effects on the lunar explosive flows are much larger than those of the terrestrial case.

Pai, S. I.↗

Investigation of hypersonic rarefied flow on a spherical nose of the AOTV

The Navier-Stokes (NS) equations were integrated numerically for investigating the flow characteristics on the forepart of the spherical nose of a space vehicle such as the AOTV or AFE by a modified Accelerated Successive Replacement (ASR) scheme under hypersonic rarefied conditions. Technical feasibility of the mathematical approach was demonstrated by computing the flowfield on a spherical nose under conditions that the AFE encounters at times t = 15 and 20 seconds after its reentry into the atmosphere. Local similar solutions for the merged layer flow along the stagnation line of the sphere were developed. These are correct to the same degree of accuracy as the NS equations. These solutions provided stagnation line boundary conditions for the domain of integration on the spherical noise. Also, a parametric study of the stagnation line solution was made with a view to understand the flow characteristics in tunnels with different ambient fluids. Analytical expressions for surface slip temperature, jump conditions, and concentration level in the presence of the real gas effects at the top of the Knudsen layer were derived and used to calculate the stagnation line flowfield with nonequilibrium dissociation and ionization. A number of graphics were drawn to illustrate the basic physics of the flowfields. The present analysis can be extended to include real gas effects and to bodies of arbitrary shapes. It can further provide boundary conditions for integrating the NS equations in the near wake region.

Jain, Amolak C.↗

Mathematical model for the Bridgman-Stockbarger crystal growing system

In a major technical breakthrough, a computer model for Bridgman-Stockbarger crystal growth was developed. The model includes melt convection, solute effects, thermal conduction in the ampule, melt, and crystal, and the determination of the curved moving crystal-melt interface. The key to the numerical method is the use of a nonuniform computational mesh which moves with the interface, so that the interface is a mesh surface. In addition, implicit methods are used for advection and diffusion of heat, concentration, and vorticity, for interface movement, and for internal gracity waves. This allows large time-steps without loss of stability or accuracy. Numerical results are presented for the interface shape, temperature distribution, and concentration distribution, in steady-state crystl growth. Solutions are presented for two test cases using water, with two different salts in solution. The two diffusivities differ by a factor of ten, and the concentrations differ by a factor of twenty.

Roberts, G. O.↗

Structured adaptive grid generation using algebraic methods

The accuracy of the numerical algorithm depends not only on the formal order of approximation but also on the distribution of grid points in the computational domain. Grid adaptation is a procedure which allows optimal grid redistribution as the solution progresses. It offers the prospect of accurate flow field simulations without the use of an excessively timely, computationally expensive, grid. Grid adaptive schemes are divided into two basic categories: differential and algebraic. The differential method is based on a variational approach where a function which contains a measure of grid smoothness, orthogonality and volume variation is minimized by using a variational principle. This approach provided a solid mathematical basis for the adaptive method, but the Euler-Lagrange equations must be solved in addition to the original governing equations. On the other hand, the algebraic method requires much less computational effort, but the grid may not be smooth. The algebraic techniques are based on devising an algorithm where the grid movement is governed by estimates of the local error in the numerical solution. This is achieved by requiring the points in the large error regions to attract other points and points in the low error region to repel other points. The development of a fast, efficient, and robust algebraic adaptive algorithm for structured flow simulation applications is presented. This development is accomplished in a three step process. The first step is to define an adaptive weighting mesh (distribution mesh) on the basis of the equidistribution law applied to the flow field solution. The second, and probably the most crucial step, is to redistribute grid points in the computational domain according to the aforementioned weighting mesh. The third and the last step is to reevaluate the flow property by an appropriate search/interpolate scheme at the new grid locations. The adaptive weighting mesh provides the information on the desired concentration of points to the grid redistribution scheme. The evaluation of the weighting mesh is accomplished by utilizing the weight function representing the solution variation and the equidistribution law. The selection of the weight function plays a key role in grid adaptation. A new weight function utilizing a properly weighted boolean sum of various flowfield characteristics is defined. The redistribution scheme is developed utilizing Non-Uniform Rational B-Splines (NURBS) representation. The application of NURBS representation results in a well distributed smooth grid by maintaining the fidelity of the geometry associated with boundary curves. Several algebraic methods are applied to smooth and/or nearly orthogonalize the grid lines. An elliptic solver is utilized to smooth the grid lines if there are grid crossings. Various computational examples of practical interest are presented to demonstrate the success of these methods.

Yang, Jiann-Cherng↗

Collision frequency of artificial satellites - The creation of a debris belt

As the number of artificial satellites in earth orbit increases, the probability of collisions between satellites also increases. Satellite collisions would produce orbiting fragments, each of which would increase the probability of further collisions, leading to the growth of a belt of debris around the earth. This process parallels certain theories concerning the growth of the asteroid belt. The debris flux in such an earth-orbiting belt could exceed the natural meteoroid flux, affecting future spacecraft designs. A mathematical model was used to predict the rate at which such a belt might form. Under certain conditions the belt could begin to form within this century and could be a significant problem during the next century. The possibility that numerous unobserved fragments already exist from spacecraft explosions would decrease this time interval. However, early implementation of specialized launch constraints and operational procedures could significantly delay the formation of the belt.

Donald J Kessler↗