Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial discretization”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Topological terms with qubit regularization and relativistic quantum circuits

Qubit regularization provides a rich framework to explore quantum field theories. The freedom to choose how the important symmetries of the theory are embedded in the qubit regularization scheme allows us to construct new lattice models with rich phase diagrams. Some of the phases can contain topological terms which lead to critical phases. In this work we introduce and study the SU(3)-F qubit regularization scheme to embed the SO(3) spin-symmetry. We argue that qubit models in this regularization scheme contain several phases including a critical phase which describes the k = 1 Wess-Zumino-Witten (WZW) conformal field theory (CFT) at long distances, and two massive phases one of which is trvially gapped and the other which breaks the lattice translation symmetry. We construct a simple space-time Euclidean lattice model with a single coupling U and study it using the Monte Carlo method. We show the model has a critical phase at small U and a trivially massive phase at large U with a first order transition separating the two. Another feature of our model is that it is symmetric under space-time rotations, which means the temporal and spatial lattice spacing are connected to each other. The unitary time evolution operator obtained by a Wick rotation of the transfer matrix of our model can help us compute the physics of the k = 1 WZW CFT in real time without the need for tuning the temporal lattice spacing to zero. We use this idea to introduce the concept of a relativistic quantum circuit on a discrete space-time lattice.

Bhattacharya, Tanmoy↗
NASA NTRSDate not supplied

Nonlinear Homogenization of Finitely Deformed Viscoelastic-Viscoplastic Composites Using Mechanics of Structure Genome

The objective of this paper is to develop a micromechanics approach to homogenizing finitely deformed viscoelastic-viscoplastic composites using the mechanics of structure genome. The incremental constitutive relation for glassy polymers, formulated in the spatial configuration, is implemented in the present approach.This involves (1) pulling-back the constitutive model to the material configuration and (2)choosing the deformation gradient tensor and the first Piola–Kirchhoff stress tensor as the strain and the stress measures during homogenization, respectively. An Euler–Newton predictor–corrector method is developed for homogenization. Each step involves formulating a variational statement using the mechanics of structure genome, discretizing the statement in a finite-dimensional space, and solving the problem using an Euler/multilevel Newton method. The present approach is demonstrated by homogenizing fiber- and particle-reinforced composites undergoing uniaxial, biaxial, or shear deformation, at different stain rates.

Multi-scale modeling, High Strain Composites, Visc↗
NASA NTRS2024

Real-Time Temperature Field State Observer and Digital Replica to Support Cyber-Physical Testing

Cyber-physical testing is a class of experimental methods in which a system is partitioned into physical and numerical components to study its overall behavior. Transfer systems are typically needed to capture relevant interactions by enforcing appropriate interface conditions between the physical and numerical portions. While cyber-physical testing is often performed using lumped-parameter systems, thermomechanical cyber-physical testing requires the use of thermal actuators as transfer systems to apply a thermal condition over a spatially distributed region or surface. Thus, the interactions between the numerical and physical components are of a distributed nature, i.e., a field rather than point values. In this paper, we develop and experimentally validate a state observer for temperature fields to enable thermomechanical cyber-physical testing. The state observer provides continuous estimates of temperature over a surface of the thermal actuator using only temperature measurements at discrete locations. The estimated temperature field is then leveraged to account for localized mechanical defects in the interface condition. The method consists of building a digital replica of one portion of the control plant. By imposing defects in the digital replica it is used to generate and output the response of a thermal load pattern that represents that of a physical system with mechanical defects.

Manuel Salmeron↗
DOE OSTI2023

Fluoride-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model

In this report we present work performed in Fiscal Year 2022 that demonstrates the modeling and simulation of a fully coupled neutronics thermal hydraulics reference plant model for a fluoride-cooled high-temperature pebble-bed reactor. The multiphysics model is developed on the Nuclear Regulatory Commission’s Comprehensive Reactor Analysis Bundle (BlueCRAB) available on the Idaho National Laboratory’s high-performance computer, which natively and seamlessly couples Griffin, Pronghorn, and the BISON Multiphysics Object-Oriented Simulation Environment based applications. Griffin provides reactor physics capabilities, including depletion to the equilibrium core, k-eigenvalue, adjoint, and transient. The unique direct equilibrium core capability in Griffin is based on a streamline methodology to spatially deplete the pebbles into burnup groups. Pronghorn solves the porous medium equations for the fluid regions and conduction in the solid regions and incorporates a fluidic diode model to simulate the transition from forced to natural convection during accident scenarios. MOOSE modules solves thermal conduction problems for the pebbles and tristructural isotropic in the pebble-bed core, thus providing the fuel and moderator spatial fields for each pebble burnup group. The neutronics feedback relies primarily on fuel, moderator, and reflector temperatures as as well as the FLiBe salt density. Here, we present results for the uncoupled equilibrium core and perform comparisons to equivalent Monte Carlo models. The power distributions and kinetic parameters obtained with Griffin are consistent with those computed with Griffin. We demonstrate a noticeable improvement with the use of discrete ordinates method (SN) transport. The coupled steady-state equilibrium core provides the initial condition for two time-dependent problems: a control rod withdrawal event and an unprotected loss of flow event. In both cases, the reactor design is self-stabilizing and the solutions are consistent with the expected physics. Although this model is prototypical regarding BlueCRAB’s capabilities, its results are consistent with published work by Kairos Power and other research entities. Significant improvements to the model are planned in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1992

