Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Weight function”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Progress in modelling fast-ion D-alpha spectra and neutral particle analyzer fluxes using FIDASIM

FIDASIM is a code that models signals produced by charge-exchange reactions between neutrals and ions (both fast and thermal) in magnetically confined plasmas. With the ion distribution function as input, the code predicts the efflux to a neutral particle analyzer diagnostic and the photon radiance of Balmer-alpha light to a fast-ion D α diagnostic, in addition to many other related quantities. A new, parallelized version of the Monte Carlo code FIDASIM has been developed in Fortran90 that is substantially faster than the original interactive data language version. Modified algorithms include more accurate treatments of the time dependent collisional-radiative equations that describe neutral energy levels, of the cloud of ‘halo’ neutrals that surround the injected neutral beam, and of finite Larmor radius effects. Enhanced physics capabilities include modelling ‘passive’ signals from cold edge neutrals, the ability to treat general three-dimensional magnetic confinement configurations, and calculations of diagnostic-specific weight functions that enable tomographic reconstructions of the fast-ion distribution function. Neutral beam attenuation, beam emission, and fast-ion birth profiles are also modelled. Finally, the new algorithms have been successfully validated against experimental data and new features have been tested through benchmarks between two independently developed versions of the code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Instability and treatments of the coupled discrete element and lattice Boltzmann method by the immersed moving boundary scheme

The immersed moving boundary (IMB) scheme has been extensively used to couple the discrete element method (DEM) with the lattice Boltzmann method (LBM). In the literature, only the formulation of IMB for lattice nodal cells covered by a single-solid particle was given. The treatment of situations where a nodal cell is covered by two or more solid particles is seldom discussed. It is found that some numerical instability can occur for such situations due to an inappropriate computation of the weighting function in the IMB formulation. This work presents an enhanced treatment that can resolve the issue and validates it using some benchmark tests. Furthermore, to avoid the extra costs associated with the treatment and simplify the complicated procedure introduced, a simplified IMB scheme is proposed. The accuracy of both enhanced and simplified IMB schemes are validated by test cases including single-particle sedimentation, two-particle drafting-kissing-tumbling phenomenon, and multiple-particle sedimentation. Then, the robustness of both schemes is examined and discussed using a specially designed flow past cylinders test. Overall, the simplified IMB scheme is proved to be robust and sufficiently accurate and simpler and more effective than the enhanced scheme.

42 ENGINEERING↗
DOE OSTI2024

Non-uniform active learning for Gaussian process models with applications to trajectory informed aerodynamic databases

The ability to non-uniformly weight the input space is desirable for many applications, and has been explored for space-filling approaches. Increased interests in linking models, such as in a digital twinning framework, increases the need for sampling emulators where they are most likely to be evaluated. In particular, here we apply non-uniform sampling methods for the construction of aerodynamic databases. This paper combines non-uniform weighting with active learning for Gaussian Processes (GPs) to develop a closed-form solution to a non-uniform active learning criterion. We accomplish this by utilizing a kernel density estimator as the weight function. We demonstrate the need and efficacy of this approach with an atmospheric entry example that accounts for both model uncertainty as well as the practical state space of the vehicle, as determined by forward modeling within the active learning loop.

42 ENGINEERING↗
DOE OSTI2020

A mass-transfer particle-tracking method for simulating transport with discontinuous diffusion coefficients

The problem of a spatially discontinuous diffusion coefficient (D(x)) is one that may be encountered in hydrogeologic systems due to natural geological features or as a consequence of numerical discretization of flow properties. To date, mass-transfer particle-tracking (MTPT) methods, a family of Lagrangian methods in which diffusion is jointly simulated by random walk and diffusive mass transfers, have been unable to solve this problem. This manuscript presents a new mass-transfer (MT) algorithm that enables MTPT methods to accurately solve the problem of discontinuous D(x). To achieve this, we derive a semi-analytical solution to the discontinuous D(x) problem by employing a predictor-corrector approach, and we use this semi-analytical solution as the weighting function in a reformulated MT algorithm. Furthermore, this semi-analytical solution is generalized for cases with multiple 1D interfaces as well as for 2D cases, including a 2 × 2 tiling of 4 subdomains that corresponds to a numerically-generated diffusion field. The solutions generated by this new mass-transfer algorithm closely agree with an analytical 1D solution or, in more complicated cases, trusted numerical results, demonstrating the success of our proposed approach.

