Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HOMOGENIZATION 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Change-point Detection and Image Segmentation for Time Series of Astrophysical Images

Many astrophysical phenomena are time-varying, in the sense that their intensity, energy spectrum, and/or the spatial distribution of the emission suddenly change. This paper develops a method for modeling a time series of images. Under the assumption that the arrival times of the photons follow a Poisson process, the data are binned into 4D grids of voxels (time, energy band, and x-y coordinates), and viewed as a time series of non-homogeneous Poisson images. The method assumes that at each time point, the corresponding multiband image stack is an unknown 3D piecewise constant function including Poisson noise. It also assumes that all image stacks between any two adjacent change points (in time domain) share the same unknown piecewise constant function. The proposed method is designed to estimate the number and the locations of all of the change points (in time domain), as well as all of the unknown piecewise constant functions between any pairs of the change points. The method applies the minimum description length principle to perform this task. A practical algorithm is also developed to solve the corresponding complicated optimization problem. Simulation experiments and applications to real data sets show that the proposed method enjoys very promising empirical properties. Applications to two real data sets, the XMM observation of a flaring star and an emerging solar coronal loop, illustrate the usage of the proposed method and the scientific insight gained from it.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1974

First-order reactant in homogeneous turbulence before the final period of decay

A method is described for studying theoretically the concentration fluctuations of a dilute contaminate undergoing a first-order chemical reaction. The method is based on Deissler's (1958) theory for homogeneous turbulence for times before the final period, and it follows the approach used by Loeffler and Deissler (1961) to study temperature fluctuations in homogeneous turbulence. Four-point correlation equations are obtained; it is assumed that terms containing fifth-order correlation are very small in comparison with those containing fourth-order correlations, and can therefore be neglected. A spectrum equation is obtained in a form which can be solved numerically, yielding the decay law for the concentration fluctuations in homogeneous turbulence for the period much before the final period of decay.

Kumar, P.↗
DOE OSTI2020

Investigation of finite-volume methods to capture shocks and turbulence spectra in compressible flows

The aim of the present paper is to provide a comparison between several finite-volume methods of different numerical accuracy: the second-order Godunov method with PPM interpolation and the high-order finite-volume WENO method. In this work, the results show that while on a smooth problem the high-order method performs better than the second-order one, when the solution contains a shock all the methods collapse to first-order accuracy. In the context of the decay of compressible homogeneous isotropic turbulence with shocklets, the actual overall order of accuracy of the methods reduces to second-order, despite the use of fifth-order reconstruction schemes at cell interfaces. Most important, results in terms of turbulent spectra are similar regardless of the numerical methods employed, except that the PPM method fails to provide an accurate representation in the high-frequency range of the spectra. It is found that this specific issue comes from the slope-limiting procedure and a novel hybrid PPM/WENO method is developed that has the ability to capture the turbulent spectra with the accuracy of a high-order method, but at the cost of the second-order Godunov method. Overall, it is shown that virtually the same physical solution can be obtained much faster by refining a simulation with the second-order method and carefully chosen numerical procedures, rather than running a coarse high-order simulation. Our results demonstrate the importance of evaluating the accuracy of a numerical method in terms of its actual spectral dissipation and dispersion properties on mixed smooth/shock cases, rather than by the theoretical formal order of convergence rate.

97 MATHEMATICS AND COMPUTING↗
NASA NTRSDate not supplied

Homogenized Ground-Based and Profile Ozone Datasets From the TOAR-II/HEGIFTOM Project: Methods and Station Trends

Within the framework of the second phase of the Tropospheric Ozone Assessment Report (TOAR-II), it was recognized that an essential first step for deriving accurate trends from the ground-based networks that monitor ozone in the free troposphere is putting the measurements on the same basis with respect to absolute references and processing methods. The relevant procedures are referred to as “harmonization” or “homogenization”. The TOAR II working group, “HEGIFTOM” (Harmonization and Evaluation of Ground-based Instruments for Free-Tropospheric Ozone Measurements), has carried out harmonization for five types of network (Figure below) instruments (mid-1990s to 2020): ozonesondes, commercial aircraft IAGOS landing/takeoff profiles, Fourier-Transform Infrared spectrometer (FTIR), tropospheric Lidar, and Brewer/Dobson Umkehr. First, we summarize the homogenization effort for each network that provides new quality-assessed ozone profile or segment (partial column) data sets, including uncertainty estimates and quality flags. Second, with the HEGIFTOM datasets forming the basis for a global assessment of tropospheric ozone column trends, results derived with various trend detection algorithms will be presented. The ultimate goal is evaluation of the consistency of the calculated trends among different techniques at selected stations and/or regions.

ozone↗
DOE OSTI2024

Fractionation of UF6 and daughter progeny in storage cylinders from external heating

