Search NASASearch

SEARCH · Search NASA

Results for “Mathematical Methods”

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 91 records · Page 5

On the reflectivity of complex mesh surfaces

Poorer than expected surface reflectivity was observed in an early Tracking and Data Relay Satellite System antenna utilizing a tricot mesh weave. This poor reflectivity was determined to be caused by inadequate electrical contact at wire crossover points. A proper mathematical and numerical approach to assess the impact of wire junctions on reflectivity performance is developed. A mathematical method is presented for computing the surface reflectivity of complex mesh configurations like those on unfurlable-type spacecraft antennas. The method is based on the Floquet mode expansion to establish an integral equation for mesh wire currents. The equation is solved using the method of moments with triangular basis functions. It is observed that it is necessary to give special attention to the junction treatment among different branches of the mesh configurations. A vector junction current approach that resulted in satisfactory solutions for the current is described. The results of numerical simulations are compared against measured data and excellent agreement is observed.

Imbriale, William A.

Atmospheric constituent density profiles from full disk solar occultation experiments

Mathematical methods are described which permit the derivation of the number of density profiles of atmospheric constituents from solar occultation measurements. The algorithm is first applied to measurements corresponding to an arbitrary solar-intensity distribution to calculate the normalized absorption profile. The application of Fourier transform to the integral equation yields a precise expression for the corresponding number density, and the solution is employed with the data given in the form of Laguerre polynomials. The algorithm is employed to calculate the results for the case of uniform distribution of solar intensity, and the results demonstrate the convergence properties of the method. The algorithm can be used to effectively model representative model-density profiles with constant and altitude-dependent scale heights.

Lumpe, J. D.

Feedback control for unsteady flow and its application to the stochastic Burgers equation

The study applies mathematical methods of control theory to the problem of control of fluid flow with the long-range objective of developing effective methods for the control of turbulent flows. Model problems are employed through the formalism and language of control theory to present the procedure of how to cast the problem of controlling turbulence into a problem in optimal control theory. Methods of calculus of variations through the adjoint state and gradient algorithms are used to present a suboptimal control and feedback procedure for stationary and time-dependent problems. Two types of controls are investigated: distributed and boundary controls. Several cases of both controls are numerically simulated to investigate the performances of the control algorithm. Most cases considered show significant reductions of the costs to be minimized. The dependence of the control algorithm on the time-descretization method is discussed.

Choi, Haecheon

Elements of the Wing Section Theory and of the Wing Theory

This report contains those results of the theory of wings and of wing sections which are of immediate practical value. They are proved and demonstrated by the use of the simple conceptions of "kinetic energy" and "momentum" only, familiar to every engineer; and not by introducing "isogonal transformations" and "vortices," which latter mathematical methods are not essential to the theory and better are used only in papers intended for mathematicians and special experts.

Munk, Max M

A resume of the advances in theoretical aeronautics made by Max M. Munk

In order to apply profitably the mathematical methods of hydrodynamics to aeronautical problems, it is necessary to make simplifications in the physical conditions of the latter. To begin with, it is allowable in many problems, as Prandtl has so successfully shown, to treat the air as having constant density and as free of viscosity. But this is not sufficient. It is also necessary to specify certain shapes for the solid bodies whose motion through the air is discussed, shapes suggested by the actual solids - airships or airfoils - it is true, but so chosen that they lead to solvable problems. In a valuable paper presented by Dr. Max M. Munk, of the National Advisory Committee for Aeronautics, Washington, to the Delft Conference in April, 1924, these necessary simplifying assumptions are discussed in detail. It is the purpose of the present paper to present in as simple a manner as possible some of the interesting results obtained by Dr. Munk's methods.

Ames, Joseph S

The Theories of Turbulence

The theory of turbulence reached its full growth at the end of the 19th century as a result of the work by Boussinesq and Reynolds. It then underwent a long period of stagnation which ended under the impulse given to it by the development of wind tunnels caused by the needs of aviation. Numerous researchers, attempted to put Reynolds' elementary statistical theory into a more precise form. During the war, some isolated scientists - von Weizsacker and Heisenberg in Germany, Kolmogoroff in Russia, Onsager in the U.S.A. - started a program of research. By a system of assumptions which make it possible to approach the structure of turbulence in well-defined limiting conditions quantitatively, they obtained a certain number of laws on the correlations and the spectrum. Since the late reports have improved the mathematical language of turbulence, it was deemed advisable to start with a detailed account of the mathematical methods applicable to turbulence, inspired at first by the work of the French school, above all for the basic principles, then the work of the foreigners, above all for the theory of the spectrum.

