SEARCH · Engineering Papers

Results for “numerical 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

NASA NTRS1996

Buoyancy-Driven Heat Transfer During Application of a Thermal Gradient for the Study of Vapor Deposition at Low Pressure Using and Ideal Gas

A mathematical model has been developed to determine heat transfer during vapor deposition of source materials under a variety of orientations relative to gravitational accelerations. The model demonstrates that convection can occur at total pressures as low as 10-2 mm Hg. Through numerical computation, using physical material parameters of air, a series of time steps demonstrates the development of flow and temperature profiles during the course of vapor deposition. These computations show that in unit gravity vapor deposition occurs by transport through a fairly complicated circulating flow pattern when applying heat to the bottom of the vessel with parallel orientation with respect to the gravity vector. The model material parameters for air predict the effect of kinematic viscosity to be of the same order as thermal diffusivity, which is the case for Prandtl number approx. 1 fluids. Qualitative agreement between experiment and the model indicates that 6-(2-methyl-4-nitroanilino)-2,4-hexadiyn-l-ol (DAMNA) at these pressures indeed approximates an ideal gas at the experiment temperatures, and may validate the use of air physical constants. It is apparent that complicated nonuniform temperature distribution in the vapor could dramatically affect the homogeneity, orientation, and quality of deposited films. The experimental test i's a qualitative comparison of film thickness using ultraviolet-visible spectroscopy on films generated in appropriately oriented vapor deposition cells. In the case where heating of the reaction vessel occurs from the top, deposition of vapor does not normally occur by convection due to a stable stratified medium. When vapor deposition occurs in vessels heated at the bottom, but oriented relative to the gravity vector between these two extremes, horizontal thermal gradients induce a complex flow pattern. In the plane parallel to the tilt axis, the flow pattern is symmetrical and opposite in direction from that where the vessel is positioned vertically. The ground-based experiments are sufficient preliminary tests of theory and should be of significant interest regarding vapor deposited films in microgravity.

Frazier, D. O.↗

NASA NTRS1996

Preliminary Work in Obtaining Site-Directed Mutants of Hen Egg White Lysozyme

Protein crystal growth studies are recognized as a critical endeavor in the field of molecular biotechnology. The scientific applications of this field include the understanding of how enzymes function and the accumulation of accurate information of atomic structures, a key factor in the process of rational drug design. NASA has committed substantial investment and resources to the field of protein crystal growth and has conducted many microgravity protein crystal growth experiments aboard shuttle flights. Crystals grown in space tend to be larger, denser and have a more perfect habit and geometry. These improved properties gained in the microgravity environment of space result largely from the reduction of solutal convection, and the elimination of sedimentation at the growing crystal surface. Shuttle experiments have yielded many large, high quality crystals that are suitable for high resolution X-ray diffraction analysis. Examples of biologically important macromolecules which have been successfully crystallized during shuttle missions include: lysozyme, isocitrate lyase, gamma-interferon, insulin, human serum albumin and canavalin. Numerous other examples are also available. In addition to obtaining high quality crystals, investigators are also interested in learning the mechanisms by which the growth events take place. Crystallization experiments indicate that for the enzyme HEWL, measured growth rates do not follow mathematical models for 2D nucleation and dislocation-led growth of tetragonal protein crystals. As has been suggested by the laboratory of Marc L. Pusey, a possible explanation for the disagreement between observation and data is that HEWL tetraconal crystals form by aggregated units of lysozyme in supersaturated solutions. Surface measurement data was shown to fit very well with a model using an octamer unit cell as the growth unit. According to this model, the aggregation pathway and subsequent crystal growth is described by: monomer < ------ > dimer < ------- > tetramer < ------ > octamer < ------ > higher order. It is believed that multimer aggregation of lysozyme occurs by interaction at specific binding sites on the surface of the protein crystals. If the presence of discrete binding sites and the aggregation hypothesis is true, then it follows that the alteration of the binding site(s) should have significant effect on the measurements obtained during growth experiments. Site-directed mutagenesis allows the specific alteration of proteins by replacement, deletion or addition of specific amino acid residues. This report outlines the approach for this strategy and the progress made thus far toward that end.

Holmes, Leonard D.↗

NASA NTRS2003

Numerical Modeling of Nanoelectronic Devices

