Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mathematics and Computing, Physics”

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

Initial Ada components evaluation

The SAIC has the responsibility for independent test and validation of the SSE. They have been using a mathematical functions library package implemented in Ada to test the SSE IV and V process. The library package consists of elementary mathematical functions and is both machine and accuracy independent. The SSE Ada components evaluation includes code complexity metrics based on Halstead's software science metrics and McCabe's measure of cyclomatic complexity. Halstead's metrics are based on the number of operators and operands on a logical unit of code and are compiled from the number of distinct operators, distinct operands, and total number of occurrences of operators and operands. These metrics give an indication of the physical size of a program in terms of operators and operands and are used diagnostically to point to potential problems. McCabe's Cyclomatic Complexity Metrics (CCM) are compiled from flow charts transformed to equivalent directed graphs. The CCM is a measure of the total number of linearly independent paths through the code's control structure. These metrics were computed for the Ada mathematical functions library using Software Automated Verification and Validation (SAVVAS), the SSE IV and V tool. A table with selected results was shown, indicating that most of these routines are of good quality. Thresholds for the Halstead measures indicate poor quality if the length metric exceeds 260 or difficulty is greater than 190. The McCabe CCM indicated a high quality of software products.

Moebes, Travis↗

Active optics simulation system

The active optics simulation system (AOSS) is a set of computer programs and associated software to be used in the development, design, and evaluation of a primary mirror control system for a large space telescope, (e.g., the tentatively proposed 3-meter telescope). The mathematical models of component subsystems and the solution of the physical processes that occur within the mirror surface control system were obtained, and based on these models AOSS simulates the behavior of the entire mirror surface control system as well as the behavior of the component subsystems. The program has a modular structure so that any subsystem module can be replaced or modified with minimum disruption of the rest of the simulation program.

Chi, C. H.↗

GENOA-PFA: Progressive Fracture in Composites Simulated Computationally

GENOA-PFA is a commercial version of the Composite Durability Structural Analysis (CODSTRAN) computer program that simulates the progression of damage ultimately leading to fracture in polymer-matrix-composite (PMC) material structures under various loading and environmental conditions. GENOA-PFA offers several capabilities not available in other programs developed for this purpose, making it preferable for use in analyzing the durability and damage tolerance of complex PMC structures in which the fiber reinforcements occur in two- and three-dimensional weaves and braids. GENOA-PFA implements a progressive-fracture methodology based on the idea that a structure fails when flaws that may initially be small (even microscopic) grow and/or coalesce to a critical dimension where the structure no longer has an adequate safety margin to avoid catastrophic global fracture. Damage is considered to progress through five stages: (1) initiation, (2) growth, (3) accumulation (coalescence of propagating flaws), (4) stable propagation (up to the critical dimension), and (5) unstable or very rapid propagation (beyond the critical dimension) to catastrophic failure. The computational simulation of progressive failure involves formal procedures for identifying the five different stages of damage and for relating the amount of damage at each stage to the overall behavior of the deteriorating structure. In GENOA-PFA, mathematical modeling of the composite physical behavior involves an integration of simulations at multiple, hierarchical scales ranging from the macroscopic (lamina, laminate, and structure) to the microscopic (fiber, matrix, and fiber/matrix interface), as shown in the figure. The code includes algorithms to simulate the progression of damage from various source defects, including (1) through-the-thickness cracks and (2) voids with edge, pocket, internal, or mixed-mode delaminations.

Murthy, Pappu L. N.↗

Multihierarchy Gaussian Process Models for Probabilistic Aerodynamic Databases using Uncertain Nominal and Off-Nominal Configuration Data

Probabilistic aerodynamic databases are a crucial component of the development lifecycle for aerospace vehicles. A key challenge when building aerodynamic databases is that most data used to construct them represent various simplifications of the real flight vehicle. For example, wind tunnel models often simplify the vehicle geometry and surface roughness characteristics, while CFD computations often make simplifications to the physics being modeled, such as fully laminar or turbulent calculations. Multifidelity data fusion models rely on a user being able to define a hierarchy of fidelity levels anchored to some "truth" data. This approach is unsatisfactory when no data can be considered to accurately reflect real flight conditions. In this work, we provide an alternative approach by presenting a consistent mathematical framework for building probabilistic aerodynamic databases in the form of a conditional probability distribution described by an ensemble of multifidelity Gaussian Processes. Instead of relying on a single hierarchy of data fidelity levels, the presented framework identifies a "nominal" configuration and potential corrections to the nominal which represent specific physical phenomena not represented in the nominal data. The nominal and correction functions themselves are constructed as multifidelity Gaussian Processes and linearly combined to form an ensemble model which fuses the uncertainties associated nominal and correction models. Results obtained using the proposed framework on a simplified Orion Crew Module wind tunnel dataset demonstrate the predictive capability of the multihierarchy framework. We further demonstrate the benefits of such a probabilistic aerodynamic database approach through function sampling and computing the conditional distributions of derived quantities, such as the trim angle of attack and aerodynamic coefficients at trim.

