Search NASA⌕ Search

SEARCH · Search NASA

Results for “MATHEMATICAL SOLUTIONS”

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 253 records · Page 14

A critical discussion of theories of flame spread across solid and liquid fuels.

Mathematical descriptions of flames spreading over liquid and solid fuels are obtained, using basic assumptions derived from observations or physical reasoning. A review of existing theories shows that they are incomplete in that they either treat an uncoupled problem of the condensed phase where the spreading rate and heat flux at the surface are given, or they merely determine the spreading rate in terms of a new vaguely defined eigenvalue. An important difference between the liquid and solid cases, due to convection, is pointed out, and it is shown that solid-fuel flame-spread theories which claim to apply to the liquid case, in reality do not apply to it. A mathematical formulation and a method of solution are presented for the phenomenon of flame spread over solid fuels with forward heat conduction in both the solid and the gas. The method uses an energy integral over the field to determine the spreading rate in terms of the basic properties of the fuel and air.

Sirignano, W. A.↗

Electrophoretic separator for purifying biologicals

Mathematical expressions were developed to describe the interrelationships between operating requirements (capabilities), cell parameters, and system constraints in terms of design criteria definition. The mathematical model was programmed for computer solution. The model was exercised to identify performance-limiting characteristics, and analyses were conducted to predict operation in space of an experiment involving separation of four components. An engineering model of a flowing electrophoretic separator was constructed. The design is directed toward verifying improvements in resolution and throughput of a thicker cell than can be used on earth.

Source record↗

Ray trajectories in a torus: An application of MACSYMA to a complex numerical computation

The study of ray trajectories of plasma waves in a torodial geometry using MACSYMA is an example of how symbolic, numerical, and graphical facilities can be used in concert to accomplish a complex computational goal. Computational features of this study which are of particular significance include: the derivation of code (i.e. writing functions to generate program fragments), the use of array functions to simplify the specification of a numerical iteration scheme, and the graphical presentation of the results. Mathematically, this study originates in the solution of a linear inhomogeneous partial differential equation in 3 dimensions by the method of characteristics. While it is possible to describe this equation compactly by using vector notation, and by specifying the spatial variation of the coefficients in terms of intermediate parameters, the transformation of the equation into a form amenable to solution is very tedious. A MACSYMA program is presented for obtaining description of the rf field structure excited by a waveguide located at the edge of a toroidal plasma confinement device.

Kulp, J. L.↗

An optimality criterion for sizing members of heated structures with temperature constraints

A thermal optimality criterion is presented for sizing members of heated structures with multiple temperature constraints. The optimality criterion is similar to an existing optimality criterion for design of mechanically loaded structures with displacement constraints. Effectiveness of the thermal optimality criterion is assessed by applying it to one- and two-dimensional thermal problems where temperatures can be controlled by varying the material distribution in the structure. Results obtained from the optimality criterion agree within 2 percent with results from a closed-form solution and with results from a mathematical programming technique. The thermal optimality criterion augments existing optimality criteria for strength and stiffness related constraints and offers the possibility of extension of optimality techniques to sizing structures with combined thermal and mechanical loading.

Rao, G. V.↗

Thermal catastrophe in the plasma sheet boundary layer

This letter presents a first step towards a substorm model including particle heating and transport in the plasma sheet boundary layer (PSBL). The heating mechanism discussed is resonant absorption of Alfven waves. For some assumed MHD perturbation incident from the tail lobes onto the plasma sheet, the local heating rate in the PSBL has the form of a resonance function of the one-fluid plasma temperature. Balancing the local heating by convective transport of the heated plasma toward the central plasma sheet, an 'equation of state" is found for the steady-state PSBL whose solution has the form of a mathematical catastrophe: at a critical value of a parameter containing the incident power flux, the local density, and the convection velocity, the equilibrium temperature jumps discontinuously. Associating this temperature increase with the abrupt onset of the substorm expansion phase, the catastrophe model indicates at least three ways in which the onset may be triggered. Several other consequences related to substorm dynamics are suggested by the simple catastrophe model.

Smith, Robert A.↗

Recursive forward dynamics for multiple robot arms moving a common task object

Recursive forward dynamics algorithms are developed for an arbitrary number of robot arms moving a commonly held object. The multiarm forward dynamics problem is to find the angular accelerations at the joints and the contact forces that the arms impart to the task object. The problem also involves finding the acceleration of this object. The multiarm forward dynamics solutions provide a thorough physical and mathematical understanding of the way several arms behave in response to a set of applied joint moments. Such an understanding simplifies and guides the subsequent control design and experimentation process. The forward dynamics algorithms also provide the necessary analytical foundation for conducting analysis and simulation studies. The multiarm algorithms are based on the filtering and smoothing approach recently advanced for single-arm dynamics, and they can be built up modularly from the single-arm algorithms. The algorithms compute recursively the joint-angle accelerations, the contact forces, and the task-object accelerations. Algorithms are also developed to evaluate in closed form the linear transformations from the active joint moments to the joint-angle accelerations, to the task-object accelerations., and to the task-object contact forces. A possible computing architecture is presented as a precursor to a more complete investigation of the computational performance of the dynamics algorithms.

Rodriguez, G.↗

The thermal catastrophe model of substorms

