Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linearized distribution”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Ideal internal kink stability in presence of plasma flow and neoclassical toroidal viscosity due to energetic particles

Abstract The influence of energetic particles (EPs) on the ideal internal kink mode, in rotating tokamak plasmas, is numerically investigated by simultaneously solving MHD-kinetic hybrid equations together with a toroidal momentum balance equation utilizing the MARS-Q code (Liu et al 2013 Phys. Plasmas 20 042503). The neoclassical toroidal viscous (NTV) torque, induced by precessional drift resonances of trapped energetic particles, acts as the momentum sink term to damp the plasma flow. Quasi-linear initial value simulations show local reduction of the flow amplitude and enhancement of the flow shear near the q = 1 rational surface ( q is the safety factor) due to EP induced NTV. Both effects in turn destabilize the internal kink mode. These numerical findings are robust against the initial linear stability of internal kink, the initial plasma flow profile, as well as the equilibrium distribution model for EPs.

Physics↗
DOE OSTI2022

An efficient method to propagate model uncertainty when inverting seismic data for time domain seismic moment tensors

SUMMARY We present a computationally efficient method to approximately propagate uncertainty when linearly inverting seismic data for point source, time variable moment tensor components. The method is based on the assumption that the data residual, given by the difference between the observed seismic data and the data predicated by a linear inversion, contains the effects of both data and model uncertainty. Our method uses a distribution of data residuals, added directly to the data, in a pseudo-Monte Carlo scheme. Using the assumption that the data residual is a stochastic process, we use the well-known Karhunen–Loève (KL) theorem to construct a distribution of data residuals, where the required basis functions are constructed using Fourier series. The Fourier series are scaled by a product of a random variable and the real-valued spectral amplitudes of the original data residual’s spectrum. Thus, the Fourier series and spectral amplitudes are eigenfunction-eigenvalue pairs used in the KL-based construction of data residual distribution. Using tests with synthetic data, we show that our method compares closely with a Finite Difference Monte Carlo (FDMC) method that we presented previously. More importantly, the method presented here is computationally several orders of magnitude faster than our previous FDMC method, and requires no a priori assumptions of model and/or data uncertainty.

Poppeliers, Christian (ORCID:0000000159526849)↗
DOE OSTI2021

Consistency of cosmic shear analyses in harmonic and real space

ABSTRACT Recent cosmic shear studies have reported discrepancies of up to 1σ on the parameter ${S_{8}=\sigma _{8}\sqrt{{\Omega _{\rm m}}/0.3}}$ between the analysis of shear power spectra and two-point correlation functions, derived from the same shear catalogues. It is not a priori clear whether the measured discrepancies are consistent with statistical fluctuations. In this paper, we investigate this issue in the context of the forthcoming analyses from the third year data of the Dark Energy Survey (DES Y3). We analyse DES Y3 mock catalogues from Gaussian simulations with a fast and accurate importance sampling pipeline. We show that the methodology for determining matching scale cuts in harmonic and real space is the key factor that contributes to the scatter between constraints derived from the two statistics. We compare the published scales cuts of the KiDS, Subaru-HSC, and DES surveys, and find that the correlation coefficients of posterior means range from over 80 per cent for our proposed cuts, down to 10 per cent for cuts used in the literature. We then study the interaction between scale cuts and systematic uncertainties arising from multiple sources: non-linear power spectrum, baryonic feedback, intrinsic alignments, uncertainties in the point spread function, and redshift distributions. We find that, given DES Y3 characteristics and proposed cuts, these uncertainties affect the two statistics similarly; the differential biases are below a third of the statistical uncertainty, with the largest biases arising from intrinsic alignment and baryonic feedback. While this work is aimed at DES Y3, the tools developed can be applied to Stage-IV surveys where statistical errors will be much smaller.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Effects of nuclear structure and quantum interference on diffractive vector meson production in ultra-peripheral nuclear collisions

Here we study diffractive vector meson production in ultra-peripheral collisions (UPCs) of heavy nuclei, utilizing a theoretical framework based on the Color Glass Condensate (CGC) formalism. We focus on Au + Au, U + U, Ru + Ru, Zr + Zr, and Pb + Pb collisions, examining the transverse momentum dependence of vector meson production cross-sections and co s(2ΔΦ)asymmetries in the decay product distributions to explore the role of nuclear geometry. The angular modulation is due to the linear polarization of the incoming photons and quantum interference effects. We extract nuclear radii and find them to be consistent with experimental data from the STAR collaboration. The amplitudes of the cos(2ΔΦ)modulation in the cross-section and the extracted radii depend on the nuclear geometry. This dependence is dominated by the geometry-dependent variation of the minimum impact parameter required for ultra-peripheral collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Gamma Decay of the 154 Sm Isovector Giant Dipole Resonance: Smekal-Raman Scattering as a Novel Probe of Nuclear Ground-State Deformation