MECHANICS, STATISTICAL

Statistical computation of tolerance limits

Based on a new theory, two computer codes were developed specifically to calculate the exact statistical tolerance limits for normal distributions within unknown means and variances for the one-sided and two-sided cases for the tolerance factor, k. The quantity k is defined equivalently in terms of the noncentral t-distribution by the probability equation. Two of the four mathematical methods employ the theory developed for the numerical simulation. Several algorithms for numerically integrating and iteratively root-solving the working equations are written to augment the program simulation. The program codes generate some tables of k's associated with the varying values of the proportion and sample size for each given probability to show accuracy obtained for small sample sizes.

Wheeler, J. T.

Boundary layer receptivity and control

Receptivity processes initiate natural instabilities in a boundary layer. The instabilities grow and eventually break down to turbulence. Consequently, receptivity questions are a critical element of the analysis of the transition process. Success in modeling the physics of receptivity processes thus has a direct bearing on technological issues of drag reduction. The means by which transitional flows can be controlled is also a major concern: questions of control are tied inevitably to those of receptivity. Adjoint systems provide a highly effective mathematical method for approaching many of the questions associated with both receptivity and control. The long term objective is to develop adjoint methods to handle increasingly complex receptivity questions, and to find systematic procedures for deducing effective control strategies. The most elementary receptivity problem is that in which a parallel boundary layer is forced by time-harmonic sources of various types. The characteristics of the response to such forcing form the building blocks for more complex receptivity mechanisms. The first objective of this year's research effort was to investigate how a parallel Blasius boundary layer responds to general direct forcing. Acoustic disturbances in the freestream can be scattered by flow non-uniformities to produce Tollmien-Schlichting waves. For example, scattering by surface roughness is known to provide an efficient receptivity path. The present effort is directed towards finding a solution by a simple adjoint analysis, because adjoint methods can be extended to more complex problems. In practice, flows are non-parallel and often three-dimensional. Compressibility may also be significant in some cases. Recent developments in the use of Parabolized Stability Equations (PSE) offer a promising possibility. By formulating and solving a set of adjoint parabolized equations, a method for mapping the efficiency with which external forcing excites the three-dimensional motions of a non-parallel boundary layer was developed. The method makes use of the same computationally efficient formulation that makes the PSE currently so appealing. In the area of flow control, adjoint systems offer a powerful insight into the effect of control forces. One of the simplest control strategies for boundary layers involves the application of localized mean wall suction.

Hill, D. C.

Predictions Of Accretion Of Ice On A Commercial Airplane

Report presents computational study of accretion of ice on wings and on horizontal tail surfaces of simplified version of Boeing 737 airplane. Calculations performed by use of LEWICE3D computer code, which incorporates panel-method mathematical model of three-dimensional potential flow and mathematical models of trajectories of drops of water entrained in flow, impingement of drops on surfaces, and transfer of heat.

Bidwell, C. S.

Dual-Use Space Technology Transfer Conference and Exhibition, Volume 1

This document contains papers presented at the Dual-Use Space Technology Transfer Conference and Exhibition held at the Johnson Space Center February 1-3, 1994. Possible technology transfers covered during the conference were in the areas of information access; innovative microwave and optical applications; materials and structures; marketing and barriers; intelligent systems; human factors and habitation; communications and data systems; business process and technology transfer; software engineering; biotechnology and advanced bioinstrumentation; communications signal processing and analysis; new ways of doing business; medical care; applications derived from control center data systems; human performance evaluation; technology transfer methods; mathematics, modeling, and simulation; propulsion; software analysis and decision tools systems/processes in human support technology; networks, control centers, and distributed systems; power; rapid development perception and vision technologies; integrated vehicle health management; automation technologies; advanced avionics; ans robotics technologies. More than 77 papers, 20 presentations, and 20 exhibits covering various disciplines were presented b experts from NASA, universities, and industry.

Kumar Krishen

Dual-Use Space Technology Transfer Conference and Exhibition

This is the second volume of papers presented at the Dual-Use Space Technology Transfer Conference and Exhibition held at the Johnson Space Center February 1-3, 1994. Possible technology transfers covered during the conference were in the areas of information access; innovative microwave and optical applications; materials and structures; marketing and barriers; intelligent systems; human factors and habitation; communications and data systems; business process and technology transfer; software engineering; biotechnology and advanced bioinstrumentation; communications signal processing and analysis; medical care; applications derived from control center data systems; human performance evaluation; technology transfer methods; mathematics, modeling, and simulation; propulsion; software analysis and decision tools; systems/processes in human support technology; networks, control centers, and distributed systems; power; rapid development; perception and vision technologies; integrated vehicle health management; automation technologies; advanced avionics; and robotics technologies.

