Search NASA⌕ Search

SEARCH · Search NASA

Results for “MATHEMATICS - BOUNDARY VALUES”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

Edward Lawrence Powers, Jr. December 30, 1915 - August 1, 2005; In Memorium: "It is only needed to say that he was a man ahead of his time, a man who led radiation research in all that he did." Michael G. Simic, Ph.D., D.Sc.; August, 2005

Larry Powers was a man of commitment: to roots, family, science, and large appetites. He grew up in Charleston, South Carolina, during the Great Depression that defined his remarkable work ethic and value for contribution. He graduated from the municipal College of Charleston in 1938 with majors in biology, chemistry, and mathematics. That same year, he married his wife Elly of Charleston, with whom he had seven daughters by the early years of his eras of discovery in 1959. He was awarded an Honorary Doctor of Letters from the College of Charleston in 1974, and retired there with Elly in 1987 as a Professor in Residence. Having obtained his Ph.D. from Johns Hopkins University in genetics and protozoology (1941), and having begun his professional career at the University of Notre Dame (Instructor and Assistant Professor, 1941-1945), Larry was often heard to comment on traditions of excellence in education and football, with adamant priority given to the first but without allowance for failure in the second.

Richmond, Robert C.↗
NASA NTRS1983

Space Shuttle propulsion parameter estimation using optimal estimation techniques, volume 1

The mathematical developments and their computer program implementation for the Space Shuttle propulsion parameter estimation project are summarized. The estimation approach chosen is the extended Kalman filtering with a modified Bryson-Frazier smoother. Its use here is motivated by the objective of obtaining better estimates than those available from filtering and to eliminate the lag associated with filtering. The estimation technique uses as the dynamical process the six degree equations-of-motion resulting in twelve state vector elements. In addition to these are mass and solid propellant burn depth as the ""system'' state elements. The ""parameter'' state elements can include aerodynamic coefficient, inertia, center-of-gravity, atmospheric wind, etc. deviations from referenced values. Propulsion parameter state elements have been included not as options just discussed but as the main parameter states to be estimated. The mathematical developments were completed for all these parameters. Since the systems dynamics and measurement processes are non-linear functions of the states, the mathematical developments are taken up almost entirely by the linearization of these equations as required by the estimation algorithms.

Source record↗
NASA NTRS2017

Cause and Cure - Deterioration in Accuracy of CFD Simulations With Use of High-Aspect-Ratio Triangular Tetrahedral Grids

Traditionally high-aspect ratio triangular/tetrahedral meshes are avoided by CFD re-searchers in the vicinity of a solid wall, as it is known to reduce the accuracy of gradient computations in those regions and also cause numerical instability. Although for certain complex geometries, the use of high-aspect ratio triangular/tetrahedral elements in the vicinity of a solid wall can be replaced by quadrilateral/prismatic elements, ability to use triangular/tetrahedral elements in such regions without any degradation in accuracy can be beneficial from a mesh generation point of view. The benefits also carry over to numerical frameworks such as the space-time conservation element and solution element (CESE), where triangular/tetrahedral elements are the mandatory building blocks. With the requirement of the CESE method in mind, a rigorous mathematical framework that clearly identities the reason behind the difficulties in use of such high-aspect ratio triangular/tetrahedral elements is presented here. As will be shown, it turns out that the degree of accuracy deterioration of gradient computation involving a triangular element is hinged on the value of its shape factor Gamma def = sq sin Alpha1 + sq sin Alpha2 + sq sin Alpha3, where Alpha1; Alpha2 and Alpha3 are the internal angles of the element. In fact, it is shown that the degree of accuracy deterioration increases monotonically as the value of Gamma decreases monotonically from its maximal value 9/4 (attained by an equilateral triangle only) to a value much less than 1 (associated with a highly obtuse triangle). By taking advantage of the fact that a high-aspect ratio triangle is not necessarily highly obtuse, and in fact it can have a shape factor whose value is close to the maximal value 9/4, a potential solution to avoid accuracy deterioration of gradient computation associated with a high-aspect ratio triangular grid is given. Also a brief discussion on the extension of the current mathematical framework to the tetrahedral-grid case along with some of the practical results of this extension is also provided. Furthermore, through the use of numerical simulations of practical viscous problems involving high-Reynolds number flows, the effectiveness of the gradient evaluation procedures within the CESE framework (that have their basis on the analysis presented here) to produce accurate and stable results on such high-aspect ratio meshes is also showcased.

