Engineering Papers⌕ Search

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.

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

At least 1,225 records · Page 68

A system identification approach for non-intrusive reduced order modeling of radiation-induced photocurrents

In this study, development of compact photocurrent models is currently dominated by analytical techniques that rely on physical assumptions to render the governing equations solvable in a closed form. Violation of these assumptions can reduce the accuracy of the models and/or limit their scope. In this paper we show that system identification of nonlinear state-space systems can serve as an alternative numerical basis for non-intrusive reduced order modeling of photocurrent effects. To that end we develop a compact gray box photocurrent model (GBPM) by using a state-space representation with a low-dimensional latent state equation that mimics a mathematical model for the response of an idealized class of devices to ionizing radiation. In so doing we obtain a model that learns the dynamics of a quantity of interest directly from its measurements without requiring snapshots of the internal device state or its discretized model, and can be inferred from very small data sets. To demonstrate the approach we train the GBPM using a small experimental data set for a Z5236 Zener diode and a small synthetic data set obtained by simulating a synthetic pn-junction device. We then compare the GBPMs with black box models trained on the same data and show that performance of the latter is limited by the size of the data set, while the former are able to achieve excellent performance in both the reproductive and the predictive regimes.

97 MATHEMATICS AND COMPUTING↗

Computer program for investigating effects of nonlinear suspension-system elastic properties on parachute inflation loads and motions

A computer program is presented by which the effects of nonlinear suspension-system elastic characteristics on parachute inflation loads and motions can be investigated. A mathematical elastic model of suspension-system geometry is coupled to the planar equations of motion of a general vehicle and canopy. Canopy geometry and aerodynamic drag characteristics and suspension-system elastic properties are tabular inputs. The equations of motion are numerically integrated by use of an equivalent fifth-order Runge-Kutta technique.

Poole, L. R.↗

Further analytical study of hybrid rocket combustion

Analytical studies of the transient and steady-state combustion processes in a hybrid rocket system are discussed. The particular system chosen consists of a gaseous oxidizer flowing within a tube of solid fuel, resulting in a heterogeneous combustion. Finite rate chemical kinetics with appropriate reaction mechanisms were incorporated in the model. A temperature dependent Arrhenius type fuel surface regression rate equation was chosen for the current study. The governing mathematical equations employed for the reacting gas phase and for the solid phase are the general, two-dimensional, time-dependent conservation equations in a cylindrical coordinate system. Keeping the simplifying assumptions to a minimum, these basic equations were programmed for numerical computation, using two implicit finite-difference schemes, the Lax-Wendroff scheme for the gas phase, and, the Crank-Nicolson scheme for the solid phase.

Hung, W. S. Y.↗

A nonlinear model for gas chromatograph systems

Fundamental engineering design techniques and concepts were studied for the optimization of a gas chromatograph-mass spectrometer chemical analysis system suitable for use on an unmanned, Martian roving vehicle. Previously developed mathematical models of the gas chromatograph are found to be inadequate for predicting peak heights and spreading for some experimental conditions and chemical systems. A modification to the existing equilibrium adsorption model is required; the Langmuir isotherm replaces the linear isotherm. The numerical technique of Crank-Nicolson was studied for use with the linear isotherm to determine the utility of the method. Modifications are made to the method eliminate unnecessary calculations which result in an overall reduction of the computation time of about 42 percent. The Langmuir isotherm is considered which takes into account the composition-dependent effects on the thermodynamic parameter, mRo.

Feinberg, M. P.↗

Scientific achievements from ten years of lunar laser ranging

In the 10 years since lunar laser ranging became a reality the need to analyze the observations has motivated improvements in several aspects of the mathematical model of earth-moon dynamics. Application of the data to improved estimates of the physical parameters of the earth-moon system has yielded significant astronomical, selenophysical, geophysical, and cosmological results. The scientific impact, both in improved theories and in numerical applications, is surveyed. The underlying physics and major difficulties are discussed, as well as the scientific results.

Mulholland, J. D.↗

Supercomputers Of The Future

