Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical Representation”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

The Challenge of Identifying Controls on Cloud Properties and Precipitation Onset for Cumulus Congestus Sampled During MC3E

Here, the controls on precipitation onset and the transition from shallow cumulus to congestus are explored using a suite of 16 large–eddy simulations based on the 25 May 2011 event from the Midlatitude Continental Convective Clouds Experiment (MC3E). The thermodynamic variables in the model are relaxed at various timescales to observationally constrained temperature and moisture profiles in order to better reproduce the observed behavior of precipitation onset and total precipitation. Three of the simulations stand out as best matching the precipitation observations and also perform well for independent comparisons of cloud fraction, precipitation area fraction, and evolution of cloud top occurrence. All three simulations exhibit a destabilization over time, which leads to a transition to deeper clouds, but the evolution of traditional stability metrics by themselves is not able to explain differences in the simulations. Conditionally sampled cloud properties (in particular, mean cloud buoyancy), however, do elicit differences among the simulations. The inability of environmental profiles alone to discern subtle differences among the simulations and the usefulness of conditionally sampled model quantities argue for hybrid observational/modeling approaches. These combined approaches enable a more complete physical understanding of cloud systems by combining observational sampling of time–varying three–dimensional meteorological quantities and cloud properties, along with detailed representation of cloud microphysical and dynamical processes from numerical models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Compressed sparse tensor based quadrature for vibrational quantum mechanics integrals

A new method for fast evaluation of high dimensional integrals arising in quantum mechanics is proposed. Here, the method is based on sparse approximation of a high dimensional function followed by a low-rank compression. In the first step, we interpret the high dimensional integrand as a tensor in a suitable tensor product space and determine its entries by a compressed sensing based algorithm using only a few function evaluations. Secondly, we implement a rank reduction strategy to compress this tensor in a suitable low-rank tensor format using standard tensor compression tools. This allows representing a high dimensional integrand function as a small sum of products of low dimensional functions. Finally, a low dimensional Gauss–Hermite quadrature rule is used to integrate this low-rank representation, thus alleviating the curse of dimensionality. Finally, numerical tests on synthetic functions, as well as on energy correction integrals for water and formaldehyde molecules demonstrate the efficiency of this method using very few function evaluations as compared to other integration strategies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Value of abstraction in performance assessment – When is a higher level of detail necessary?

In this study, different approaches in performance assessment (PA) of the long-term safety of a repository for radioactive waste were examined. This investigation was carried out as part of the DECOVALEX-2023 project, an international collaborative effort for research and model comparison. One specific task of the DECOVALEX-2023 project was the Salt Performance Assessment Modelling task (Salt PA), which aimed at comparing various models and methods employed in the performance assessment of deep geological repositories in salt. In the context of the Salt PA task, three distinct teams from SNL (United States), Quintessa Ltd (United Kingdom), and GRS (Germany) examined the consequences of employing different levels of abstractions when modelling the repository's geometry and implementing various features and processes, using the example of a simple hypothetical repository structure in domal salt. Each team applied their own tools: PFLOTRAN (SNL), QPAC (Quintessa) and LOPOS (GRS). These differ essentially regarding numerical concept and degree of detail in the representation of the underlying physical processes. The discussion focused on when simplifications can be appropriately applied and what consequences result from them. Furthermore, it was explored when and if a higher level of fidelity in geometry or physical processes is required.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Ballooning theory for micro-tearing mode in tokamak

This paper aims to investigate the impact of magnetic drift on the linear micro-tearing mode by using a kinetic approach to derive a reduced two-field eigen system in real space. Here, since the magnetic drift in real space has derivatives, it is more efficient to solve the mode equations in a Fourier-ballooning representation using the two-dimensional (2D) ballooning transform. The lowest-order eigen system in the Fourier-ballooning representation consists of two integral equations, which are numerically solved using the finite difference method for both eigenvalues and wave functions. The main results will be presented through graphical eigenvalue scans for each parameter. Furthermore, we present a graphical comparison between the predictions of the ballooning theory and GENE gyrokinetic code simulation in the pedestal region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

PSInet: a new global water potential network