𝛾 decays of the isovector giant dipole resonance (IVGDR) of the deformed nucleus 154 Sm were measured using 2$^{+}_{1}$-Smekal-Raman and elastic scattering of linearly polarized, quasimonochromatic photon beams. The two scattering processes were disentangled through their distinct angular distributions. Their branching ratio and cross sections were determined at six excitation energies covering the 154 Sm IVGDR. Both agree with the predictions of the geometrical model for the IVGDR and confirm 𝛾 decay as an observable sensitive to the structure of the resonance. Consequently, the data place strong constraints on the nuclear shape, including the degree of triaxiality. The derived 154 Sm shape parameters 𝛽 = 0.2925⁢(25) and 𝛾 = 5.0⁢(15)° agree well with other measurements and recent Monte Carlo shell-model calculations.

150 ≤ A ≤ 189↗
DOE OSTI2020

Analytic marginalization of N(z) uncertainties in tomographic galaxy surveys

In this paper, we present a new method to marginalize over uncertainties in redshift distributions, N(z), within tomographic cosmological analyses applicable to current and upcoming photometric galaxy surveys. We allow for arbitrary deviations from the best-guess N(z) governed by a general covariance matrix describing the uncertainty in our knowledge of redshift distributions. In principle, this is marginalization over hundreds or thousands of new parameters describing potential deviations as a function of redshift and tomographic bin. However, by linearly expanding the theory predictions around a fiducial model, this marginalization can be performed analytically, resulting in a modified data covariance matrix that effectively downweights the modes of the data vector that are more sensitive to redshift distribution variations. We showcase this method by applying it to the galaxy clustering measurements from the Hyper Suprime-Cam first data release. We illustrate how to marginalize over sample variance of the calibration sample and a large general systematic uncertainty in photometric estimation methods, and explore the impact of priors imposing smoothness in the redshift distributions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Analyzing the Validity of Brazilian Testing Using Digital Image Correlation and Numerical Simulation Techniques

Characterizing the mechanical behavior of rocks plays a crucial role to optimize the fracturing process in unconventional reservoirs. However, due to the intrinsic anisotropy and heterogeneity in unconventional resources, fracture process prediction remains the most significant challenge for sustainable and economic hydrocarbon production. During the deformation tracking under compression, deploying conventional methods (strain gauge, extensometer, etc.) is insufficient to measure the deformation since the physical attachment of the device is restricted to the size of the sample, monitoring limited point-wise deformation, producing difficulties in data retrieval, and a tendency to lose track in failure points, etc. Where conventional methods are limited, the application of digital image correlation (DIC) provides detailed and additional information of strain evolution and fracture patterns under loading. DIC is an image-based optical method that records an object with a camera and monitors the random contrast speckle pattern painted on the facing surface of the specimen. To overcome the existing limitations, this paper presents numerical modeling of Brazilian disc tests under quasi-static conditions to understand the full-field deformation behaviors and finally, it is validated by DIC. As the direct tensile test has limitations in sample preparation and test execution, the Brazilian testing principle is commonly used to evaluate indirectly the tensile strength of rocks. The two-dimensional numerical model was built to predict the stress distribution and full-field deformation on Brazilian disc under compression based on the assumptions of a homogenous, isotropic and linear elastic material. The uniaxial compression test was conducted using the DIC technique to determine the elastic properties of Spider Berea sandstone, which were used as inputs for the simulation model. The model was verified by the analytical solution and compared with the digital image correlation. The numerical simulation results showed that the solutions matched reasonably with the analytical solutions where the maximum deviation of stress distribution was obtained as 14.59%. The strain evolution (normal and shear strains) and displacements along the central horizontal and vertical planes were investigated in three distinguishable percentages of peak loads (20%, 40%, and 90%) to understand the deformation behaviors in rock. The simulation results demonstrated that the strain evolution contours consistently matched with DIC generated contours with a reasonable agreement. The changes in displacement along the central horizontal and vertical planes showed that numerical simulation and DIC generated experimental results were repeatable and matched closely. In terms of validation, Brazilian testing to measure the indirect tensile strength of rocks is still an issue of debate. The numerical model of fracture propagation supported by digital image correlation from this study can be used to explain the fracturing process in the homogeneous material and can be extended to non-homogeneous cases by incorporating heterogeneity, which is essential for rock mechanics field applications.

