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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Pure quark and gluon observables in collinear drop

We construct a class of pure quark and gluon observables by using the collinear drop grooming technique. The construction is based on linear combinations of multiple cumulative distributions of the jet mass in collinear drop, whose specific weights are fully predicted perturbatively. This yields observables which obtain their values purely from quarks (or purely from gluons) in a wide region of phase space. We demonstrate this by showing that these observables are effective in two phase space regions, one dominated by perturbative resummation and one dominated by nonperturbative effects. The nonperturbative effects are included using shape functions which only appear as a common factor in the linear combinations constructed. We test this construction using a numerical analysis with next-to-leading logarithmic resummation and various shape function models, as well as analyzing these observables with PYTHIA and VINCIA. Choices for the collinear drop parameters are optimized for experimental use.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Phase Demodulation by Frequency Chirping in Coherence Microwave Photonic Interferometry

This paper presents a signal processing method to demodulate the optical interference phase of cascaded individual optical fiber intrinsic Fabry-Perot interferometric (IFPI) sensors in a coherent microwave-photonic interferometry (CMPI) distributed sensing system. The new method utilizes the chirp effect of electro-optic modulator (EOM) to create a quasi-quadrature optical interference phase shift between two adjacent pulses which correspond to two adjacent reflection points in the time domain. The phase shift can be controlled by adjusting the bias voltage that is applied to the EOM. The interference phase is calculated by elliptically fitting the phase shift. The interference phase change is proportional to the optical path difference (OPD) change of the interferometer, and the sign can be used to differentiate the increase or decrease of the OPD. The method is demonstrated for distributed strain sensing, showing good linearity, high resolution and large dynamic range.

42 ENGINEERING↗
DOE OSTI2023

Bayesian Monte Carlo Evaluation of Imperfect (n, 233 U) Data and Model

Conventional nuclear data evaluation methods using generalized linear least squares make the following assumptions: prior and posterior probability distribution functions (PDFs) of all model parameters and data are normal (Gaussian); the linear approximation is sufficiently accurate to minimize the cost function (even for nonlinear models); the model (e.g., of neutron cross section) and experimental data (including covariance data) are without defect and prior PDFs of parameters and measured data are known perfectly. Neglect of covariance between model parameters and measured data in conventional evaluations contributes to imperfections. These assumptions are inherent to the generalized linear least squares minimization method commonly used for resolved resonance region neutron cross section evaluations but are often not justified due to the presence of non-normal PDFs, nonlinear models (e.g., R-matrix formalism), and inherent imperfections in data and models (e.g., imperfect covariance data). Here, these assumptions are removed in a mathematical framework of Bayes’ theorem, which is implemented using the Metropolis-Hastings Monte Carlo method. Most importantly, new parameters are introduced to parameterize discrepancies between the theoretical model and measured data to quantify judgement about discrepancies or imperfections in a reproducible manner. An evaluation of 233U in the eV region using the ENDF-B/VIII.0 library and transmission data (Guber et al.) is presented, and posterior parameters are compared to those obtained by conventional evaluation methods. This example illustrates the effects of removing the most harmful assumption: that of model-data perfection.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

A Search for Polarized Thermal Emission from Directly Imaged Exoplanets and Brown Dwarf Companions to Nearby Stars

Aerosols in the atmospheres of cloudy gas giant exoplanets and brown dwarfs scatter and polarize these objects’ thermal emission. If such an object has an oblate shape or nonuniform cloud distribution, the net degree of linear polarization can show an increase ranging from several tenths of a percent to a few percent. Modern high-contrast imaging polarimeters are now poised to detect such low-polarization signals, opening up a new window into the rotational velocities and cloud properties of substellar companions to nearby stars. In this paper, we present the results of a near-IR survey searching for linearly polarized thermal emission from a sample of two planetary-mass companions and five brown dwarf companions using GPI and SPHERE-IRDIS. We probe the subpercent linear polarization regime that typifies polarized free-floating brown dwarfs and place limits on each object’s degree of linear polarization. We relate our upper limits on each target’s degree of linear polarization to its rotation rate, and place our results in the context of rotation rates measured using high-resolution spectroscopy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Nonlinear kinetic simulation study of the ion–ion streaming instability in single- and multi-ion species plasmas