Nanoelectronic Modeling 3-D (NEMO 3-D) is a computer program for numerical modeling of the electronic structure properties of a semiconductor device that is embodied in a crystal containing as many as 16 million atoms in an arbitrary configuration and that has overall dimensions of the order of tens of nanometers. The underlying mathematical model represents the quantummechanical behavior of the device resolved to the atomistic level of granularity. The system of electrons in the device is represented by a sparse Hamiltonian matrix that contains hundreds of millions of terms. NEMO 3-D solves the matrix equation on a Beowulf-class cluster computer, by use of a parallel-processing matrix vector multiplication algorithm coupled to a Lanczos and/or Rayleigh-Ritz algorithm that solves for eigenvalues. In a recent update of NEMO 3-D, a new strain treatment, parameterized for bulk material properties of GaAs and InAs, was developed for two tight-binding submodels. The utility of the NEMO 3-D was demonstrated in an atomistic analysis of the effects of disorder in alloys and, in particular, in bulk In(x)Ga(l-x)As and in In0.6Ga0.4As quantum dots.

Klimeck, Gerhard↗

NASA NTRS1976

A mathematical analysis of the theory of interplanetary scintillation in the weak scattering approximation

A simplified analytical technique is presented for modeling the interplanetary scintillation of radio sources of finite angular size with a power-law electron-density-fluctuation power spectrum. The simplification results from the representation of the scintillation spectrum in confluent hypergeometric functions. The approximations presented allow fast numerical evaluation of a spectrum for a weakly scattering but extended medium with less than 10% error over the entire spectrum. Parameters describing anisotropic electron irregularities as well as anisotropic source structure are included, and the dependence of the spectrum normalization on the scales of the medium is derived explicitly. The parametric description of the domains of convergence of the approximate expansions also provides a simple conceptualization of the relative contributions of the scattered radiation along the line of sight to the observed spectrum. This is particularly useful for sources of finite angular size. This technique is applied to previously published observations.

Mitchell, D. G.↗

NASA NTRS1983

Development of a distributed-parameter mathematical model for simulation of cryogenic wind tunnels

A one-dimensional distributed-parameter dynamic model of a cryogenic wind tunnel was developed which accounts for internal and external heat transfer, viscous momentum losses, and slotted-test-section dynamics. Boundary conditions imposed by liquid-nitrogen injection, gas venting, and the tunnel fan were included. A time-dependent numerical solution to the resultant set of partial differential equations was obtained on a CDC CYBER 203 vector-processing digital computer at a usable computational rate. Preliminary computational studies were performed by using parameters of the Langley 0.3-Meter Transonic Cryogenic Tunnel. Studies were performed by using parameters from the National Transonic Facility (NTF). The NTF wind-tunnel model was used in the design of control loops for Mach number, total temperature, and total pressure and for determining interactions between the control loops. It was employed in the application of optimal linear-regulator theory and eigenvalue-placement techniques to develop Mach number control laws.

Tripp, J. S.↗

NASA NTRS2000

SSME Bearing and Seal Tester Data Compilation, Analysis and Reporting; and Refinement of the Cryogenic Bearing Analysis Mathematical Model

SRS and NASA/MSFC have developed software with unique capabilities to couple bearing kinematic modeling with high fidelity thermal modeling. The core thermomechanical modeling software was developed by SRS and others in the late 1980's and early 1990's under various different contractual efforts. SRS originally developed software that enabled SHABERTH (Shaft Bearing Thermal Model) and SINDA (Systems Improved Numerical Differencing Analyzer) to exchange data and autonomously allowing bearing component temperature effects to propagate into the steady state bearing mechanical model. A separate contract was issued in 1990 to create a personal computer version of the software. At that time SRS performed major improvements to the code. Both SHABERTH and SINDA were independently ported to the PC and compiled. SRS them integrated the two programs into a single program that was named SINSHA. This was a major code improvement.

Moore, James↗

NASA NTRS2000

SSME Bearing and Seal Tester Data Compilation, Analysis, and Reporting; and Refinement of the Cryogenic Bearing Analysis Mathematical Model

Moore, James↗

NASA NTRS1977

Design of experiment for earth rotation and baseline parameter determination from very long baseline interferometry

The possibility of recovering earth rotation and network geometry (baseline) parameters are emphasized. The numerical simulated experiments performed are set up in an environment where station coordinates vary with respect to inertial space according to a simulated earth rotation model similar to the actual but unknown rotation of the earth. The basic technique of VLBI and its mathematical model are presented. The parametrization of earth rotation chosen is described and the resulting model is linearized. A simple analysis of the geometry of the observations leads to some useful hints on achieving maximum sensitivity of the observations with respect to the parameters considered. The basic philosophy for the simulation of data and their analysis through standard least squares adjustment techniques is presented. A number of characteristic network designs based on present and candidate station locations are chosen. The results of the simulations for each design are presented together with a summary of the conclusions.

Dermanis, A.↗

NASA NTRS1977

Similarity laws of lunar and terrestrial volcanic flows

