Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical techniques”

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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Comparison and investigation of production of specific nuclides by different fission yield data

Generally, passive gamma spectroscopy (PGS) is known as a useful technique for nondestructive measurement of burnup and nuclide production in normal spent nuclear fuel. However, PGS can be used not only for normal spent fuel, but also for molten fuel generated in severe accidents involving core meltdowns. In such cases, it is estimated by the radioactivity ratios of specific nuclides. At the time of the decommissioning of the Three Mile Island's unit 2 reactor, fission products (FPs) such as Cs{sup 134}, Cs{sup 137}, Ce{sup 144}, Eu{sup 154}, Ru{sup 106}, and Sb{sup 125} were concerned. The radioactivity ratios are estimated by numerical simulations. However, there are various uncertainties in numerical simulations. In this study, we evaluated the effect of the difference in fission yield data. We compared the productions of concerned nuclides of Three Mile Island's unit 2 reactor by using Tokyo-tech and JENDL FP Fission Yields Data File 2011(FPY- 2011). Tokyo-tech is newly created by Tsubakihara et al. at Tokyo Institute of Technology. FPY-2011 is currently the latest fission yield data in Japan. In the concerned nuclides, the relative difference (Tokyo-tech/FPY-2011 -1) in the production of Sb{sup 125} is 4.5% and that of Eu{sup 154} is -5.6%. The relative difference in the production of Sb{sup 125} is found to be caused by the difference in the independent yields of Sn{sup 125}, In{sup 125}, Cd{sup 125}, Ag{sup 125}, and Pd{sup 125}. The relative difference in Eu{sup 154} is found to be caused by the relative difference in the production of Sm{sup 152} and the relative difference in the cumulative yield of Eu{sup 153}.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

HAIMOS Ensemble Forecasts for Intra-day and Day- Ahead GHI, DNI and Ramps

The objective of this research is to develop a hybrid physics-based/data-driven forecast model to improve direct normal and global horizontal irradiance (DNI and GHI) prediction for horizons ranging from 1 to 72 hours. Project objectives also address key gaps in state-of-the-art solar forecasting: accurate probabilistic solar forecasts and the forecasting of large irradiance ramps (ramp onset and magnitude). The proposed model ensembles Numerical Weather Prediction (NWP) forecasts, determinist physics-based algorithms, and new-generation cloud cover products (high-resolution rapid refresh satellite images and Large Eddy Simulations). The result is the Hybrid Adaptive Input Model Objective Selection (HAIMOS) ensemble model. HAIMOS blends state of the art machine learning methodologies with physics-based models for cloud cover and cloud optical depth forecasts. The technical activities followed a two-pronged strategy. First, the preprocessing of data, the selection of inputs to the nonlinear approximators, the type of approximator and objective functions, and post-processing ensembling techniques included in HAIMOS were all optimized adaptively to find the best model for a specific goal (reduce DNI/GHI forecast error, improve the prediction of ramp onset, etc.). Second, a large effort was put in improving cloud identification and the forecast of cloud cover and cloud optical depth. To this end, new-generation cloud parametrization products were developed in this work. These include improved algorithms to assist in cloud identification, cloud classification and cloud parametrization from satellite images – three key factors in the accuracy of 1 to 6-hours irradiance forecasts and prediction of ramp onset. Furthermore, we also included cloud information extracted from high resolution rapid refresh satellite images (GOES-16) and Large Eddy Simulations (LES). LES was used to model the atmosphere in detail over locations of interest and produce cloud optical depth forecasts. Once these data streams were validated, they were used as input data to the HAIMOS forecast. The model was developed using data from several climatologically distinct locations with potential for high solar penetration. In the last year of the project, we conducted a validation campaign according to the guidelines stipulated by the Topic Area 1 project as described in the FOA. This effort brings, for the first time, proven machine-learning methodologies for generating state-of-the-art solar forecasts interweaved with detailed physics-based models for cloud detection, and cloud optical depth forecasts. HAIMOS will generate accurate irradiance probabilistic forecast to assist in reducing solar generation prediction error. Globally optimized solar forecast models are more likely to impact solar energy stakeholders. The goal of this project was to increase the state-of-the-art forecast skill from their present values of 10 to 35%. At the end of the project, we achieved between 30% and 50% forecast skill across a wide range of horizons for both GHI and DNI.