Abstract Given the pressing challenges posed by climate change, it is crucial to develop a deeper understanding of the impacts of escalating drought and heat stress on terrestrial ecosystems and the vital services they offer. Soil and plant water potential play a pivotal role in governing the dynamics of water within ecosystems and exert direct control over plant function and mortality risk during periods of ecological stress. However, existing observations of water potential suffer from significant limitations, including their sporadic and discontinuous nature, inconsistent representation of relevant spatio-temporal scales and numerous methodological challenges. These limitations hinder the comprehensive and synthetic research needed to enhance our conceptual understanding and predictive models of plant function and survival under limited moisture availability. In this article, we present PSInet (PSI—for the Greek letter Ψ used to denote water potential), a novel collaborative network of researchers and data, designed to bridge the current critical information gap in water potential data. The primary objectives of PSInet are as follows. (i) Establishing the first openly accessible global database for time series of plant and soil water potential measurements, while providing important linkages with other relevant observation networks. (ii) Fostering an inclusive and diverse collaborative environment for all scientists studying water potential in various stages of their careers. (iii) Standardizing methodologies, processing and interpretation of water potential data through the engagement of a global community of scientists, facilitated by the dissemination of standardized protocols, best practices and early career training opportunities. (iv) Facilitating the use of the PSInet database for synthesizing knowledge and addressing prominent gaps in our understanding of plants’ physiological responses to various environmental stressors. The PSInet initiative is integral to meeting the fundamental research challenge of discerning which plant species will thrive and which will be vulnerable in a world undergoing rapid warming and increasing aridification.

Forestry↗
DOE OSTI2022

Gaussian processes meet NeuralODEs: a Bayesian framework for learning the dynamics of partially observed systems from scarce and noisy data

We present a machine learning framework (GP-NODE) for Bayesian model discovery from partial, noisy and irregular observations of nonlinear dynamical systems. The proposed method takes advantage of differentiable programming to propagate gradient information through ordinary differential equation solvers and perform Bayesian inference with respect to unknown model parameters using Hamiltonian Monte Carlo sampling and Gaussian Process priors over the observed system states. This allows us to exploit temporal correlations in the observed data, and efficiently infer posterior distributions over plausible models with quantified uncertainty. The use of the Finnish Horseshoe as a sparsity-promoting prior for free model parameters also enables the discovery of parsimonious representations for the latent dynamics. A series of numerical studies is presented to demonstrate the effectiveness of the proposed GP-NODE method including predator–prey systems, systems biology and a 50-dimensional human motion dynamical system. This article is part of the theme issue ‘Data-driven prediction in dynamical systems’.

Science & Technology - Other Topics↗
DOE OSTI2021

Steepest-descent algorithm for simulating plasma-wave caustics via metaplectic geometrical optics

The design and optimization of radiofrequency-wave systems for fusion applications is often performed using ray-tracing codes, which rely on the geometrical-optics (GO) approximation. However, GO fails at wave cutoffs and caustics. To accurately model the wave behavior in these regions, more advanced and computationally expensive “full-wave” simulations are typically used, but this is not strictly necessary. A new generalized formulation called metaplectic geometrical optics (MGO) has been proposed that reinstates GO near caustics. The MGO framework yields an integral representation of the wavefield that must be evaluated numerically in general. We present an algorithm for computing these integrals using Gauss-Freud quadrature along the steepest-descent contours. Benchmarking is performed on the standard Airy problem, for which the exact solution is known analytically. Furthermore, the numerical MGO solution provided by the new algorithm agrees remarkably well with the exact solution and significantly improves on previously derived analytical approximations of the MGO integral.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Time and Frequency Domain Methods for Basis Selection in Random Linear Dynamical Systems

Polynomial chaos methods have been extensively used to analyze systems in uncertainty quantification. Furthermore, several approaches exist to determine a low-dimensional approximation (or sparse approximation) for some quantity of interest in a model, where just a few orthogonal basis polynomials are required. In this work, we consider linear dynamical systems consisting of ordinary differential equations with random variables. The aim of this paper is to explore methods for producing low-dimensional approximations of the quantity of interest further. We investigate two numerical techniques to compute a low-dimensional representation, which both fit the approximation to a set of samples in the time domain. On the one hand, a frequency domain analysis of a stochastic Galerkin system yields the selection of the basis polynomials. It follows a linear least squares problem. On the other hand, a sparse minimization yields the choice of the basis polynomials by information from the time domain only. An orthogonal matching pursuit produces an approximate solution of the minimization problem. Finally, we compare the two approaches using a test example from a mechanical application.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A beginner's guide to non-abelian iPEPS for correlated fermions