Equations are presented for the physical mechanisms involved in the resonant absorption of Alfven waves in the plasma sheet boundary layer. It is shown that energy absorbed by the plasma-sheet particles is a function of the central plasma sheet temperature. The heating curve, when coupled with convective transport, yielded an equation of state for the steady state plasma sheet, whose solution has the form of a mathematical catastrophe. The master equation includes dynamic terms describing the transition across the thermal catastrophe, making it possible to evaluate the time scale for the catastrophe to occur.

Goertz, C. K.↗

Recursive forward dynamics for multiple robot arms moving a common task object

Recursive forward dynamics algorithms are developed and presented for an arbitrary number of robot arms moving a commonly held object. The multiarm forward dynamics problem is to find the angular accelerations at the joints and the contact forces that the arms impart to the task object. The problem also involves finding the acceleration of this object. The multiarm forward dynamics solutions provide a thorough physical and mathematical understanding of the way several arms behave in response to a set of applied joint moments. Such an understanding simplifies and guides the subsequent control design and experimentation process. The forward dynamics algorithms also provide the necessary analytical foundation for conducting analysis and simulation studies. The multiarm algorithms are based on the filtering and smoothing approach recently advanced for single-arm dynamics, and they can be built up modularly from the single-arm algorithms. The algorithms compute recursively the joint angle accelerations, the contact forces, and the task-object accelerations. Algorithms are also developed to evaluate in closed form the linear transformations from the active joint moments to the joint angle accelerations, to the task object accelerations, and to the task-object contact forces. A possible computing architecture is presented as a precursor to a more complete investigation of the computational performance of the dynamics algorithms.

Rodriguez, Guillermo↗

Simulations of cloudy, gaseous galactic disks

Aspects of the simulation of cloudy, gaseous galactic disks are discussed. The motivation for conducting such studies is addressed, and the formulation of the fundamental physical problem in mathematical terms is considered. The computational solution study of the formulated physical problem is examined, and the results and interpretation of computational studies are discussed.

Roberts, William W., Jr.↗

The Influence of Finite-size Sources in Acousto-ultrasonics

This work explores the effects that the finite normal axisymmetric traction loading of an infinite isotropic plate has on wave propagation in acousto-ultrasonics (AU), in which guided waves are created using two normal incidence transducers. Although the work also addresses the effects of the transducer pressure distribution and pulse shape, this thesis concentrates on two main questions: how does the transducer's diameter control the phase velocity and frequency spectrum of the response, and how does the plate thickness relate to the plate's excitability? The mathematics of the time-harmonic solution and the physical principles and the practical considerations for AU wave generation are explained. Transient sources are modeled by the linear superposition of the time-harmonic solutions found using the Hankel transform and they are then compared to experimental data to provide insight into the relation between the size of the transducer and the preferred phase velocity.

Pavlakovic, Brian N.↗

The Application of a Massively Parallel Computer to the Simulation of Electrical Wave Propagation Phenomena in the Heart Muscle Using Simplified Models

The simulation of heart arrhythmia and fibrillation are very important and challenging tasks. The solution of these problems using sophisticated mathematical models is beyond the capabilities of modern super computers. To overcome these difficulties it is proposed to break the whole simulation problem into two tightly coupled stages: generation of the action potential using sophisticated models. and propagation of the action potential using simplified models. The well known simplified models are compared and modified to bring the rate of depolarization and action potential duration restitution closer to reality. The modified method of lines is used to parallelize the computational process. The conditions for the appearance of 2D spiral waves after the application of a premature beat and the subsequent traveling of the spiral wave inside the simulated tissue are studied.

Karpoukhin, Mikhii G.↗

2D and 3D Method of Characteristic Tools for Complex Nozzle Development

This report details the development of a 2D and 3D Method of Characteristic (MOC) tool for the design of complex nozzle geometries. These tools are GUI driven and can be run on most Windows-based platforms. The report provides a user's manual for these tools as well as explains the mathematical algorithms used in the MOC solutions.

Rice, Tharen↗

Pulse Phase Thermography for Defect Detection and Visualization

Pulse Phase Thermography (PPT) has been reported as a novel powerful technique of the thermal NDE. It employs application of the Discrete Fourier Transform (DFT) to thermal images obtained following flash heating of the front surface of a specimen. The computed phasegrams are excellent for defect visualization in a wide range of materials. This is in part due to their low sensitivity to uneven heating. This work is an attempt to analyze advantages and limitations of PPT. Results of application of the DFT to simulated temperature decays are presented. The temperature evolution on a surface has been simulated based on an analytical solution of the one-dimensional heat diffusion problem. A more sophisticated study has been done for different sizes of defects using numerical solution of the three-dimensional mathematical model. Capabilities of PPT for in-depth scanning and for monitoring of the material loss are discussed. The recommendations for the practical application of the PPT are presented. Experimental results obtained following these recommendations are reported.

Marinetti, Sergio↗

A Study into Validating a Coupled Method of Characteristics and Direct Simulation Monte Carlo Method Against Empirical Data

The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.

Direct Simulation Monte Carlo Analysis Code↗

A Study into Validating A Coupled Method of Characteristics And Direct Simulation Monte Carlo Method Against Empirical Data

The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.

DAC↗