Abstract Uranium hexafluoride (UF 6 ) is a significant concern for material accountancy and verification in the international safeguards community. Verification of the contents of UF 6 cylinders is generally attempted with gamma spectroscopy but the current methods assume a uniform, homogeneous UF 6 mass distribution within the cylinder. In this work, it was found experimentally and confirmed via modeling, that under an external heat load (the sun), the UF 6 and its daughter products undergo fractionation in the cylinder. This fractionation of the UF 6 and daughter products can cause an errant measurement of the enrichment of the cylinder when using the current verification methods.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Towards dynamical low-rank approximation for neutrino kinetic equations. Part I: Analysis of an idealized relaxation model

Dynamical low-rank approximation (DLRA) is an emerging tool for reducing computational costs and provides memory savings when solving high-dimensional problems. Here, in this work, we propose and analyze a semi-implicit dynamical low-rank discontinuous Galerkin (DLR-DG) method for the space homogeneous kinetic equation with a relaxation operator, modeling the emission and absorption of particles by a background medium. Both DLRA and the discontinuous Galerkin (DG) scheme can be formulated as Galerkin equations. To ensure their consistency, a weighted DLRA is introduced so that the resulting DLR-DG solution is a solution to the fully discrete DG scheme in a subspace of the standard DG solution space. Similar to the standard DG method, we show that the proposed DLR-DG method is well-posed. We also identify conditions such that the DLR-DG solution converges to the equilibrium. Numerical results are presented to demonstrate the theoretical findings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Simulation of gas mixture dynamics in a pipeline network using explicit staggered-grid discretization

Here we develop an explicit staggered finite difference discretization scheme for simulating the transport of highly heterogeneous gas mixtures through pipeline networks. This study is motivated by the proposed blending of hydrogen into natural gas pipelines to reduce end use carbon emissions while using existing pipeline systems throughout their planned lifetimes. Our computational method accommodates an arbitrary number of constituent gases with very different physical properties that may be injected into a network with significant spatiotemporal variation. In this setting, the gas flow physics are highly location- and time- dependent, so that local composition and nodal mixing must be accounted for. The resulting conservation laws are formulated in terms of pressure, partial densities and flows, and volumetric and mass fractions of the constituents. We include non-ideal equations of state that employ linear approximations of gas compressibility factors, so that the pressure dynamics propagate locally according to a variable wave speed that depends on mixture composition and density. We derive compatibility relationships for network edge boundary values that are more complex than for a homogeneous gas. The simulation method is evaluated on initial boundary value problems for a single pipe and a small network, is cross-validated with a lumped element simulation, and used to demonstrate a local monitoring and control policy for maintaining allowable concentration levels.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Particulate Fuel Modeling of MC 2 -3 using Iterative Local Spatial Self-shielding Method

We report a new spatial self-shielding method for particulate fuels has been developed based on disadvantage factors and implemented in the MC 2 -3 code. This method named the iterative local spatial self-shielding (ILSS) method considers the shadowing effect of randomly distributed particles on spatial self-shielding in particles through a homogenized composition region added outside the particle of interest at the center. The self-shielded cross sections of the central particle are determined iteratively since they are used in determining the cross sections of the homogenized composition region. The ILSS method was verified for infinite stochastic medium problems of single and multiple types of particles, VHTR unit cell problems, and HTTR assembly problems. The verification test results show that the ILSS method accurately predicts the stochastic particle shadowing effect and reaction rates in particles, whereas the regular array model and the stochastic collision probability method underpredict the particle shadowing effect and overestimate reaction rates in particles.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1989

'Nose method' of calculating critical cooling rates for glass formation

The use of the so-called 'nose method' for computing critical cooling rates for glass formation is examined and compared with other methods, presenting data for the glass-forming systems SiO2, GeO2, and P2O5. It is shown that, for homogeneous crystallization, the nose-method will give an overestimate of Rc, a conclusion which was drawn after assessing the enfluence of a range of values for the parameters which control crystal growth and nucleation. The paper also proposes an alternative simple procedure (termed the 'cutoff method') for computing critical cooling rates from T-T-T diagrams, which was shown in the SiO2 and GeO2 systems to be superior to the nose method.

Weinberg, Michael C.↗
DOE OSTI2022

Effect of Grain Boundary Misorientation on Spall Strength in Ta via Shock-Free Simulations with Relatively Few Atoms

A suite of 37 molecular dynamics simulations is conducted at two system sizes to systematically characterize the role of grain boundary (GB) misorientation on spall strength in pure BCC tantalum (Ta). The systems studied consist of bicrystals with a single [110] symmetric tilt grain boundary. Two loading conditions are compared: (i) homogeneous extension under uniaxial strain simulated in this study and (ii) piston/flyer impact of sample, which induces heterogeneous deformation via shockwave propagation along the length of the sample. The piston/flyer impact is taken from the literature and run on the same set of GB misorientation angles using LAMMPS. The major finding here is that both methods result in similar spall strength predictions, but the homogeneous extension method generally requires two to three orders of magnitude fewer atoms and similar reductions in computational costs. Spall strength results systematically overpredict using this method, by about 10% for the dataset three orders of magnitude smaller than piston/flyer simulations, and 5% for the dataset two orders of magnitude smaller. Lastly, the effect of system size and pre-compression magnitude on spall strength is systematically characterized.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Planar Systems for Quantum Information