42 ENGINEERING↗
DOE OSTI2024

Pan-tropical daily L-band microwave land surface emissivity retrieval from GNSS-R observations

The uncertainties in microwave land surface emissivity (MLSE) measurements have long limited the use of spaceborne microwave radiometer data. As an emerging observation method, Global Navigation Satellite System Reflectometry (GNSS-R) has demonstrated great potential in several land and ocean applications. In this study, a method for obtaining daily MLSE dataset in the pan-tropical region from Cyclone GNSS (CYGNSS) observations is presented and evaluated. The CYGNSS observations are first aggregated into the Equal-Area-Scalable-Earth (EASE) 2.0 36 km grid by a combined weight function of distance, time, and signal-to-noise ratio variance. Then, the method employs a pixel-by-pixel regression algorithm to conduct the daily MLSE retrieval using reference emissivity derived from the Soil Moisture Active Passive (SMAP) brightness temperature. The CYGNSS MLSE shows good agreement with SMAP MLSE, delivering an overall root-mean-square error (RMSE) of 0.022 and 0.017 for horizontal and vertical polarization, respectively, during the training set spanning the whole year of 2018. Furthermore, on the test set from January 2019 to May 2019, the RMSE values amounted to 0.030 and 0.023 for horizontal and vertical polarization, respectively. Temperature records from the International Soil Moisture Network are employed to calculate the emissivity and for in-situ validation, which yield an RMSE of 0.034 and 0.026 for the two polarizations, respectively. The proposed algorithm provides an encouraging approach to obtain accurate daily MLSE dataset for microwave remote sensing. Compared to the SMAP MLSE, the CYGNSS MLSE has a remarkable improvement of 86% in temporal resolution, greatly complementing the existing microwave emissivity datasets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Models of bars − II. Exponential profiles

ABSTRACT We present a new model for galactic bars with exponentially falling major axis luminosity profiles and Gaussian cross-sections. This is based on the linear superposition of Gaussian potential–density pairs with an exponential weight function, using an extension of the method originally introduced by Long & Murali. We compute the density, potential, and forces, using Gaussian quadrature. These quantities are given as explicit functions of position. There are three independent scaled bar parameters that can be varied continuously to produce bespoke bars of a given mass and shape. We categorize the effective potential by splitting a reduced parameter space into six regions. Unusually, we find bars with three stable Lagrange points on the major axis are possible. Our model reveals a variety of unexpected orbital structure, including a bifurcating x1 orbit coexisting with a stable x4 orbit. Propeller orbits are found to play a dominant role in the orbital structure, and we find striking similarities between our bar configuration and the model of Kaufmann & Contopoulos. We find a candidate orbital family, sired from the propeller orbits, that may be responsible for the observed high-velocity peaks in the Milky Way’s bar. As a cross-check, we inspect, for the first time, the proper motions of stars in the high-velocity peaks, which also match our suggested orbital family well. This work adds to the increasing body of evidence that real galactic bars may be supported at least partly by propeller orbits rather than solely by elliptical-like orbits of the x1 family.

McGough, D. P.↗
DOE OSTI2021

Mitigating the effects of undersampling in weak lensing shear estimation with metacalibration

