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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Research on knowledge representation, machine learning, and knowledge acquisition

Research in knowledge representation, machine learning, and knowledge acquisition performed at Knowledge Systems Lab. is summarized. The major goal of the research was to develop flexible, effective methods for representing the qualitative knowledge necessary for solving large problems that require symbolic reasoning as well as numerical computation. The research focused on integrating different representation methods to describe different kinds of knowledge more effectively than any one method can alone. In particular, emphasis was placed on representing and using spatial information about three dimensional objects and constraints on the arrangement of these objects in space. Another major theme is the development of robust machine learning programs that can be integrated with a variety of intelligent systems. To achieve this goal, learning methods were designed, implemented and experimented within several different problem solving environments.

Buchanan, Bruce G.↗
DOE OSTI2018

Multidimensional Analogs of the Cauchy–Riemann System and Representations of Their Solutions via Harmonic Functions

At present, there are numerous multidimensional generalizations of holomorphic vectors. The most general of these is the four-dimensional generalization of the Cauchy–Riemann system. In the present work, by introducing two quaternion functions and the notion of quaternion differentiation, we obtain, for the first time, a five-dimensional generalization of holomorphic vectors. By using the representation of holomorphic vectors via the quaternion harmonic function and its derivatives, we consider the Riemann–Hilbert problem and one problem in a layer. A new solution of the Riemann–Hilbert problem in the five-dimensional half space is obtained.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2004

Kalman Filter Estimation of Spinning Spacecraft Attitude using Markley Variables

There are several different ways to represent spacecraft attitude and its time rate of change. For spinning or momentum-biased spacecraft, one particular representation has been put forward as a superior parameterization for numerical integration. Markley has demonstrated that these new variables have fewer rapidly varying elements for spinning spacecraft than other commonly used representations and provide advantages when integrating the equations of motion. The current work demonstrates how a Kalman filter can be devised to estimate the attitude using these new variables. The seven Markley variables are subject to one constraint condition, making the error covariance matrix singular. The filter design presented here explicitly accounts for this constraint by using a six-component error state in the filter update step. The reduced dimension error state is unconstrained and its covariance matrix is nonsingular.

Sedlak, Joseph E.↗
DOE OSTI2020

A Numerical Study of Hypoxia in Chesapeake Bay Using an Unstructured Grid Model: Validation and Sensitivity to Bathymetry Representation

A three-dimensional unstructured-grid hydrodynamic and water quality model (Semi-implicit Cross-scale Hydroscience Intergrated System Model-Integrated Compartment Model) is applied successfully for Chesapeake Bay. The model is validated with observations of salinity, chlorophyll-a, dissolved oxygen, nutrients, and phytoplankton productions from the year 1991 to 1995 for the mainstem and some major tributaries, based on multiple model skill scores. Model experiments are conducted to test the importance of having (1) an accurate representation of bathymetry to correctly predict hypoxia and other processes and (2) a high-resolution model grid for tributaries to correctly simulate water quality variables. Further, comparison with the model experiment results with bathymetry smoothing indicates that bathymetry smoothing, as commonly used for many systems, changes the stratification and lateral circulation pattern, resulting in more salt intrusion into shallow water regions, and an increase in the freshwater age. Consequently, a model with bathymetry smoothing can lead to an unrealistic prediction of the distribution of hypoxia and phytoplankton production. Local grid refinement shows significant improvement of model simulations on local stratification and water quality variables. Overall, the use of high-resolution unstructured grid model leads to a faithful representation of the complex geometry, and thus a seamless cross-scale capability for simulating water quality processes in the Bay including tributaries and tidal creeks.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1989

Polarization filtering of SAR data

A theoretical analysis of polarization filtering for the bistatic case is developed for optimum discrimination between two types of targets. The resulting method is half analytical and half numerical. Because it is based on the Stokes matrix representation, the targets of interest can be extended targets. The scattered field from such targets is partially polarized. This method is then applied to the monostatic case with numerical examples relying on the JPL (Jet Propulsion Laboratory) full-polarimetric L-band radar data. A matched filter to maximize the power ratio between urban and natural targets is developed. The results show that the same filter is optimal for both ocean and forest targets as natural targets.