This project aims to develop two‐dimensional (2D) moiré materials as a quantum simulator to implement model Hamiltonians and their phase diagrams. Progress in quantum information science (QIS) requires the development of advanced quantum materials systems. The rich family of layered van der Waals materials and their heterostructures present opportunities to create previously unrealized types of applications for QIS. Specifically, when two layers of van der Waals materials are overlaid with a small twist angle or/and lattice mismatch, a moiré superlattice with a period of about ten nanometers is formed. This provides a periodic trapping potential for electrons. Electrons can tunnel between the traps and repel each other by their mutual Coulomb interactions. The platform of 2D moiré materials provides many attractive features, including tunability of length and energy scales, charge density, and even lattice symmetry. It presents new possibilities for realizing quantum simulation of the many-body physics in a solid-state platform. This integrated team of six investigators seeks to develop relevant theoretical treatments to link ab-initio studies of 2D moiré materials to model Hamiltonians and to evaluate correlated phases predicted by these model Hamiltonians in the relevant regimes. On the experimental side, the team aims to develop methods to realize a homogeneous and highly controlled potential landscape for the electrons and to initiate, protect, and measure their quantum many-body states.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Influences of Alkali Metal Cation Interactions on N2 Coordination to Iron(I) Beta-Diketiminate Complexes in Electrolyte Solutions

N2 functionalization reactions often use stable, isolable N2-bound complexes and low temperatures or highly activated reagents. However, little work has focused on systems where N2 binds weakly, and the N2 binding equilibrium can be modulated. Here, we describe analysis of surprising trends in alkali metal cation solvation on the formation of an anionic iron(I) beta-diketiminate complex and the subsequent N2 binding equilibrium. Variable temperature UV-vis spectroscopy shows that increased ionic strength in the coordinating solvent tetrahydrofuran promotes N2 coordination. Additionally, less solvated alkali metal cations in noncoordinating 2-methyltetrahydrofuran (MeTHF) stabilize N2 coordination, likely through the formation of contact ion pairs in solution that have enhanced N2 binding. Ionic strength has little effect on the energetics of N2 coordination to these proposed contact ion pairs in MeTHF. When the electrochemical analysis is done at –78 °C, there is an irreversible reduction of iron(II) and an anodically shifted oxidation, which is attributed to the formation of the spectroscopically observed N2 complex. Finally, bulk electrolysis is used to verify the electrochemical formation of the N2-bound anionic iron(I) complex. These studies demonstrate a new strategy for promoting N2 coordination by changing solution properties, which is relevant to homogeneous electrochemical N2 reduction methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1974

Development and creation of conditions for the radiation of biological objects with high energy protons

Depth distribution of an absorbed proton radiation dose was studied by both calculation and experimental methods on a cylindrical, homogeneous plexiglass phantom that represented the geometric form of a dog. The basis of the heterogeneous phantom was the natural skeleton of a dog, muscle tissue was simulated by a mixture of paraffin and silica filler, and lung tissue by gelatine capsules.

Popov, V. P.↗
NASA NTRS1976

Transfer matrix approach to layered systems with axial symmetry

The stress and displacement distribution in a layered medium is found by means of transfer matrices. The surface loading exhibits axial symmetry, and each layer is of infinite extent in the horizontal direction, of constant depth, and is considered to be linearly elastic, homogeneous, and isotropic. The method developed has the built-in advantage of enforcing interface continuity conditions automatically. Its application to layered composites shows the flexibility with which it predicts the local as well as the global response of the medium.

Bahar, L. Y.↗
NASA NTRS1985

Application of singular value decomposition to structural dynamics systems with constraints

Singular value decomposition is used to construct a coordinate transformation for a linear dynamic system subject to linear, homogeneous constraint equations. The method is compared with two commonly used methods, namely classical Gaussian elimination and Walton-Steeves approach. Although the classical method requires fewer numerical operations, the singular value decomposition method is more accurate and convenient in eliminating the dependent coordinates. Numerical examples are presented to demonstrate the application of the method.

Juang, J.-N.↗
NASA NTRS1983

Robot arm dynamic model reduction for control

General methods are described by which the mathematical complexities of explicit and exact state equations of robot arms can be reduced to a simplified and compact state equation representation without introducing significant errors into the robot arm dynamic model. The model reduction methods are based on homogeneous coordinates and on the Langrangian algorithm for robot arm dynamics, and utilize matrix, vector and numeric analysis techniques. The derivation of differential vector representation of centripetal and Coriolis forces which has not yet been established in the literature is presented.

Bejczy, A. K.↗