Luminance-model-based DCT quantization for color image compression

A model is developed to approximate visibility thresholds for discrete cosine transform (DCT) coefficient quantization error based on the peak-to-peak luminance of the error image. Experimentally measured visibility thresholds for R, G, and B DCT basis functions can be predicted by a simple luminance-based detection model. This model allows DCT coefficient quantization matrices to be designed for display conditions other than those of the experimental measurements: other display luminances, other veiling luminances, and other spatial frequencies (different pixel spacings, viewing distances, and aspect ratios).

Ahumada, Albert J., Jr.↗
DOE OSTI2022

Computational Methods to Characterize Panel Loading Conditions for Accelerated Testing

Panel cracking and degradation due to wind loading are known to have a detrimental effect on power output. Recreating these damaging conditions in a controlled experimental setting requires an understanding of the loads being generated at the panel surface as a function of both wind speed, wind direction, and panel orientation. To better understand these relationships, a computational fluid dynamics (CFD) simulation package was constructed using an open-source Python library for solving partial differential equations with the finite element method. This CFD package allows the simulated wind speed and panel orientation of a multi-panel array to be easily changed and provides the capability to measure the traction forces at discrete points along the sun-facing and ground-facing surface of an interior panel residing in the wake of one or more upstream panels. A parameter exploration was performed in which the panel angle was varied in even increments from -70 deg to 70 deg and the wind speed was varied from 2 to 30 m/s. During post processing, the measured traction along the panel surface was averaged spatially and interpreted as a time-varying signal, where further processing of this signal yielded the root mean square amplitude and a characteristic frequency associated with the loading. This study found that higher wind speeds are generally associated with increased amplitude of loading and that the panel orientation angle can significantly exacerbate or mitigate this loading. These outcomes are presented along with current work on higher-fidelity verification simulations and recommendations for performing accelerated experimental testing.

loading↗
NASA NTRS1996

Numerical Modeling of Spray Combustion with an Unstructured-Grid Method

The present unstructured-grid method follows strictly the basic finite volume forms of the conservation laws of the governing equations for the entire flow domain. High-order spatially accurate formulation has been employed for the numerical solutions of the Navier-Stokes equations. A two-equation k-epsilon turbulence model is also incorporated in the unstructured-grid solver. The convergence of the resulted linear algebraic equation is accelerated with preconditioned Conjugate Gradient method. A statistical spray combustion model has been incorporated into the present unstructured-grid solver. In this model, spray is represented by discrete particles, rather than by continuous distributions. A finite number of computational particles are used to predict a sample of total population of particles. Particle trajectories are integrated using their momentum and motion equations and particles exchange mass, momentum and energy with the gas within the computational cell in which they are located. The interaction calculations are performed simultaneously and eliminate global iteration for the two-phase momentum exchange. A transient spray flame in a high pressure combustion chamber is predicted and then the solution of liquid-fuel combusting flow with a rotating cup atomizer is presented and compared with the experimental data. The major conclusion of this investigation is that the unstructured-grid method can be employed to study very complicated flow fields of turbulent spray combustion. Grid adaptation can be easily achieved in any flow domain such as droplet evaporation and combustion zone. Future applications of the present model can be found in the full three-dimensional study of flow fields of gas turbine and liquid propulsion engine combustion chambers with multi-injectors.

Shang, H. M.↗
NASA NTRS2021

Intergenerational Inequities in Exposure to Climate Extremes: Young Generations Are Severely Threatened By Climate Change

