Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mathematical physics 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 181 records · Page 10

Recent advances in the solution of three-dimensional flow over wings with leading edge vortex separation

Recent advances in a panel method for the solution of three-dimensional flow about wing and wing-body combinations with leading-edge vortex separation are presented. These advances were achieved as part of an ultimately successful assault on two shortcomings of the method, namely convergence failures in seemingly random cases, and overprediction of lift coefficient for high aspect-ratio wings. Advances include the implementation of improved panel numerics for the purpose of eliminating the highly non-linear effects of ring vortices around doublet panel edges, and the development of a least squares procedure for damping vortex sheet geometry update instabilities. A variety of cases generated by the computer program implementing the method are presented. These cases are of two types. The first type consists of numerical studies, which verify the underlying mathematical assumptions of the method and moreover show that the results are strongly invariant with respect to such user dependent input as wing panel layout, initial sheet shape, sheet rollup, etc. The second type consists of cases run for the purpose of comparing computed results with experimental data, and these comparisons verify the underlying physical assumptions made by the method.

Johnson, F. T.↗

NASA Langley's Approach to the Sandia's Structural Dynamics Challenge Problem

The objective of this challenge is to develop a data-based probabilistic model of uncertainty to predict the behavior of subsystems (payloads) by themselves and while coupled to a primary (target) system. Although this type of analysis is routinely performed and representative of issues faced in real-world system design and integration, there are still several key technical challenges that must be addressed when analyzing uncertain interconnected systems. For example, one key technical challenge is related to the fact that there is limited data on target configurations. Moreover, it is typical to have multiple data sets from experiments conducted at the subsystem level, but often samples sizes are not sufficient to compute high confidence statistics. In this challenge problem additional constraints are placed as ground rules for the participants. One such rule is that mathematical models of the subsystem are limited to linear approximations of the nonlinear physics of the problem at hand. Also, participants are constrained to use these models and the multiple data sets to make predictions about the target system response under completely different input conditions. Our approach involved initially the screening of several different methods. Three of the ones considered are presented herein. The first one is based on the transformation of the modal data to an orthogonal space where the mean and covariance of the data are matched by the model. The other two approaches worked solutions in physical space where the uncertain parameter set is made of masses, stiffnesses and damping coefficients; one matches confidence intervals of low order moments of the statistics via optimization while the second one uses a Kernel density estimation approach. The paper will touch on all the approaches, lessons learned, validation 1 metrics and their comparison, data quantity restriction, and assumptions/limitations of each approach. Keywords: Probabilistic modeling, model validation, uncertainty quantification, kernel density

Horta, Lucas G.↗

Surface Oscillations of An Electromagnetically Levitated Droplet

Oscillating droplets are of interest in a number of disciplines. A practical application is the oscillating drop method, which is a technique for measuring surface tension and viscosity of liquid metals. It is especially suited to undercooled and highly reactive metals, because it is usually performed by electromagnetic levitation. The natural oscillation frequency a f the droplets is related to the surface tension of the material, and the decay of oscillations is related to its viscosity. The fluid flow inside the droplet must be laminar in order for this technique to yield good results. Because no experimental method has yet been developed to visualize flow in electromagnetically-levitated oscillating metal droplets, mathematical modeling is required to determine whether or not turbulence occurs. Three mathematical models of the flow: (1) assuming laminar conditions, (2) using the turbulence model, and (3) using the RNG turbulence model, respectively, are compared and contrasted to determine the physical characteristics of the flow. It is concluded that the RNG model is the best suited for describing this problem. The goal of the presented work was to characterize internal flow in an oscillating droplet of liquid metal, and to verify the accuracy of the characterization by comparing calculated surface tension and viscosity values to experimentally-obtained results.

Berry, S.↗

Introduction to COFFE: The Next-Generation HPCMP CREATE-AV CFD Solver

