Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mathematical 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 1,387 records · Page 77

An overview of methods to identify and manage uncertainty for modelling problems in the water-environment-agriculture cross-sector

Uncertainty pervades the representation of systems in the water–environment–agriculture cross-sector. Successful methods to address uncertainties have largely focused on standard mathematical formulations of biophysical processes in a single sector, such as partial or ordinary differential equations. More attention to integrated models of such systems is warranted. Model components representing the different sectors of an integrated model can have less standard, and different, formulations to one another, as well as different levels of epistemic knowledge and data informativeness. Thus, uncertainty is not only pervasive but also crosses boundaries and propagates between system components. Uncertainty assessment (UA) cries out for more eclectic treatment in these circumstances, some of it being more qualitative and empirical. Here in this paper, we discuss the various sources of uncertainty in such a cross-sectoral setting and ways to assess and manage them. We have outlined a fast-growing set of methodologies, particularly in the computational mathematics literature on uncertainty quantification (UQ), that seem highly pertinent for uncertainty assessment. There appears to be considerable scope for advancing UA by integrating relevant UQ techniques into cross-sectoral problem applications. Of course this will entail considerable collaboration between domain specialists who often take first ownership of the problem and computational methods experts.

97 MATHEMATICS AND COMPUTING↗

Spacecraft attitude determination using the earth's magnetic field

A method is presented by which the attitude of a low-Earth orbiting spacecraft may be determined using a vector magnetometer, a digital Sun sensor, and a mathematical model of the Earth's magnetic field. The method is currently being implemented for the Solar Maximum Mission spacecraft (as a backup for the failing star trackers) as a way to determine roll gyro drift.

Simpson, David G.↗

An improved iterative HDG approach for partial differential equations

Here we propose and analyze an iterative high-order hybridized discontinuous Galerkin (iHDG) discretization for linear partial differential equations. We improve our previous work [45] in several directions: 1) the improved iHDG approach converges in a finite number of iterations for the scalar transport equation; 2) it is unconditionally convergent for both the linearized shallow water system and the convection–diffusion equation; 3) it has improved stability and convergence rates; 4) we uncover a relationship between the number of iterations and time stepsize, solution order, meshsize and the equation parameters. This allows us to choose the time stepsize such that the number of iterations is approximately independent of the solution order and the meshsize; and 5) we provide both strong and weak scalings of the improved iHDG approach up to 16,384 cores. A connection between iHDG and time integration methods such as parareal and implicit/explicit methods are discussed. Extensive numerical results for linear (and nonlinear) PDEs are presented to verify the theoretical findings.

97 MATHEMATICS AND COMPUTING↗

Application of HCMM data to regional geologic analysis for mineral and energy resource evaluation

Experimentation with several potentially promising techniques led to the selection of a fairly simple scheme for registration of data from the HCMM using an affine transformation. A method based solely on remote sensing was developed to estimate those meteorological effects which are required for thermal inertia mapping. It assumes that the atmospheric fluxes are spatially invariant and that the solar, sky and sensible heat fluxes can be approximated by a simple mathematical form. Coefficients are determined from least squares method by fitting observational data to the thermal model.

Watson, K.↗

A Weak Galerkin Method for the Reissner–Mindlin Plate in Primary Form

We developed a new finite element method for the Reissner–Mindlin equations in its primary form by using the weak Galerkin approach. Like other weak Galerkin finite element methods, this one is highly flexible and robust by allowing the use of discontinuous approximating functions on arbitrary shape of polygons and, at the same time, is parameter independent on its stability and convergence. Furthermore, error estimates of optimal order in mesh size h are established for the corresponding weak Galerkin approximations. Numerical experiments are conducted for verifying the convergence theory, as well as suggesting some superconvergence and a uniform convergence of the method with respect to the plate thickness.

97 MATHEMATICS AND COMPUTING↗

Fourth Order Finite Difference Methods for the Wave Equation with Mesh Refinement Interfaces

In this work, we analyze two types of summation-by-parts finite difference operators for approximating the second derivative with variable coefficient. The first type uses ghost points, while the second type does not use any ghost points. A previously unexplored relation between the two types of summation-by-parts operators is investigated. By combining them we develop a new fourth order accurate finite difference discretization with hanging nodes on the mesh refinement interface. We take the model problem as the two-dimensional acoustic wave equation in second order form in terms of acoustic pressure, and we prove energy stability for the proposed method. Compared to previous approaches using ghost points, the proposed method leads to a smaller system of linear equations that needs to be solved for the ghost point values. Another attractive feature of the proposed method is that the explicit time step does not need to be reduced relative to the corresponding periodic problem. Numerical experiments, both for smoothly varying and discontinuous material properties, demonstrate that the proposed method converges to fourth order accuracy. A detailed comparison of the accuracy and the time-step restriction with the simultaneous-approximation-term penalty method is also presented. (An erratum is attached.)