ABSTRACT metacalibration is a state-of-the-art technique for measuring weak gravitational lensing shear from well-sampled galaxy images. We investigate the accuracy of shear measured with metacalibration from fitting elliptical Gaussians to undersampled galaxy images. In this case, metacalibration introduces aliasing effects leading to an ensemble multiplicative shear bias about 0.01 for Euclid and even larger for the Roman Space Telescope, well exceeding the missions’ requirements. We find that this aliasing bias can be mitigated by computing shapes from weighted moments with wider Gaussians as weight functions, thereby trading bias for a slight increase in variance of the measurements. We show that this approach is robust to the point-spread function in consideration and meets the stringent requirements of Euclid for galaxies with moderate to high signal-to-noise ratios. We therefore advocate metacalibration as a viable shear measurement option for weak lensing from upcoming space missions.

Kannawadi, Arun↗
DOE OSTI2022

Simulations for x-ray imaging of wave-packet dynamics

Previous work on imaging wave-packet dynamics with x-ray scattering revealed that the scattering patterns deviate substantially from the notion of instantaneous momentum density of the wave packet. Here we show that scattering patterns can provide clear insights into the electron wave-packet dynamics if the final state of the scattered electron and the scattered photon momentum are determined simultaneously. The scattering probability is shown to be proportional to the modulus square of the Fourier transform of the instantaneous electronic spatial wave function weighted by the final state of the electron. Several cases for the choice of final state of the electron are explored. First, the case where the final state can be measured up to a given principal quantum number n and orbital angular momentum l are presented. Next, the case where the final states can only be determined up to a given energy is discussed. Lastly, the case of an initial wave packet consisting of a large amount of a known stationary state and a small amount of an unknown stationary state is examined. The scattering profile is used to determine the properties of the unknown state in the wave packet.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Modeling Metallic Fuel using Peridynamics

Based on available modeling and simulation capabilities of peridynamics module in MOOSE framework for oxide fuel, the overall goal of this project is to further develop the peridynamics capabilities for modeling metallic fuel. It includes two major tasks: 1) develop validated scheme to handle the shape tensor singularity due to insufficient active neighbors of a material particle in the peridynamic correspondence model for fracture problems, and 2) develop failure modeling scheme including failure criterion for metallic fuels. Before the peridynamics can be applied to model metallic fuel, the formulation instability of the peridynamic correspondence model should be addressed. The PI first worked on developing new stabilization method to improve the performance of the peridynamic correspondence model and reduce the possibility of getting a singular shape tensor while applying the model for fracture problems. The new stabilization scheme uses bond-associated weight function rather than bond-associated horizon. Compared to bondassociated horizon stabilized method, this new stabilization scheme has better performance with improved prediction accuracy and reduced free surface effect. Using this newly developed stabilization, materials models from BISON can be directly used in peridynamics for metallic fuels, such as fission rate and burnup dependent creep and swell models. Publication of this work is under preparation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Group Structure Machine Learning Proposal

Nuclear data is the linchpin underwriting several fundamental capabilities and mission needs at LANL. New techniques such as machine learning can be brought to bear to solve old problems such as multigroup cross-section accuracy. In neutron transport, generating multigroup cross sections is a complex and arcane task, but a crucial one, as accurate solutions require appropriate cross sections. There are two key challenges when generating multigroup cross sections: (1) choosing an accurate weight function, and (2) choosing appropriate energy boundaries. Often, energy boundaries are chosen using “expert judgment” that is not documented and is difficult to replicate. The long-standing Los Alamos 30-group structure has been in use since at least 1969 and is still in use today. Simplistic attempts over the years since to improve on the 30-group structure have been met with limited success. Machine learning algorithms would enable the selection of appropriate, problem-dependent group boundaries without an inordinate investment of scientist time. We will develop workflows and tools to enable these improved group boundary choices, which will reduce uncertainty and increase predictive capability of neutron-transport applications at LANL.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Mapped Moments to a Cartesian Grid (MMCG) Value-Added Product Report

