SEARCH · Engineering Papers
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.
A genetic algorithm is used for design of infrared filters and in the understanding of the material structure of a resonant tunneling diode. These two components are examples of microdevices and nanodevices that can be numerically simulated using fundamental mathematical and physical models.
Mullite bond coats filled with silicon particles have recently been proposed for the next generation of environmental barrier coatings (EBCs). The sacrificial oxidation of the silicon particles eliminates some of the environmental durability issues associated with traditional EBC designs. To investigate this novel design approach, a mathematical formulation and resulting numerical method have been developed for simulating the oxidation of silicon particles embedded in a mullite matrix. The focus is on the oxidation of the silicon in a dry oxygen environment where the oxidant is presumably oxygen. Oxygen transport from the exterior environment through the mullite is assumed to occur via two lattice defect diffusion mechanisms. As a result, the rate of oxygen transport is dictated by the self-diffusivity of the lattice defects. The numerical method is implemented using the COMSOL Multiphysics® program (COMSOL, Inc.). The method is applied to simulate the oxidation of silicon particles in a mullite pellet in recent oxidation experiments. Numerical results compare favorably with the oxidation behavior observed in the oxidation experiments. The numerical method is also used to investigate the effects of various features of the mullite-silicon particle system on the oxidation behavior.
Explore the source record for details and available documents.
Mathematical model includes nonlinear effects. In development simulates growth of silicon ribbon from melt. Takes account of entire temperature and stress history of ribbon. Numerical simulations performed with new model helps in search for temperature distribution, pulling speed, and other conditions favoring growth of wide, flat, relatively defect-free silicon ribbons for solar photovoltaic cells at economically attractive, high production rates. Also applicable to materials other than silicon.
Here, we offer a modern interpretation of Lagrangian hydrodynamics as employed in Lagrangian simulations of compressible fluid flow. Our main result is to show that artificial viscosity, traditionally viewed as a numerical artifice to control unphysical oscillations in flows with shocks, actually represents a physical process and is necessary to derive accurate simulations in any compressible flow. We begin by reviewing the origins of two numerical devices, artificial viscosity and finite-volume methods. We proceed to construct a mathematical (PDE) model that incorporates those numerics and in which a new length scale, the observer, arises representing the discretization. Associated with that length scale, there are new inviscid fluxes that are the artificial viscosity as first formulated by Richtmyer and an artificial heat flux postulated by Noh but typically not included in Lagrangian codes. We discuss the connection of our results to bivelocity hydrodynamics. We conclude with some speculation as to the direction of future developments in multidimensional Lagrangian codes as computers get faster and have larger memories.
A numerical model of water and heat flow in geologic media was developed, verified, and tested. The hydraulic parameters (transmittivity and storativity) and the location of a linear hydrologic barrier were simulated and compared with results from field experiments involving two injection-storage-recovery cycles. For both cycles, the initial simulated and observed temperatures agree (55c).
Numerous results were obtained relevant to remote sensing, variational objective analysis, and data assimilation. A list of publications relevant in whole or in part is attached. The principal investigator gave many invited lectures, disseminating the results to the meteorological community as well as the statistical community. A list of invited lectures at meetings is attached, as well as a list of departmental colloquia at various universities and institutes.
Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet and environmentally friendly aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for aircraft noise reduction and are characterized by a frequency-dependent impedance. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic wave scattering by geometric bodies treated with acoustic liners. This work considers using either an impedance or an admittance (inverse of impedance) boundary condition to allow for acoustic scattering problems to be modeled with geometries consisting of both un-lined and lined surfaces. Three acoustic liner models are discussed: the Extended Helmholtz Resonator Model, the Three-Parameter Impedance Model, and the Broadband Impedance Model. In both the Helmholtz and Three-Parameter models, liner impedance is specified at a given frequency, whereas the Broadband model allows for the investigation of multiple frequencies simultaneously. The impedance and admittance boundary conditions for acoustic liners are derived for each model and coupled with a time-domain boundary integral equation. The scattering solution is obtained iteratively using a boundary element method with constant spatial and third-order temporal basis functions. Time-domain boundary integral equations are unfortunately prone to numerical instabilities due to resonant frequencies resulting from non-trivial solutions in the interior domain. When reformulated with the Burton-Miller method, the instabilities are eliminated. Using a Burton-Miller reformulation, the stability of the boundary element method assuming a liner boundary condition is assessed using eigenvalue analysis. The stability of each liner model is discussed, and it is shown that the Three-Parameter and Broadband models are sufficient for modeling an acoustic liner on the surface of scattering bodies. The Helmholtz model demonstrates strict limitations for stability, whereas the Three-Parameter and Broadband models are stable for most cases. Also included in this work is an assessment of the spatial accuracy of the time-domain boundary element method with respect to the surface element basis functions, as well as a performance study of the numerical algorithm.
Mathematical models of wave propagation in three-fluid plasma and numerical solutions to dispersion and propagation properties
A mathematical theory and an accompanying numerical scheme have been developed for predicting the oxidation behavior of C/SiC composite structures. The theory is derived from the mechanics of the flow of ideal gases through a porous solid. Within the mathematical formulation, two diffusion mechanisms are possible: (1) the relative diffusion of one species with respect to the mixture, which is concentration gradient driven and (2) the diffusion associated with the average velocity of the gas mixture, which is total gas pressure gradient driven. The result of the theoretical formulation is a set of two coupled nonlinear differential equations written in terms of the oxidant and oxide partial pressures. The differential equations must be solved simultaneously to obtain the partial vapor pressures of the oxidant and oxides as a function of space and time. The local rate of carbon oxidation is determined as a function of space and time using the map of the local oxidant partial vapor pressure along with the Arrhenius rate equation. The nonlinear differential equations are cast into matrix equations by applying the Bubnov-Galerkin weighted residual method, allowing for the solution of the differential equations numerically. The end result is a numerical scheme capable of determining the variation of the local carbon oxidation rates as a function of space and time for any arbitrary C/SiC composite structures.
A user's manual for a one dimensional thermal model to predict the temperature profiles of a deep body of water for any number of annual cycles is presented. The model is essentially a set of partial differential equations which are solved by finite difference methods using a high speed digital computer. The model features the effects of area change with depth, nonlinear interaction of wind generated turbulence and buoyancy, adsorption of radiative heat flux below the surface, thermal discharges, and the effects of vertical convection caused by discharge. The main assumption in the formulation is horizontal homogeneity. The environmental impact of thermal discharges from power plants is emphasized. Although the model is applicable to most lakes, a specific site (Lake Keowee, S.C.) application is described in detail. The programs are written in FORTRAN 5.
Abstract The cluster expansion method (CEM) is a widely used lattice-based technique in the study of multicomponent alloys. Despite its prevalent use, a clear understanding of expansion terms is lacking. We present a modern mathematical formalism of the CEM and introduce thecluster decomposition—a unique and basis-independent decomposition for functions of the atomic configuration in a crystal. We identify the cluster decomposition as an invariant ANOVA decomposition; and demonstrate how functional analysis of variance and sensitivity analysis can be used to interpret interactions among species. Furthermore, we show how the mathematical structure of the cluster decomposition enables numerical evaluation that scales with the number of clusters and is independent of the number of species. Overall, our work enables rigorous interpretations of interactions among species, provides opportunities to explore parameter estimation beyond linear regression, introduces a numerical efficient implementation, and enables analysis of cluster expansions based on established mathematical and statistical principles.
The numerical calculation is described of the steady-state flow of electrons in an axisymmetric, spherical, electrostatic collector for a range of boundary conditions. The trajectory equations of motion are solved alternately with Poisson's equation for the potential field until convergence is achieved. A direct (noniterative) numerical technique is used to obtain the solution to Poisson's equation. Space charge effects are included for initial current densities as large as 100 A/sq cm. Ways of dealing successfully with the difficulties associated with these high densities are discussed. A description of the mathematical model, a discussion of numerical techniques, results from two typical runs, and the FORTRAN computer program are included.
With the rise of Large Language Model (LLM) generated code, including in domains like scientific computing, ensuring not only syntactical, but also mathematical correctness, has become a critical task. Traditional formal methods approaches often struggle with the ambiguity of floating-point code, and full symbolic execution is extremely costly and limited. We propose a chain-of-reasoning approach that iteratively lifts basic semantics from code into the SPIRAL system and then establishes numerical equivalency to the desired mathematical operation. Here, we leverage the ample mathematical knowledge already formalized in SPIRAL to enable the system to recognize not just different implementations of the same algorithm but fully separate approaches to solving the given problem. The chain establishes tight error bounds on the output of given code with respect to the true continuous solution it approximates, quantifying all sources of error. We demonstrate this approach by establishing the correctness of a pseudospectral solver for a simple 1-dimensional Poisson problem.
Flow in fractured porous media occurs in the earth’s subsurface, in biological tissues, and in man-made materials. Fractures have a dominating influence on flow processes, and the last decade has seen an extensive development of models and numerical methods that explicitly account for their presence. To support these developments, four benchmark cases for single-phase flow in three-dimensional fractured porous media are presented. Furthermore, the cases are specifically designed to test the methods’ capabilities in handling various complexities common to the geometrical structures of fracture networks. Based on an open call for participation, results obtained with 17 numerical methods were collected. This paper presents the underlying mathematical model, an overview of the features of the participating numerical methods, and their performance in solving the benchmark cases.
We consider the radiative transfer problem in a plane-parallel slab of thermal electrons in the presence of an ultra-strong magnetic field (B approximately greater than B(sub c) approx. = 4.4 x 10(exp 13) G). Under these conditions, the magnetic field behaves like a birefringent medium for the propagating photons, and the electromagnetic radiation is split into two polarization modes, ordinary and extraordinary, that have different cross-sections. When the optical depth of the slab is large, the ordinary-mode photons are strongly Comptonized and the photon field is dominated by an isotropic component. Aims. The radiative transfer problem in strong magnetic fields presents many mathematical issues and analytical or numerical solutions can be obtained only under some given approximations. We investigate this problem both from the analytical and numerical point of view, provide a test of the previous analytical estimates, and extend these results with numerical techniques. Methods. We consider here the case of low temperature black-body photons propagating in a sub-relativistic temperature plasma, which allows us to deal with a semi-Fokker-Planck approximation of the radiative transfer equation. The problem can then be treated with the variable separation method, and we use a numerical technique to find solutions to the eigenvalue problem in the case of a singular kernel of the space operator. The singularity of the space kernel is the result of the strong angular dependence of the electron cross-section in the presence of a strong magnetic field. Results. We provide the numerical solution obtained for eigenvalues and eigenfunctions of the space operator, and the emerging Comptonization spectrum of the ordinary-mode photons for any eigenvalue of the space equation and for energies significantly lesser than the cyclotron energy, which is on the order of MeV for the intensity of the magnetic field here considered. Conclusions. We derived the specific intensity of the ordinary photons, under the approximation of large angle and large optical depth. These assumptions allow the equation to be treated using a diffusion-like approximation.
Some of the recent work at NASA Langley on transition, turbulence, and reacting flows is summarized. Much of this effort has been motivated by outstanding technological problems in high-speed flow. A class of numerical algorithms for these problems has been developed and a variety of physical problems have been simulated. Descriptions are provided of the basic mathematical models, the nature of the numerical methods, and some of the recent simulations.
A numerical study for confined, axisymmetrical, turbulent diffusion flames is presented. Local mean gas properties are predicted by solving the appropriate conservation equations in the finite difference form with the corresponding boundary conditions. The k-epsilon two-equation turbulence model is employed to describe the turbulent nature of the flow. A two-step kinetic model is assumed to govern the reaction mechanism. The finite reaction rate is the smaller of an Arrhenius type of reaction rate and a modified version of eddy-breakup model. Reasonable agreement is observed between calculations and measurements, but to obtain better agreement, more work is needed on improvements of the above mathematical models. However, the present numerical study offers an improvement in the analysis and design of the gas turbine combustors.