Dubois, Pascale C.↗
NASA NTRS1976

Comparisons of several aerodynamic methods for application to dynamic loads analyses

The results of a study are presented in which the applicability at subsonic speeds of several aerodynamic methods for predicting dynamic gust loads on aircraft, including active control systems, was examined and compared. These aerodynamic methods varied from steady state to an advanced unsteady aerodynamic formulation. Brief descriptions of the structural and aerodynamic representations and of the motion and load equations are presented. Comparisons of numerical results achieved using the various aerodynamic methods are shown in detail. From these results, aerodynamic representations for dynamic gust analyses are identified. It was concluded that several aerodynamic methods are satisfactory for dynamic gust analyses of configurations having either controls fixed or active control systems that primarily affect the low frequency rigid body aircraft response.

Kroll, R. I.↗
DOE OSTI2015

Short‐term time step convergence in a climate model

Abstract This paper evaluates the numerical convergence of very short (1 h) simulations carried out with a spectral‐element (SE) configuration of the Community Atmosphere Model version 5 (CAM5). While the horizontal grid spacing is fixed at approximately 110 km, the process‐coupling time step is varied between 1800 and 1 s to reveal the convergence rate with respect to the temporal resolution. Special attention is paid to the behavior of the parameterized subgrid‐scale physics. First, a dynamical core test with reduced dynamics time steps is presented. The results demonstrate that the experimental setup is able to correctly assess the convergence rate of the discrete solutions to the adiabatic equations of atmospheric motion. Second, results from full‐physics CAM5 simulations with reduced physics and dynamics time steps are discussed. It is shown that the convergence rate is 0.4—considerably slower than the expected rate of 1.0. Sensitivity experiments indicate that, among the various subgrid‐scale physical parameterizations, the stratiform cloud schemes are associated with the largest time‐stepping errors, and are the primary cause of slow time step convergence. While the details of our findings are model specific, the general test procedure is applicable to any atmospheric general circulation model. The need for more accurate numerical treatments of physical parameterizations, especially the representation of stratiform clouds, is likely common in many models. The suggested test technique can help quantify the time‐stepping errors and identify the related model sensitivities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

The best models of metabolism

Biochemical systems are among of the oldest application areas of mathematical modeling. Spanning a time period of over one hundred years, the repertoire of options for structuring a model and for formulating reactions has been constantly growing, and yet, it is still unclear whether or to what degree some models are better than others and how the modeler is to choose among them. In fact, the variety of options has become overwhelming and difficult to maneuver for novices and experts alike. This review outlines the metabolic model design process and discusses the numerous choices for modeling frameworks and mathematical representations. It tries to be inclusive, even though it cannot be complete, and introduces the various modeling options in a manner that is as unbiased as that is feasible. However, the review does end with personal recommendations for the choices of default models. WIREs Syst Biol Med 2017, 9:e1391. doi: 10.1002/wsbm.1391 This article is categorized under: Analytical and Computational Methods > Dynamical Methods Models of Systems Properties and Processes > Mechanistic Models Biological Mechanisms > Metabolism

Voit, Eberhard O.↗
DOE OSTI2021

A unifying modeling abstraction for infinite-dimensional optimization

Infinite-dimensional optimization (InfiniteOpt) problems involve modeling components (variables, objectives, and constraints) that are functions defined over infinite-dimensional domains. Examples include continuous-time dynamic optimization (time is an infinite domain and components are functions of time), PDE optimization problems (space and time are infinite domains and components are functions of space-time), as well as stochastic and semi-infinite optimization (random space is an infinite domain and components are a function of such random space). InfiniteOpt problems also arise from combinations of these problem classes (e.g., stochastic PDE optimization). Given the infinite-dimensional nature of objectives and constraints, one often needs to define appropriate quantities (measures) to properly pose the problem. Moreover, InfiniteOpt problems often need to be transformed into a finite dimensional representation so that they can be solved numerically. Here in this work, we present a unifying abstraction that facilitates the modeling, analysis, and solution of InfiniteOpt problems. The proposed abstraction enables a general treatment of infinite-dimensional domains and provides a measure-centric paradigm to handle associated variables, objectives, and constraints. This abstraction allows us to transfer techniques across disciplines and with this identify new, interesting, and useful modeling paradigms (e.g., event constraints and risk measures defined over time domains). Our abstraction serves as the backbone of an intuitive Julia-based modeling package that we call InfiniteOpt.jl. We demonstrate the developments using diverse case studies arising in engineering.