97 MATHEMATICS AND COMPUTING↗

GASP Advancements With Symbolic Computations, Optimization, and Decoupled Numerical Methods

NASA is exploring advancements in hybrid-electric propulsion concepts for aircraft design. A new analysis capability has been deployed at NASA’s Ames Research Center using a new mathematical modeling framework in Python that leverages symbolic representation of mathematical expressions, optimization-based problem formulations, and standalone numerical methods. This capability, called Gascon, has been verified against NASA’s General Aviation Synthesis Program (GASP) results for several vehicles that exercise the analysis capability for a variety of aircraft classes and engine models, including both turbofans and turboprops. We believe this capability will set the stage for more rapid development of novel aircraft models with more flexible design assessments using optimization, and more traceable analysis for reproducibility. Additionally, a novel Pareto optimization-based design space exploration that shows the trade between range flown and fuel burn by manipulating the thrust split during cruise for a true parallel hybrid aircraft was developed by leveraging the unique features of Gascon.

aircraft design↗

Nuclear Forensics Analysis with Missing and Uncertain Data

We have applied a new imputation-based method for analyzing incomplete data, called Monte Carlo Bayesian Database Generation (MCBDG), to the Spent Fuel Isotopic Composition (SFCOMPO) database. About 60% of the entries are absent for SFCOMPO. The method estimates missing values of a property from a probability distribution created from the existing data for the property, and then generates multiple instances of the completed database for training a machine learning algorithm. Uncertainty in the data is represented by an empirical or an assumed error distribution. The method makes few assumptions about the underlying data, and compares favorably against results obtained by replacing missing information with constant values.

97 MATHEMATICS AND COMPUTING↗

Questions of Simplifying the Methods of Testing Photoconverters Based on Cu(In,Ga)Se{sub 2}. Part 1: A Basis for Mathematical Modeling of the Current–Voltage Characteristics

An analysis of possible errors of laboratory studies of current–voltage (I–V) characteristics for solar cells based on Cu(In,Ga)Se{sub 2} (CIGS) using different lamps is provided, and the possibility of simplifying the procedure of testing the cells and panels by calculations and modeling is presented. Formulas are presented to develop mathematical models of the correlation between the quantum efficiency and the spectrum of the used lamp to obtain (I–V) characteristics of solar cells based on CIGS, which simplifies the testing procedure and data collection.

14 SOLAR ENERGY↗

Multilateration - A nondegenerate method of obtaining station coordinates and satellite ephemerides

A technique for the determination of three-dimensional station coordinates and satellite ephemerides is developed which is based on the principle of multilateration. The method makes use of a system of six ground stations which simultaneously measure the slant range between each station and one or two satellites. It is demonstrated that a minimum of six stations is required in order to yield a system which will be free of mathematical degeneracies. It will be seen that the method of multilateration is not dependent upon the position of the satellite or any other dynamical considerations in the equations used to determine the relative station coordinates. In fact, the satellite coordinates are obtained as a direct by-product of the method. Numerical results are presented which indicate that the method of multilateration can determine the relative three-dimensional station coordinates with an accuracy that is limited only by the hardware measurement system. If a highly accurate laser ranging system is used, then accuracies in the 1-cm range can be expected.

Ong, K. M.↗

A general framework for substructuring-based domain decomposition methods for models having nonlocal interactions [minus appendix B]

A rigorous mathematical framework is provided for a substructuring-based domain-decomposition approach for nonlocal problems that feature interactions between points separated by a finite distance. Here, by substructuring it is meant that a traditional geometric configuration for local partial differential equation problems is used in which a computational domain is subdivided into non-overlapping subdomains. In the nonlocal setting, this approach is substructuring-based in the sense that those subdomains interact with neighboring domains over interface regions having finite volume, in contrast to the local PDE setting in which interfaces are lower dimensional manifolds separating abutting subdomains Key results include the equivalence between the global, single-domain nonlocal problem and its multi-domain reformulation, both at the continuous and discrete levels. These results provide the rigorous foundation necessary for the development of efficient solution strategies for nonlocal domain-decomposition methods.

97 MATHEMATICS AND COMPUTING↗

The arbitrary order mixed mimetic finite difference method for the diffusion equation