Gaussian Processes↗

An efficient code for the simulation of nonhydrostatic stratified flow over obstacles

The physical model and computational procedure of the code is described in detail. The code is validated in tests against a variety of known analytical solutions from the literature and is also compared against actual mountain wave observations. The code will receive as initial input either mathematically idealized or discrete observational data. The form of the obstacle or mountain is arbitrary.

Pihos, G. G.↗

A magnetic hysteresis model

The Passive Aerodynamically Stabilized Magnetically Damped Satellite (PAMS) will be deployed from the Space Shuttle and used as a target for a Shuttle-mounted laser. It will be a cylindrical satellite with several corner cube reflectors on the ends. The center of mass of the cylinder will be near one end, and aerodynamic torques will tend to align the axis of the cylinder with the spacecraft velocity vector. Magnetic hysteresis rods will be used to provide passive despin and oscillation-damping torques on the cylinder. The behavior of the hysteresis rods depends critically on the 'B/H' curves for the combination of materials and rod length-to-diameter ratio ('l-over-d'). These curves are qualitatively described in most Physics textbooks in terms of major and minor 'hysteresis loops'. Mathematical modeling of the functional relationship between B and H is very difficult. In this paper, the physics involved is not addressed, but an algorithm is developed which provides a close approximation to empirically determined data with a few simple equations suitable for use in computer simulations.

Flatley, Thomas W.↗

Buoyancy-Driven Heat Transfer During Application of a Thermal Gradient for the Study of Vapor Deposition at Low Pressure Using and Ideal Gas

A mathematical model has been developed to determine heat transfer during vapor deposition of source materials under a variety of orientations relative to gravitational accelerations. The model demonstrates that convection can occur at total pressures as low as 10-2 mm Hg. Through numerical computation, using physical material parameters of air, a series of time steps demonstrates the development of flow and temperature profiles during the course of vapor deposition. These computations show that in unit gravity vapor deposition occurs by transport through a fairly complicated circulating flow pattern when applying heat to the bottom of the vessel with parallel orientation with respect to the gravity vector. The model material parameters for air predict the effect of kinematic viscosity to be of the same order as thermal diffusivity, which is the case for Prandtl number approx. 1 fluids. Qualitative agreement between experiment and the model indicates that 6-(2-methyl-4-nitroanilino)-2,4-hexadiyn-l-ol (DAMNA) at these pressures indeed approximates an ideal gas at the experiment temperatures, and may validate the use of air physical constants. It is apparent that complicated nonuniform temperature distribution in the vapor could dramatically affect the homogeneity, orientation, and quality of deposited films. The experimental test i's a qualitative comparison of film thickness using ultraviolet-visible spectroscopy on films generated in appropriately oriented vapor deposition cells. In the case where heating of the reaction vessel occurs from the top, deposition of vapor does not normally occur by convection due to a stable stratified medium. When vapor deposition occurs in vessels heated at the bottom, but oriented relative to the gravity vector between these two extremes, horizontal thermal gradients induce a complex flow pattern. In the plane parallel to the tilt axis, the flow pattern is symmetrical and opposite in direction from that where the vessel is positioned vertically. The ground-based experiments are sufficient preliminary tests of theory and should be of significant interest regarding vapor deposited films in microgravity.

Frazier, D. O.↗

Analyzing Space-Based Interferometric Measurements of Stars and Network Measurements of Gamma-Ray Bursts

Since the announcement of the discovery of sources of bursts of gamma-ray radiation in 1973, hundreds more reports of such bursts have now been published. Numerous artificial satellites have been equipped with gamma-ray detectors including the very successful Compton Gamma Ray Observatory BATSE instrument. Unfortunately, we have made no progress in identifying the source(s) of this high energy radiation. We suspected that this was a consequence of the method used to define gamma-ray burst source "error boxes." An alternative procedure to compute gamma-ray burst source positions, with a purely physical underpinning, was proposed in 1988 by Taff. Since then we have also made significant progress in understanding the analytical nature of the triangulation problem and in computing actual gamma-ray burst positions and their corresponding error boxes. For the former, we can now mathematically illustrate the crucial role of the area occupied by the detectors, while for the latter, the Atteia et al. (1987) catalog has been completely re-reduced. There are very few discrepancies in locations between our results and those of the customary "time difference of arrival" procedure. Thus, we have numerically demonstrated that the end result, for the positions, of these two very different-looking procedures is the same. Finally, for the first time, we provide a sample of realistic "error boxes" whose non-simple shapes vividly portray the difficulty of burst source localization.