Objective analysis (OA) is a method of mapping unstructured data to a structured grid. In the context of scanning radar data, OA is used to interpolate data in antenna coordinates (range, azimuth, and elevation) onto a regularly spaced Cartesian grid (Trapp and Doswell 2000). The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s Mapped Moments to a Cartesian Grid (MMCG) Value-Added Product (VAP) uses the Python ARM Radar Toolkit (Py-ART), a data model-driven interactive architecture for working with weather radar data, to map the data to a Cartesian grid (Helmus and Collis 2016). MMCG, with Py-ART built in, has the ability to take radar data in antenna coordinates and map the gates to a Cartesian grid using inverse distance weight functions such as Cressman (square) and Barnes (exponential), but also can filter the data during the interpolation. MMCG also allows arbitrary formulations for the radius of influence, which are matched to particular radar scanning strategies. This creates a complex parameter space for optimizing the retention of storm structure detail while minimizing artifacts. MMCG takes data processed with ARM’s Corrected Precipitation Radar Moments in Antenna Coordinates (CMAC) VAP and maps it to a Cartesian grid as the output product. A variety of fields that have been mapped to the Cartesian grid are then saved to plots to complement each grid file.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Transport @ LANL: An Overview [Slides]

Deterministic neutron transport requires discretization of the Boltzmann transport equation in space, angle, energy, and time. Discretization in energy is typically performed with the multigroup approximation, which tries to preserve continuous energy cross section information by using a weighted average over discrete energy groups. There are two important requirements for accurate multigroup cross sections: appropriate energy group boundaries and an accurate weight function. The former can be selected with appropriate optimizers; the latter either requires flux solutions from comparable simulations or expert knowledge.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Towards Realistic and High Fidelity Models for Nuclear Reactor Power Synthesis Simulation with Self-Powered Neutron Detectors

As presented in this report, a weighting function–based inferencing method is being applied to synthesize the power distribution in next-generation and university research reactors based on simulated self power neutron detector (SPND) responses. The overall goal is to assess the impacts of sensor uncertainty and true power distribution perturbations on the error in the synthesized power distribution. Regarding sensor uncertainty, the NuScale Small Modular Reactor (SMR) and the Westinghouse AP1000 serve as testbeds for analyzing the impact of varying the sensor uncertainty, as well as varying the number of sensors per sensor string in the reactor core. The reactor models are informed by Monte Carlo N-Particle (MCNP) neutron flux tallies. For the NuScale SMR and Westinghouse AP1000, the SPND response functions (i.e., the response of the SPNDs to individual segments of fuel) were determined homogeneously. Regarding an analysis of power distribution perturbation detection, the Texas A&M Testing, Research, Isotopes, General Atomics Reactor (TAMU TRIGA) reactor was used as a demonstration case with one particular arrangement of SPNDs; the response functions for this reactor model were determined heterogeneously, making this a uniquely high-fidelity demonstration of perturbation detection. Finally, SPND models generated in the Geometry and Tracking 4 (Geant4) code have been generated and tested for comparison with traditionally implemented analytical SPND models, with the intent for Geant4 integration with the full methodological framework. SPND current outputs as a function of distance from some fuel assembly segment in the NuScale SMR are compared with the analytically determined currents. Results from the sensor uncertainty simulations for the NuScale SMR and AP1000 indicate that the average error in the inferred power distribution on the fuel assembly segment level is reasonably low, being slightly less than the random uncertainty applied to all respective SPNDs in both cores. For example, if all SPNDs in the core have a random uncertainty of 5%, then the corresponding fuel assembly segment level error (i.e. difference between the true and inferred local power) is ~2–3%. However, the maximum error in the inferred power distribution on the fuel assembly segment level can be considerably high (>15%) when SPND random uncertainties start to exceed ~3%. In general, the average and maximum errors in the inferred power distribution were slightly higher in the AP1000 as opposed to the NuScale SMR for the sensor string configurations considered herein. Another result determined from analysis of the sensor uncertainty simulations was that increasing the number of SPNDs per string does not clearly reduce inferred power distribution error and can in fact make the error large in some cases; however, this assessment may skewed due to imposed iteration limits. Results from the perturbation detection demonstration using the high-fidelity TAMU TRIGA model indicate that, given the arrangement of 17 SPND strings and 4 SPNDs per string considered herein, there is a clear, provable ability to infer a localized Gaussian-type peak perturbation in the 3D power distribution. Such a perturbation was detected with an average fuel assembly segment level error of 0.19%, and the general visualization of the detected perturbation clearly indicates that the magnitude and shape were appropriately resolved. Finally, the electrical current output generated by the Geant4 modeled SPND indicates significant magnitude differences than the analytically modeled SPND, demonstrating the need for accurate SPND models which account for finite sensor geometry effects to inform the power synthesis work described herein.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Quantized Constant-Q Gabor Atoms for Sparse Binary Representations of Cyber-Physical Signatures