HPCMP CREATE-AV Conservative Field Finite Element (COFFE) is a modular, extensible, robust numerical solver for the Navier-Stokes equations that invokes modularity and extensibility from its first principles. COFFE implores a flexible, class-based hierarchy that provides a modular approach consisting of discretization, physics, parallelization, and linear algebra components. These components are developed with modern software engineering principles to ensure ease of uptake from a user's or developer's perspective. The Streamwise Upwind/Petrov-Galerkin (SU/PG) method is utilized to discretize the compressible Reynolds-Averaged Navier-Stokes (RANS) equations tightly coupled with a variety of turbulence models. The mathematics and the philosophy of the methodology that makes up COFFE are presented.

Glasby, Ryan S.↗

Plastic scintillator block as photon beam monitor for EGRET calibration

The EGRET (Energetic Gamma Ray Experiment Telescope) detector has been calibrated at SLAC (Stanford Linear Accelerator) and, to a lesser degree, at the MIT Bates Linear Accelerator Center. To monitor the photon beams for the calibration, a plastic scintillator block, 5 cm x 5 cm in cross section, 15 cm in length, and viewed by a single photomultiplier tube, was used for the entire beam energy range of 15 MeV to 10 GeV. The design operation, and method of analysis of the beam intensity are presented. A mathematical framework has been developed to treat the general case of a beam with multiphoton beam pulses and with a background component. A procedure to deal with the fluctuations of the beam intensity over a data-taking period was also developed. The photon beam monitor is physically sturdy, electronically steady, simple to construct, and easy to operate. Its major merits lie in its sheer simplicity of construction and operation and in the wide energy range it can cover.

Lin, Y. C.↗

Asymptotic solution for solar electric low thrust orbit raising with eclipse penalty

This paper describes an analytic solution for the solar electric low thrust orbit raising problem derived by application of the two-variable asymptotic expansion method. Included in the formulation are the effects of (1) high initial eccentricity, (2) large change of inclination, (3) variable mass, and (4) earth shadow penalty including the rotation of the shadow axis relative to the osculating orbit. An advantage of this mathematical strategy is that the effects of various steering policies can be evaluated readily from a physical point of view.

Flandro, G. A.↗

Mathematical Modeling of an Oscillating Droplet

Oscillating droplets are of interest in a number of disciplines. A practical application is the oscillating drop method, which is a technique for measuring surface tension and viscosity of liquid metals. It is especially suited to undercooled and highly reactive metals, because it is performed by electromagnetic levitation. The natural oscillation frequency of the droplets is related to the surface tension of the material, and the decay of oscillations is related to its viscosity. The fluid flow inside the droplet must be laminar in order for this technique to yield good results. Because no experimental method has yet been developed to visualize flow in electromagnetically-levitated oscillating metal droplets, mathematical modeling is required to determine whether or not turbulence occurs. Three mathematical models of the flow: (1) assuming laminar conditions, (2) using the k-epsilon turbulence model, and (3) using the RNG turbulence model, respectively, are compared and contrasted to determine the physical characteristics of the flow. It is concluded that the RNG model is the best suited for describing this problem. The goal of the presented work was to characterize internal flow in an oscillating droplet of liquid metal, and to verify the accuracy of the characterization by comparing calculated surface tension and viscosity.

Berry, S.↗

Cumulative reports and publications through December 31, 1991

A reports and publications list is given from the Institute for Computer Applications in Science and Engineering (ICASE) through December 31, 1991. The major categories of the current ICASE research program are; numerical methods, control and parameter identification problems, computational problems in engineering and the physical sciences, and computer systems and software. Since ICASE reports are intended to be preprints of articles that will appear in journals or conference proceedings, the published reference is included when available.

Mathematics↗

Heat pipe dynamic behavior