Taff, L. G.↗

Modeling Magnetic Properties in EZTB

A software module that calculates magnetic properties of a semiconducting material has been written for incorporation into, and execution within, the Easy (Modular) Tight-Binding (EZTB) software infrastructure. [EZTB is designed to model the electronic structures of semiconductor devices ranging from bulk semiconductors, to quantum wells, quantum wires, and quantum dots. EZTB implements an empirical tight-binding mathematical model of the underlying physics.] This module can model the effect of a magnetic field applied along any direction and does not require any adjustment of model parameters. The module has thus far been applied to study the performances of silicon-based quantum computers in the presence of magnetic fields and of miscut angles in quantum wells. The module is expected to assist experimentalists in fabricating a spin qubit in a Si/SiGe quantum dot. This software can be executed in almost any Unix operating system, utilizes parallel computing, can be run as a Web-portal application program. The module has been validated by comparison of its predictions with experimental data available in the literature.

Lee, Seungwon↗

Convergence of an iterative procedure for large-scale static analysis of structural components

The paper proves convergence of an iterative procedure for calculating the deflections of built-up component structures which can be represented as consisting of a dominant, relatively stiff primary structure and a less stiff secondary structure, which may be composed of one or more substructures that are not connected to one another but are all connected to the primary structure. The iteration consists in estimating the deformation of the primary structure in the absence of the secondary structure on the assumption that all mechanical loads are applied directly to the primary structure. The j-th iterate primary structure deflections at the interface are imposed on the secondary structure, and the boundary loads required to produce these deflections are computed. The cycle is completed by applying the interface reaction to the primary structure and computing its updated deflections. It is shown that the mathematical condition for convergence of this procedure is that the maximum eigenvalue of the equation relating primary-structure deflection to imposed secondary-structure deflection be less than unity, which is shown to correspond with the physical requirement that the secondary structure be more flexible at the interface boundary.

Austin, F.↗

Telescience Support Center Data System Software

The Telescience Support Center (TSC) team has developed a databasedriven, increment-specific Data Require - ment Document (DRD) generation tool that automates much of the work required for generating and formatting the DRD. It creates a database to load the required changes to configure the TSC data system, thus eliminating a substantial amount of labor in database entry and formatting. The TSC database contains the TSC systems configuration, along with the experimental data, in which human physiological data must be de-commutated in real time. The data for each experiment also must be cataloged and archived for future retrieval. TSC software provides tools and resources for ground operation and data distribution to remote users consisting of PIs (principal investigators), bio-medical engineers, scientists, engineers, payload specialists, and computer scientists. Operations support is provided for computer systems access, detailed networking, and mathematical and computational problems of the International Space Station telemetry data. User training is provided for on-site staff and biomedical researchers and other remote personnel in the usage of the space-bound services via the Internet, which enables significant resource savings for the physical facility along with the time savings versus traveling to NASA sites. The software used in support of the TSC could easily be adapted to other Control Center applications. This would include not only other NASA payload monitoring facilities, but also other types of control activities, such as monitoring and control of the electric grid, chemical, or nuclear plant processes, air traffic control, and the like.

Rahman, Hasan↗

Prediction of Ablation Rates from Solid Surfaces Exposed to High Temperature Gas Flow

