Analysis of lower atmospheric data for diffusion studies
Mathematical diffusion models and meteorological inputs for estimating toxic fuel hazards at Kennedy Space Center
SEARCH · Search NASA
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.
Mathematical diffusion models and meteorological inputs for estimating toxic fuel hazards at Kennedy Space Center
An analogy between the mathematical modeling of transonic potential flow and the flow in a cavitating bearing is described. Based on the similarities, characteristics of the cavitated region and jump conditions across the film reformation and rupture fronts are developed using the method of weak solutions. The mathematical analogy is extended by utilizing a few computational concepts of transonic flow to numerically model the cavitating bearing. Methods of shock fitting and shock capturing are discussed. Various procedures used in transonic flow computations are adapted to bearing cavitation applications, for example, type differencing, grid transformation, an approximate factorization technique, and Newton's iteration method. These concepts have proved to be successful and have vastly improved the efficiency of numerical modeling of cavitated bearings.
An analogy between the mathematical modeling of transonic potential flow and the flow in a cavitating bearing is described. Based on the similarities, characteristics of the cavitated region and jump conditions across the film reformation and rupture fronts are developed using the method of weak solutions. The mathematical analogy is extended by utilizing a few computational concepts of transonic flow to numerically model the cavitating bearing. Methods of shock fitting and shock capturing are discussed. Various procedures used in transonic flow computations are adapted to bearing cavitation applications, for example, type differencing, grid transformation, an approximate factorization technique, and Newton's iteration method. These concepts have proved to be successful and have vastly improved the efficiency of numerical modeling of cavitated bearings.
Report presents study of finite-element mathematical modeling as used in analyzing stresses and strains at joints between thin, shell-like components (e.g., ducts) and thicker components (e.g., flanges or engine blocks). First approach uses global/local model to evaluate system. Provides correct total response and correct representation of stresses away from any discontinuities. Second approach involves development of special transition finite elements to model transitions between shells and thicker structural components.
Physical and mathematical formulae and concepts are fundamental elements of scientific and engineering software. These classical equations and methods are time tested, universally accepted, and relatively unambiguous. The existence of this classical ontology suggests an ideal problem for automated comprehension. This problem is further motivated by the pervasive use of scientific code and high code development costs. To investigate code comprehension in this classical knowledge domain, a research prototype has been developed. The prototype incorporates scientific domain knowledge to recognize code properties (including units, physical, and mathematical quantity). Also, the procedure implements programming language semantics to propagate these properties through the code. This prototype's ability to elucidate code and detect errors will be demonstrated with state of the art scientific codes.
This document walks a user through a demonstration of working with PDV digitizer data using python. This demonstration and included suggested exercises will be used at the 2024 PDV workshop hands-on session as an example and skill-development training session. The tutorial allows the user to generate synthetic but realistic PDV waveform data and visualize/inspect the results using spectrograms and waveform viewing tools.
Cyber-physical systems (CPS) are engineered systems that rely on the smooth integration of computational algorithms and physical elements. This integration presents new challenges for verifying that systems will behave as expected. The goal of this presentation is to present current challenges and potential solutions for the formal verification of cyber-physical systems. For cyber systems, formal methods refer to systematically rigorous mathematical techniques employed in the specification, development, analysis, and verification of both software and hardware systems. Recent advancements in computer science have yielded sophisticated tools specifically designed to address challenges associated with formal methods in complex systems. These tools leverage various foundational concepts such as logic, formal languages, program semantics, type systems, type theory, and automata theory. A notable achievement in the application of formal methods is the seL4 microkernel, claimed to be the first general-purpose operating-system kernel to be verified. Its proof implies the absence of bugs and guarantees that the kernel meets specifications. For physical systems, dynamic and control theory has a history of using rigorous analytic techniques to prove functional correctness. Lyapunov, optimal, classical, modern, and robust control theories all provide rigorous mathematical methods both to analyze system performance and to design controller that can be guaranteed to meet certain objectives. Recent computational techniques like level set theory and reachability analysis provide assertions that a system's state will avoid unsafe regions. Even though success has been independently achieved for cyber systems and physical systems, the integration of such systems creates new challenges. In particular, there is an obvious discrepancy between finite-state machines and infinite-state systems, resulting in different approaches for modeling and analyzing these system. While it is possible to simulate hybrid systems, this provides only a demonstration of a performance and not proof. For hybrid systems, current formal methods and system analysis approaches typically require a workarounds to work on hybrid systems like CPS. This paper will outline the state of the art and limits of current practice for formally verifying CPS and will identify possible research directions that require attention.
An important consideration in the global local finite-element method (GLFEM) is the availability of global functions for the given problem. The role and mathematical requirements of these global functions in a GLFEM analysis of localized stress states in prismatic structures are discussed. A method is described for determining these global functions. Underlying this method are theorems due to Toupin and Knowles on strain energy decay rates, which are related to a quantitative expression of Saint-Venant's principle. It is mentioned that a mathematically complete set of global functions can be generated, so that any arbitrary interface condition between the finite element and global subregions can be represented. Convergence to the true behavior can be achieved with increasing global functions and finite-element degrees of freedom. Specific attention is devoted to mathematically two-dimensional and three-dimensional prismatic structures. Comments are offered on the GLFEM analysis of NASA flat panel with a discontinuous stiffener. Methods for determining global functions for other effects are also indicated, such as steady-state dynamics and bodies under initial stress.
The relation between localization of interference fringes in classical and holographic interferometry is reviewed and an application of holographic interferometry is considered for which the object is a transparent medium with nonhomogeneous refractive index. The technique is based on the analysis of the optical path length change of the object wave as it propagates through a transparent medium. Phase shifts due to variations of the speed of light within the medium give rise to an interference pattern. The resulting interferogram can be used to determine the physical properties of the medium or transparent object. Such properties include the mass density of fluids, electron densities of plasmas, the temperature of fluids, the chemical species concentration of fluids, and the state of stress in solids. The optical wave used can be either a simple plane or spherical wave, or it may be a complicated spatial wave scattered by a diffusing screen. The mathematical theory on the formation and analysis of localized fringes, the general theoretical concepts used, and a computer code for analysis are included along with the inversion of fringe order data.
Not Available
This talk explores how advanced visualization technologies are transforming analytical reasoning and knowledge discovery in energy research, drawing on recent work at the National Laboratory of the Rockies' Computational Science Center. Through a series of scientific case studies, we demonstrate how immersive and high-resolution visualization environments enable scientists and engineers to identify previously unseen patterns and features - insights that often remain hidden in traditional desktop-based analysis. By embedding richer information into interactive analytics tools, these approaches support the exploration of complex, multivariate parameter spaces, where interaction itself catalyzes understanding. Beyond capability, we emphasize the critical role of visualization design grounded in perception and cognition, showing how visual encodings directly influence analytical outcomes. Spanning applications from materials science to integrated energy systems, these visualization approaches accelerate innovation and improve decision-making by enabling deeper, more reliable insight into increasingly complex energy data.
The NASTRAN thermal analyzer (NTA) which performs large-scale unified thermo-structural analyses with the NASTRAN (NASA structural analysis) computer program is described. The mathematical similitude between these two distinct disciplines of thermal and structure is examined. It serves as the theoretical basis upon which the implementation of the thermal capability in NASTRAN was accomplished. The program structure, the functional flow, the solution algorithms, the organization of an input data deck and the solution capabilities of NTA are summarized. Emphasis is placed on the interface of the unified approach in thermo-structural analyses where stresses, deflections, vibrations and bucklings induced by the effect of temperature change are of concern. Attentions are also directed to the preprocessor and post processors. As a specially designed preprocessor, the VIEW program is capable of generating exchange factors which can be output, at user's option, in formats compatible with that required by NTA. Two post processors that serve specific objectives are included. They are the thermal variance analysis and the graphical displaying capability of temperatures in color or black and white.
The structural analysis methods research has several goals. One goal is to develop analysis methods that are general. This goal of generality leads naturally to finite-element methods, but the research will also include other structural analysis methods. Another goal is that the methods be amenable to error analysis; that is, given a physical problem and a mathematical model of that problem, an analyst would like to know the probable error in predicting a given response quantity. The ultimate objective is to specify the error tolerances and to use automated logic to adjust the mathematical model or solution strategy to obtain that accuracy. A third goal is to develop structural analysis methods that can exploit parallel processing computers. The structural analysis methods research will focus initially on three types of problems: local/global nonlinear stress analysis, nonlinear transient dynamics, and tire modeling.
This paper uses a mathematical framework to identify the interdependencies of key variables in supersonic transportation demand analysis. The existing quantitative models for supersonic transportation demand analysis are compared for consistency in modeling the interdependencies. Unlike other quantitative models, the Low-Boom Systems Analysis Model (LBSAM2) can propagate important design characteristics of a supersonic transport aircraft concept to the level of economic metrics (such as the number of future supersonic passengers), with a consistent coupling of the market demand analysis, detailed mission analysis, and low-boom constraint. This allows the use of the detailed demand analysis results from LBSAM2 to maximize the economic viability of a supersonic transport aircraft by finding favorable system-level trades between weight, range, fuel burn, and assumed sonic boom ground noise limit for supersonic overland flight. In this paper, LBSAM2 is integrated with conceptual low-boom design to improve the economic viability of low-boom supersonic transport aircraft. A brief discussion of uncertainties in the LBSAM2 analysis is also included, focusing on their impacts on the relative economic advantages between low-boom concepts.
This paper uses a mathematical framework to identify the interdependencies of key variables in supersonic transportation demand analysis. The existing quantitative models for supersonic transportation demand analysis are compared for consistency in modeling the interdependencies. Unlike other quantitative models, the Low-Boom Systems Analysis Model (LBSAM2) can propagate important design characteristics of a supersonic transport aircraft concept to the level of economic metrics (such as the number of future supersonic passengers), with a consistent coupling of the market demand analysis, detailed mission analysis, and low-boom constraint. This allows the use of the detailed demand analysis results from LBSAM2 to maximize the economic viability of a supersonic transport aircraft by finding favorable system-level trades between weight, range, fuel burn, and assumed sonic boom ground noise limit for supersonic overland flight. In this paper, LBSAM2 is integrated with conceptual low-boom design to improve the economic viability of low-boom supersonic transport aircraft. A brief discussion of uncertainties in the LBSAM2 analysis is also included, focusing on their impacts on the relative economic advantages between low-boom concepts.
A general survey of the progress made in the areas of mathematical modelling of the system dynamics, structural analysis, development of control algorithms, and simulation of environmental disturbances is presented. The use of graph theory techniques is employed to examine the effects of inherent damping associated with LSST systems on the number and locations of the required control actuators. A mathematical model of the forces and moments induced on a flexible orbiting beam due to solar radiation pressure is developed and typical steady state open loop responses obtained for the case when rotations and vibrations are limited to occur within the orbit plane. A preliminary controls analysis based on a truncated (13 mode) finite element model of the 122m. Hoop/Column antenna indicates that a minimum of six appropriately placed actuators is required for controllability. An algorithm to evaluate the coefficients which describe coupling between the rigid rotational and flexible modes and also intramodal coupling was developed and numerical evaluation based on the finite element model of Hoop/Column system is currently in progress.
Early space radiation shield code development relied on Monte Carlo methods and made important contributions to the space program. Monte Carlo methods have resorted to restricted one-dimensional problems leading to imperfect representation of appropriate boundary conditions. Even so, intensive computational requirements resulted and shield evaluation was made near the end of the design process. Resolving shielding issues usually had a negative impact on the design. Improved spacecraft shield design requires early entry of radiation constraints into the design process to maximize performance and minimize costs. As a result, we have been investigating high-speed computational procedures to allow shield analysis from the preliminary concept to the final design. For the last few decades, we have pursued deterministic solutions of the Boltzmann equation allowing field mapping within the International Space Station (ISS) in tens of minutes using standard Finite Element Method (FEM) geometry common to engineering design methods. A single ray trace in such geometry requires 14 milliseconds and limits application of Monte Carlo methods to such engineering models. A potential means of improving the Monte Carlo efficiency in coupling to spacecraft geometry is given.
This study presents a mathematical formulation developed for aerodynamic sensitivity coefficients based on a discretized form of the compressible 2D Euler equations. A brief motivating introduction to the aerodynamic sensitivity analysis and the reasons behind an integrated flow/sensitivity analysis for design algorithms are presented. The finite difference approach and the quasi-analytical approach are used to determine the aerodynamic sensitivity coefficients. A new flow prediction concept, which is an outcome of the direct method in the quasi-analytical approach, is developed and illustrated with an example. Surface pressure coefficient distributions of a nozzle-afterbody configuration obtained from the predicted flowfield solution are compared successfully with their corresponding values obtained from a flowfield analysis code and the experimental data.