58 GEOSCIENCES↗
DOE OSTI2023

Bead–Spring Simulation of Ionomer Melts—Studying the Effects of Chain-Length and Associating Group Fraction on Equilibrium Structure and Extensional Flow Behavior

Ionomers are associative polymers with diverse applications ranging from selective membranes and high-performance adhesives to abrasion- and chemical-resistant coatings, insulation layers, vacuum packaging, and foamed sheets. Within equilibrium melt, the ionic or associating groups are known to form thermally reversible, associative clusters whose presence can significantly affect the system’s mechanical, viscoelastic, and transport properties. It is, thus, of great interest to understand how to control such clusters’ size distribution, shape, and stability through the designed choice of polymer architecture and the ionic groups’ fraction, arrangement, and interaction strength. In this work, we represent linear associating polymers using a Kremer–Grest type bead–spring model and perform large-scale MD simulations to explore the effect of polymer chain-length (l) and fraction (f s ) of randomly placed associating groups on the size distribution and stability of formed clusters. We consider different chain-lengths (below and above entanglement), varying fractions of associating groups (represented by ‘sticky’ beads) between 5 and 20%, and a fixed sticky–sticky nonbond interaction strength of four times that between regular non-associating beads. For all melts containing associating groups the equilibrium structure factor S(q) displays a signature ionomer peak at low wave vector q whose intensity increases with increasing f s and l. The average cluster size Nc increases with f s . However, the effect of chain-length on Nc appears to be pronounced only at higher values of f s . Under extensional flows, the computed stress (and viscosity) is higher at higher f s and l regardless of strain rate. Beyond a critical strain rate, we observe fragmentation of the associative clusters, which has interesting effects on the stress/viscous response.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Flexible dynamic boundary microgrid operation considering network and load unbalances

Flexible microgrids with dynamic boundaries have recently been introduced in the literature. With the ability to reconfigure the topology of the microgrids dynamically through remotely controlled switches, flexible microgrids with dynamic boundaries can further improve the resiliency and energy efficiency of microgrids with distributed energy resources (DERs). This paper focuses on the optimal operation considering one of the predominant characteristics of microgrids and distribution systems – unbalanced networks and loads. In existing literature, balanced modeling of microgrids is more common due to its attractive simplicity. The three-phase power unbalance has not been considered as a constraint on the generation units in a microgrid. Further, negative sequence constraints have also been neglected. In this article, we propose a set of constraints that is specifically related to the capabilities of inverter interfaced resources to supply unbalanced current/power when the microgrid is islanded from the main distribution grid. We incorporate the new set of constraints into two optimization formulations leveraging two convex relaxations of the three-phase power flow equations: mixed-integer linear programming (MILP) and mixed-integer semidefinite programming (MISDP) that optimize the dispatch of controllable switches and DERs in the microgrid. The algorithms are then extended to networked microgrids with grid-forming sources. We test the algorithms on a realistic community microgrid model in Puerto Rico as well as standardized IEEE distribution test feeders. The testing results demonstrate the performance of the proposed algorithms. The MILP is fast and scalable, and the MISDP enforces the negative sequence voltage constraints.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

A tri-level distribution locational marginal price-based demand response framework

Here, in this paper, we propose a tri-level, nested, two-stage price-based demand response (PBDR) framework that considers distribution locational marginal price (DLMP) as DR enabler between load-serving entities (LSE), demand response providers (DRPs), and customers in the day-ahead distribution market. It enables LSE and customer interactions by using multiple DRPs, positioned in-between, and independently optimizes their objectives. The problem is formulated using linear power flow with approximated power losses and its application in DLMP as DR pricing. The tri-level problem is solved using a nested reformulation & decomposition (R&D) method and tested on the real Indian-108 bus distribution system under various dynamic pricings. Further, the temporal–spatial variations in DLMPs are assessed using fairness criteria. Numerical analyses demonstrate that DLMP applications can effectively improve economic efficiency, and transparency in DR programs valuation with a favorable fairness margin. The results show that DLMP as DR pricing signal induces (0-2) % variation in DLMP for DR participation up to 10 %. Further, it gives over 90 % fairness over temporal–spatial variation for all the customers.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Rapid subsurface analysis of frequency-domain thermoreflectance images with K-means clustering