CESE Method↗
DOE OSTI2024

An efficient quantum circuit for block encoding a pairing Hamiltonian

We present an efficient quantum circuit for block encoding a pairing Hamiltonian often studied in nuclear physics. Our block encoding scheme does not require mapping the creation and annihilation operators to the Pauli operators and representing the Hamiltonian as a linear combination of unitaries. Instead, we show how to encode the Hamiltonian directly using controlled swap operations. We analyze the gate complexity of the block encoding circuit and show that it scales polynomially with respect to the number of qubits required to represent a quantum state associated with the pairing Hamiltonian. We also show how the block encoding circuit can be combined with the quantum singular value transformation to construct an efficient quantum circuit for approximating the density of states of a pairing Hamiltonian. The techniques presented can be extended to encode more general second-quantized Hamiltonians.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2000

Robust Fault Detection for Aircraft Using Mixed Structured Singular Value Theory and Fuzzy Logic

The purpose of fault detection is to identify when a fault or failure has occurred in a system such as an aircraft or expendable launch vehicle. The faults may occur in sensors, actuators, structural components, etc. One of the primary approaches to model-based fault detection relies on analytical redundancy. That is the output of a computer-based model (actually a state estimator) is compared with the sensor measurements of the actual system to determine when a fault has occurred. Unfortunately, the state estimator is based on an idealized mathematical description of the underlying plant that is never totally accurate. As a result of these modeling errors, false alarms can occur. This research uses mixed structured singular value theory, a relatively recent and powerful robustness analysis tool, to develop robust estimators and demonstrates the use of these estimators in fault detection. To allow qualitative human experience to be effectively incorporated into the detection process fuzzy logic is used to predict the seriousness of the fault that has occurred.

Collins, Emmanuel G.↗
NASA NTRS1981

Shaded computer graphic techniques for visualizing and interpreting analytic fluid flow models

Mathematical models which predict the behavior of fluid flow in different experiments are simulated using digital computers. The simulations predict values of parameters of the fluid flow (pressure, temperature and velocity vector) at many points in the fluid. Visualization of the spatial variation in the value of these parameters is important to comprehend and check the data generated, to identify the regions of interest in the flow, and for effectively communicating information about the flow to others. The state of the art imaging techniques developed in the field of three dimensional shaded computer graphics is applied to visualization of fluid flow. Use of an imaging technique known as 'SCAN' for visualizing fluid flow, is studied and the results are presented.

Parke, F. I.↗
NASA NTRS1986

A simple nonlinear joint model

Hertzian contact theory is applied to a butt joint with specially mismatched bearing surfaces to devise a simple mathematical model of nonlinear axial force-displacement behavior in jointed members. Normalized tangent stiffness-force plots, for several values of a joint imperfection parameter, are presented for the sample case of solid structural members of circular section. The results illustrate the potential problem of high joint compliance at low axial-force levels, as well as the generally desirable stiffening and linearizing effects of preload. A nonlinear oscillator problem based on the static model is also formulated and solved to illustrate the effect of amplitude on natural frequency. As expected, natural frequency is low when amplitude is small. The results call attention to the important roles that tight tolerances and preload are expected to play in the design and fabrication of deployable and erectable truss-type space structures.

Fichter, W. B.↗
NASA NTRS1990

Low order dynamical models for mixing layers

The goal of this project was to investigate the feasibility of developing useful low order dynamical systems to allow one to determine cheaply the effect of various initial conditions on the evolution of mixing layers. A mathematically rigorous approach for a time evolving base flow as found in a mixing layer seems to be the singular value decomposition of the linear operator transforming some initial perturbations into the solution at a given later time. Unfortunately, this is very costly, if at all feasible, and only deals with the linear evolution of a disturbance. A different approach based on weakly nonlinear analysis was attempted, but numerous difficulties were encountered. Several formulations were proposed which may be ill-posed. The reasons for this possible failure of the attempted approach is related to the fundamental limitations of weakly nonlinear analysis. The clarification of these issues is attempted in this report.

Riahi, D. N.↗
NASA NTRS2005

Function Based Risk Assessment: Mapping Function to Likelihood

The concept of function offers a high potential for thinking and reasoning about designs as well as providing a common thread for relating together other design information. This paper focuses specifically on the relation between function and risk by presenting a mathematical mapping from product function to risk likelihood. This risk information is composed of design parameters, failure modes, and likelihood values. A spacecraft orientation subsystem, subsystem used to guide science instruments, and a Bell 206 rotorcraft are used to test the mapping which continues research on these products relating function to failure. Finally, a case study is presented in which the risk element likelihood is calculated for a fuel cell which is in the conceptual design phase at NASA JPL.

