Search NASA⌕ Search

SEARCH · Search NASA

Results for “Applied Mathematics”

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 397 records · Page 22

Adiabatic invariants in stellar dynamics. 1: Basic concepts

The adiabatic criterion, widely used in astronomical dynamics, is based on the harmonic oscillator. It asserts that the change in action under a slowly varying perturbation is exponentially small. Recent mathematical results that precisely define the conditions for invariance show that this model does not apply in general. In particular, a slowly varying perturbation may cause significant evolution stellar dynamical systems even if its time scale is longer than any internal orbital time scale. This additional 'heating' may have serious implications for the evolution of star clusters and dwarf galaxies which are subject to long-term environmental forces. The mathematical developments leading to these results are reviewed, and the conditions for applicability to and further implications for stellar systems are discussed. Companion papers present a computational method for a general time-dependent disturbance and detailed example.

Weinberg, Martin D.↗

Applying Formal Methods to Safety-Critical Systems

How do you know a proof is correct? Traditionally, mathematical proofs are socially verified – at least one human, following a set of implicit rules of natural language and logic, determines if the proof is believable. If the proof becomes overly tedious and/or is essential to some safety- or mission-critical application, it becomes necessary to determine the soundness to a higher standard. 'Formal methods' refer to mathematically rigorous techniques and tools that enable specification, design, and verification of hardware and software systems. The specification used in formal methods are statements in a mathematical logic while the formal verifications are deductions in that logic. Formal methods can be difficult or time/resource intensive, but offer a higher level of assurance than standard verification through testing or handwritten proofs. This talk will introduce formal methods, motivated by applications of interest to NASA, including uncrewed aircraft operations in the national airspace, urban air environments, and wildfire areas. The audience will be given a crash course in mechanically verified proofs in the Prototype Verification System (PVS), an interactive theorem prover.

Formal Methods↗

Polynomial Chaos Surrogate Construction for Random Fields with Parametric Uncertainty

Engineering and applied science rely on computational experiments to rigorously study physical systems. The mathematical models used to probe these systems are highly complex, and sampling-intensive studies often require prohibitively many simulations for acceptable accuracy. Surrogate models provide a means of circumventing the high computational expense of sampling such complex models. In particular, polynomial chaos expansions (PCEs) have been successfully used for uncertainty quantification studies of deterministic models where the dominant source of uncertainty is parametric. We discuss an extension to conventional PCE surrogate modeling to enable surrogate construction for stochastic computational models that have intrinsic noise in addition to parametric uncertainty. We develop a PCE surrogate on a joint space of intrinsic and parametric uncertainty, enabled by Rosenblatt transformations, which are evaluated via kernel density estimation of the associated conditional cumulative distributions. Furthermore, we extend the construction to random field data via the Karhunen–Loève expansion. We then take advantage of closed-form solutions for computing PCE Sobol indices to perform a global sensitivity analysis of the model which quantifies the intrinsic noise contribution to the overall model output variance. Additionally, the resulting joint PCE is generative in the sense that it allows generating random realizations at any input parameter setting that are statistically approximately equivalent to realizations from the underlying stochastic model. The method is demonstrated on a chemical catalysis example model and a synthetic example controlled by a parameter that enables a switch from unimodal to bimodal response distributions.

97 MATHEMATICS AND COMPUTING↗

Turbulence modeling methods for the compressible Navier-Stokes equations

Turbulence modeling methods for the compressible Navier-Stokes equations, including several zero- and two-equation eddy-viscosity models, are described and applied. Advantages and disadvantages of the models are discussed with respect to mathematical simplicity, conformity with physical theory, and numerical compatibility with methods. A new two-equation model is introduced which shows advantages over other two-equation models with regard to numerical compatibility and the ability to predict low-Reynolds-number transitional phenomena. Calculations of various transonic airfoil flows are compared with experimental results. A new implicit upwind-differencing method is used which enhances numerical stability and accuracy, and leads to rapidly convergent steady-state solutions.