14 SOLAR ENERGY↗
DOE OSTI2026

Chemical and Radiological Compatibility Testing of 3D Printed Materials

To dramatically increase the adaptability, performance, and safety of processes in support of the Defense Waste Processing Facility (DWPF), Savannah River National Lab (SRNL) plans to perform chemical and radiological compatibility testing on a wide variety of 3D printed materials of interest. The 3D printing process provides numerous strategic operational benefits such as rapid prototyping of complex designs and geometry specific to the needs of the nuclear waste disposition process, as well as on-demand rapid prototyping and iteration with materials that aren’t as accessible through traditional manufacturing methods. Reaction chemistry in simulated waste batches can be matched closely to its radioactive counterpart, but glass reactor vessels have limitations. Vessel geometry can play a big factor in mixing transport limitations, process chemistry, and degradation reaction kinetics. In addition, additive manufacturing allows for much more detailed vessel design than traditional alternatives. Waste processing techniques in DWPF also encounter extreme chemical environments including high pH, strong acids, abrasive slurries, and significant irradiation. To meet these challenges, a matrix of various polymer, ceramic, and metal additive manufacturing materials have been exposed to a suite of chemical environments of interest as well as radioactive dose (such as gamma radiation from 60 Co) to properly test their durability under these conditions. Mass change has been monitored over a period of up to a week in these conditions, as well as added characterization for surface modification through Scanning Electron Microscopy/Electron Dispersive X-ray analysis (SEM/EDX). Further chemical characterization has been monitored through Fourier-Transform InfraRed Spectroscopy (FTIR), with planned investigation via thermal and tensile strength degradation. While the direct product of this research is identification of material(s) that can withstand specific hazardous environments encountered by the mercury water wash tank in DWPF process simulation experiments, the reference base of materials will be used for many other nuclear processes in the pursuit of rapidly developed, cost-efficient, and highly specific devices for environmental remediation and much more.