Report evaluates supercomputer needs of five key disciplines: turbulence physics, aerodynamics, aerothermodynamics, chemistry, and mathematical modeling of human vision. Predicts these fields will require computer speed greater than 10(Sup 18) floating-point operations per second (FLOP's) and memory capacity greater than 10(Sup 15) words. Also, new parallel computer architectures and new structured numerical methods will make necessary speed and capacity available.

Peterson, Victor L.↗

Analyzing Impacts Of Meteoroids And Debris On Spacecraft

BUMPERII is modular computer program implementing numerical solution technique to calculate probability of no penetration (PNP) of spacecraft by orbiting artificial debris or impacts of meteoroids. Equation used to calculate PNP based on Poisson-distribution model applicable to similar analysis of smaller craft, but reflects more rigorous mathematical modeling of size, shape, orientation, and other characteristics of spacecraft that affect impacts, and necessary for PNP analysis of such larger structures as components of space station. Written in FORTRAN 77.

Bjorkman, Mike↗

The obliquity of Mars and 'climate friction'

A mathematical theory is presented which explains the increase of the mean obliquity of Mars over geologic time due to the so called 'climate friction' (i.e., the climatic changes associated with obliquity oscillations of Mars). The theory is compared with a 10 m.y. numerical integration of the equations performed for a hypothetically large amount of climate friction for two cases of the obliquity oscillations: (1) a single sinusoid and (2) a sum of three sinusoids. The theory and numerics agree for both cases within about 12 percent on the size of the secular increase in obliquity. One possible mechanism of climate friction investigated is 'postglacial rebound' on Mars. According to this theory, giant polar caps form when the obliquity is low, and slowly squeeze out an equatorial bulge. When the obliquity is high, the caps disappear, but the bulge takes some time to collapse, due to mantle viscosity, causing it to oscillate out of phase with the obliquity oscillations. This causes a secular increase in the average obliquity.

Rubincam, David P.↗

The Story of PI

The early history and the uses of the mathematical notation - pi - are presented through both film footage and computer animation in this 'Project Mathematics' series video. Pi comes from the first letter in the Greek word for perimeter. Archimedes, and early Greek mathematician, formulated the equations for the computation of a circle's area using pi and was the first person to seriously approximate pi numerically, although only to a few decimal places. By 1985, pi had been approximated to over one billion decimal places and was found to have no repeating pattern. One use of pi is the application of its approximation calculation as an analytical tool for determining the accuracy of supercomputers and software designs.

Apostol, Tom M.↗

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.↗

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.↗

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↗

From Machine Learning to Machine Reasoning: A Model-based Approach to Analyze Equipment Reliability Data

In current nuclear power plants (NPPs) a large amount of condition-based data which can be used to assess and monitor component health and performance. Assessing component health from such data can be performed with a large variety of methods. While the analysis of numeric data can be performed with several methods, the extraction of information from textual data remains a challenge. Currently employed natural language processing (NLP) methods do not really provide quantitative information that might be contained in IRs. In addition, the integration of numeric and textual data to identify possible causal relationships between data elements is still an unresolved challenge. This paper presents an approach to extract information from textual (e.g., incident or maintenance reports) and numeric data that relies on model based system engineer (MBSE) models. MBSE are diagrams designed to represent system and component dependencies (from both a form and functional point of view). In our approach, MBSE models emulate system engineer knowledge about component/system architecture. NLP methods are employed to perform syntactic and semantic analyses. Syntactic analysis analyzes the grammatical structure of a sentence while semantic analysis is designed to analyze the logic structure of a sentence. An innovative element of our approach is that semantic analysis uses MBSE models to identify links between textual elements. Similarly, numeric data is directly linked to elements of the MBSE models in order to map which functions are being monitored.

97 - MATHEMATICS AND COMPUTING↗

Bayesian Spatially Varying Multi-Regularization Image Deblurring

Many scientific experiments such as those found in astronomy, geology, microbiology, and X-ray radiography require the use of high-energy instruments to capture images. Since blur and noise are inevitably present in any imaging system, the images must be \deblurred" to extract the full information content. Mathematically, image deblurring is an ill-posed inverse problem that requires regularization. The regularization, in turn, has a large effect on the deblurred image: different regularization strengths, and types, lead to drastically different reconstructions. Moreover, many images contain a mixture of smooth and sharp features which suggests the use of multi-regularization, i.e., varying the type of regularization (e.g. Tikhonov or total variation) across the image. We address these issues by formulating the image deblurring problem within a hierarchical Bayesian framework in which we spatially adapt the strength of the regularization and also vary the regularization type across the image. In this way, the image itself, along with corresponding regularization strength at each pixel, are described jointly by a posterior distribution which we can sample by Markov chain Monte Carlo (MCMC) methods. We illustrate our techniques on simplified test problems and apply them to high-energy X-ray images taken at the Nevada National Security Site. Numerical tests show that our new method is robustly applicable and increases the quality of the image reconstruction when compared to other (Bayesian) methods.

97 MATHEMATICS AND COMPUTING↗

A New Semistructured Algebraic Multigrid Method

Multigrid methods are well suited to large massively parallel computer architectures because they are mathematically optimal and display good parallelization properties. Since current architecture trends are favoring regular compute patterns to achieve high performance, the ability to express structure has become much more important. The hypre software library provides high-performance multigrid preconditioners and solvers through conceptual interfaces, including a semistructured interface that describes matrices primarily in terms of stencils and logically structured grids. This paper presents a new semistructured algebraic multigrid (SSAMG) method built on this interface. The numerical convergence and performance of a CPU implementation of this method are evaluated for a set of semistructured problems. In conclusion, SSAMG achieves significantly better setup times than hypre’s unstructured AMG solvers and comparable convergence. In addition, the new method is capable of solving more complex problems than hypre’s structured solvers.

97 MATHEMATICS AND COMPUTING↗

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.↗

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.↗

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↗