Tumer, Irem Y.↗
NASA NTRS2008

Thermal Capacitance (Slug) Calorimeter Theory Including Heat Losses and Other Decaying Processes

A mathematical model, termed the Slug Loss Model, has been developed for describing thermal capacitance (slug) calorimeter behavior when heat losses and other decaying processes are not negligible. This model results in the temperature time slope taking the mathematical form of exponential decay. When data is found to fit well to this model, it allows a heat flux value to be calculated that corrects for the losses and may be a better estimate of the cold wall fully catalytic heat flux, as is desired in arc jet testing. The model was applied to the data from a copper slug calorimeter inserted during a particularly severe high heating rate arc jet run to illustrate its use. The Slug Loss Model gave a cold wall heat flux 15% higher than the value of 2,250 W/sq cm obtained from the conventional approach to processing the data (where no correction is made for losses). For comparison, a Finite Element Analysis (FEA) model was created and applied to the same data, where conduction heat losses from the slug were simulated. The heat flux determined by the FEA model was found to be in close agreement with the heat flux determined by the Slug Loss Model.

Hightower, T. Mark↗
NASA NTRS2015

Software Validation via Model Animation

This paper explores a new approach to validating software implementations that have been produced from formally-verified algorithms. Although visual inspection gives some confidence that the implementations faithfully reflect the formal models, it does not provide complete assurance that the software is correct. The proposed approach, which is based on animation of formal specifications, compares the outputs computed by the software implementations on a given suite of input values to the outputs computed by the formal models on the same inputs, and determines if they are equal up to a given tolerance. The approach is illustrated on a prototype air traffic management system that computes simple kinematic trajectories for aircraft. Proofs for the mathematical models of the system's algorithms are carried out in the Prototype Verification System (PVS). The animation tool PVSio is used to evaluate the formal models on a set of randomly generated test cases. Output values computed by PVSio are compared against output values computed by the actual software. This comparison improves the assurance that the translation from formal models to code is faithful and that, for example, floating point errors do not greatly affect correctness and safety properties.

Dutle, Aaron M.↗
DOE OSTI2024

A computational analysis of effective R-values of buried ducts – the dynamic performance of buried ducts

Here, this paper evaluates the thermal performance of ducts partially or fully buried in loose-fill attic insulation. The overall thermal resistance between the ducts and the attic is referred to as an effective R-value. This paper shows a strong dependency of assumed attic temperature on the effective R-value. Based on the results, the effective R-value can be about twice as much with an attic temperature of 130 °F [54.4 °C], compared to when the attic temperature is 80 °F [26.7 °C]. Thus, this paper provides a polynomial regression equation based on a large set of simulations to determine the effective R-value of buried ducts depending on attic temperature and whether the HVAC system runs in cooling or heating mode. Further, the work presented in this paper investigated the potential impact of convective airflow within the attic insulation, particularly around the ducts. The analysis was based on computational fluid dynamics (CFD) and indicated that convectional forces are presented around the exterior surface of the ducts, but with negligible impact on the overall heat balance between the duct and the attic space.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Flexibility Options: A Proposed Product for Managing Imbalance Risk

The presence of variable renewable energy resources with uncertain outputs in day-ahead electricity markets results in additional balancing needs in real-time. Addressing those needs cost-effectively and reliably within a competitive market with unbundled products is challenging as both the demand for and the availability of flexibility depends on day-ahead energy schedules. Existing approaches for reserve procurement usually rely either on oversimplified demand curves that do not consider how system conditions that particular day affect the value of flexibility, or on bilateral trading of hedging instruments that are not co-optimized with day-ahead schedules. This article proposes a new product, ‘Flexibility Options', to address these two limitations. The demand for this product is endogenously determined in the day-ahead market and it is met cost-effectively by considering real-time supply curves for product providers, which are co-optimized with the energy supply. As we illustrate with numerical examples and mathematical analysis, the product addresses the hedging needs of participants with imbalances cost-effectively, provides a less intermittent revenue stream for participants with flexible outputs, promotes value-driven pricing of flexibility, and ensures that the system operator is revenue-neutral. This article provides a comprehensive design that can be further tested and applied in large-scale systems.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2026