Wilson, Nathan W. [Savannah River National Laborat↗
DOE OSTI2023

Customizable adaptive regularization techniques for B-spline modeling

B-spline models are a powerful way to represent scientific data sets with a functional approximation. However, these models can suffer from spurious oscillations when the data to be approximated are not uniformly distributed. Model regularization (i.e., smoothing) has traditionally been used to minimize these oscillations; unfortunately, it is sometimes impossible to sufficiently remove unwanted artifacts without smoothing away key features of the data set. In this article, we present a method of model regularization that preserves significant features of a data set while minimizing artificial oscillations. Our method varies the strength of a smoothing parameter throughout the domain automatically, removing artifacts in poorly-constrained regions while leaving other regions unchanged. Further, the proposed method selectively incorporates regularization terms based on first and second derivatives to maintain model accuracy while minimizing numerical artifacts. The behavior of our method is validated on a collection of two- and three-dimensional data sets produced by scientific simulations. In addition, a key tuning parameter is highlighted and the effects of this parameter are presented in detail. This paper is an extension of our previous conference paper at the 2022 International Conference on Computational Science (ICCS) (Lenz et al., 2022).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Accelerating Multivariate Functional Approximation Computation with Domain Decomposition Techniques⋆

Modeling large datasets through Multivariate Functional Approximations (MFA) provide an elegant way to handle many visualization and scientific analysis workflows. The process necessitates scalable data partitioning methods to compute MFA representations efficiently without compromising the accuracy or continuity of the reconstructed solution. We propose a domain -decomposed method for computing the MFA with B -spline bases, which reduces the total work per task and uses a restricted Additive Schwarz (RAS) method to converge the control point data degrees -of -freedom along subdomain boundaries. We provide an in-depth analysis of the parallel approach with domain decomposition solvers, aiming to minimize local subdomain error residuals and recover high -order continuity at subdomain interfaces with appropriate choices of knot overlaps. The communication cost, determined by the overlap regions in the RAS implementation, is optimized to recover the numerical error profile of the single subdomain case. Our proposed method stands in contrast to previous methods, which typically only recover either C 0 or at best C 1 continuity for arbitrary B -spline degree expansions, or those that require post -processing to blend discontinuities in the reconstructed data. We demonstrate the effectiveness of our approach using analytical and real -world datasets in 1D, 2D, and 3D through both strong and weak scaling studies. The performance results indicate that the overall cost of computing the approximation is directly proportional to the underlying nearest -neighbor communication implementation, and is only weakly dependent on the overlap region size that determines the size of the messages. This finding underscores the efficiency and scalability of our proposed method, making it a promising solution for handling large datasets in scientific workflows.

additive Schwarz solvers↗
DOE OSTI2021

Motif magnetism and quantum many-body scars

In this work, we generally expect quantum systems to thermalize and satisfy the eigenstate thermalization hypothesis (ETH), which states that finite-energy-density eigenstates are thermal. However, some systems, such as many-body localized systems and systems with quantum many-body scars, violate ETH and have high-energy athermal eigenstates. In systems with scars, most eigenstates thermalize, but a few atypical scar states do not. Scar states can give rise to a periodic revival when time-evolving particular initial product states, which can be detected experimentally. Recently, a family of spin Hamiltonians was found with magnetically ordered three-colored eigenstates that are quantum many-body scars [Lee et al., Phys. Rev. B 101, 241111(R) (2020)]. These models can be realized in any lattice that can be tiled by triangles, such as the triangular or kagome lattices, and have been shown to have close connections to the physics of quantum spin liquids in the Heisenberg kagome antiferromagnet. In this paper, we introduce a generalized family of n-colored Hamiltonians with “spiral colored” eigenstates made from n-spin motifs such as polygons or polyhedra. We show how these models can be realized in many different lattice geometries and provide numerical evidence that they can exhibit quantum many-body scars with periodic revivals that can be observed by time-evolving simple product states. The simple structure of these Hamiltonians makes them promising candidates for future experimental studies of quantum many-body scars.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Reinforcement learning based automated history matching for improved hydrocarbon production forecast

History matching aims to find a numerical reservoir model that can be used to predict the reservoir performance. An engineer and model calibration (data inversion) method are required to adjust various parameters/properties of the numerical model in order to match the reservoir production history. In this study, we develop deep neural networks within the reinforcement learning framework to achieve automated history matching that will reduce engineers’ efforts, human bias, automatically and intelligently explore the parameter space, and remove the need of large set of labeled training data. To that end, a fast-marching-based reservoir simulator is encapsulated as an environment for the proposed reinforcement learning. The deep neural-network-based learning agent interacts with the reservoir simulator within reinforcement learning framework to achieve the automated history matching. Reinforcement learning techniques, such as discrete Deep Q Network and continuous Deep Deterministic Policy Gradients, are used toth, used to train the learning agents. The continuous actions enable the Deep Deterministic Policy Gradients to explore more states at each iteration in a a learning episode; consequently, a better history matching is achieved using this algorithm as compared to Deep Q Network. For simplified dual-target composite reservoir models, the best history-matching performances of the discrete and continuous learning methods in terms of normalized root mean square errors are 0.0447 and 0.0038, respectively. Furthermore, our study shows that continuous action space achieved by the deep deterministic policy gradient drastically outperforms deep Q network.

42 ENGINEERING↗
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 OSTI2024

Numerical investigation of equivalence ratio effects on Flex-Fuel mixing controlled combustion enabled by prechamber ignition

There is interest in adapting renewable and low-carbon intensity fuels to heavy-duty engines to help displace criteria and greenhouse gas emissions associated with diesel combustion. Low-carbon fuels have inherently low cetane numbers and thus pose a significant challenge when considering direct substitution in diesel combustion systems. Use of an actively fueled prechamber as an ignition source to initiate mixing-controlled combustion (PC-MCC) of low-cetane fuels is becoming an attractive combustion mode to alleviate the identified reactivity deficit. This work focuses on the utilization of an active prechamber to facilitate diesel-like combustion of bioethanol-gasoline fuel blends in a heavy-duty engine. Recent results in this emerging technology suggest that ignition quality is uniquely coupled to prechamber equivalence ratio. In this numerical study, the fundamental implications of prechamber equivalence ratio on the ignition performance of direct injected fuels ranging from E10 to E100 are investigated using CONVERGE. Parametric studies of the prechamber operating strategy were assessed at diesel-like conditions to characterize the performance trends relative to a diesel baseline at the same boundary conditions. Simulation results indicate that PC-MCC is flex-fuel capable and achieves diesel-like ignition qualities and combustion processes for all fuels considered under stoichiometric and rich prechamber conditions. To characterize the equivalence ratio trends observed, a novel inflow boundary modeling technique was utilized to prescribe turbulent jets in place of the prechamber where selective speciation of the jet composition was conducted to isolate the ignition contributions of radical and reacting species in lean, stoichiometric, and rich jets. The inflow boundary modeling suggests that excess fuel and combustion intermediates present in jets produced from stratified and fuel-rich prechamber operation promotes higher jet temperatures and as such are superior ignition sources. Relative to a lean prechamber jet, the peak temperature of a reacting jet from a rich prechamber was up to 600 K hotter at fixed distances from the orifice exit. Furthermore, radicals also demonstrated an influence on the ignition process, but the combustion mode was identified to be thermally dominant.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Rotor blade imbalance fault detection for variable-speed marine current turbines via generator power signal analysis

Marine hydrokinetic (MHK) turbines extract renewable energy from oceanic environments. However, due to the harsh conditions that these turbines operate in, system performance naturally degrades over time. Thus, ensuring efficient condition-based maintenance is imperative towards guaranteeing reliable operation and reduced costs for marine hydrokinetic power. This work proposes a novel framework aimed at identifying and classifying the severity of rotor blade pitch imbalance faults experienced by marine current turbines (MCTs). In the framework, a Continuous Morlet Wavelet Transform (CMWT) is first utilized to acquire the wavelet coefficients encompassed within the 1P frequency range of the turbine's rotor shaft. From these coefficients, several statistical indices are tabulated into a six-dimensional feature space. Next, Principle Component Analysis (PCA) is employed on the resulting feature space for dimensionality reduction, and then the application of a K-Nearest Neighbor (KNN) machine learning algorithm is utilized for fault detection and severity classification. The effectiveness of the proposed framework is validated using a high-fidelity MCT numerical simulation platform, where results demonstrate that the presence of a pitch imbalance fault can be accurately detected 100% of the time and correctly classified based upon severity more than 97% of the time.

42 ENGINEERING↗
DOE OSTI2024

Inverse-designed metasurfaces with facile fabrication parameters

Abstract Optical metasurfaces are planar nanostructured devices that are industrially attractive in part because they utilize high-throughput microelectronic fabrication techniques for implementation. It is therefore critical to develop design paradigms that can balance the realization of highly efficient wavefront responses together with device manufacturability. We introduce a gradient-based design framework for freeform metasurfaces in which nanoscale elements are explicitly constrained to feature basic shapes, nearly uniform feature sizes, and exceptionally low aspect ratios. In spite of the apparent uniformity of the metasurface geometric features, the devices are able to utilize nonlocal near-field optical coupling to achieve highly efficient and extreme wavefront scattering beyond conventional design methodologies. Utilizing this approach, we design facile high-numerical-aperture devices such as beam deflectors and large-area metalenses capable of diffraction-limited focusing. We anticipate that these concepts can facilitate the design and integration of metasurfaces into monolithic optical systems.

Optics↗
DOE OSTI2020

Consistent boundary conditions for PDE filter regularization in topology optimization

Abstract Design variables in density-based topology optimization are typically regularized using filtering techniques. In many cases, such as stress optimization, where details at the boundaries are crucially important, the filtering in the vicinity of the design domain boundary needs special attention. One well-known technique, often referred to as “padding,” is to extend the design domain with extra layers of elements to mitigate artificial boundary effects. We discuss an alternative to the padding procedure in the context of PDE filtering. To motivate this augmented PDE filter, we make use of the potential form of the PDE filter which allows us to add penalty terms with a clear physical interpretation. The major advantages of the proposed augmentation compared with the conventional padding is the simplicity of the implementation and the possibility to tune the boundary properties using a scalar parameter. Analytical results in 1D and numerical results in 2D and 3D confirm the suitability of this approach for large-scale topology optimization.

42 ENGINEERING↗
DOE OSTI2022

Significant improvement of conversion efficiency by passivation of low-angle grain boundaries in flexible low-cost single-crystal-like GaAs thin-film solar cells directly deposited on metal tape

We develop an easily implementable yet effective process to achieve high power conversion efficiency in newly-developed flexible single-crystal-like GaAs thin-film photovoltaic solar cells on cheap metal tape. We study various techniques on the fabricated devices to passivate the defect states associated with dangling bonds at low-angle grain boundaries (LA-GBs) existing in the single-crystal-like GaAs material. We report that sulfur by trioctylphosphine sulfide (TOP:S) solution treatment can partially but effectively passivate the defect states at the LA-GBs to increase all the solar cell performance parameters of open-circuit voltage VOC (33.2%), short-circuit current JSC (19.5%), and fill factor FF (25.8%), resulting in the improvement of the power conversion efficiency more than double. The material and device characterizations and numerical modeling are performed to propose a mechanism, also suggesting that the power conversion efficiency higher than 20%, which is an important target milestone in emerging solar cell technologies, can be achieved. This work presents the first effective passivation of the bulk defects in GaAs at the device level and can also offer a solution for other photonic and electronic devices using non-single-crystalline GaAs films for flexible electronics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Quantifying thermal barrier coating delamination through luminescence modeling

Thermal barrier coatings (TBCs) are widely used in the hot sections of gas turbine systems as they are remarkably efficient in insulating the underlying superalloys, leading to higher operating temperatures and therefore enhancing the performance of the engines. However, this benefit is only possible if the integrity of the TBC under aggressive thermo-mechanical environments is ensured. Delamination is a common but hard to detect failure mode. We present, in this work, supported by experimentation, the implementation of a modeling approach applying the Kubelka-Munk theory to provide numerical quantification of luminescence contrast and intensity due to top coat delamination in TBCs. The method relies on the drastic change in reflectivity when a delamination forms, exploiting it for high-contrast luminescence mapping. Two distinct TBC configurations containing an erbium-doped yttria-stabilized zirconia (YSZ:Er 3+ ) layer for delamination sensing were used to validate this model. A delamination zone induced by Rockwell indentation was successfully tracked by measuring an increase of the intensity of the erbium emission line at 562 nm. Luminescence-based methods for delamination detection can provide a revolutionary non-invasive technique, with potential for both off-line and on-line engine monitoring.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Low-cost sensor-based damage localization for large-area monitoring of FRP composites

Abstract In recent years, there has been growing interest in self-sensing structural materials across research and industry sectors. Detecting and locating structural damage typically requires numerous sensors wired to a data acquisition (DAQ) circuit, rendering implementation impractical in real structures. This paper proposes an innovative, cost-effective sensor network for damage detection and localization in fiber-reinforced polymer composites. The innovation encompasses three key elements: (1) utilizing carbon fiber tows within the composite as piezoresistive sensors, eliminating the need for additional foreign sensor devices; (2) introducing a novel sensor layout wherein sensor tow branches with varied resistance values are connected in parallel, reducing the number of connections to the DAQ circuit and cutting manufacturing costs significantly; (3) developing a practical sensor terminal fabrication technique to minimize manufacturing expenses. The proposed design methodology for the branch resistance values is first validated using a demonstration panel. Subsequently, the overall strategy is assessed by conducting impact tests on carbon and glass fiber-reinforced composite specimens. Results validate the sensor’s ability to accurately detect and locate structural damage.

Demo, Luke B. (ORCID:0000000208806522)↗
DOE OSTI2020

Digital-analog quantum computation

Digital quantum computing paradigm offers highly desirable features such as universality, scalability, and quantum error correction. However, physical resource requirements to implement useful error-corrected quantum algorithms are prohibitive in the current era of NISQ devices. As an alternative path to performing universal quantum computation, within the NISQ era limitations, we propose to merge digital single-qubit operations with analog multiqubit entangling blocks in an approach we call digital-analog quantum computing (DAQC). Along these lines, although the techniques may be extended to any resource, we propose to use unitaries generated by the ubiquitous Ising Hamiltonian for the analog entangling block and we prove its universal character. Here, we construct explicit DAQC protocols for efficient simulations of arbitrary inhomogeneous Ising, two-body, and M-body spin Hamiltonian dynamics by means of single-qubit gates and a fixed homogeneous Ising Hamiltonian. Additionally, we compare a sequential approach where the interactions are switched on and off (stepwise DAQC) with an always-on multiqubit interaction interspersed by fast single-qubit pulses (banged DAQC). Finally, we perform numerical tests comparing purely digital schemes with DAQC protocols, showing a remarkably better performance of the latter. The proposed DAQC approach combines the robustness of analog quantum computing with the flexibility of digital methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Bayesian optimization algorithms for accelerator physics

Accelerator physics relies on numerical algorithms to solve optimization problems in online accelerator control and tasks such as experimental design and model calibration in simulations. The effectiveness of optimization algorithms in discovering ideal solutions for complex challenges with limited resources often determines the problem complexity these methods can address. The accelerator physics community has recognized the advantages of Bayesian optimization algorithms, which leverage statistical surrogate models of objective functions to effectively address complex optimization challenges, especially in the presence of noise during accelerator operation and in resource-intensive physics simulations. In this review article, we offer a conceptual overview of applying Bayesian optimization techniques toward solving optimization problems in accelerator physics. We begin by providing a straightforward explanation of the essential components that make up Bayesian optimization techniques. We then give an overview of current and previous work applying and modifying these techniques to solve accelerator physics challenges. Finally, we explore practical implementation strategies for Bayesian optimization algorithms to maximize their performance, enabling users to effectively address complex optimization challenges in real-time beam control and accelerator design. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Challenges and opportunities in verification and validation of low temperature plasma simulations and experiments

This paper describes the verification and validation (V&V) framework developed for the stochastic Particle-in-Cell, Direct Simulation Monte Carlo code Aleph. An ideal framework for V&V from the viewpoint of the authors is described where a physics problem is defined, and relevant physics models and parameters to the defined problem are assessed and captured in a Phenomena Identification and Ranking Table (PIRT). Furthermore, numerous V&V examples guided by the PIRT for a simple gas discharge are shown to demonstrate the V&V process applied to a real-world simulation tool with the overall goal to demonstrably increase the confidence in the results for the simulation tool and its predictive capability. Although many examples are provided here to demonstrate elements of the framework, the primary goal of this work is to introduce this framework and not to provide a fully complete implementation, which would be a much longer document. Comparisons and contrasts are made to more usual approaches to V&V, and techniques new to the low-temperature plasma community are introduced. Specific challenges relating to the sufficiency of available data (e.g., cross sections), the limits of ad hoc validation approaches, the additional difficulty of utilizing a stochastic simulation tool, and the extreme cost of formal validation are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