Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mathematical 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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Algebraic grid generation

Algebraic grid generation is the direct expression of a physical coordinate system as a function of a uniform grid in a rectangular computational coordinate system. Algebraic grid generation is based on mathematical interpolation and is presented in general terms of multivariate transfinite interpolation. The multisurface method and the two-boundary technique are described as univariate procedures that can be applied within the context of transfinite interpolation. A technique for grid clustering is described. Problems that are commonly encountered in three-dimensional grid generation are discussed and approaches for dealing with complex physical domains using multiple computational grid blocks are presented.

Smith, Robert E.↗
NASA NTRS2011

Project Management Using Modern Guidance, Navigation and Control Theory

Implementing guidance, navigation, and control (GN&C) theory principles and applying them to the human element of project management and control is not a new concept. As both the literature on the subject and the real-world applications are neither readily available nor comprehensive with regard to how such principles might be applied, this paper has been written to educate the project manager on the "laws of physics" of his or her project (not to teach a GN&C engineer how to become a project manager) and to provide an intuitive, mathematical explanation as to the control and behavior of projects. This paper will also address how the fundamental principles of modern GN&C were applied to the National Aeronautics and Space Administration's (NASA) Constellation Program (CxP) space suit project, ensuring the project was managed within cost, schedule, and budget. A project that is akin to a physical system can be modeled and managed using the same over arching principles of GN&C that would be used if that project were a complex vehicle, a complex system(s), or complex software with time-varying processes (at times nonlinear) containing multiple data inputs of varying accuracy and a range of operating points. The classic GN&C theory approach could thus be applied to small, well-defined projects; yet when working with larger, multiyear projects necessitating multiple organizational structures, numerous external influences, and a multitude of diverse resources, modern GN&C principles are required to model and manage the project. The fundamental principles of a GN&C system incorporate these basic concepts: State, Behavior, Feedback Control, Navigation, Guidance and Planning Logic systems. The State of a system defines the aspects of the system that can change over time; e.g., position, velocity, acceleration, coordinate-based attitude, and temperature, etc. The Behavior of the system focuses more on what changes are possible within the system; this is denoted in the state of the system. The behavior of a system, as captured in the system modeling, when properly done will aid in accurately predicting future system performance. The Feedback Control system understands the state and behavior of the system and uses feedback to adjust control inputs into the system. The feedback, which is the right arm of the Control system, allows change to be affected in the overall system; it therefore is important to not only correctly identify the system feedback inputs, but also the system response to the feedback inputs. The Navigation system takes multiple data inputs and based on a priori knowledge of the inputs, develops a statistically based weighting of the inputs and measurements to determine the system's state. Guidance and Planning Logic of the system, complete with an understanding of where the system is (provided by the Navigation system), will in turn determine where the system needs to be and how to get it there. With any system/project, it is critical that the objective of the system/project be clearly defined -- not only to plan but to measure performance and to aid in guiding the system or the project. The system principles discussed above, which can be and have been applied to the current CxP space suit development project, can also be mapped to real-world constituents, thus allowing project managers to apply systems theories that are well defined in engineering and mathematics to a discipline (i.e., Project Management) that historically has been based in personal experience and intuition. This mapping of GN&C theory to Project Management will, in turn, permit a direct, methodical approach to Project Management, planning and control providing a tool to help predict (and guide) performance and an understanding of the project constraints, how the project can be controlled, and the impacts to external influences and inputs. This approach, to a project manager, flows down to the three bottom-line variables of cost, schedule, and scope ando the needed control of these three variables to successfully perform and complete a project.

Hill, Terry R.↗
NASA NTRSDate not supplied

A Physics-Based Work-Energy Formulation for Real-Time Trajectory Guidance of A Lunar Lander