Under continued global warming, extreme events such as heat waves will continue to rise in frequency, intensity, duration, and spatial extent over the next decades. Younger generations are therefore expected to face more such events across their lifetimes compared with older generations. This raises important issues of solidarity and fairness across generations that have fueled a surge of climate protests led by young people in recent years and that underpin issues of intergenerational equity raised in recent climate litigation. However, the standard scientific paradigm is to assess climate change in discrete time windows or at discrete levels of warming, a “period” approach that inhibits quantification of how much more extreme events a particular generation will experience over its lifetime compared with another. By developing a “cohort” perspective to quantify changes in lifetime exposure to climate extremes and compare across generations (see the first figure), we estimate that children born in 2020 will experience a two- to sevenfold increase in extreme events, particularly heat waves, compared with people born in 1960, under current climate policy pledges. Our results highlight a severe threat to the safety of young generations and call for drastic emission reductions to safeguard their future. Meteorological extremes, hazards, or climate change impacts are mostly studied as they evolve over time under varying emission scenarios and socioeconomic pathways. For example, applying a heat wave indicator (see table S1) to four bias-adjusted global climate models indicates that the land area annually affected by such heat waves will increase from ~15% around 2020 to ~22% by 2100 under a scenario compatible with limiting global warming to 1.5°C, and to ~46% under a scenario in line with current emission reduction pledges (see the first figure). Recent studies extended this approach, studying aspects of climate change as a function of global mean temperature (GMT) increments, highlighting the scenario-independence of several extreme event indicators but remaining, in essence, a comparison of time windows. By contrast, we performed a birth cohort analysis by combining a collection of multimodel extreme event projections with country-scale life expectancy information, gridded population data, and future global temperature trajectories from the Intergovernmental Panel on Climate Change (IPCC) Special Report on Global Warming of 1.5°C (see supplementary materials). By integrating the exposure of an average person in a country or region to extreme events across their lifetime, we encapsulate spatiotemporal changes in climate hazards, population density, cohort size, and life expectancy (see the first figure).

global warming↗
DOE OSTI2025

Co-designing Spectral Transformation Oracles with Hybrid Oscillator-Qubit Quantum Processors: From Algorithms to Compilation

We co-design a family of quantum eigenvalue transformation oracles that can be efficiently implemented on hybrid discrete- or continuous-variable (qubit or qumode) hardware. To illustrate the oracle’s representation-theoretic power and near-term experimental accessibility, we encode a Gaussian imaginary time-evolution spectral filter. As a result, we define a continuous linear combination of unitaries block encoding. Due to the ancillary qumode’s infinite-dimensional nature, continuous-variable qumodes constitute a powerful compilation tool for encoding continuous spectral functions without discretization errors while minimizing resource requirements. We then focus on the ubiquitous task of preparing eigenstates in quantum spin models. For completeness, we provide an end-to-end compilation which expresses high-level oracles in terms of an experimentally realizable instruction set architecture in both 1D and 2D. Finally, we examine the leading-order effects of physical errors and highlight open research directions. Our algorithms scale linearly with the spatial extent of the target system and are applicable to both near-term and large-scale quantum processors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1973

Surface electrical properties experiment

The surface electrical properties (SEP) experiment was used to explore the subsurface material of the Apollo 17 landing site by means of electromagnetic radiation. The experiment was designed to detect electrical layering, discrete scattering bodies, and the possible presence of water. From the analysis of the data, it was expected that values of the electrical properties (dielectric constant and loss tangent) of lunar material in situ would be obtained. The basic principle of the SEP experiment is interferometry. This principle involves only the interference of two or more waves to produce an interference pattern. The inversion of the interference pattern in terms of the spatial distribution of the electrical properties is the basic aim of the experiment (fig. 15-1). The experiment is most easily understood in terms of a single dipole antenna for radiating electromagnetic energy and a loop receiver for measuring the magnitudes of the fields.

Simmons, G.↗
NASA NTRS1990

Modeling bidirectional reflectance of forests and woodlands using Boolean models and geometric optics

Geometric-optical discrete-element mathematical models for forest canopies have been developed using the Boolean logic and models of Serra. The geometric-optical approach is considered to be particularly well suited to describing the bidirectional reflectance of forest woodland canopies, where the concentration of leaf material within crowns and the resulting between-tree gaps make plane-parallel, radiative-transfer models inappropriate. The approach leads to invertible formulations, in which the spatial and directional variance provides the means for remote estimation of tree crown size, shape, and total cover from remotedly sensed imagery.

Strahler, Alan H.↗
NASA NTRS1992

Combustion noise and combustion instabilities in propulsion systems