The nonlinear evolution of the ion–ion streaming instability (IISI) is studied using numerical techniques novel to this problem that afford direct insight into the evolution of the particle distributions of each species. Here, during the linear phase of the instability, we demonstrate quantitative agreement with linear kinetic theory. Subsequently, the electrostatic field generated by the IISI causes ring-like velocity distributions of ions to form that are both heated and slowed to varying degrees relative to their initial flows. Due to variation in the trapping conditions for ion species of differing charge-to-mass ratio, when flows of multiple species interact, the nonlinear evolution of each species can be starkly different: we show a case where a lighter ion species is completely stopped by a heavier ion species via the IISI alone (i.e., without collisions) and, for the first time, demonstrate how the IISI can introduce a relative flow between ion species that initially have the same flow velocities, thereby separating them.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Impacts of Experimentally Obtained Harmonic Spectrums of Residential Appliances on Distribution Feeder: Preprint

This paper focuses on the harmonic load flow study for a set of residential loads and the impacts of the harmonic contents on the distribution feeder. Owing to the increased use of the power electronic based appliances and energy efficient home equipment like the lighting loads, the percentage of non-linear loads has increased in the distribution system which can impact the system performance, loss, and stability. Thus a comprehensive knowledge of their power quality and harmonic analysis is essential to improve the load models, voltage stability assessment, determining possible interactions at harmonic frequencies, protection planning, and the effect of the system impedance. For this study, the household appliances include the lighting, power electronics, resistive, and motor loads with a nominal supply voltage of 120 V single-phase or 240 V split-phase at a nominal frequency of 60 Hz. The harmonic spectrum as obtained from the experimental data is used to solve the harmonic load flow for a detailed secondary network including the distribution transformer extracted from a real distribution feeder's model. Additionally, the impact of the voltage harmonic on a three phase test load is presented through experimental results. The harmonic data of the transformer terminal voltage as obtained from the former load flow study is scaled to generate the source voltage for the three phase configuration of the grid simulator to make the test setup very similar to a real scenario.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

A Particle Method for the Multispecies Landau Equation

Abstract The multispecies Landau collision operator describes the two-particle, small scattering angle or grazing collisions in a plasma made up of different species of particles such as electrons and ions. Recently, a structure preserving deterministic particle method (Carrillo et al. in J. Comput. Phys. 7:100066, 2020) has been developed for the single species spatially homogeneous Landau equation. This method relies on a regularization of the Landau collision operator so that an approximate solution, which is a linear combination of Dirac delta distributions, is well-defined. Based on a weak form of the regularized Landau equation, the time dependent locations of the Dirac delta functions satisfy a system of ordinary differential equations. In this work, we extend this particle method to the multispecies case, and examine its conservation of mass, momentum, and energy, and decay of entropy properties. We show that the equilibrium distribution of the regularized multispecies Landau equation is a Maxwellian distribution, and state a critical condition on the regularization parameters that guarantees a species independent equilibrium temperature. A convergence study comparing an exact multispecies Bobylev-Krook-Wu (BKW) solution to the particle solution shows approximately 2nd order accuracy. Important physical properties such as conservation, decay of entropy, and equilibrium distribution of the particle method are demonstrated with several numerical examples.

Mathematics↗
DOE OSTI2024

Operando visualization of porous metal additive manufacturing with foaming agents through high-speed x-ray imaging

Porous metals find extensive applications in soundproofing, filtration, catalysis, and energy-absorbing structures, thanks to their unique internal pore structure and high specific strength. In recent years, there has been an increasing interest in fabricating porous metals using additive manufacturing (AM), leveraging its unique advantages, including improved design freedom, spatial material control, and cost-effective small-batch production. In this study, we conducted pioneering operando visualization of AM porous metal using a laser powder bed fusion (L-PBF) setup combined with a high-speed synchrotron x-ray imaging system. Single track printing experiments using Ti6Al4V (Ti64) combined with titanium hydride (TiH 2 ) and sodium carbonate (Na 2 CO 3 ) as foaming agents, with varying mixing ratios were performed under different processing conditions. Here. the results elucidate the dynamic development of porosity formation. The average pore size is significantly influenced by the particle size of foaming agents when pore coalescence is absent. For all foaming agent content tested in the current study, the number of pores is found to be more sensitive to changes in laser power than in laser scanning speed. Increasing linear energy density (increasing laser power or reducing laser scanning speed) promotes the foaming agent activation thereby porosity formation. However, high linear energy density skews pore distribution towards the surface despite forming deeper melt pools. In addition, the impact of additional factors including foaming agent's laser absorptivity and decomposition kinetics with respect to AM time scales should be carefully considered to avoid ineffective activation of foaming agents during the AM of porous metals.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Hydrogen Adsorption, Reactivity, and Catalysis on Colloidal Iron Carbide Nanoparticles