Infinite projected entangled pair states (iPEPS) have emerged as a powerful tool for studying interacting two-dimensional fermionic systems. In this review, we discuss the iPEPS construction and some basic properties of this tensor network (TN) ansatz. Special focus is put on (i) a gentle introduction of the diagrammatic TN representations forming the basis for deriving the complex numerical algorithm, and (ii) the technical advance of fully exploiting non-abelian symmetries for fermionic iPEPS treatments of multi-band lattice models. The exploitation of non-abelian symmetries substantially increases the performance of the algorithm, enabling the treatment of fermionic systems up to a bond dimension D=24 D = 24 on a square lattice. A variety of complex two-dimensional (2D) models thus become numerically accessible. Here, we present first promising results for two types of multi-band Hubbard models, one with 2 2 bands of spinful fermions of \mathrm{SU}(2)_\mathrm{spin} \otimes \mathrm{SU}(2)_\mathrm{orb} S U ( 2 ) s p i n ⊗ S U ( 2 ) o r b symmetry, the other with 3 3 flavors of spinless fermions of \mathrm{SU}(3)_\mathrm{flavor} S U ( 3 ) f l a v o r symmetry.

Bruognolo, Benedikt↗
NASA NTRS1974

Analysis of three-point-bend test for materials with unequal tension and compression properties

An analysis capability is described for the three-point-bend test applicable to materials of linear but unequal tensile and compressive stress-strain relations. The capability consists of numerous equations of simple form and their graphical representation. Procedures are described to examine the local stress concentrations and failure modes initiation. Examples are given to illustrate the usefulness and ease of application of the capability. Comparisons are made with materials which have equal tensile and compressive properties. The results indicate possible underestimates for flexural modulus or strength ranging from 25 to 50 percent greater than values predicted when accounting for unequal properties. The capability can also be used to reduce test data from three-point-bending tests, extract material properties useful in design from these test data, select test specimen dimensions, and size structural members.

Chamis, C. C.↗
NASA NTRS1974

Nonlinear cellular motions in Poiseuille channel flow

In a number of nonlinear solutions to the equations of channel flow the velocity is decomposed into a mean part plus a nonlinear disturbance. The idea that nonlinear effects place a limitation on the amplitudes of the disturbance flow is considered. In the reported investigation the disturbance flow is represented by drastically limited Fourier expansions in the downstream coordinate. The resulting equations are solved numerically with high accuracy to obtain a good representation of the cross-stream structure of the solution. The results of the investigation show that indeed the nonlinear terms always limit the amplitude of the disturbance flow in this approximation.

Zahn, J.-P.↗
NASA NTRS1980

Vector analogues of the Maggi-Rubinowicz theory of edge diffraction

The Maggi-Rubinowicz technique for scalar and electromagnetic fields is interpreted as a transformation of an integral over an open surface to a line integral around its rim. Maggi-Rubinowicz analogues are found for several vector physical optics representations. For diffraction from a circular aperture, a numerical comparison between these formulations shows the two methods are in agreement. To circumvent certain convergence difficulties in the Maggi-Rubinowicz integrals that occur as the observer approaches the shadow boundary, a variable mesh integration is used. For the examples considered, where the ratio of the aperture diameter to wavelength is about ten, the Maggi-Rubinowicz formulation yields an 8 to 10 fold decrease in computation time relative to the physical optics formulation.

Meneghini, R.↗
NASA NTRS1981

On a fourth order accurate implicit finite difference scheme for hyperbolic conservation laws. II - Five-point schemes

This paper presents a family of two-level five-point implicit schemes for the solution of one-dimensional systems of hyperbolic conservation laws, which generalized the Crank-Nicholson scheme to fourth order accuracy (4-4) in both time and space. These 4-4 schemes are nondissipative and unconditionally stable. Special attention is given to the system of linear equations associated with these 4-4 implicit schemes. The regularity of this system is analyzed and efficiency of solution-algorithms is examined. A two-datum representation of these 4-4 implicit schemes brings about a compactification of the stencil to three mesh points at each time-level. This compact two-datum representation is particularly useful in deriving boundary treatments. Numerical results are presented to illustrate some properties of the proposed scheme.

Harten, A.↗
NASA NTRS1982

The tailward magnetopause field beyond 10 RE