A mathematical model and a solution algorithm is developed to study the physics of high temperature heat transfer and material ablation and identify the problems associated with the flow of hydrogen gas at very high temperatures and velocities through pipes and various components of Nuclear Thermal Rocket (NTR) motors. Ablation and melting can be experienced when the inner solid surface of the cooling channels and the diverging-converging nozzle of a Nuclear Thermal Rocket (NTR) motor is exposed to hydrogen gas flow at temperatures around 2500 degrees Kelvin and pressures around 3.4 MPa. In the experiments conducted on typical NTR motors developed in 1960s, degradation of the cooling channel material (cracking in the nuclear fuel element cladding) and in some instances melting of the core was observed. This paper presents the results of a preliminary study based on two types of physics based mathematical models that were developed to simulate the thermal-hydrodynamic conditions that lead to ablation of the solid surface of a stainless steel pipe exposed to high temperature hydrogen gas near sonic velocities. One of the proposed models is one-dimensional and assumes the gas flow to be unsteady, compressible and viscous. An in-house computer code was developed to solve the conservations equations of this model using a second-order accurate finite-difference technique. The second model assumes the flow to be three-dimensional, unsteady, compressible and viscous. A commercial CFD code (Fluent) was used to solve the later model equations. Both models assume the thermodynamic and transport properties of the hydrogen gas to be temperature dependent. In the solution algorithm developed for this study, the unsteady temperature of the pipe is determined from the heat equation for the solid. The solid-gas interface temperature is determined from an energy balance at the interface which includes heat transfer from or to the interface by conduction, convection, radiation, and ablation. Two different ablation models are proposed to determine the heat loss from the solid surface due to the ablation of the solid material. Both of them are physics based. Various numerical simulations were carried out using both models to predict the temperature distribution in the solid and in the gas flow, and then predict the ablation rates at a typical NTR motor hydrogen gas temperature and pressure. Solid mass loss rate per foot of a pipe was also calculated from these predictions. The results are presented for fully developed turbulent flow conditions in a sample SS pipe with a 6 inch diameter.

Akyuzlu, Kazim M.↗

User's manual for the REEDM (Rocket Exhaust Effluent Diffusion Model) computer program

The REEDM computer program predicts concentrations, dosages, and depositions downwind from normal and abnormal launches of rocket vehicles at NASA's Kennedy Space Center. The atmospheric dispersion models, cloud-rise models, and other formulas used in the REEDM model are described mathematically Vehicle and source parameters, other pertinent physical properties of the rocket exhaust cloud, and meteorological layering techniques are presented as well as user's instructions for REEDM. Worked example problems are included.

Bjorklund, J. R.↗

Linear Calibration of Radiographic Mineral Density Using Video-Digitizing Methods

Radiographic images can provide quantitative as well as qualitative information if they are subjected to densitometric analysis. Using modern video-digitizing techniques, such densitometry can be readily accomplished using relatively inexpensive computer systems. However, such analyses are made more difficult by the fact that the density values read from the radiograph have a complex, nonlinear relationship to bone mineral content. This article derives the relationship between these variables from the nature of the intermediate physical processes, and presents a simple mathematical method for obtaining a linear calibration function using a step wedge or other standard.

Martin, R. Bruce↗

Linear Calibration of Radiographic Mineral Density Using Video-Digitizing Methods

Radiographic images can provide quantitative as well as qualitative information if they are subjected to densitometric analysis. Using modem video-digitizing techniques, such densitometry can be readily accomplished using relatively inexpensive computer systems. However, such analyses are made more difficult by the fact that the density values read from the radiograph have a complex, nonlinear relationship to bone mineral content. This article derives the relationship between these variables from the nature of the intermediate physical processes, and presents a simple mathematical method for obtaining a linear calibration function using a step wedge or other standard.

Martin, R. Bruce↗

A Hierarchical Approach to Fracture Mechanics

Recent research conducted under NASA LaRC's Creativity and Innovation Program has led to the development of an initial approach for a hierarchical fracture mechanics. This methodology unites failure mechanisms occurring at different length scales and provides a framework for a physics-based theory of fracture. At the nanoscale, parametric molecular dynamic simulations are used to compute the energy associated with atomic level failure mechanisms. This information is used in a mesoscale percolation model of defect coalescence to obtain statistics of fracture paths and energies through Monte Carlo simulations. The mathematical structure of predicted crack paths is described using concepts of fractal geometry. The non-integer fractal dimension relates geometric and energy measures between meso- and macroscales. For illustration, a fractal-based continuum strain energy release rate is derived for inter- and transgranular fracture in polycrystalline metals.

Saether, Erik↗

Development of a three-dimensional time-dependent flow field model

A three-dimensional, time-dependent mathematical model to represent Mobile Bay was developed. Computer programs were developed which numerically solve the appropriate conservation equations for predicting bay and estuary flow fields. The model is useful for analyzing the dispersion of sea water into fresh water and the transport of sediment, and for relating field and physical model data.

Farmer, R. C.↗

Research in applied mathematics, numerical analysis, and computer science

Research conducted at the Institute for Computer Applications in Science and Engineering (ICASE) in applied mathematics, numerical analysis, and computer science is summarized and abstracts of published reports are presented. The major categories of the ICASE research program are: (1) numerical methods, with particular emphasis on the development and analysis of basic numerical algorithms; (2) control and parameter identification; (3) computational problems in engineering and the physical sciences, particularly fluid dynamics, acoustics, and structural analysis; and (4) computer systems and software, especially vector and parallel computers.

Source record↗