Coakley, T. J.↗

An Integrated Analysis of the Physiological Effects of Space Flight: Executive Summary

A large array of models were applied in a unified manner to solve problems in space flight physiology. Mathematical simulation was used as an alternative way of looking at physiological systems and maximizing the yield from previous space flight experiments. A medical data analysis system was created which consist of an automated data base, a computerized biostatistical and data analysis system, and a set of simulation models of physiological systems. Five basic models were employed: (1) a pulsatile cardiovascular model; (2) a respiratory model; (3) a thermoregulatory model; (4) a circulatory, fluid, and electrolyte balance model; and (5) an erythropoiesis regulatory model. Algorithms were provided to perform routine statistical tests, multivariate analysis, nonlinear regression analysis, and autocorrelation analysis. Special purpose programs were prepared for rank correlation, factor analysis, and the integration of the metabolic balance data.

Leonard, J. I.↗

Maximum dynamic responses using matched filter theory and random process theory

This paper describes and illustrates two ways of performing time-correlated gust-load calculations. The first is based on Matched Filter Theory; the second on Random Process Theory. The two yield theoretically identical results and both employ novel applications of the theories and unconventional interpretations of the intermediate and final results. Both approaches are computationally fast and are general enough to be applied to dynamic-response problems other than gust loads. A brief mathematical development and example calculations using both Matched Filter Theory and Random Process Theory are presented.

Pototzky, Anthony S.↗

Convection-driven tectonics on Venus

An analysis is presented of convective stress coupling to an elastic lithosphere as applied to Venus. Theoretical solutions are introduced for the response of a mathematically thick elastic plate overlying a Newtonian viscous medium with an exponential depth dependence of viscosity, and a Green's function solution is obtained for the viscous flow driven by a harmonic density distribution at a specified depth. An elastic-plastic analysis is carried out for the deformation of a model Venus lithosphere. The results predict that dynamic uplift of Venusian topography must be accompanied by extensive brittle failure and viscous flow in the lithosphere.

Phillips, Roger J.↗

Analysis of two-dimensional ducts with sudden area changes using the boundary element method

In this paper, the boundary element method (BEM) is applied to two-dimensional acoustic ducts with sudden area changes. The mathematical basis of BEM along with the technique of handling singularities is presented in a systematic manner. The computer code is verified by comparison with the exact solutions and a number of illustrative examples are included.

Pates, Carl S., III↗

Theoretical studies of the physics of the solar atmosphere

Significant advances in our theoretical basis for understanding several physical processes related to dynamical phenomena on the sun were achieved. We have advanced a new model for spicules and fibrils. We have provided a simple physical view of resonance absorption of MHD surface waves; this allowed an approximate mathematical procedure for obtaining a wealth of new analytical results which we applied to coronal heating and p-mode absorption at magnetic regions. We provided the first comprehensive models for the heating and acceleration of the transition region, corona, and solar wind. We provided a new view of viscosity under coronal conditions. We provided new insights into Alfven wave propagation in the solar atmosphere. And recently we have begun work in a new direction: parametric instabilities of Alfven waves.

Hollweg, Joseph V.↗

On the distribution of pitch angles in external galactic spirals NGC 1232 and NGC 5457

A numerical method, originally developed to analyze the morphology of global and local structure in prototype galaxies, is modified for analyzing observed disk-shape galaxies. Two digitized spiral galaxies NGC 1232 and NGC 5457 with varying degrees of contrast between arm and interarm regions are analyzed. A synergism of partitioning methods and a geometric mean least-squares regression algorithm serves to isolate local arm segments, spurs, feathers, and secondary features and to measure their pitch angles and lengths. The global arms are actually highly disjointed, with arm segments frequently revealing pitch angles between 30 and 50 deg, certainly greater than those of the parent arms. Prominent spurs tend to exhibit a much greater pitch angle. The automated mathematical algorithm is shown to have negligible numerical biasing and could be applied to any number of spiral galaxies manifesting flocculent structure, either prototype or observed, and could possibly be used as a tool for classification of multiple-armed-type galaxies.