Throughout the years, many researchers have calculated and optimized trajectory solutions for lunar landing systems by employing sophisticated mathematical methods, that include: Hamilton’s Principle of Variation, Pontryagin’s maximum principle, and well known convex-optimization techniques among others. Many of these approaches typically require expensive computational resources to achieve convergence in the solution. In an effort to reduce complexity and the computational load required to obtain real-time guidance commands, a simple physics-based work-energy approach has been formulated. This approach is based on the dissipation of the mechanical energy of the vehicle to its final desired energy state required to achieve a safe landing. The rocket engine(s) employed during landing (among other maneuvers) dissipates mechanical energy by both doing work against the velocity vector of the vehicle (thus defining the trajectory path), and by jettisoning mass. Therefore, by solving the energy dissipation problem at every step of the maneuver, a much simpler formulation that naturally and quickly attains convergence is obtained. This formulation is not limited to approach, landing, and divert maneuvers, but in principle it can be employed during de-orbiting, braking burn, ascent, as well as orbit insertion.

Guidance↗
NASA NTRSDate not supplied

A Physics-Based Work-Energy Formulation for Real-Time Trajectory Guidance of a Lunar Lander

Throughout the years, many researchers have calculated and optimized trajectory solutions for lunar landing systems by employing sophisticated mathematical methods, that include: Hamilton’s Principle of Variation, Pontryagin’s maximum principle, and well known convex-optimization techniques among others. Many of these approaches typically require expensive computational resources to achieve convergence in the solution. In an effort to reduce complexity and the computational load required to obtain real-time guidance commands, a simple physics-based work-energy approach has been formulated. This approach is based on the dissipation of the mechanical energy of the vehicle to its final desired energy state required to achieve a safe landing. The rocket engine(s) employed during landing (among other maneuvers) dissipates mechanical energy by both doing work against the velocity vector of the vehicle (thus defining the trajectory path), and by jettisoning mass. Therefore, by solving the energy dissipation problem at every step of the maneuver, a much simpler formulation that naturally and quickly attains convergence is obtained. This formulation is not limited to approach, landing, and divert maneuvers, but in principle it can be employed during de-orbiting, braking burn, ascent, as well as orbit insertion.

Guidance↗
NASA NTRS2009

A Verification-Driven Approach to Traceability and Documentation for Auto-Generated Mathematical Software

Model-based development and automated code generation are increasingly used for production code in safety-critical applications, but since code generators are typically not qualified, the generated code must still be fully tested, reviewed, and certified. This is particularly arduous for mathematical and control engineering software which requires reviewers to trace subtle details of textbook formulas and algorithms to the code, and to match requirements (e.g., physical units or coordinate frames) not represented explicitly in models or code. Both tasks are complicated by the often opaque nature of auto-generated code. We address these problems by developing a verification-driven approach to traceability and documentation. We apply the AUTOCERT verification system to identify and then verify mathematical concepts in the code, based on a mathematical domain theory, and then use these verified traceability links between concepts, code, and verification conditions to construct a natural language report that provides a high-level structured argument explaining why and how the code uses the assumptions and complies with the requirements. We have applied our approach to generate review documents for several sub-systems of NASA s Project Constellation.

Denney, Ewen W.↗
NASA NTRS1973

Active cavity radiometer.

The active cavity radiometer (ACR) is a pyrheliometer that accurately defines the absolute radiation scale. The physics of the pyrheliometric method and the ACR approach to this method are presented in detail. A mathematical abstraction of the method is generated through a quasi-equilibrium analysis of the power balance of the ACR's cavity detector. An error analysis is carried out on the quasi-equilibrium equation to determine the uncertainties of ACR measurements relative to the absolute radiation scale. The uncertainty of ACR measurements as a function of irradiance level is presented in graphical form.

Willson, R. C.↗
NASA NTRS1980

An improved panel method for the solution of three-dimensional leading-edge vortex flows. Volume 1: Theory document