Here, we propose an arbitrary-order accurate mimetic finite difference (MFD) method for the approximation of diffusion problems in mixed form on unstructured polygonal and polyhedral meshes. As usual in the mimetic numerical technology, the method satisfies local consistency and stability conditions, which determines the accuracy and the well-posedness of the resulting approximation. The method also requires the definition of a high-order discrete divergence operator that is the discrete analog of the divergence operator and is acting on the degrees of freedom. The new family of mimetic methods is proved theoretically to be convergent and optimal error estimates for flux and scalar variable are derived from the convergence analysis. A numerical experiment confirms the high-order accuracy of the method in solving diffusion problems with variable diffusion tensor. It is worth mentioning that the approximation of the scalar variable presents a superconvergence effect.

97 MATHEMATICS AND COMPUTING↗

Hermite Methods for the Scalar Wave Equation

Arbitrary order dissipative and conservative Hermite methods for the scalar wave equation are presented here. Both methods use $(m+1)^d$ degrees of freedom per node for the displacement in $d$ dimensions; the dissipative and conservative methods achieve orders of accuracy $(2m-1)$ and $2m$, respectively. Stability and error analyses as well as implementation strategies for accelerators are also given.

97 MATHEMATICS AND COMPUTING↗

Prediction of countercurrent flow limitation and its uncertainty in a PWR hot leg

In previous studies, we did small-scale air-water experiments for the hot leg and pressurizer surge line in a PWR and three-dimensional (3D) numerical simulations for PWR conditions to evaluate countercurrent flow limitation (CCFL) there, and we developed a one-dimensional (1D) prediction method to generalize prediction of CCFL in a piping system such as the hot leg and pressurizer surge line. In this study, we evaluated computation errors of the 1D prediction method in comparison with CCFL data in the PWR hot leg models with the diameter D and length L of D = 50-650 mm and L/D = 4.5-9.5, and then evaluated uncertainties of the 1D prediction method in comparison with 3D numerical simulations for PWR conditions with D = 750 mm and L/D = 8.6. We found the computation errors were within the difference of dC = ±0.05 for the CCFL constant of C = 0.63 in the correlation with the Wallis parameters, but the computation errors include experimental errors and they are not uncertainties for the PWR conditions. The total uncertainty of the 1D computations for the actual PWR hot leg (including effects of the diameter and fluid properties) was estimated to be dC = ±0.035. The CCFL correlation and its uncertainty of dC = ±0.035 derived from the 1D prediction method were very close to those previously obtained from CCFL data and 3D simulations (dC = ±0.03). (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Feature isolation and quantification of evolving datasets

Identifying and isolating features is an important part of visualization and a crucial step for the analysis and understanding of large time-dependent data sets (either from observation or simulation). In this proposal, we address these concerns, namely the investigation and implementation of basic 2D and 3D feature based methods to enhance current visualization techniques and provide the building blocks for automatic feature recognition, tracking, and correlation. These methods incorporate ideas from scientific visualization, computer vision, image processing, and mathematical morphology. Our focus is in the area of fluid dynamics, and we show the applicability of these methods to the quantification and tracking of three-dimensional vortex and turbulence bursts.

Source record↗

Green's function methods for multiphysics simulations (Final Report)

Green's functions are important tools for analyzing mathematical properties of partial differential equations (PDEs), and for numerically solving PDEs, especially when equations for the same operator but with multiple right hand sides need to be solved simultaneously. This proposal aims at developing efficient and accurate numerical methods for computing Green's functions, which can be used to tackle a challenging question in multiphysics simulation of DOE-mission science problems: how to couple quantum physics with classical physics. The key mathematical difficulty is properly formulate a "boundary condition" for the region described by quantum physics, and conventional approaches often model such boundary conditions in an empirical way. The proposed Green's function methods use the Dirichlet-to-Neumann map to formulate a boundary condition that is non-empirical and can couple the quantum and classical regions in an in principle exact way. The key ingredient of the new methods is to construct the Dirichlet-to-Neumann map in an efficient, accurate and versatile manner. The new methods have provably low complexity and are ideally suited for massively parallel and emerging many-core computational systems.

97 MATHEMATICS AND COMPUTING↗

On OM-decomposable sets

We introduce and study the family of sets in a finite-dimensional Euclidean space which can be written as the Minkowski sum of an open, bounded, and convex set and a closed convex cone. We establish several properties of the class of such sets, called OM-decomposable, some of which extend related properties holding for the class of Motzkin decomposable sets (i.e., those for which the convex and bounded set in the decomposition is requested to be closed, hence compact).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