42 ENGINEERING↗
DOE OSTI2020

A Model for Turbulence Spectra in the Equilibrium Range of the Stable Atmospheric Boundary Layer

Stratification can cause turbulence spectra to deviate from Kolmogorov's isotropic -5/3 power law scaling in the universal equilibrium range at high Reynolds numbers. However, a consensus has not been reached with regard to the exact shape of the spectra. Here we propose a shape of the turbulent kinetic energy and temperature spectra in horizontal wavenumber for the equilibrium range that consists of three regimes at small Froude number: the buoyancy subrange, a transition region, and the isotropic inertial subrange through dimensional analysis and substantial revision of previous theoretical approximation. These spectral regimes are confirmed by various observations in the atmospheric boundary layer. The representation of the transition region in direct numerical simulations will require large-scale separation between the Dougherty-Ozmidov scale and the Kolmogorov scale for strongly stratified turbulence at high Reynolds numbers, which is still challenging computationally. In addition, we suggest that the failure of Monin-Obukhov similarity theory in the very stable atmospheric boundary layer is due to the fact that it does not consider the buoyancy scale that characterizes the transition region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Multiscale Modeling Framework Using Element‐Based Galerkin Methods for Moist Atmospheric Limited‐Area Simulations

This paper presents a multiscale modeling framework (MMF) to model moist atmospheric limited-area weather. The MMF resolves large-scale convection using a coarse grid while simultaneously resolving local features through numerous fine local grids and coupling them seamlessly. Both large- and small-scale processes are modeled using the compressible Navier-Stokes equations within the Nonhydrostatic Unified Model of the Atmosphere (NUMA), and are discretized using a continuous element-based Galerkin method (spectral elements) with high-order basis functions. Consequently, the large-scale and small-scale models share the same dynamical core but have the flexibility to be adjusted individually. The proposed MMF method is tested in 2D and 3D idealized limited-area weather problems involving storm clouds produced by squall line and supercell simulations. Numerical results from the MMF showed enhanced representation of cloud processes compared to the coarse model.