Supported iron carbide particles have long served as catalysts for CO hydrogenation (the Fischer Tropsch synthesis, FTS) and continue to be attractive. Despite this, little is known about their chemistry. Reported here is a colloidal Fe x C nanoparticle (NP) model system that allows direct observation of surface hydrogen and CO, as well as quantification of the surface H. Dodecylamine-capped Fe x C NPs (DDA-Fe x C NPs) were synthesized through solution-phase carburization of Fe NPs and form stable colloids in low-polarity organic solvents. Treatment of these colloids with H 2 or D 2 produced highly hydrogenated materials, and FTIR spectra of DDA-Fe x C-D n showed that most of the D binds to carbides, with at least four distinct v(C-D) modes. The surface C-H(D) bonds were all reactive, transferring hydrogen to alkenes and other reagents in solution. Titration and ICP measurements showed a 0.17:1 ratio of added H:total Fe, or roughly. 40 H per 1.8 nm DDA-Fe x C NP. Conversely, CO was preferentially bound to surface Fe sites, with FTIR spectra showing a single broad v(CO) that shifted with CO coverage or co-adsorption of H 2 . The DDA-Fe x C NPs were active catalysts for both olefin hydrogenation and the FTS, under mild conditions and without catalyst pre-treatment. The CO hydrogenation reactions yielded a broad distribution of long-chain linear paraffins and olefins. Though quantitative comparisons with typical FTS results are not possible because of our use of sealed batch reactors and other factors, the observations of high catalytic reactivity demonstrate the relevance of this model system to iron-carbide catalysis. DFT calculations on model slab surfaces with varying iron carbide stoichiometries revealed that the thermodynamically preferred surface adsorption sites are C for H ads and Fe for CO ads . A variety of binding sites and binding energies were found for each adsorbate. We are unaware of previous studies indicating that a diverse array of C-H bonds is the primary source of reactive H on iron carbides. Experimentally, the diversity of *C-H sites was evident in reactions with H-atom donors and abstractors of different strengths, from both the reaction stoichiometries and IR spectra. The different surface–H binding energies correlate with the v(C-D) stretching frequencies. Furthermore, these insights into complex iron carbide surfaces and catalysis could assist catalyst design, and they showcase the importance of stoichiometric studies of reaction intermediates.

08 HYDROGEN↗
DOE OSTI2021

Bringing discrete-time Langevin splitting methods into agreement with thermodynamics

In light of the recently published complete set of statistically correct Grønbech–Jensen (GJ) methods for discrete-time thermodynamics, we revise a differential operator splitting method for the Langevin equation in order to comply with the basic GJ thermodynamic sampling features, namely, the Boltzmann distribution and Einstein diffusion, in linear systems. This revision, which is based on the introduction of time scaling along with flexibility of a discrete-time velocity attenuation parameter, provides a direct link between the ABO splitting formalism and the GJ methods. This link brings about the conclusion that any GJ method has at least weak second order accuracy in the applied time step. It further helps identify a novel half-step velocity, which simultaneously produces both correct kinetic statistics and correct transport measures for any of the statistically sound GJ methods. Explicit algorithmic expressions are given for the integration of the new half-step velocity into the GJ set of methods. Finally, numerical simulations, including quantum-based molecular dynamics (QMD) using the QMD suite Los Alamos Transferable Tight-Binding for Energetics, highlight the discussed properties of the algorithms as well as exhibit the direct application of robust, time-step-independent stochastic integrators to QMD.

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

Beyond Fisher forecasting for cosmology

The planning and design of future experiments rely heavily on forecasting to assess the potential scientific value provided by a hypothetical set of measurements. The Fisher information matrix, due to its convenient properties and low computational cost, provides an especially useful forecasting tool. However, the Fisher matrix only provides a reasonable approximation to the true likelihood when data are nearly Gaussian distributed and observables have nearly linear dependence on the parameters of interest. Also, Fisher forecasting techniques alone cannot be used to assess their own validity. Thorough sampling of the exact or mock likelihood can definitively determine whether a Fisher forecast is valid, though such sampling is often prohibitively expensive. Here we propose a simple test, based on the Derivative Approximation for likelihoods (DALI) technique, to determine whether the Fisher matrix provides a good approximation to the exact likelihood. We show that the Fisher matrix becomes a poor approximation to the true likelihood in regions where two-dimensional slices of level surfaces of the DALI approximation to the likelihood differ from two-dimensional slices of level surfaces of the Fisher approximation to the likelihood. We demonstrate that our method accurately predicts situations in which the Fisher approximation deviates from the true likelihood for various cosmological models and several data combinations, with only a modest increase in computational cost compared to standard Fisher forecasts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Multiparticle-hole excitations in nuclei near N = Z = 20: $^{41}$K