The vapor flow in a heat pipe was mathematically modeled and the equations governing the transient behavior of the core were solved numerically. The modeled vapor flow is transient, axisymmetric (or two-dimensional) compressible viscous flow in a closed chamber. The two methods of solution are described. The more promising method failed (a mixed Galerkin finite difference method) whereas a more common finite difference method was successful. Preliminary results are presented showing that multi-dimensional flows need to be treated. A model of the liquid phase of a high temperature heat pipe was developed. The model is intended to be coupled to a vapor phase model for the complete solution of the heat pipe problem. The mathematical equations are formulated consistent with physical processes while allowing a computationally efficient solution. The model simulates time dependent characteristics of concern to the liquid phase including input phase change, output heat fluxes, liquid temperatures, container temperatures, liquid velocities, and liquid pressure. Preliminary results were obtained for two heat pipe startup cases. The heat pipe studied used lithium as the working fluid and an annular wick configuration. Recommendations for implementation based on the results obtained are presented. Experimental studies were initiated using a rectangular heat pipe. Both twin beam laser holography and laser Doppler anemometry were investigated. Preliminary experiments were completed and results are reported.

Issacci, F.↗

Hidden Statistics of Schroedinger Equation

Work was carried out in determination of the mathematical origin of randomness in quantum mechanics and creating a hidden statistics of Schr dinger equation; i.e., to expose the transitional stochastic process as a "bridge" to the quantum world. The governing equations of hidden statistics would preserve such properties of quantum physics as superposition, entanglement, and direct-product decomposability while allowing one to measure its state variables using classical methods.

Zak, Michail↗

Early Estimation of Solar Activity Cycle: Potential Capability and Limits

The variable solar magnetic activity known as the 11-year solar cycle has the longest history of solar observations. These cycles dramatically affect conditions in the heliosphere and the Earth's space environment. Our current understanding of the physical processes that make up global solar dynamics and the dynamo that generates the magnetic fields is sketchy, resulting in unrealistic descriptions in theoretical and numerical models of the solar cycles. The absence of long-term observations of solar interior dynamics and photospheric magnetic fields hinders development of accurate dynamo models and their calibration. In such situations, mathematical data assimilation methods provide an optimal approach for combining the available observational data and their uncertainties with theoretical models in order to estimate the state of the solar dynamo and predict future cycles. In this presentation, we will discuss the implementation and performance of an Ensemble Kalman Filter data assimilation method based on the Parker migratory dynamo model, complemented by the equation of magnetic helicity conservation and longterm sunspot data series. This approach has allowed us to reproduce the general properties of solar cycles and has already demonstrated a good predictive capability for the current cycle, 24. We will discuss further development of this approach, which includes a more sophisticated dynamo model, synoptic magnetogram data, and employs the DART Data Assimilation Research Testbed.

Kitiashvili, Irina N.↗

Flight test planning and parameter extraction for rotorcraft system identification

The present study is concerned with the mathematical modelling of aircraft dynamics on the basis of an investigation conducted with the aid of the Rotor System Research Aircraft (RSRA). The particular characteristics of RSRA make it possible to investigate aircraft properties which cannot be readily studied elsewhere, for example in the wind tunnel. The considered experiment had mainly the objective to develop an improved understanding of the physics of rotor flapping dynamics and rotor loads in maneuvers. The employed approach is based on a utilization of parameter identification methodology (PID) with application to helicopters. A better understanding of the contribution of the main rotor to the overall aircraft forces and moments is also to be obtained. Attention is given to the mathematical model of a rotorcraft system, an integrated identification method, flight data processing, and the identification of RSRA mathematical models.

Wang, J. C.↗

Three-dimensional solidification and melting using magnetic field control

A new two-fluid mathematical model for fully three dimensional steady solidification under the influence of an arbitrary acceleration vector and with or without an arbitrary externally applied steady magnetic field have been formulated and integrated numerically. The model includes Joule heating and allows for separate temperature dependent physical properties within the melt and the solid. Latent heat of phase change during melting/solidification was incorporated using an enthalpy method. Mushy region was automatically captured by varying viscosity orders of magnitude between liquidus and solidus temperature. Computational results were obtained for silicon melt solidification in a parallelepiped container cooled from above and from a side. The results confirm that the magnetic field has a profound influence on the solidifying melt flow field thus changing convective heat transfer through the boundaries and the amount and shape of the solid accrued. This suggests that development of a quick-response algorithm for active control of three dimensional solidification is feasible since it would require low strength magnetic fields.

Dulikravich, George S.↗

Correlation of finite element free vibration predictions using random vibration test data