The geomagnetic field in the near-earth magnetotail is made up of three major contributions related to the dipole field, the magnetotail current system, and the magnetopause current system exclusive of the return currents in the magnetotail current system. The present investigation is concerned with the third contribution. The calculation of the magnetospheric field is discussed, taking into account the calculation of the field inside the magnetopause on the basis of the numerical integration of the Biot Savart integral. The representation of the field for distances greater than 10 earth radii on the dark side of the magnetosphere is considered. The field internal to the magnetopause is expressed as the negative gradient of a scalar potential which is expressed as a sum of solutions to Laplace's equation in cylindrical coordinates. The origin of coordinates is the earth.

Beard, D. B.↗
NASA NTRS1984

A parametric study of the effect of inlet lip shape upon the radiated sound field

Far field sound radiation predictions for four different axisymmetric inlet lips excited by different tangential acoustic modes at several cut-off ratios are presented. These results were obtained by numerical integration of a special cylindrically symmetric integral representation of the external solutions of the Helmholtz equation which yields unique solutions at all wave numbers. The paper presents plots which detail the dependence of the relative SPL (dB) in the field upon the engine inlet lip shape, the modal input, and the cut-off ratio. Examination of these data indicate that: (1) as the inlet lip becomes larger the predominant acoustic radiation peak in the field becomes narrower and moves towards the centerline of the inlet; (2) as the order of the tangential acoustic mode of the driver increases the radiated sound peak again becomes narrower but moves away from the inlet centerline; and (3) as the cut-off ratio is increased for a specific tangential acoustic mode the predominant radiation peak becomes narrower and moves towards the centerline of the inlet.

Meyer, W. L.↗
NASA NTRS1985

Medium range forecasting: The experience of ECMWF

The results of five years of operational medium-range forecasts at ECMWF is discussed. It is shown that a considerable improvement in predictive skill has taken place, resulting in more accurate forecasts and an extension in time of useful predictive skill. For the Northern Hemisphere extratropics, the improvement in skill during the winter season is about three days for the 500-mb geopotential compared to the very first experiments in the early 1970s. Three are large geographical and temporal variations in the performance of the model, as well as large-scale errors of a systematic nature. Numerical experiments have indicated great sensitivity to the representation of orography and to tropical forcing.

Bengtsson, L.↗
NASA NTRS1990

Perspectives on the use of rule-based control

Issues regarding the application of artificial intelligence techniques to real-time control are discussed. Advantages associated with knowledge-based programming are discussed. A proposed rule-based control technique is summarized and applied to the problem of automated aircraft emergency procedure execution. Although emergency procedures are by definition predominately procedural, their numerous evaluation and decision points make a declarative representation of the knowledge they encode highly attractive, resulting in an organized and easily maintained software hierarchy. Simulation results demonstrate that real-time performance can be obtained using a microprocessor-based controller. It is concluded that a rule-based control system design approach may prove more useful than conventional methods under certain circumstances, and that declarative rules with embedded procedural code provide a sound basis for the construction of complex, yet economical, control systems.

Handelman, David A.↗
NASA NTRS2014

Near-Real Time Satellite-Retrieved Cloud and Surface Properties for Weather and Aviation Safety Applications

Cloud properties determined from satellite imager radiances provide a valuable source of information for nowcasting and weather forecasting. In recent years, it has been shown that assimilation of cloud top temperature, optical depth, and total water path can increase the accuracies of weather analyses and forecasts. Aircraft icing conditions can be accurately diagnosed in near-­‐real time (NRT) retrievals of cloud effective particle size, phase, and water path, providing valuable data for pilots. NRT retrievals of surface skin temperature can also be assimilated in numerical weather prediction models to provide more accurate representations of solar heating and longwave cooling at the surface, where convective initiation. These and other applications are being exploited more frequently as the value of NRT cloud data become recognized. At NASA Langley, cloud properties and surface skin temperature are being retrieved in near-­‐real time globally from both geostationary (GEO) and low-­‐earth orbiting (LEO) satellite imagers for weather model assimilation and nowcasting for hazards such as aircraft icing. Cloud data from GEO satellites over North America are disseminated through NCEP, while those data and global LEO and GEO retrievals are disseminated from a Langley website. This paper presents an overview of the various available datasets, provides examples of their application, and discusses the use of the various datasets downstream. Future challenges and areas of improvement are also presented.

Minnis, Patrick↗