Krishen, Kumar

Progress in The Semantic Analysis of Scientific Code

This paper concerns a procedure that analyzes aspects of the meaning or semantics of scientific and engineering code. This procedure involves taking a user's existing code, adding semantic declarations for some primitive variables, and parsing this annotated code using multiple, independent expert parsers. These semantic parsers encode domain knowledge and recognize formulae in different disciplines including physics, numerical methods, mathematics, and geometry. The parsers will automatically recognize and document some static, semantic concepts and help locate some program semantic errors. These techniques may apply to a wider range of scientific codes. If so, the techniques could reduce the time, risk, and effort required to develop and modify scientific codes.

Stewart, Mark

An Experiment in Scientific Program Understanding

This paper concerns a procedure that analyzes aspects of the meaning or semantics of scientific and engineering code. This procedure involves taking a user's existing code, adding semantic declarations for some primitive variables, and parsing this annotated code using multiple, independent expert parsers. These semantic parsers encode domain knowledge and recognize formulae in different disciplines including physics, numerical methods, mathematics, and geometry. The parsers will automatically recognize and document some static, semantic concepts and help locate some program semantic errors. Results are shown for three intensively studied codes and seven blind test cases; all test cases are state of the art scientific codes. These techniques may apply to a wider range of scientific codes. If so, the techniques could reduce the time, risk, and effort required to develop and modify scientific codes.

Stewart, Mark E. M.

OTIS 3.2 Software Released

Trajectory, mission, and vehicle engineers concern themselves with finding the best way for an object to get from one place to another. These engineers rely upon special software to assist them in this. For a number of years, many engineers have used the OTIS program for this assistance. With OTIS, an engineer can fully optimize trajectories for airplanes, launch vehicles like the space shuttle, interplanetary spacecraft, and orbital transfer vehicles. OTIS provides four modes of operation, with each mode providing successively stronger optimization capability. The most powerful mode uses a mathematical method called implicit integration to solve what engineers and mathematicians call the optimal control problem. OTIS 3.2, which was developed at the NASA Glenn Research Center, is the latest release of this industry workhorse and features new capabilities for parameter optimization and mission design. OTIS stands for Optimal Control by Implicit Simulation, and it is implicit integration that makes OTIS so powerful at solving trajectory optimization problems. Why is this so important? The optimization process not only determines how to get from point A to point B, but it can also determine how to do this with the least amount of propellant, with the lightest starting weight, or in the fastest time possible while avoiding certain obstacles along the way. There are numerous conditions that engineers can use to define optimal, or best. OTIS provides a framework for defining the starting and ending points of the trajectory (point A and point B), the constraints on the trajectory (requirements like "avoid these regions where obstacles occur"), and what is being optimized (e.g., minimize propellant). The implicit integration method can find solutions to very complicated problems when there is not a lot of information available about what the optimal trajectory might be. The method was first developed for solving two-point boundary value problems and was adapted for use in OTIS. Implicit integration usually allows OTIS to find solutions to problems much faster than programs that use explicit integration and parametric methods. Consequently, OTIS is best suited to solving very complicated and highly constrained problems.

Riehl, John P.

The Impact of ARM on Climate Modeling

Climate models are among humanity's most ambitious and elaborate creations. They are designed to simulate the interactions of the atmosphere, ocean, land surface, and cryosphere on time scales far beyond the limits of deterministic predictability, and including the effects of time-dependent external forcings. The processes involved include radiative transfer, fluid dynamics, microphysics, and some aspects of geochemistry, biology, and ecology. The models explicitly simulate processes on spatial scales ranging from the circumference of the Earth down to one hundred kilometers or smaller, and implicitly include the effects of processes on even smaller scales down to a micron or so. The atmospheric component of a climate model can be called an atmospheric global circulation model (AGCM). In an AGCM, calculations are done on a three-dimensional grid, which in some of today's climate models consists of several million grid cells. For each grid cell, about a dozen variables are time-stepped as the model integrates forward from its initial conditions. These so-called prognostic variables have special importance because they are the only things that a model remembers from one time step to the next; everything else is recreated on each time step by starting from the prognostic variables and the boundary conditions. The prognostic variables typically include information about the mass of dry air, the temperature, the wind components, water vapor, various condensed-water species, and at least a few chemical species such as ozone. A good way to understand how climate models work is to consider the lengthy and complex process used to develop one. Lets imagine that a new AGCM is to be created, starting from a blank piece of paper. The model may be intended for a particular class of applications, e.g., high-resolution simulations on time scales of a few decades. Before a single line of code is written, the conceptual foundation of the model must be designed through a creative envisioning that starts from the intended application and is based on current understanding of how the atmosphere works and the inventory of mathematical methods available.