This experimental study of high-spin structure near the N = Z = 20 region focuses on 41 K and also reports three newly observed γ transitions in 41 Ca from the same reaction. High-spin states were populated using the 26 Mg( 18 O, p2nγ) 41 K and 26 Mg( 18 O, 3nγ) 41 Ca reactions at a beam energy of 50 MeV at the Florida State University John D. Fox Superconducting Linear Accelerator Laboratory, employing the FSU high-purity germanium detector array. The level scheme of 41 K was extended to 12325 keV, possibly up to J π = 25/2 − or 27/2 + , by means of 25 new γ-ray transitions, and that of 41 Ca to 9916 keV. Linear polarization and angular-distribution measurements were used to provide spin and parity information for several states in the 41 K decay scheme. The results are compared with spsdpf cross-shell shell-model calculations using the FSU interaction. Configurations involving zero or one nucleon promoted from the sd to the fp shell reproduce the energies of many known states reasonably well, while multi-particle excitations reveal a more complex interplay of single-particle motion, collective pairing, and deformation, posing an interesting challenge for future theoretical work.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

A two-level GPU-accelerated incomplete LU preconditioner for general sparse linear systems

This paper presents a parallel preconditioning approach based on incomplete LU (ILU) factorizations in the framework of Domain Decomposition (DD) for general sparse linear systems. We focus on distributed memory parallel architectures, specifically, those that are equipped with graphic processing units (GPUs). In addition to block-Jacobi, we present general purpose two-level ILU Schur complement-based approaches, where different strategies are presented to solve the coarse-level reduced system. These strategies are combined with modified ILU methods in the construction of the coarse-level operator, in order to effectively remove smooth errors by targeting an algebraically smooth vector. We leverage available GPU-based sparse matrix kernels to accelerate the setup and the solve phases of the proposed ILU preconditioner. We evaluate the efficiency of the proposed methods as a smoother for algebraic multigrid (AMG) and as a preconditioner for Krylov subspace methods on challenging anisotropic diffusion problems and a collection of general sparse matrices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Evaluating Space Object Conjunction Probabilities Using Characteristic Function Inversion

This report discusses an approach to computing the probability of a conjunction between two space objects in the short-term encounter scenario. A conjunction is defined here as an event where the miss distance between the objects is less than some specified value. The scenario assumptions are that the motion of the objects is linear, their positions are Gaussian distributed, and their velocities are known and constant. Under these assumptions, the squared-miss distance is shown to have the generalized chi-square distribution. An established statistical technique called characteristic function inversion is employed to evaluate the distribution and obtain conjunction probabilities. The method is closely related to a recent approach based on moment generating function inversion, and a qualitative comparison of the approaches is provided. Last, the method is tested on several benchmark test cases where it agrees with numerical integration on the cases with conjunction probabilities above 10 –12 . However, the exact probability in these cases is usually not needed and this probability can be bounded above using an independent Gaussian approximation. Overall, the report shows how to compute conjunction probabilities using a standard statistical method, though numerical integration seems to perform equally well.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

An entropy-based debiasing approach to quantifying experimental coverage for novel applications of interest in the nuclear community

This manuscript proposes a novel information-theoretic approach to the quantification of experimental relevance, i.e., coverage, to achieve optimal data assimilation results for nuclear engineering applications. Specifically, this work posits the need for a new metric, called coverage (q C ) of an application’s quantity of interest, i.e., eigenvalue or power peaking for an advanced reactor concept, defined herein as the theoretically maximum achievable reduction in the quantity’s uncertainty given measurements from a pool of experiments in a manner that is independent of the data assimilation procedure employed. Currently, reduction in a quantity’s uncertainty is strongly biased by the underlying assumptions of the assimilation procedure to account for the under-determined nature of such problems and the similarity criterion employed to identify relevant experiments. To address this challenge, this work has developed a coverage metric, q C , based on mutual information, which establishes a new conceptual framework for assessing coverage, one that is independent of the model parameters and responses degree of variations in both the experimental and application domains, i.e., linear vs non-linear, and their prior uncertainty distributions, i.e., Gaussian vs. non-Gaussian. The q C is an entropic measure capable of addressing coverage for general nonlinear problems with non-Gaussian uncertainties and inclusive of the measurement uncertainties from multiple experiments. Numerical experiments from manufactured analytical problems as well as a set of benchmarks from the ICSBEP handbook are employed to demonstrate its theoretical and practical performance as compared to the c k -based experiment selection methodology, commonly employed in the neutronic community. The manuscript then employs other well-known adaptations to existing data assimilation methodologies for nonlinear and non-Gaussian problems capable of achieving the coverage posited by q C .