This paper is concerned with some aspects of non-linear behavior of unsteady motions in combustion chambers. The emphasis is on conditions under which organized oscillations having discrete frequencies may exist in the presence of random motions. In order to treat the two types of motions together, and particularly to investigate coupling between noise and combustion instabilities, the unsteady field is represented as a synthesis of acoustic modes having time-varying amplitudes. Each of the amplitudes are written as the sum of two parts, one associated with the random field and the remainder representing the organized oscillations. After spatial averaging, the general problem is reduced to solution of a set of second-order ordinary differential equations whose structure depends on the sorts of nonlinear processes accounted for. This formulation accommodates any physical process; in particular, terms are included to represent noise sources, although only limited modeling is discussed. Our results suggest that random sources of noise have only small effects on combustion instabilities and seem not to be a cause of unstable motions. However, the coupling between the two sorts of unsteady motions may be important as an essential process in a proposed scheme for noise control. It is now a familiar observation that many nonlinear deterministic systems are capable of exhibiting apparently random motions called 'chaos.' This is a particularly interesting possibility for systems which also executed non-deterministic random motions. In combustion chambers, a nonlinear deterministic system (acoustical motions) exists in the presence of noise produced by flow separation, turbulent motions, and energy released by combustion processes. The last part of the paper is directed to the matter of discovering whether or not chaotic motions exist in combustion systems. Analysis has not progressed sufficiently far to answer the question. We report here recent results of processing data taken in one combustor to determine the dimensions of any attractors in the motions. No evidence has been found for chaos in the strict sense, but the method seems to be an important means of investigating the nonlinear behavior of combustion systems.

F. E. C. Culick↗
DOE OSTI2024

Data from "A Bayesian Record Linkage Approach to Applications in Tree Demography Using Overlapping LiDAR Scans"

Processed LiDAR data and environmental covariates from 2015 and 2019 LiDAR scans in the Vicinity of Snodgrass Mountain (Western Colorado, USA), in a geographic subset used in primary analysis for the research paper.This package contains LiDAR-derived canopy height maps for 2015 and 2019, crown polygons derived from the height maps using a segmentation algorithm, and environmental covariates supporting the model of forest growth. Source datasets include August 2015 and August 2019 discrete-return LiDAR point clouds collected by Quantum Geospatial for terrain mapping purposes on behalf of the Colorado Hazard Mapping Program and the Colorado Water Conservation Board. Both datasets adhere to the USGS QL2 quality standard. The point cloud data were processed using the R package lidR to generate a canopy height model representing maximum vegetation height above the ground surface, using a pit-free algorithm.This dataset was compiled to assess how spatial patterns of tree growth in montane and subalpine forests are influenced by water and energy availability. Understanding these growth patterns can provide insight into forest dynamics in the Southern Rocky Mountains under changing climatic conditions.This dataset contains .tif, .csv, and .txt files. This dataset additionally includes a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; and a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1997

External Boundary Conditions for Three-Dimensional Problems of Computational Aerodynamics