atmospheric circulation

Cryovolcanic Emplacement of Domes on Europa

Here we explore the hypothesis that certain domes on Europa may have been produced by the extrusion of viscous cryolavas. A new mathematical method for the emplacement and relaxation of viscous lava domes is presented and applied to putative cryovolcanic domes on Europa. A similarity solution approach is applied to the governing equation for fluid flow in a cylindrical geometry, and dome relaxation is explored assuming a volume of cryolava has been rapidly emplaced onto the surface. Nonphysical sin- gularities inherent in previous models for dome relaxation have been eliminated, and cryolava cooling is represented by a time-variable viscosity. We find that at the onset of relaxation, bulk kinematic viscosities may lie in the range between 10(exp 3) and 10(exp 6) sq m/s, while the actual fluid lava viscosity may be much lower. Plausible relaxation times to form the domes, which are linked to bulk cryolava rheology, are found to range from 3.6 days to 7.5 years. We find that cooling of the cryolava, while dominated by conduction through an icy skin, should not prevent fluids from advancing and relaxing to form domes within the timescales considered. Determining the range of emplacement conditions for putative cryolava domes will shed light on Europa's resurfacing history. In addition, the rheologies and compositions of erupted cryolavas have implications for subsurface cryomagma ascent and local surface stress conditions on Europa.

volcanism

Provably Correct Floating-Point Implementation of a Point-In-Polygon Algorithm

The problem of determining whether or not a point lies inside a given polygon occurs in many applications. In air traffic management concepts, a correct solution to the point-in-polygon problem is critical to geofencing systems for Unmanned Aerial Vehicles and in weather avoidance applications. Many mathematical methods can be used to solve the point-in-polygon problem. Unfortunately, a straightforward floating- point implementation of these methods can lead to incorrect results due to round-off errors. In particular, these errors may cause the control flow of the program to diverge with respect to the ideal real-number algorithm. This divergence potentially results in an incorrect point-in- polygon determination even when the point is far from the edges of the polygon. This paper presents a provably correct implementation of a point-in-polygon method that is based on the computation of the winding number. This implementation is mechanically generated from a source- to-source transformation of the ideal real-number specification of the algorithm. The correctness of this implementation is formally verified within the Frama-C analyzer, where the proof obligations are discharged using the Prototype Verification System (PVS).

Moscato, Mariano M.

Decomposition of Electron Ionization Mass Spectra for Space Application Using a Monte-Carlo approach

Rationale: Quadrupole mass spectrometers equipped with an electron ionization (EI) sources have been widely used in space exploration to investigate the composition of planetary surfaces and atmospheres. However, the complexity of the samples and the minimal calibration for the fragmentation of molecules in the ionization chambers have prevented the deconvolution of the majority of the mass spectra obtained at different targets, thus limiting the determination of the exact composition of the samples analyzed. We propose a Monte‐Carlo approach to solve this issue mathematically. Methods: We decomposed simulated mass spectra of mixtures acquired with unit resolving power mass spectrometers and EI sources into the sum of the single components fragmentation patterns weighted by their relative concentration using interior‐point least‐square fitting. To fit compounds with poorly known fragmentation patterns, we used a Monte‐Carlo method to vary the intensity of individual fragment ions. We then decomposed the spectrum thousands of times to obtain a statistical distribution. Results: By performing the deconvolution on a mixture of seven different molecules with interfering fragmentation patterns (H2O, O2, CH4, Ar, N2, C2H4, and C2H6) we show that this approach retrieves the mixing ratio of the individual components more accurately than regular mass spectra decomposition methods that rely on fragmentation patterns from general databases. It also provides the probability density function for each species's mixing ratio. Conclusions: By removing the solution degeneracy in the decomposition of mass spectra, the method described herein could significantly increase the scientific retrieval from archived space flight mass spectrometry data, where calibration of the ionization source is no longer an option.

Thomas Gautier