Kang, Soonpil [Naval Postgraduate School, Monterey↗
DOE OSTI2016

Perspective: How good is DFT for water?

Kohn-Sham density functional theory (DFT) has become established as an indispensable tool for investigating aqueous systems of all kinds, including those important in chemistry, surface science, biology, and the earth sciences. Nevertheless, many widely used approximations for the exchange-correlation (XC) functional describe the properties of pure water systems with an accuracy that is not fully satisfactory. The explicit inclusion of dispersion interactions generally improves the description, but there remain large disagreements between the predictions of different dispersion-inclusive methods. We present here a review of DFT work on water clusters, ice structures, and liquid water, with the aim of elucidating how the strengths and weaknesses of different XC approximations manifest themselves across this variety of water systems. Our review highlights the crucial role of dispersion in describing the delicate balance between compact and extended structures of many different water systems, including the liquid. By referring to a wide range of published work, we argue that the correct description of exchange-overlap interactions is also extremely important, so that the choice of semi-local or hybrid functional employed in dispersion-inclusive methods is crucial. The origins and consequences of beyond-2-body errors of approximate XC functionals are noted, and we also discuss the substantial differences between different representations of dispersion. We propose a simple numerical scoring system that rates the performance of different XC functionals in describing water systems, and we suggest possible future developments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Beatification: Flattening the Poisson bracket for two-dimensional fluid and plasma theories

A perturbative method called beatification is presented for a class of two-dimensional fluid and plasma theories. The Hamiltonian systems considered, namely, the Euler, Vlasov-Poisson, Hasegawa-Mima, and modified Hasegawa-Mima equations, are naturally described in terms of noncanonical variables. The beatification procedure amounts to finding the correct transformation that removes the explicit variable dependence from a noncanonical Poisson bracket and replaces it with a fixed dependence on a chosen state in the phase space. As such, beatification is a major step toward casting the Hamiltonian system in its canonical form, thus enabling or facilitating the use of analytical and numerical techniques that require or favor a representation in terms of canonical, or beatified, Hamiltonian variables.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Detection of correlated fragments in a sequence of images by superimposed Fourier holograms

The problem of detecting correlated fragments in a sequence of images recorded by the superimposing holograms within the Fourier holography scheme with angular multiplication of a spatially modulated reference beam is considered. The approach to the solution of this problem is based on the properties of the variance of the image sum. It is shown that this problem can be solved by providing a constant distance between the signal and reference images when recording superimposed holograms and a partial mutual correlatedness of reference images. The detection efficiency is analysed from the point of view of estimated image data capacity, the degree of mutual correlation of reference images, and the hologram recording conditions. The results of a numerical experiment under the most complicated conditions (representation of images by realisations of homogeneous random fields) confirm the theoretical conclusions. (holographic memory)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Non-coherent continuum scattering as a line polarization mechanism

Line scattering polarization can be strongly affected by Rayleigh scattering at neutral hydrogen and Thomson scattering at free electrons. Often a depolarization of the continuum results, but the Doppler redistribution produced by the continuum scatterers, which are light (hence, fast), induces more complex interactions between the polarization in spectral lines and in the continuum. Here we formulate and solve the radiative transfer problem of scattering line polarization with non-coherent continuum scattering consistently. The problem is formulated within the spherical tensor representation of atomic and light polarization. The numerical method of solution is a generalization of the Accelerated Lambda Iteration that is applied to both the atomic system and the radiation field. We show that the redistribution of the spectral line radiation due to the non-coherence of the continuum scattering may modify the shape of the emergent fractional linear polarization patterns significantly, even yielding polarization signals above the continuum level in intrinsically unpolarizable lines.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Cross slip of extended dislocations in face-centered cubic metals through phase-field modeling

Cross slip is a dislocation mechanism that significantly impacts the mechanical behavior of engineering alloys. Here, in this work, we advance a 3D phase-field dislocation dynamics (PFDD) mesoscale technique to simulate cross slip across a broad range of face-centered cubic (FCC) metals. The formulation incorporates elastic anisotropy, an FCC numerical grid, and a high-fidelity representation of the entire γ -surface from density functional theory for eight FCC metals and no adjustable parameters or rules. The relaxed core structures under zero stress for all metals are predicted to extend in plane. The analytical model for stacking fault width agrees well with the PFDD result under the assumption of elastic isotropy but overestimates it under elastic anisotropy, when the degree of anisotropy is large. The dynamic simulations are designed to elucidate the material parameters that influence the propensity for cross slip. Whether cross slip occurs under a non-Schmid stress or to bypass a hard obstacle, the critical stress to cross slip scales strongly with the anisotropic energy coefficient for a screw dislocation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Spectral-partitioned Kohn-Sham density functional theory

Here we introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Vegetation Demographics in Earth System Models: a review of progress and priorities

Numerous current efforts seek to improve the representation of ecosystem ecology and vegetation demographic processes within Earth System Models (ESMs). Furthermore, these developments are widely viewed as an important step in developing greater realism in predictions of future ecosystem states and fluxes. Increased realism, however, leads to increased model complexity, with new features raising a suite of ecological questions that require empirical constraints. We review the developments that permit the representation of plant demographics in ESMs, and identify issues raised by these developments that highlight important gaps in ecological understanding. These issues inevitably translate into uncertainty in model projections but also allow models to be applied to new processes and questions concerning the dynamics of real-world ecosystems. We also argue that stronger and more innovative connections to data, across the range of scales considered, are required to address these gaps in understanding. The development of first-generation land surface models as a unifying framework for ecophysiological understanding stimulated much research into plant physiological traits and gas exchange. Constraining predictions at ecologically relevant spatial and temporal scales will require a similar investment of effort and intensified inter-disciplinary communication.

54 ENVIRONMENTAL SCIENCES↗