Russell, William S.↗

Statistically generated weighted curve fit of residual functions for modal analysis of structures

A statistically generated weighting function for a second-order polynomial curve fit of residual functions has been developed. The residual flexibility test method, from which a residual function is generated, is a procedure for modal testing large structures in an external constraint-free environment to measure the effects of higher order modes and interface stiffness. This test method is applicable to structures with distinct degree-of-freedom interfaces to other system components. A theoretical residual function in the displacement/force domain has the characteristics of a relatively flat line in the lower frequencies and a slight upward curvature in the higher frequency range. In the test residual function, the above-mentioned characteristics can be seen in the data, but due to the present limitations in the modal parameter evaluation (natural frequencies and mode shapes) of test data, the residual function has regions of ragged data. A second order polynomial curve fit is required to obtain the residual flexibility term. A weighting function of the data is generated by examining the variances between neighboring data points. From a weighted second-order polynomial curve fit, an accurate residual flexibility value can be obtained. The residual flexibility value and free-free modes from testing are used to improve a mathematical model of the structure. The residual flexibility modal test method is applied to a straight beam with a trunnion appendage and a space shuttle payload pallet simulator.

Bookout, P. S.↗

Vibration Analysis of a Split Path Gearbox

Split path gearboxes can be attractive alternatives to the common planetary designs for rotorcraft, but because they have seen little use, they are relatively high risk designs. To help reduce the risk of fielding a rotorcraft with a split path gearbox, the vibration and dynamic characteristics of such a gearbox were studied. A mathematical model was developed by using the Lagrangian method, and it was applied to study the effect of three design variables on the natural frequencies and vibration energy of the gearbox. The first design variable, shaft angle, had little influence on the natural frequencies. The second variable, mesh phasing, had a strong effect on the levels of vibration energy, with phase angles of 0 deg and 180 deg producing low vibration levels. The third design variable, the stiffness of the shafts connecting the spur gears to the helical pinions, strongly influenced the natural frequencies of some of the vibration modes, including two of the dominant modes. We found that, to achieve the lowest level of vibration energy, the natural frequencies of these two dominant modes should be less than those of the main excitation sources.

Krantz, Timothy L.↗

Prediction of Fatigue Crack Growth Using Regularized Numerical Models

Though it is known in the engineering community that successful analyses rest upon the proper balance of (1) theoretical analysis of mathematical models, (2) physical experimentation and (3) computational simulation, this balance is currently handled in sometimes unwieldy and inefficient manner. It is proposed to investigate and develop rigorous and computationally efficient method to effectively combine all available information, from both experimental measurements and mathematical models, in the emulation of physical systems. This will be specifically applied to fatigue crack growth in metallic structures of interest to NASA.

Meade, Andrew J.↗

Investigation of Cloud Interior Structure Using a Tomographic Technique Based on Radon Transform. Application to the Measurements Made by the Research Scanning Polarimeter During CAMP2Ex

The Research Scanning Polarimeter (RSP) is an airborne along-track scanner measuring the polarized and total reflectances in 9 spectral channels. The RSP's uniquely high angular resolution coupled with the high frequency of measurements allows for characterization of liquid-water cloud droplet sizes using the rainbow structure observed in the polarized reflectance measurements. The RSP observations also provide for derivation of 2D field of extinction coefficient inside the cloud using a tomographic technique. This technique utilizes the family of cloud shapes derived using "cutout" technique and corresponding to a number of thresholds in total reflectance (separating bright cloud from its darker background). These shapes can be interpreted as level curves of an abstract "reflectance density", which is then used for derivation of the directional cloud optical thickness (dCOT) tomogram. The latter is a collection of dCOTs parameterized by the angles and offsets of the corresponding view rays (chords) relative to the cloud center. After this, the inverse Radon Transform (the mathematical basis of the X-ray com-puter tomography) is applied to the dCOT tomogram yielding 2D spatial distribution of the extinction coefficient. The later can be converted into droplet number concentration using the droplet size profiles derived from the RSP's polarized reflectance measurements. After tests on synthetic data this technique was applied to real RSP measurements from NASA's Cloud, Aerosol and Monsoon Processes Philippines Experiment (CAMP2Ex) conducted in the vicinity of the Philippines during the Southwest Monsoon (August–September 2019). We have investigated the interiors of a number of clouds observed during CAMP2Ex focusing on Cu and CuCg (Tcu) cases. Our retrievals were routinely compared with the correlative measurements by lidar (HSRL-2) and cloud radar (APR-3) that were deployed on the same airborne platform (NASA P-3B) during this field experiment.