Operator-level quantum acceleration of non-logconcave sampling

Sampling from probability distributions of the form 𝝈 ∝ e −𝜷V , where V is a continuous potential, is a fundamental task across physics, chemistry, biology, computer science, and statistics. However, when V is nonconvex, the resulting distribution becomes non-logconcave, and classical methods such as Langevin dynamics often exhibit poor performance. We introduce a quantum algorithm that provably accelerates a broad class of continuous-time sampling dynamics. For Langevin dynamics, our method encodes the target Gibbs measure into the amplitudes of aquantum state, identified as the kernel of a block matrix derived from a factorization of the Witten Laplacian operator. This connection enables Gibbs sampling via singular value thresholding and yields up to a quartic quantum speedup over best-knownclassical Langevin-based methods in the non-logconcave setting. Building on this framework, we further develop the first quantum algorithm that accelerates replica exchange Langevin diffusion, a widely used method for sampling from complex, rugged energy landscapes.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

Structure function calculations for Ostwald Ripening processes

A program for computing the structure function for configurations involved in Ostwald Ripening was written. The basic algorithms are derived from a mathematical analysis of a two-dimensional model system developed by Bortz, et. al. (1974). While it is expected that the values form the computer simulations will reflect Ostwald Ripening, at this point the program is still being tested. Some preliminary runs seem to justify the expectations.

Hassan, Razi A.↗
NASA NTRS2005

Specific Impulse and Mass Flow Rate Error

Specific impulse is defined in words in many ways. Very early in any text on rocket propulsion a phrase similar to .specific impulse is the thrust force per unit propellant weight flow per second. will be found.(2) It is only after seeing the mathematics written down does the definition mean something physically to scientists and engineers responsible for either measuring it or using someone.s value for it.

Gregory, Don A.↗
NASA NTRS2006

Hard and Soft Constraints in Reliability-Based Design Optimization

This paper proposes a framework for the analysis and design optimization of models subject to parametric uncertainty where design requirements in the form of inequality constraints are present. Emphasis is given to uncertainty models prescribed by norm bounded perturbations from a nominal parameter value and by sets of componentwise bounded uncertain variables. These models, which often arise in engineering problems, allow for a sharp mathematical manipulation. Constraints can be implemented in the hard sense, i.e., constraints must be satisfied for all parameter realizations in the uncertainty model, and in the soft sense, i.e., constraints can be violated by some realizations of the uncertain parameter. In regard to hard constraints, this methodology allows (i) to determine if a hard constraint can be satisfied for a given uncertainty model and constraint structure, (ii) to generate conclusive, formally verifiable reliability assessments that allow for unprejudiced comparisons of competing design alternatives and (iii) to identify the critical combination of uncertain parameters leading to constraint violations. In regard to soft constraints, the methodology allows the designer (i) to use probabilistic uncertainty models, (ii) to calculate upper bounds to the probability of constraint violation, and (iii) to efficiently estimate failure probabilities via a hybrid method. This method integrates the upper bounds, for which closed form expressions are derived, along with conditional sampling. In addition, an l(sub infinity) formulation for the efficient manipulation of hyper-rectangular sets is also proposed.

Crespo, L.uis G.↗
NASA NTRS2010

Kepler Mission's Focal Plane Characterization Models Implementation

The Kepler Mission photometer is an unusually complex array of CCDs. A large number of time-varying instrumental and systemic effects must be modeled and removed from the Kepler pixel data to produce light curves of sufficiently high quality for the mission to be successful in its planet-finding objective. After the launch of the spacecraft, many of these effects are difficult to remeasure frequently, and various interpolations over a small number of sample measurements must be used to determine the correct value of a given effect at different points in time. A library of software modules, called Focal Plane Characterization (FC) Models, is the element of the Kepler Science Data Pipeline (hereafter "pipeline") that handles this. FC, or products generated by FC, are used by nearly every element of the SOC processing chain. FC includes Java components: database persistence classes, operations classes, model classes, and data importers; and MATLAB code: model classes, interpolation methods, and wrapper functions. These classes, their interactions, and the database tables they represent, are discussed. This paper describes how these data and the FC software work together to provide the pipeline with the correct values to remove non-photometric effects caused by the photometer and its electronics from the Kepler light curves. The interpolation mathematics is reviewed, as well as the special case of the sky-to-pixel,pixel-to-sky coordinate transformation code, which incorporates a compound model that is unique in the SOC software.

mission↗