Increased data acquisition by uncalibrated, heterogeneous digital sensor systems such as smartphones present new challenges. Binary metrics are proposed for the quantification of cyber-physical signal characteristics and features, and a standardized constant-Q variation of the Gabor atom is developed for use with wavelet transforms. Two different continuous wavelet transform (CWT) reconstruction formulas are presented and tested under different signal to noise ratio (SNR) conditions. A sparse superposition of Nth order Gabor atoms worked well against a synthetic blast transient using the wavelet entropy and an entropy-like parametrization of the SNR as the CWT coefficient-weighting functions. The proposed methods should be well suited for sparse feature extraction and dictionary-based machine learning across multiple sensor modalities.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗
DOE OSTI2023

Hybrid Temperature-Reactivity PID Controllers for Nuclear Thermal Propulsion Startup

In this work, a novel temperature-reactivity controller is introduced to rotate control drums based on a chamber temperature signal. It consists of two components: (1) a reactivity-driven proportional controller converting the temperature demand into a reactivity signal, and (2) a temperature-driven proportional-derivative controller which relies on a derivative term to anticipate sudden changes in demand and limit temperature overshoot. A weighting function is proposed to continuously transition from one component to the other.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Hybrid Temperature-Reactivity PID Controllers for Nuclear Thermal Propulsion Startup

In this work, a novel temperature-reactivity controller is introduced to rotate control drums based on a chamber temperature signal. It consists of two components: (1) a reactivity-driven proportional controller converting the temperature demand into a reactivity signal, and (2) a temperature-driven proportional-derivative controller which relies on a derivative term to anticipate sudden changes in demand and limit temperature overshoot. A weighting function is proposed to continuously transition from one component to the other.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Efficient generalized boundary detection

Fast, efficient, and robust compression-based methods for detecting boundaries in arbitrary datasets, including sequences (1D datasets), are desired. The methods, each employing three simple algorithms, approximate the information distance between two adjacent sliding windows within a dataset. One of the algorithms calculates an initial ordered list of subsequences; while a second algorithm updates the ordered list of subsequences by dropping a first entry and appending a last entry rather than calculating completely new ordered lists with each iteration. Large values in the distance metric are indicative of boundary locations. A smoothed z-score or a wavelet-based algorithm may then be used to locate peaks in the distance metric, thereby identifying boundary locations. An adaptive version of the method employs a collection of window sizes and corresponding weighting functions, making it more amenable to real datasets with unknown, complex, and changing structures.

Ting, Christina↗
DOE OSTI2026

Adaptive spectra-to-exposure conversion using ridge regularized polynomial response models

Real-time gamma spectra-to-exposure conversion in aerial and ground monitoring commonly relies on calibration-derived, detector- or system-specific conversion coefficients that are assumed to generalize across operational environments. In practice, deployment specific differences in spectral composition and transport conditions can introduce systematic bias relative to reference instruments, motivating methods that adapt coefficients using minimal field supervision while explicitly limiting overfitting. In this work, we present a conservative coefficient adaptation framework that updates a baseline polynomial energy-weighting function using ridge-regularized regression, with leave-one-out cross-validation (LOOCV) used to select the regularization strength. The findings support ridge-constrained minimal-supervision adaptation as a practical mechanism to suppress site-specific bias without destabilizing a calibration-derived baseline.

61 RADIATION PROTECTION AND DOSIMETRY↗