Finite element analysis is regularly used during the engineering cycle of mechanical systems to predict the response to static, thermal, and dynamic loads. The finite element model (FEM) used to represent the system is often correlated with physical test results to determine the validity of analytical results provided. Results from dynamic testing provide one means for performing this correlation. One of the most common methods of measuring accuracy is by classical modal testing, whereby vibratory mode shapes are compared to mode shapes provided by finite element analysis. The degree of correlation between the test and analytical mode shapes can be shown mathematically using the cross orthogonality check. A great deal of time and effort can be exhausted in generating the set of test acquired mode shapes needed for the cross orthogonality check. In most situations response data from vibration tests are digitally processed to generate the mode shapes from a combination of modal parameters, forcing functions, and recorded response data. An alternate method is proposed in which the same correlation of analytical and test acquired mode shapes can be achieved without conducting the modal survey. Instead a procedure is detailed in which a minimum of test information, specifically the acceleration response data from a random vibration test, is used to generate a set of equivalent local accelerations to be applied to the reduced analytical model at discrete points corresponding to the test measurement locations. The static solution of the analytical model then produces a set of deformations that once normalized can be used to represent the test acquired mode shapes in the cross orthogonality relation. The method proposed has been shown to provide accurate results for both a simple analytical model as well as a complex space flight structure.

Chambers, Jeffrey A.↗

Lightning Effects in the Payload Changeout Room

Analytical and empirical studies have been performed to provide better understanding of the electromagnetic environment inside the Payload Changeout Room and Orbiter payload bay resulting from lightning strikes to the launch pad lightning protection system. The analytical studies consisted of physical and mathematical modeling of the pad structure and the Payload Changeout Room. Empirical testing was performed using a lightning simulator to simulate controlled (8 kA) lightning strikes to the catenary wire lightning protection system. In addition to the analyses and testing listed above, an analysis of the configuration with the vehicle present was conducted, in lieu of testing, by the Finite Difference, Time Domain method.

Thomas, Garland L.↗

Unsteady Aerodynamic Modeling in Roll for the NASA Generic Transport Model

Reducing the impact of loss-of-control conditions on commercial transport aircraft is a primary goal of the NASA Aviation Safety Program. One aspect in developing the supporting technologies is to improve the aerodynamic models that represent these adverse conditions. Aerodynamic models appropriate for loss of control conditions require a more general mathematical representation to predict nonlinear unsteady behaviors. In this paper, a more general mathematical model is proposed for the subscale NASA Generic Transport Model (GTM) that covers both low and high angles of attack. Particular attention is devoted to the stall region where full-scale transports have demonstrated a tendency for roll instability. The complete aerodynamic model was estimated from dynamic wind-tunnel data. Advanced computational methods are used to improve understanding and visualize the flow physics within the region where roll instability is a factor.

Murphy, Patrick C.↗

On ionospheric aerodynamics

This paper presents theoretical methods to determine the gas dynamic and the electrostatic effects due to the interaction caused by a rapidly moving body in the ionosphere. The principles of the methods are derived from the kinetic theory of collision-free plasma. It is shown that the collective behavior of the collision-free plasma makes it possible to use the fluid approach to treat the problems of ionospheric aerodynamics. Various solutions to the system of fluid and field equations that have direct bearing on the ionospheric aerodynamics are presented and discussed. Physical significances of the mathematical results are stressed. Some outstanding unsolved problems in ionospheric aerodynamics are elaborated and discussed.

Liu, V. C.↗

Uncertainty of in-flight thrust determination

Methods for estimating the measurement error or uncertainty of in-flight thrust determination in aircraft employing conventional turbofan/turbojet engines are reviewed. While the term 'in-flight thrust determination' is used synonymously with 'in-flight thrust measurement', in-flight thrust is not directly measured but is determined or calculated using mathematical modeling relationships between in-flight thrust and various direct measurements of physical quantities. The in-flight thrust determination process incorporates both ground testing and flight testing. The present text is divided into the following categories: measurement uncertainty methodoogy and in-flight thrust measurent processes.

Abernethy, Robert B.↗