An improved panel method for the solution of three dimensional flow and wing and wing-body combinations with leading edge vortex separation is presented. The method employs a three dimensional inviscid flow model in which the configuration, the rolled-up vortex sheets, and the wake are represented by quadratic doublet distributions. The strength of the singularity distribution as well as shape and position of the vortex spirals are computed in an iterative fashion starting with an assumed initial sheet geometry. The method calculates forces and moments as well as detail surface pressure distributions. Improvements include the implementation of improved panel numerics for the purpose of elimination the highly nonlinear effects of ring vortices around double panel edges, and the development of a least squares procedure for damping vortex sheet geometry update instabilities. A complete description of the method is included. A variety of cases generated by the computer program implementing the method are presented which verify the mathematical assumptions of the method and which compare computed results with experimental data to verify the underlying physical assumptions made by the method.

Johnson, F. T.↗
NASA NTRS1991

Multiple methods integration for structural mechanics analysis and design

A new research area of multiple methods integration is proposed for joining diverse methods of structural mechanics analysis which interact with one another. Three categories of multiple methods are defined: those in which a physical interface are well defined; those in which a physical interface is not well-defined, but selected; and those in which the interface is a mathematical transformation. Two fundamental integration procedures are presented that can be extended to integrate various methods (e.g., finite elements, Rayleigh Ritz, Galerkin, and integral methods) with one another. Since the finite element method will likely be the major method to be integrated, its enhanced robustness under element distortion is also examined and a new robust shell element is demonstrated.

Housner, J. M.↗
NASA NTRS2019

High-Performance Computing Optimization for Aladyn – Adaptive Neural Network Molecular Dynamics Mini-Application