A mathematical model for the terrestrial and lunar volcanic flows is proposed. This mathematical model, which is one-dimensional, steady duct flow of a mixture of a gas and small solid particles (rock), has been analyzed in detail. The similarities and the differences of the essential features of lunar and terrestrial volcanic flows are determined. Numerical results for the equilibrium two-phase flows of lunar and terrestrial volcanics under similar conditions are presented. The main results of the theoretical model are: the lunar crater is much larger than the corresponding terrestrial crater; the exit velocity from the lunar volcanic flow may be higher than the lunar escape velocity but the exit velocity of terrestrial volcanic flow is much less than that of the linear case (this result confirms the hypothesis that Australian tektites came from the moon, as a stream of rock and gas of extremely high speed), and the thermal effects on the lunar volcanic flow are much larger than those in the terrestrial case.

Pai, S. I.↗

NASA NTRS1993

Mathematical Relationships Between Two Sets of Laser Anemometer Measurements for Resolving the Total Velocity Vector

The mathematical relations between the measured velocity fields for the same compressor rotor flow field resolved by two fringe type laser anemometers at different observational locations are developed in this report. The relations allow the two sets of velocity measurements to be combined to produce a total velocity vector field for the compressor rotor. This report presents the derivation of the mathematical relations, beginning with the specification of the coordinate systems and the velocity projections in those coordinate systems. The vector projections are then transformed into a common coordinate system. The transformed vector coordinates are then combined to determine the total velocity vector. A numerical example showing the solution procedure is included.

Owen, Albert K.↗

NASA NTRS1994

Acoustics of laminar boundary layers breakdown

Boundary layer flow transition has long been suggested as a potential noise source in both marine (sonar-dome self noise) and aeronautical (aircraft cabin noise) applications, owing to the highly transient nature of process. The design of effective noise control strategies relies upon a clear understanding of the source mechanisms associated with the unsteady flow dynamics during transition. Due to formidable mathematical difficulties, theoretical predictions either are limited to early linear and weakly nonlinear stages of transition, or employ acoustic analogy theories based on approximate source field data, often in the form of empirical correlation. In the present work, an approach which combines direct numerical simulation of the source field with the Lighthill acoustic analogy is utilized. This approach takes advantage of the recent advancement in computational capabilities to obtain detailed information about the flow-induced acoustic sources. The transitional boundary layer flow is computed by solving the incompressible Navier-Stokes equations without model assumptions, thus allowing a direct evaluation of the pseudosound as well as source functions, including the Lighthill stress tensor and the wall shear stress. The latter are used for calculating the radiated pressure field based on the Curle-Powell solution of the Lighthill equation. This procedure allows a quantitative assessment of noise source mechanisms and the associated radiation characteristics during transition from primary instability up to the laminar breakdown stage. In particular, one is interested in comparing the roles played by the fluctuating volume Reynolds stress and the wall-shear-stresses, and in identifying specific flow processes and structures that are effective noise generators.

Wang, Meng↗

NASA NTRS2005

Spinning Spacecraft Attitude Estimation Using Markley Variables: Filter Implementation And Results