Bayesian data assimilation↗
DOE OSTI2021

Supercooled Liquid Water Detection Capabilities from Ka-Band Doppler Profiling Radars: Moment-Based Algorithm Formulation and Assessment

The occurrence of supercooled liquid water in mixed-phase cloud (MPC) affects their cloud microphysical and radiative properties. The prevalence of MPCs in the mid- and high latitudes translates these effects to significant contributions to Earth’s radiative balance and hydrological cycle. The current study develops and assesses a radar-only, moment-based phase partition technique for the demarcation of supercooled liquid water volumes in arctic, MPC conditions. The study utilizes observations from the Ka band profiling radar, the collocated high spectral resolution lidar, and ambient temperature profiles from radio sounding deployments following a statistical analysis of 5.5 years of data (January 2014–May 2019) from the Atmospheric Radiation Measurement observatory at the North Slope of Alaska. The ice/liquid phase partition occurs via a per-pixel, neighborhood-dependent algorithm based on the premise that the partitioning can be deduced by examining the mean values of locally sampled probability distributions of radar-based observables and then compare those against the means of climatologically derived, per-phase probability distributions. Analyzed radar observables include linear depolarization ratio (LDR), spectral width, and vertical gradients of reflectivity factor and radial velocity corrected for vertical air motion. Results highlight that the optimal supercooled liquid water detection skill levels are realized for the radar variable combination of spectral width and reflectivity vertical gradient, suggesting that radar-based polarimetry, in the absence of full LDR spectra, is not as critical as Doppler capabilities. The cloud phase masking technique is proven particularly reliable when applied to cloud tops with an Equitable Threat Score (ETS) of 65%; the detection of embedded supercooled layers remains much more uncertain (ETS = 27%).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Stochastic Averaging of Radiative Transfer Coefficients for Relativistic Electrons

Abstract Synchrotron emissivities, absorptivities, and Faraday rotation and conversion coefficients are needed in modeling a variety of astrophysical sources, including Event Horizon Telescope (EHT) sources. We develop a method for estimating transfer coefficients that exploits their linear dependence on the electron distribution function, decomposing the distribution function into a sum of parts each of whose emissivity can be calculated easily. We refer to this procedure as stochastic averaging and apply it in two contexts. First, we use it to estimate the emissivity of an isotropicκdistribution function with a high-energy cutoff. The resulting coefficients can be evaluated efficiently enough to be used directly in ray-tracing calculations, and we provide an example calculation. Second, we use stochastic averaging to assess the effect of subgrid turbulence on the volume-averaged emissivity and along the way provide a prescription for a turbulent emissivity. We find that for parameters appropriate to EHT sources turbulence reduces the emissivity slightly. In the infrared, turbulence can dramatically increase the emissivity.

Astronomy & Astrophysics↗
DOE OSTI2021

Efficient Automated Driving Strategies Leveraging Anticipation and Optimal Control

Automated vehicles and advanced driver assistance systems bring computation, sensing, and communication technologies that exceed human abilities in some ways. For example, automated vehicles may sense a panorama all at once, do not suffer from human impairments and distractions, and could wirelessly communicate precise data with neighboring vehicles. Prototype and commercial deployments have demonstrated the capability to relieve human operators of some driving tasks up to and including fully autonomous taxi rides in some areas. The ultimate impact of this technology’s large-scale market penetration on energy efficiency remains unclear, with potential negative factors like road use by empty vehicles competing with positive ones like automatic eco-driving. Fundamentally enabled by historic and look-ahead data, this dissertation addresses the use of automated driving and driver assistance to optimize vehicle motion for energy efficiency. Facets of this problem include car following, co-optimized acceleration and lane change planning, and collaborative multi-agent guidance. Optimal control, especially model predictive control, is used extensively to improve energy efficiency while maintaining safe and timely driving via constraints. Techniques including chance constraints and mixed integer programming help overcome uncertainty and non-convexity challenges. Extensions of these techniques to tractor trailers on sloping roads are provided by making use of linear parameter-varying models. To approach the wheel-input energy eco-driving problem over generally shaped sloping roads with the computational potential for closed-loop implementation, a linear programming formulation is constructed. Distributed and collaborative techniques that enable connected and automated vehicles to accommodate their neighbors in traffic are also explored and compared to centralized control. Using simulations and vehicle-in-the-loop car following experiments, the proposed algorithms are benchmarked against others that do not make use of look-ahead information.

Dollar, Robert Austin↗