K-means clustering analysis is applied to frequency-domain thermoreflectance (FDTR) hyperspectral image data to rapidly screen the spatial distribution of thermophysical properties at material interfaces. Performing FDTR while raster scanning a sample consisting of 8.6 μm of doped-silicon (Si) bonded to a doped-Si substrate identifies spatial variation in the subsurface bond quality. Routine thermal analysis at select pixels quantifies this variation in bond quality and allows assignment of bonded, partially bonded, and unbonded regions. Performing this same routine thermal analysis across the entire map, however, becomes too computationally demanding for rapid screening of bond quality. To address this, K-means clustering was used to reduce the dimensionality of the dataset from more than 20 000 pixel spectra to just K = 3 component spectra. The three component spectra were then used to express every pixel in the image through a least-squares minimized linear combination providing continuous interpolation between the components across spatially varying features, e.g., bonded to unbonded transition regions. Fitting the component spectra to the thermal model, thermal properties for each K cluster are extracted and then distributed according to the weighting established by the regressed linear combination. Thermophysical property maps are then constructed and capture significant variation in bond quality over 25 μm length scales. The use of K-means clustering to achieve these thermal property maps results in a 74-fold speed improvement over explicit fitting of every pixel.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Hybrid Simulation and Quasi-linear Theory of Bi-Kappa Proton Instabilities

The quasi-steady states of collisionless plasmas in space (e.g., in the solar wind and planetary environments) are governed by the interactions of charged particles with wave fluctuations. These interactions are responsible not only for the dissipation of plasma waves but also for their excitation. The present analysis focuses on two instabilities, mirror and electromagnetic ion cyclotron instabilities, associated with the same proton temperature anisotropy T ⊥ > T ∥ (where ⊥, ∥ are directions defined with respect to the local magnetic field vector). Theories relying on standard Maxwellian models fail to link these two instabilities (i.e., predicted thresholds) to the proton quasi-stable anisotropies measured in situ in a completely satisfactory manner. Here we revisit these instabilities by modeling protons with the generalized bi-Kappa (bi-κ power-law) distribution, and by a comparative analysis of a 2D hybrid simulation with the velocity-moment-based quasi-linear (QL) theory. It is shown that the two methods feature qualitative and, even to some extent, quantitative agreement. The reduced QL analysis based upon the assumption of a time-dependent bi-Kappa model thus becomes a valuable theoretical approach that can be incorporated into the present studies of solar wind dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Robust scalable initialization for Bayesian variational inference with multi-modal Laplace approximations

Predictive modeling typically relies on Bayesian model calibration to provide uncertainty quantification. Variational inference utilizing fully independent (“mean-field”) Gaussian distributions are often used as approximate probability density functions. This simplification is attractive since the number of variational parameters grows only linearly with the number of unknown model parameters. However, the resulting diagonal covariance structure and unimodal behavior can be too restrictive to provide useful approximations of intractable Bayesian posteriors that exhibit highly non-Gaussian behavior, including multimodality. High-fidelity surrogate posteriors for these problems can be obtained by considering the family of Gaussian mixtures. Gaussian mixtures are capable of capturing multiple modes and approximating any distribution to an arbitrary degree of accuracy, while maintaining some analytical tractability. Unfortunately, variational inference using Gaussian mixtures with full-covariance structures suffers from a quadratic growth in variational parameters with the number of model parameters. The existence of multiple local minima due to strong nonconvex trends in the loss functions often associated with variational inference present additional complications, These challenges motivate the need for robust initialization procedures to improve the performance and computational scalability of variational inference with mixture models. In this work, we propose a method for constructing an initial Gaussian mixture model approximation that can be used to warm-start the iterative solvers for variational inference. The procedure begins with a global optimization stage in model parameter space. In this step, local gradient-based optimization, globalized through multistart, is used to determine a set of local maxima, which we take to approximate the mixture component centers. Around each mode, a local Gaussian approximation is constructed via the Laplace approximation. Finally, the mixture weights are determined through constrained least squares regression. The robustness and scalability of the proposed methodology is demonstrated through application to an ensemble of synthetic tests using high-dimensional, multimodal probability density functions. Here, the practical aspects of the approach are demonstrated with inversion problems in structural dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Goal-oriented a-posteriori estimation of model error as an aid to parameter estimation