Attitude estimation is often more difficult for spinning spacecraft than for three-axis stabilized platforms due to the need to follow rapidly-varying state vector elements and the lack of three-axis rate measurements from gyros. The estimation problem simplifies when torques are negligible and nutation has damped out, but the general case requires a sequential filter with dynamics propagation. This paper describes the implementation and test results for an extended Kalman filter for spinning spacecraft attitude and rate estimation based on a novel set of variables suggested in a paper by Markley [AAS93-3301 (referred to hereafter as Markley variables). Markley has demonstrated that the new set of variables provides a superior parameterization for numerical integration of the attitude dynamics for spinning or momentum-biased spacecraft. The advantage is that the Markley variables have fewer rapidly-varying elements than other representations such as the attitude quaternion and rate vector. A filter based on these variables was expected to show improved performance due to the more accurate numerical state propagation. However, for a variety of test cases, it has been found that the new filter, as currently implemented, does not perform significantly better than a quaternion-based filter that was developed and tested in parallel. This paper reviews the mathematical background for a filter based on Markley variables. It also describes some features of the implementation and presents test results. The test cases are based on a mission using magnetometer and Sun sensor data and gyro measurements on two axes normal to the spin axis. The orbit and attitude scenarios and spacecraft parameters are modeled after one of the THEMIS (Time History of Events and Macroscale Interactions during Substorms) probes. Several tests are presented that demonstrate the filter accuracy and convergence properties. The tests include torque-free motion with various nutation angles, large constant-torque attitude slews, sensor misalignments, large initial attitude and rate errors, and cases with low data frequency. It is found that the convergence is rapid, the radius of convergence is large, and the results are reasonably accurate even in the presence of unmodeled perturbations.

Sedlak, Joseph E.↗

NASA NTRS1987

Theoretical performance of hydrogen-bromine rechargeable SPE fuel cell

A mathematical model was formulated to describe the performance of a hydrogen-bromine fuel cell. Porous electrode theory was applied to the carbon felt flow-by electrode and was coupled to theory describing the solid polymer electrolyte (SPE) system. Parametric studies using the numerical solution to this model were performed to determine the effect of kinetic, mass transfer, and design parameters on the performance of the fuel cell. The results indicate that the cell performance is most sensitive to the transport properties of the SPE membrane. The model was also shown to be a useful tool for scale-up studies.

Savinell, Robert F.↗

NASA NTRS1988

Theoretical performance of hydrogen-bromine rechargeable SPE fuel cell

Savinell, R. F.↗

NASA NTRS1990

Identification of aerodynamic models for maneuvering aircraft

Due to the requirement of increased performance and maneuverability, the flight envelope of a modern fighter is frequently extended to the high angle-of-attack regime. Vehicles maneuvering in this regime are subjected to nonlinear aerodynamic loads. The nonlinearities are due mainly to three-dimensional separated flow and concentrated vortex flow that occur at large angles of attack. Accurate prediction of these nonlinear airloads is of great importance in the analysis of a vehicle's flight motion and in the design of its flight control system. A satisfactory evaluation of the performance envelope of the aircraft may require a large number of coupled computations, one for each change in initial conditions. To avoid the disadvantage of solving the coupled flow-field equations and aircraft's motion equations, an alternate approach is to use a mathematical modeling to describe the steady and unsteady aerodynamics for the aircraft equations of motion. Aerodynamic forces and moments acting on a rapidly maneuvering aircraft are, in general, nonlinear functions of motion variables, their time rate of change, and the history of maneuvering. A numerical method was developed to analyze the nonlinear and time-dependent aerodynamic response to establish the generalized indicial function in terms of motion variables and their time rates of change.

Chin, Suei↗

NASA NTRS1993

An analytic study of a two-phase laminar airfoil in simulated heavy rain

A mathematical model for a two-phase flow laminar airfoil in simulated heavy rain has been established. The set of non-linear partial differential equations has been converted into a set of finite difference equations; appropriate initial and boundary conditions are provided. The numerical results are compared with the experimental measurements. They show good agreement in quality.

Hsu, Yu-Kao↗

NASA NTRS2023

A Theoretical and Computational Revisit of Conversions Between Whitham’s F-Function and Equivalent Area

This paper proves mathematically that the integral transforms between Whitham’s F-function and equivalent area are the inverse transforms of each other if and only if the slope of the equivalent area at the origin is zero. This mathematical fact contradicts the accepted unconditional inverse relation between Whitham’s F-function and equivalent area in the sonic boom research literature. Piecewise linear approximations of an F-function and of the second derivative of an equivalent area are used to derive numerical formulas for conversions between Whitham’s F-function and equivalent area. Numerical results are included to show convergence of the numerical conversions as the maximum length of the segments for piecewise linear approximations goes to zero. These numerical conversions are approximately the inverse transforms of each other when the second derivative of an equivalent area is continuous and the slope of the equivalent area at the origin is zero.

Whitham's F-function↗

NASA NTRS2003

Numerical Simulation of Protoplanetary Vortices

The fluid dynamics within a protoplanetary disk has been attracting the attention of many researchers for a few decades. Previous works include, to list only a few among many others, the well-known prescription of Shakura & Sunyaev, the convective and instability study of Stone & Balbus and Hawley et al., the Rossby wave approach of Lovelace et al., as well as a recent work by Klahr & Bodenheimer, which attempted to identify turbulent flow within the disk. The disk is commonly understood to be a thin gas disk rotating around a central star with differential rotation (the Keplerian velocity), and the central quest remains as how the flow behavior deviates (albeit by a small amount) from a strong balance established between gravitational and centrifugal forces, transfers mass and momentum inward, and eventually forms planetesimals and planets. In earlier works we have briefly described the possible physical processes involved in the disk; we have proposed the existence of long-lasting, coherent vortices as an efficient agent for mass and momentum transport. In particular, Barranco et al. provided a general mathematical framework that is suitable for the asymptotic regime of the disk; Barranco & Marcus (2000) addressed a proposed vortex-dust interaction mechanism which might lead to planetesimal formation; and Lin et al. (2002), as inspired by general geophysical vortex dynamics, proposed basic mechanisms by which vortices can transport mass and angular momentum. The current work follows up on our previous effort. We shall focus on the detailed numerical implementation of our problem. We have developed a parallel, pseudo-spectral code to simulate the full three-dimensional vortex dynamics in a stably-stratified, differentially rotating frame, which represents the environment of the disk. Our simulation is validated with full diagnostics and comparisons, and we present our results on a family of three-dimensional, coherent equilibrium vortices.

Lin, H.↗