This report provides a description and performance evaluation of the optimization techniques for high performance computing (HPC) implementation of the open source Computational Materials mini-application Aladyn (https://github.com/nasa/aladyn). Aladyn is a basic molecular dynamics code written in FORTRAN 2003, which is designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to efficiently reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of the more expensive quantum mechanics-based calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory (DFT) method. While achieving orders of magnitude faster computational performance than DFT, the ANN-based approach was still very computationally demanding compared to the conventional approach of using empirically fitted energy functions. After its initial development, Aladyn was evaluated and optimized by experts at the NASA Advanced Supercomputing (NAS) division to exploit modern supercomputer architectures. The code has been optimized for execution on multicore central processing units (CPUs), including Intel® Skylake microarchitecture, and on graphic accelerators, such as Nvidia® V100 graphic processing units (GPUs), using Open Multi-Processing (OpenMP) and Open Accelerators (OpenACC) programming interfaces. The optimization achieved a speedup of 4.7 times the baseline version on CPU performance and an additional 2.4 times on CPU+GPU performance. Atomistic computer simulations are a fundamental tool in materials research to model material properties form physics-based first principles. Atomic interaction, governed by Quantum Mechanics (QM) require sophisticated and highly computationally demanding mathematical models to calculate [1]. Classical methods use approximate functional forms, empirically fitted through a set of variable parameters to emulate atomic energies as direct functions of atomic coordinates [2]. While empirical potentials are computationally much simpler, allowing simulations of large-scale systems of up to a trillion (1012) atoms [3], they are substantially less accurate compared to quantum calculations and applicable only to very specific atomic configurations or predefined crystallographic phases. A recently suggested approach is to use heuristic machine learning methods [4], such as those based on Adaptive Neural Networks (ANNs) to predict atomic energies, after being trained on a sufficiently large database of QM-calculated structures [5,6]. This approach reduces significantly the computational complexity, allowing for simulations of orders of magnitude larger systems compared to QM-based methods without compromising accuracy. Still, compared to classical methods using empirical energy functions, ANN methods remain two- to three orders of magnitude more computationally demanding. Hence, the computational cost of simulations, together with the need for extensive training of ANNs, still makes the practical implementation of ANN-based methods quite challenging. The purpose of the Aladyn mini-application software [7], available as open source at https://github.com/nasa/aladyn, is to be a testbed for exploring possible optimization strategies to develop highly scalable parallel algorithms for ANN-based atomistic simulations. Aladyn is aimed at utilizing the architecture of the high-end modern highperformance computing (HPC) hardware based on multicore central processing units (CPUs) equipped with graphic processing unit (GPU) accelerators. Specifically, the goal is to optimize the performance on a single HPC compute node, before implementing scaling to multi-node parallelization using message passing interface (MPI). At the same time, the open source code of Aladyn can serve as a training model for students and professors in academia.

Yamakov, Vesselin I.↗
NASA NTRS1989

Interaction of surface radiation with convection in crystal growth by physical vapor transport

Growth of single crystals from vapor in closed ampoules is governed by an intricate interplay between mass, momentum and heat transfer processes. The objective of this study is to examine and isolate the effects of surface radiation heat transfer on the vapor transport process using a mathematical model. The model consists of a set of coupled nonlinear partial differential equations for conservation of mass, momentum, energy and species, and the integrodifferential equations which represent radiative exchange. It depends on five important physical parameters. These are Grashof number, Prandtl number, Schmidt number, aspect ratio and the radiation-conduction number. The effects of these dimensionless groupings are systematically investigated. From the cases examined, it is concluded that surface radiation can change the flow structure appreciably. This is especially true in microgravity environment where radiation competes primarily with conduction in modifying the thermal profiles. The numerical results also show that in the presence of radiation, the top heating configuration (source on top) is no longer stable and that near the growing crystal, radiation-induced vortices can introduce significant nonuniformities in the growth flux.

Kassemi, Mohammad↗
NASA NTRS2003

Mathematical Modeling of Food Supply for Long Term Space Missions Using Advanced Life Support

A habitat for long duration missions which utilizes Advanced Life Support (ALS), the Bioregenerative Planetary Life Support Systems Test Complex (BIO-Plex), is currently being built at JSC. In this system all consumables will be recycled and reused. In support of this effort, a menu is being planned utilizing ALS crops that will meet nutritional and psychological requirements. The need exists in the food system to identify specific physical quantities that define life support systems from an analysis and modeling perspective. Once these quantities are defined, they need to be fed into a mathematical model that takes into consideration other systems in the BIO-Plex. This model, if successful, will be used to understand the impacts of changes in the food system on the other systems and vice versa. The Equivalent System Mass (ESM) metric has been used to describe systems and subsystems, including the food system options, in terms of the single parameter, mass. There is concern that this approach might not adequately address the important issues of food quality and psychological impact on crew morale of a supply of fiesh food items. In fact, the mass of food can also depend on the quality of the food. This summer faculty fellow project will involve creating an appropriate mathematical model for the food plan developed by the Food Processing System for BIO-Plex. The desired outcome of this work will be a quantitative model that can be applied to the various options of supplying food on long-term space missions.

Cruthirds, John E.↗
NASA NTRS1991

One-way nesting for a primitive equation ocean model

Prognostic numerical models for atmospheric and oceanic circulations require initial fields, boundary conditions, and forcing functions in addition to a consistent set of partial differential equations, including a state relation and equations expressing conservation of mass, momentum, and energy. Depending on the horizontal domain to be modeled, the horizontal boundary conditions are either physically obvious or extremely difficult to specify consistently. If the entire atmosphere is modeled, periodic horizontal boundary conditions are appropriate. On the other hand, the physical horizontal boundaries on the entire ocean are solid walls. Obviously, the normal velocity at a solid wall is zero while the specification of the tangential velocity depends on the mathematical treatment of the horizontal viscous terms. Limitations imposed by computer capacity and cost, as well as research interests, have led to the use of limited area models to study flows in the atmosphere and ocean. The limited area models do not have physical horizontal boundaries, merely numerical ones. Correctly determining these open boundary conditions for limited-area numerical models has both intrigued and frustrated numerical modelers for decades. One common approach is to use the closed or solid wall boundary conditions for a limited-area model. The argument given for this approach is that the boundary conditions affect flow near the walls but that none of these effects are propagated into the interior. Therefore, one chooses a big enough domain that the central region of interest is not corrupted by the boundary flow. Research in progress to model the North Atlantic circulation vividly illustrates the pitfalls of this approach. Two model runs are compared: (1) the southern boundary at 20S between latitudes 0 and 40W is artificially closed; and (2) the same boundary is specified as open with an inward transport of 15 Sv (determined from a global model with the same physics) uniformly spread across the boundary. A comparison of both runs is presented.

Blake, D. W.↗
NASA NTRS1969

Flow properties of suspensions rich in solids

Mathematical evaluation of flow properties of fluids carrying high concentrations of solids in suspension relates suspension viscosity to physical properties of the solids and liquids, and provides a means for predicting flow behavior. A technique for calculating a suspensions flow rates is applicable to the design of pipelines.

Armstrong, W. P.↗
NASA NTRS1973

A representation of Jacchia's thermospheric models in spherical harmonics

The Jacchia models are represented in terms of spherical harmonic functions. This representation has the advantages of ease of comparison with theoretical and other observational models and data, mathematical analyticity and relative simplicity. The symmetry properties of the models are emphasized by this representation and some physical characteristics like the increase of the amplitude of the diurnal density variation with decreasing solar activity become more apparent.

Blum, P.↗
NASA NTRS1977

Orbit-orbit resonances among natural satellites

A qualitative physical description of the orbital resonance mechanism based on a previously published (Greenberg et al. 1972; Greenberg, 1973) mathematical analysis is presented. The fundamental principles of the resonance mechanism are illustrated by an examination of the Titan-Hyperion case. More complex resonances, including Enceladus-Dione, Mimas-Tethys are considered. Resonance phenomena involving the rings of Saturn, the inner three satellites of Uranus, and the inner three Galilean satellites of Jupiter are also discussed. The role of analyses of resonance phenomena in developing theories of satellite formation and evolution is examined.

Greenberg, R.↗
NASA NTRS1980

Cost decomposition of linear systems with application to model reduction

A means is provided to assess the value or 'cst' of each component of a large scale system, when the total cost is a quadratic function. Such a 'cost decomposition' of the system has several important uses. When the components represent physical subsystems which can fail, the 'component cost' is useful in failure mode analysis. When the components represent mathematical equations which may be truncated, the 'component cost' becomes a criterion for model truncation. In this latter event component costs provide a mechanism by which the specific control objectives dictate which components should be retained in the model reduction process. This information can be valuable in model reduction and decentralized control problems.

Skelton, R. E.↗
NASA NTRS1988

Theoretical and Experimental Studies of the Transonic Flow Field and Associated Boundary Conditions near a Longitudinally-Slotted Wind-Tunnel Wall

A theoretical examination of the slotted-wall flow field is conducted to determine the appropriate wall pressure drop (or boundary condition) equation. This analysis improves the understanding of the fluid physics of these types of flow fields and helps in evaluating the uncertainties and limitations existing in previous mathematical developments. It is shown that the resulting slotted-wall boundary condition contains contributions from the airfoil-induced streamline curvature and the non-linear, quadratic, slot crossflow in addition to an often neglected linear term which results from viscous shearing in the slot. Existing and newly acquired experimental data are examined in the light of this formulation and theoretical developments.

Everhart, Joel Lee↗
NASA NTRS1994

Active control of instabilities in laminar boundary-layer flow. Part 1: An overview

This paper (the first in a series) focuses on using active-control methods to maintain laminar flow in a region of the flow in which the natural instabilities, if left unattended, lead to turbulent flow. The authors review previous studies that examine wave cancellation (currently the most prominent method) and solve the unsteady, nonlinear Navier-Stokes equations to evaluate this method of controlling instabilities. It is definitely shown that instabilities are controlled by the linear summation of waves (i.e., wave cancellation). Although a mathematically complete method for controlling arbitrary instabilities has been developed (but not yet tested), the review, duplication, and physical explanation of previous studies are important steps for providing an independent verification of those studies, for establishing a framework for subsequent work which will involve automated transition control, and for detailing the phenomena by which the automated studies can be used to expand knowledge of flow control.

Joslin, Ronald D.↗