In this work, a Bayesian model calibration framework is presented that utilizes goal-oriented a-posterior error estimates in quantities of interest (QoIs) for classes of high-fidelity models characterized by PDEs. It is shown that for a large class of computational models, it is possible to develop a computationally inexpensive procedure for calibrating parameters of high-fidelity models of physical events when the parameters of low-fidelity (surrogate) models are known with acceptable accuracy. The main ingredients in the proposed model calibration scheme are goal-oriented a-posteriori estimates of error in QoIs computed using a so-called lower fidelity model compared to those of an uncalibrated higher fidelity model. The estimates of error in QoIs are used to define likelihood functions in Bayesian inversion analysis. A standard Bayesian approach is employed to compute the posterior distribution of model parameters of high-fidelity models. As applications, parameters in a quasi-linear second-order elliptic boundary-value problem (BVP) are calibrated using a second-order linear elliptic BVP. In a second application, parameters of a tumor growth model involving nonlinear time-dependent PDEs are calibrated using a lower fidelity linear tumor growth model with known parameter values.

A-posterior estimates↗
DOE OSTI2025

High-order limiting methods using maximum principle bounds derived from the Boltzmann equation I: Euler equations

The use of limiting methods for high-order numerical approximations of hyperbolic conservation laws generally requires defining an admissible region/bounds for the solution. In this work, we present a novel approach for computing solution bounds and limiting for the Euler equations through the kinetic representation provided by the Boltzmann equation, which allows for extending limiters designed for linear advection directly to the Euler equations. Given an arbitrary set of solution values to compute bounds over (e.g., numerical stencil) and a desired linear advection limiter, the proposed approach yields an analytic expression for the admissible region of particle distribution function values, which may be numerically integrated to yield a set of bounds for the density, momentum, and total energy. Further, these solution bounds are shown to preserve positivity of density/pressure/internal energy and, when paired with a limiting technique, can robustly resolve strong discontinuities while recovering high-order accuracy in smooth regions without any ad hoc corrections (e.g., relaxing the bounds). This approach is demonstrated in the context of an explicit unstructured high-order discontinuous Galerkin/flux reconstruction scheme for a variety of difficult problems in gas dynamics, including cases with extreme shocks and shock-vortex interactions. Furthermore, this work presents a foundation for limiting techniques for more complex macroscopic governing equations that can be derived from an underlying kinetic representation for which admissible solution bounds are not well-understood.

42 ENGINEERING↗
DOE OSTI2020

Analysis of hematite attrition in a grid jet apparatus

Particulate attrition is of interest for novel carbon-capture processes such as chemical looping combustion because the makeup cost of oxygen carrier is a significant portion of operating cost. As such, models to study and predict attrition of various oxygen carriers in fluidized bed systems are being developed. One of the regions of concern in fluidized bed systems is the high-velocity jet region near gas distributors in a fluid bed. This work studies the attrition of hematite particles using a modified ASTM apparatus to measure the particle size distribution throughout the experiment. Bed weight and gas velocity were varied. Hematite particles above the corresponding threshold value had a severe variation of particle size distributions which decreased with time. Weight fractions of the sieves over time were fit to a linear, time-variant population balance model to offer insight into particle attrition. The first-order rate constant was modified as a decaying exponential.

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

Asymmetries in ionization of atomic superposition states by ultrashort laser pulses

Abstract Progress in ultrafast science allows for probing quantum superposition states with ultrashort laser pulses in the new regime where several linear and nonlinear ionization pathways compete. Interferences of pathways can be observed in the photoelectron angular distribution and in the past they have been analyzed for atoms and molecules in a single quantum state via anisotropy and asymmetry parameters. Those conventional parameters, however, do not provide comprehensive tools for probing superposition states in the emerging research area of bright and ultrashort light sources, such as free-electron lasers and high-order harmonic generation. We propose a new set of generalized asymmetry parameters which are sensitive to interference effects in the photoionization and the interplay of competing pathways as the laser pulse duration is shortened and the laser intensity is increased. The relevance of the parameters is demonstrated using results of state-of-the-art numerical solutions of the time-dependent Schrödinger equation for ionization of helium atom and neon atom.

74 ATOMIC AND MOLECULAR PHYSICS↗