SEARCH · Engineering Papers
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.
Abstracts of research in the areas of biological rhythms, body fluids, botany, endrocrinology, enzymology, exobiology, genetics, human performance, immunology, life support systems, mathematical modeling, and numerous other topics related to space and life sciences are given.
Visualization software is widely used in scientific and engineering research. But computed visualizations can be very misleading, and the errors are easy to miss. We feel that the software producing the visualizations must be thoroughly evaluated and the evaluation process as well as the results must be made available. Testing and evaluation of visualization software is not a trivial problem. Several methods used in testing other software are helpful, but these methods are (apparently) often not used. When they are used, the description and results are generally not available to the end user. Additional evaluation methods specific to visualization must also be developed. We present several useful approaches to evaluation, ranging from numerical analysis of mathematical portions of algorithms to measurement of human performance while using visualization systems. Along with this brief survey, we present arguments for the importance of evaluations and discussions of appropriate use of some methods.
In this paper, we report the results of our recent research on the application of a multiprocessor Cray T916 supercomputer in modeling super-thermal electron transport in the earth's magnetic field. In general, this mathematical model requires numerical solution of a system of partial differential equations. The code we use for this model is moderately vectorized. By using Amdahl's Law for vector processors, it can be verified that the code is about 60% vectorized on a Cray computer. Speedup factors on the order of 2.5 were obtained compared to the unvectorized code. In the following sections, we discuss the methodology of improving the code. In addition to our goal of optimizing the code for solution on the Cray computer, we had the goal of scalability in mind. Scalability combines the concepts of portabilty with near-linear speedup. Specifically, a scalable program is one whose performance is portable across many different architectures with differing numbers of processors for many different problem sizes. Though we have access to a Cray at this time, the goal was to also have code which would run well on a variety of architectures.
In this chapter an analysis of the behavior of an arbitrary (perhaps massive) collective of computational processes in terms of an associated "world" utility function is presented We concentrate on the situation where each process in the collective can be viewed as though it were striving to maximize its own private utility function. For such situations the central design issue is how to initialize/update the collective's structure, and in particular the private utility functions, so as to induce the overall collective to behave in a way that has large values of the world utility. Traditional "team game" approaches to this problem simply set each private utility function equal to the world utility function. The "Collective Intelligence" (COIN) framework is a semi-formal set of heuristics that recently have been used to construct private utility. functions that in many experiments have resulted in world utility values up to orders of magnitude superior to that ensuing from use of the team game utility. In this paper we introduce a formal mathematics for analyzing and designing collectives. We also use this mathematics to suggest new private utilities that should outperform the COIN heuristics in certain kinds of domains. In accompanying work we use that mathematics to explain previous experimental results concerning the superiority of COIN heuristics. In that accompanying work we also use the mathematics to make numerical predictions, some of which we then test. In this way these two papers establish the study of collectives as a proper science, involving theory, explanation of old experiments, prediction concerning new experiments, and engineering insights.
The Transforming Growth Factor beta (TGFbeta) signaling pathway is a prominent regulatory signaling pathway controlling various important cellular processes. It can be induced by several factors, including ionizing radiation. It is regulated by Smads in a negative feedback loop through promoting increases in the regulatory Smads in the cell nucleus, and subsequent expression of inhibitory Smad, Smad7 to form a ubiquitin ligase with Smurf targeting active TGF receptors for degradation. In this work, we proposed a mathematical model to study the radiation-induced Smad-regulated TGF signaling pathway. By modularization, we are able to analyze each module (subsystem) and recover the nonlinear dynamics of the entire network system. Meanwhile the excitability, a common feature observed in the biological systems, along the TGF signaling pathway is discussed by mathematical analysis and numerical simulation.
Numerical simulations to assess the effectiveness of Generalized Predictive Control (GPC) for active control of dynamic systems having rigid-body modes are presented. GPC is a linear, time-invariant, multi-input/multi-output predictive control method that uses an ARX model to characterize the system and to design the controller. Although the method can accommodate both embedded (implicit) and explicit feedforward paths for incorporation of disturbance effects, only the case of embedded feedforward in which the disturbances are assumed to be unknown is considered here. Results from numerical simulations using mathematical models of both a free-free three-degree-of-freedom mass-spring-dashpot system and the XV-15 tiltrotor research aircraft are presented. In regulation mode operation, which calls for zero system response in the presence of disturbances, the simulations showed reductions of nearly 100%. In tracking mode operations, where the system is commanded to follow a specified path, the GPC controllers produced the desired responses, even in the presence of disturbances.
Introduction Before filling a propellant tank on the ground or in Space, the transfer line between the donor and receiver tanks must be cooled down preferably by sacrificing a minimum amount of the cryogenic fluid. The cryogenic line chill-down process involves a transition between different flow boiling regimes, namely, film boiling, transition film boiling, and nucleate boiling which are complex and may be quite gravity-dependent. Capturing these boiling phenomena and predicting the transition between them in a CFD framework is new and challenging both for 1g and microgravity applications. Materials & Methods The present work addresses this challenge by employing a two-phase Eulerian approach in the context of a homogeneous fluid mixture together with the Lee phase change model to capture the film boiling regime of the chill-down process using ANSYS Fluent®. The nucleate boiling regime is predicted by incorporating an in-house developed sub-grid model that accounts for bubble nucleation, bubble growth, bubble departure diameter, and their shedding frequency. The sub-grid model is encoded and implemented into Fluent via a user-defined function for the wall-fluid heat flux calculations. The mathematical formulation and numerical implementation of the CFD model are described in detail. The coupled CFD-Subgrid model is validated against published experimental data for liquid nitrogen chill-down of a heated stainless-steel pipe in 1g. Results Numerical simulation results show good agreements between the CFD predictions of the wall temperature evolution, rewetting temperature, and transition between film and nucleate boiling, with the experimental measurements published by Darr et al [2] for several different LN2 flowrates in the vertical pipe orientation. The CFD predictions for the wall temperature distribution indicate a rapid quenching of the wall at two upstream and downstream temperature sensing locations as compared to the experimental measurement. The only tuning parameter in the CFD model is the Lee mass transfer coefficient. The CFD Model predicts the Liedenfrost rewetting temperature in close agreement with the experiment. This marks a transition between stable and transitionary flow boiling regimes. The CFD-predicted boiling curve for the downstream sensor location is also compared against its experimental counterpart and indicates that the model is able to predict all the key temperature and heat flux parameters during the transitions from stable to transitionary film boiling to nucleate boiling in close agreement with the experiment. A sequence of predicted volume fraction, and temperature contours depicting these transitions will be presented.
In this paper we briefly describe a combined symbolic and numeric approach for solving mathematical models on parallel computers. An experimental software system, PIER, is being developed in Common Lisp to synthesize computationally intensive and domain formulation dependent phases of finite element analysis (FEA) solution methods. Quantities for domain formulation like shape functions, element stiffness matrices, etc., are automatically derived using symbolic mathematical computations. The problem specific information and derived formulae are then used to generate (parallel) numerical code for FEA solution steps. A constructive approach to specify a numerical program design is taken. The code generator compiles application oriented input specifications into (parallel) FORTRAN77 routines with the help of built-in knowledge of the particular problem, numerical solution methods and the target computer.
When establishing the pedigree of a simulation tool, code verification is used to ensure that the implemented numerical algorithm is a faithful representation of its underlying mathematical model. During this process, numerical results on various meshes are systematically compared to a reference analytic solution. The selection of analytic solutions can be a laborious process, as it is difficult to establish adequate code confidence without performing redundant work. In this report we address this issue by applying a physics-based process that establishes a set of reference problems. In this process, code simulation options are categorized and systematically tested, which ensures that gaps in testing are easily identified and addressed. The resulting problems are primarily intended for code verification analysis but may also be useful for comparison to other simulation codes, troubleshooting activities, or training exercises. The process is used to select fifteen code verification problems relevant for the one-dimensional steady-state heat conduction equation. These problems are applicable to a wide variety of simulation tools, but, in this work, a demonstration is performed using the finite element-based nuclear fuel performance code BISON. Convergence to the analytic solution at the theoretical rate is quantified for a selection of the problems, which establishes a baseline pedigree for the code. Not only can this standard set of conduction solutions be used for verification of other codes, but also the physics-based process for selecting problems can be utilized to quantify and expand testing for any simulation tool.
The reduction of the separation of the leading and following aircrafts is desirable to enhance the airport capacity provided that there is a physics-based operational model applicable to all regions of the flight domain (out of ground effect, OGE; near ground effect, NGE; and in ground effect, IGE) and that the quality of the quantitative input from the measurements of the prevailing atmospheric conditions and the quality of the total airport operations regarding the safety and the sound interpretation of the prevailing conditions match the quality of the analysis and numerical simulations. In the absence of an analytical solution, the physics of the flow is best expressed by a mathematical model based on numerical simulations, field and laboratory experiments, and heuristic reasoning. This report deals with the creation of a sound physics-based real-time IGE model of the aircraft wake vortices subjected to crosswind, stratification and shear.
Mathematical model and computer program for numerical solution of earth orbit differential equations of motion
Supercomputing systems are used for a wide range of computationally demanding tasks in many fields of science and engineering. They play a key role in numerical simulation, in which mathematical models are computed in order to simulate the behavior of physical systems. Scientists and engineers that use supercomputers for numerical simulation often have their productivity limited by the need to manually organize and manage extremely large amounts of data that are often produced and consumed by the software programs run on these systems. Recognizing these limitations, Kitware Inc. (Clifton Park, NY) and SLAC National Accelerator Laboratory (Menlo Park, CA) are developing an advanced software platform that can reduce the cognitive overhead required by knowledge workers when using supercomputers for numerical simulation. Phase I of the project is complete and includes the development of new capabilities for organizing simulation project files, improvements to the user interface and overall usability, and deployment of a “middle tier” server to sit between user desktop machines and supercomputers to offload much of the data management workload. The project also developed prototype software for executing sequences of numerical simulations, and a prototype for migrating supercomputing software to cloud-based computing systems to provide a potential alternative to supercomputers with different logistical and price-to-performance tradeoffs.
Two orbit computation methods were used: (1) numerical method- The solution to the satellite differential equations were solved in a step-by-step manner, using a mathematical algorithm taken from numerical analysis; and (2) analytical method - The solution was expressed by explicit functions of the independent variable. Analytical drag modules, tesseral terms initialization module, second order and long period terms module, and verification testing of the ASOP program were also considered.
Abstract. A common task in Lagrangian oceanography is to calculate a large number of drifter trajectories from a velocity field precalculated with an ocean model. Mathematically, this is simply numerical integration of an ordinary differential equation (ODE), for which a wide range of different methods exist. However, the discrete nature of the modelled ocean currents requires interpolation of the velocity field in both space and time, and the choice of interpolation scheme has implications for the accuracy and efficiency of the different numerical ODE methods. We investigate trajectory calculation in modelled ocean currents with 800 m, 4 km, and 20 km horizontal resolution, in combination with linear, cubic and quintic spline interpolation. We use fixed-step Runge–Kutta integrators of orders 1–4, as well as three variable-step Runge–Kutta methods (Bogacki–Shampine 3(2), Dormand–Prince 5(4) and 8(7)). Additionally, we design and test modified special-purpose variants of the three variable-step integrators, which are better able to handle discontinuous derivatives in an interpolated velocity field. Our results show that the optimal choice of ODE integrator depends on the resolution of the ocean model, the degree of interpolation, and the desired accuracy. For cubic interpolation, the commonly used Dormand–Prince 5(4) is rarely the most efficient choice. We find that in many cases, our special-purpose integrators can improve accuracy by many orders of magnitude over their standard counterparts, with no increase in computational effort. Equivalently, the special-purpose integrators can provide the same accuracy as standard methods at a reduced computational cost. The best results are seen for coarser resolutions (4 and 20 km), thus the special-purpose integrators are particularly advantageous for research using regional to global ocean models to compute large numbers of trajectories. Our results are also applicable to trajectory computations on data from atmospheric models.
Explore the source record for details and available documents.
Analytical drafting curves provide explicit mathematical expressions for any numerical data that appears in the form of graphical plots. The curves each have a reference coordinate axis system indicated on the curve as well as the mathematical equation from which the curve was generated.
The mathematical description and implementation of the statistical estimation procedure known as the Houston integrated spatial/spectral estimator (HISSE) is discussed. HISSE is based on a normal mixture model and is designed to take advantage of spectral and spatial information of LANDSAT data pixels, utilizing the initial classification and clustering information provided by the AMOEBA algorithm. The HISSE calculates parametric estimates of class proportions which reduce the error inherent in estimates derived from typical classify and count procedures common to nonparametric clustering algorithms. It also singles out spatial groupings of pixels which are most suitable for labeling classes. These calculations are designed to aid the analyst/interpreter in labeling patches with a crop class label. Finally, HISSE's initial performance on an actual LANDSAT agricultural ground truth data set is reported.
A numerical simulation of the actuation system for the propulsion control valve (PCV) of the NASA Langley Aircraft Landing Dynamics Facility was developed during the preliminary design of the PCV and used throughout the entire project. The simulation is based on a predictive model of the PCV which is used to evaluate and design the actuation system. The PCV controls a 1.7 million-pound thrust water jet used in propelling a 108,000-pound test carriage. The PCV can open and close in 0.300 second and deliver over 9,000 gallons of water per sec at pressures up to 3150 psi. The numerical simulation results are used to predict transient performance and valve opening characteristics, specify the hydraulic control system, define transient loadings on components, and evaluate failure modes. The mathematical model used for numerically simulating the mechanical fluid power system is described, and numerical results are demonstrated for a typical opening and closing cycle of the PCV. A summary is then given on how the model is used in the design process.