Research Scanning Polarimeter↗

Robust Containment Queries over Collections of Trimmed NURBS Surfaces via Generalized Winding Numbers

Here, we propose a containment query that is robust to the watertightness of regions bound by trimmed NURBS surfaces, as this property is difficult to guarantee for in-the-wild CAD models. Containment is determined through the generalized winding number (GWN), a mathematical construction that is indifferent to the arrangement of surfaces in the shape. Applying contemporary techniques for the 3D GWN to trimmed NURBS surfaces requires some form of geometric discretization, introducing computational inefficiency to the algorithm and even risking containment misclassifications near the surface. In contrast, our proposed method leverages properties of the 3D solid angle to solve the relevant surface integral using a boundary formulation with rapidly converging adaptive quadrature. Batches of queries are further accelerated by memoizing (i.e., caching and reusing) quadrature node positions and tangents as they are evaluated. We demonstrate that our GWN method is robust to complex trimming geometry in a CAD model, and is accurate up to arbitrary precision at arbitrary distances from the surface. The derived containment query is therefore robust to model non-watertightness while respecting all curved features of the input shape.

97 MATHEMATICS AND COMPUTING↗

Application of Monte Carlo techniques to optimization of high-energy beam transport in a stochastic environment

An algorithm employing a modified sequential random perturbation, or creeping random search, was applied to the problem of optimizing the parameters of a high-energy beam transport system. The stochastic solution of the mathematical model for first-order magnetic-field expansion allows the inclusion of state-variable constraints, and the inclusion of parameter constraints allowed by the method of algorithm application eliminates the possibility of infeasible solutions. The mathematical model and the algorithm were programmed for a real-time simulation facility; thus, two important features are provided to the beam designer: (1) a strong degree of man-machine communication (even to the extent of bypassing the algorithm and applying analog-matching techniques), and (2) extensive graphics for displaying information concerning both algorithm operation and transport-system behavior. Chromatic aberration was also included in the mathematical model and in the optimization process. Results presented show this method as yielding better solutions (in terms of resolutions) to the particular problem than those of a standard analog program as well as demonstrating flexibility, in terms of elements, constraints, and chromatic aberration, allowed by user interaction with both the algorithm and the stochastic model. Example of slit usage and a limited comparison of predicted results and actual results obtained with a 600 MeV cyclotron are given.

Parrish, R. V.↗

Mathematical modeling of the aerodynamics of high-angle-of-attack maneuvers

This paper is a review of the current state of aerodynamic mathematical modeling for aircraft motions at high angles of attack. The mathematical model serves to define a set of characteristic motions from whose known aerodynamic responses the aerodynamic response to an arbitrary high angle-of-attack flight maneuver can be predicted. Means are explored of obtaining stability parameter information in terms of the characteristic motions, whether by wind-tunnel experiments, computational methods, or by parameter-identification methods applied to flight-test data. A rationale is presented for selecting and verifying the aerodynamic mathematical model at the lowest necessary level of complexity. Experimental results describing the wing-rock phenomenon are shown to be accommodated within the most recent mathematical model by admitting the existence of aerodynamic hysteresis in the steady-state variation of the rolling moment with roll angle. Interpretation of the experimental results in terms of bifurcation theory reveals the general conditions under which aerodynamic hysteresis must exist.

Schiff, L. B.↗