We consider an unbounded steady-state flow of viscous fluid over a three-dimensional finite body or configuration of bodies. For the purpose of solving this flow problem numerically, we discretize the governing equations (Navier-Stokes) on a finite-difference grid. The grid obviously cannot stretch from the body up to infinity, because the number of the discrete variables in that case would not be finite. Therefore, prior to the discretization we truncate the original unbounded flow domain by introducing some artificial computational boundary at a finite distance of the body. Typically, the artificial boundary is introduced in a natural way as the external boundary of the domain covered by the grid. The flow problem formulated only on the finite computational domain rather than on the original infinite domain is clearly subdefinite unless some artificial boundary conditions (ABC's) are specified at the external computational boundary. Similarly, the discretized flow problem is subdefinite (i.e., lacks equations with respect to unknowns) unless a special closing procedure is implemented at this artificial boundary. The closing procedure in the discrete case is called the ABC's as well. In this paper, we present an innovative approach to constructing highly accurate ABC's for three-dimensional flow computations. The approach extends our previous technique developed for the two-dimensional case; it employs the finite-difference counterparts to Calderon's pseudodifferential boundary projections calculated in the framework of the difference potentials method (DPM) by Ryaben'kii. The resulting ABC's appear spatially nonlocal but particularly easy to implement along with the existing solvers. The new boundary conditions have been successfully combined with the NASA-developed production code TLNS3D and used for the analysis of wing-shaped configurations in subsonic (including incompressible limit) and transonic flow regimes. As demonstrated by the computational experiments and comparisons with the standard (local) methods, the DPM-based ABC's allow one to greatly reduce the size of the computational domain while still maintaining high accuracy of the numerical solution. Moreover, they may provide for a noticeable increase of the convergence rate of multigrid iterations.

Tsynkov, Semyon V.↗
NASA NTRS2019

Entropy Stable Method for the Euler Equations Revisited: Central Differencing via Entropy Splitting and SBP

The two decades old high order central differencing via entropy splitting and summation-by-parts (SBP) difference boundary closure of Ols- son & Oliger (1994), Gerritsen & Olsson (1996), and Yee et al. (2000) is revisited. The entropy splitting is a form of skew-symmetric splitting of the nonlinear Euler flux derivatives. Central differencing applied to the entropy splitting form of the Euler flux derivatives together with SBP difference operators will, hereafter, be referred to as entropy split schemes. This study is prompted by the recent growing interest in numerical methods for which a discrete entropy conservation law holds, a discrete global entropy conservation can be proved and/or the numerical method possesses a stable entropy in the framework of SBP difference operators and L2-energy norm estimate. The objective of this paper is to recast the entropy split scheme as the re- cent definition of an entropy stable method for central differencing with SBP operators for both periodic and non-periodic boundary conditions for non- linear Euler equations. Standard high order spatial central differencing as well as high order central spatial DRP (dispersion relation preserving) spatial differencing is part of the entropy stable methodology framework. Long time integration of 2D and 3D test cases is included to show the comparison of this efficient entropy stable method with the Tadmor-type of entropy conservative methods. Studies also include the comparison among the three skew-symmetric splittings on their nonlinear stability and accuracy performance without added numerical dissipations for smooth flows. These are, namely, entropy splitting, Ducros et al. splitting and the Kennedy & Grub- ber splitting.

Sjogreen, Bjorn↗
NASA NTRS1983

An investigation of coronal active region loop structures using AS&E rocket X-ray images

Simultaneous high spatial resolution observations at 6 cm in soft X-rays, in photospheric magnetograms, and in optical filtergrams were used to compare the most intense sources of centimetric emission in two active regions to coronal loops, sunspots, chromospheric structures, and photospheric magnetic fields. Results show that the majority of the bright microwave components are not associated with sunspots or X-ray emission. A nonthermal mechanism appears necessary to explain the brightest microwave components, discrete regions of continuous particle acceleration may be common in active regions. Studies of the plasma parameters of selected loops imply that the radio emission is consistent with gyro-resonance absorption at the third and fourth harmonic, at least from part of each loop. Results are presented for: (1) X-ray and microwave observations of active regions; (2) comparison of coronal holes observed in soft X-rays and Hel 10830 A spectrosheliograms; and (3) the reappearance of polar coronal holes and the evolution of the solar magnetic field.

Webb, D. F.↗
NASA NTRS1993

A multi-dimensional kinetic-based upwind solver for the Euler equations

A multidimensional kinetic fluctuation-splitting scheme has been developed for the Euler equations. The scheme is based on an N-scheme discretization of the Boltzmann equation at the kinetic level for triangulated Cartesian meshes with a diagonal-adaptive strategy. The resulting Euler scheme is a cell-vertex fluctuation-splitting scheme where fluctuations in the conserved-variable vector Q are obtained as moments of the fluctuation in the Maxwellian velocity distribution function at the kinetic level. Encouraging preliminary results have been obtained for perfect gases on Cartesian meshes with first-order spatial accuracy. The present approach represents an improvement to the well-established dimensionally-split upwind schemes.

Eppard, W. M.↗
NASA NTRS2016

Organic Analysis in Miller Range 090657 and Buckley Island 10933 CR2 Chondrites: Part 1 In-Situ Observation of Carbonaceous Material

Primitive carbonaceous chondrites contain a wide variety of organic material, ranging from soluble discrete molecules to insoluble unstructured kerogen-like component as well as structured nano-globules of macromolecular carbon. The relationship between the soluble organic molecules, macromolecular organic material, and host minerals are poorly understood. Due to the differences in extractability of soluble and insoluble organic materials, the analysis methods for each differ and are often performed independently. The combination of soluble and insoluble analyses, when performed concurrently, can provide a wider understanding on spatial distribution, and elemental, structural and isotopic composition of organic material in primitive meteorites. Furthermore, they can provide broader perspective on how extraterrestrial organic ma-terials potentially contributed to the synthesis of life's essential compounds such as amino acids, sugar acids, activated phosphates and nucleobases.

Cao, T.↗