Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical integration”

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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Proton-conducting oxides for energy conversion and storage

Proton-conducting oxides are a class of solid-state ion-conducting ceramic materials that demonstrate significant hydrogen ion (proton) conductivity at intermediate temperatures (e.g., 300–700 °C). They are garnering significant attention due to several unique characteristics that distinguish them from both higher temperature oxygen ion conducting oxides and lower temperature proton-conducting polymers. By enabling proton-mediated electrochemistry under both dry and wet environments at moderate temperatures, protonic ceramics provide unique opportunities to enhance or synergize a diverse range of complementary electrochemical and thermochemical processes. Because of this potential, significant efforts have been devoted to advancing numerous energy-related applications using these materials. This review aims to comprehensively summarize these applications and analyze the most up-to-date and future developments of proton-conducting oxides. We aim to bring together this diverse subject matter by integrating the fundamentals of proton-conducting oxides with application-oriented insights. We begin with a historical roadmap, followed by a basic overview of the materials, theories and fundamentals, and fabrication and processing technologies underlying the field. The central section of our review summarizes major applications and developments of proton-conducting ceramics, ranging from maturing applications approaching commercialization to embryonic technologies just now emerging from the lab. These include protonic ceramic fuel cells, protonic ceramic electrolysis cells, reversible protonic ceramic electrochemical cells, protonic ceramic membrane reactors, and protonic ceramic electrochemical reactors. For each application, we analyze both the prospects and challenges and offer recommendations for future research directions so that tomorrow's researchers can continue to advance the development and commercialization of these fascinating materials.

Duan, Chuancheng (ORCID:0000000218261415)↗
DOE OSTI2021

Decentralized Failure-Tolerant Optimization of Electric Vehicle Charging

We present a decentralized failure-tolerant algorithm for optimizing electric vehicle (EV) charging, using charging stations as computing agents. The algorithm is based on the alternating direction method of multipliers (ADMM) and it has the following features: (i) It handles capacity, peak demand, and ancillary services coupling constraints. (ii) It does not require a central agent collecting information and performing coordination (e.g. an aggregator), instead all agents exchange information and computations are carried out in a fully decentralized fashion. (iii) It can withstand the failure of any number of computing agents, as long as the remaining computing agents are in a connected communications network. We construct this algorithm by reformulating the optimal EV charging problem in a decomposable form, amenable to ADMM, and then developing efficient decentralized solution methods for the subproblems dealing with coupling constraints. We conduct numerical experiments on industry-scale synthetic EV charging datasets, with up to 1,152 charging stations, using a high performance computing cluster. The experiments demonstrate that the proposed algorithm can solve the optimal EV charging problem fast enough to permit the integration of EV charging with real-time electricity markets, even in the presence of failures.

42 ENGINEERING↗
DOE OSTI2023

gpde

The R package gpde is an R implementation of a probabilistic differential equation solver. The package has a basic implementation of the Chkrebtii (2013) probabilistic ODE solver and an extension using the connection between Kalman Filtering for SDEs and Gaussian Processes from work by Simo Sarkka. The package also implements the Integrated Wiener covariance function from Schober 2014 to obtain probabilistic ODE solutions with the Runge-Kutta solution as the mean result. This collection of tools has been built into an R package for easy use at LANL and in the probabilistic numerics community.

Grosskopf, Michael↗
DOE OSTI2025

Designing a Framework for Solving Multiobjective Simulation Optimization Problems

Multiobjective simulation optimization (MOSO) problems are optimization problems with multiple conflicting objectives, where evaluation of at least one of the objectives depends on a black-box numerical code or real-world experiment, which we refer to as a simulation. Whereas an extensive body of research is dedicated to developing new algorithms and methods for solving these and related problems, it is challenging and time-consuming to integrate these techniques into real-world production-ready solvers. This is partly because of the diversity and complexity of modern state-of-the-art MOSO algorithms and methods and partly because of the complexity and specificity of many real-world problems and their corresponding computing environments. The complexity of this problem is only compounded when introducing potentially complex and/or domain-specific surrogate-modeling techniques, problem formulations, design spaces, and data acquisition functions. Here, this paper carefully surveys the current state of the art in MOSO algorithms, techniques, and solvers, as well as problem types and computational environments where MOSO is commonly applied. We then present several key challenges in the design of a parallel multiobjective simulation optimization framework (ParMOO) and how they have been addressed. Finally, we provide two case studies demonstrating how customized ParMOO solvers can be quickly built and deployed to solve real-world MOSO problems.

engineering design optimization↗
DOE OSTI2026

Telomere-to-telomere assemblies of chromosome 10 reveal complex adaptive variation of 3-ketoacyl-CoA-synthases in Populus trichocarpa likely driven by Helitrons

The model woody plant Populus trichocarpa displays an atypical alkene-diverse wax cuticle likely driven by copy number variation (CNV) of 3-ketoacyl-CoA synthases ( KCS ), which has been difficult to confirm with short-read assemblies. Long-read sequencing enables the development of telomere-to-telomere resources to detect cryptic variation, including CNVs, which are currently missed. Integrating this information can improve genomic prediction for breeding and provide insights into the evolutionary basis of important traits. Our analysis of 78 long-read haplotypes from chromosome 10 identified more than twice as many KCS genes as previously reported, and numerous intragenic non-synonymous substitutions. Random Forest predictive models highlighted the importance of Potri.010G079500 in producing very long chain alkenes; however, its absence did not predict previously reported alkene-deficient phenotypes. Instead, alkene levels are best predicted by the combinations of KCS copies. Additionally, amino acid substitutions clustered around ligand and donor binding pockets, suggesting they contribute to differing wax cuticle composition. Finally, each KCS gene and copy was linked to a Helitron transposon. A phylogenetic analysis suggests Helitrons are the evolutionary mechanism for generating KCS tandem arrays. Long-read generated telomere-to-telomere assemblies of P. trichocarpa chromosome 10 revealed large-effect loci critical to genetic studies that are unattainable from short-reads. This new resource produced novel insights into genome structure and function, and a novel mechanism for generating tandem gene duplication. Our results highlight that, given current challenges in annotation and assembly, detailed and focused long-read sequences are key to interpreting complex genomic regions that contain tandem copy number variants.

09 BIOMASS FUELS↗
DOE OSTI2023

Modeling of standing wave RF cavities for tracking through multi-pass energy recovery linac

Short bunches, high current and multiple linac pass are all characteristics of Energy Recovery Linacs (ERLs), which may result in collective effects. They in turn, may affect the beam, degrading its quality, or even yield to instabilities causing a beam loss. To study and mitigate these effects one needs a numerical simulation code, that can take into account both the collective effects, as well as, particular ERL features, such as a multi-turn design that does not reach a steady state or the multiple passages of the beam through Radio-Frequency (RF) cavities at different energies. CODAL [1], a code developed by SOLEIL in collaboration with IJCLab, enables such studies. It is a 6 dimensional (6D) tracking code applying 'kicks' based on the integration of the local Hamiltonian for each element of the lattice. It is also capable of simulating space charge, wakefields and coherent synchrotron radiation. However, to correctly take into account the ERL dynamics, an upgrade had to be made to include the effect of a standing wave RF cavity in 6D. In this paper, we will concentrate on the implementation and benchmarking (with DESY’s tracking code ASTRA [2]) of both the longitudinal and the transverse models (by J.B. Rosenzweig and L. Serafini [3]), which we use to carry out tracking of fully analytical 6D RF cavity.

Bogacz, A.↗
DOE OSTI2021

Integrating computation and experiment to investigate photoelectrodes for solar water splitting at the microscopic scale.

CONSPECTUS: Photoelectrochemical water-splitting is a promising and sustainable way to store the energy of the sun in chemical bonds and use it to produce hydrogen gas, a clean fuel. The key components in photoelectrochemical cells (PECs) are photoelectrodes, including a photocathode that reduces water to hydrogen gas and a photoanode that oxidizes water to oxygen gas. Materials used in photoelectrodes for PECs must effectively absorb sunlight, yield photogenerated carriers, and exhibit electronic properties that enable the efficient shuttling of carriers to the surface to participate in relevant water-splitting reactions. Discovering and understanding the key characteristics of optimal photoelectrode materials is paramount to the realization of PEC technologies. Oxide-based photoelectrodes can satisfy many of these materials requirements, including stability in aqueous environments, band edges with reasonable alignment with the redox potentials for water splitting, and ease of synthesis. However, oxide photoelectrodes generally suffer from poor charge transport properties and considerable bulk electron-hole separation, and they have relatively large band gaps. Numerous strategies have been proposed to improve these aspects and understand how these improvements are reflected in the photoelectrochemical performance. Unfortunately, the structural and compositional complexity of multinary oxides accompanied by the inherent complexity of photoelectrochemical processes makes it challenging to understand the individual effects of composition, structure, and defects in the bulk and on the surface on a material's photoelectrochemical properties. The integration of experiment and theory has great potential to increase our atomic-level understanding of structure-composition-property relationships in oxide photoelectrodes. In this Account, we describe how integrating experiment and theory is beneficial for achieving scientific insights at the microscopic scale. We highlight studies focused on understanding the role of (i) bulk composition via solid-state solutions, intercalation, and comparison with isoelectronic compounds, (ii) dopants for both the anion and cation and their interactions with oxygen vacancies, and (iii) surface/interface structure in the photocurrent generation and photoelectrochemical performance in oxide photoelectrodes. In each instance, we outline strategies and considerations for integrating experiment and theory and describe how this integration led to valuable insights and new directions in uncovering structure-composition-property relationships. Our aim is to demonstrate the unique value of combining experiment and theory in studying photoelectrodes and to encourage the continued effort to bring experiment and theory in closer step with each other.

Wang, Wennie↗
DOE OSTI2024

Strong Surface-Enhanced Coherent Phonon Generation in van der Waals Materials

Terahertz (THz) coherent phonons have emerged as promising candidates for the next generation of high-speed, low-energy information carriers in atomically thin phononic or phonon-integrated on-chip devices. However, effectively manipulating THz coherent phonons remains a significant challenge. Here, in this study, we investigated THz coherent phonon generation in exfoliated van der Waals (vdW) flakes of Fe 3 GeTe 2 , Fe 5 GeTe 2 , and FePS 3 . We successfully generated the THz A 1g coherent phonon mode in these vdW flakes. An innovative approach involved partially exfoliating vdW flakes on a gold substrate and partially on a silicon (Si) substrate to compare the THz coherent phonon generation between both sides. Interestingly, we observed a significantly enhanced THz coherent phonon in the vdW/gold area compared with that in the vdW/Si area. Frequency-domain Raman mapping across the vdW flakes corroborated these findings. Numerical simulations further indicated a stronger enhanced surface field in vdW/gold structures than in vdW/Si structures. Consequently, we attribute the observed enhancement in THz coherent phonon generation to the increased surface field on the gold substrate. This enhancement was consistent across the three different vdW materials studied, suggesting the universality of this strategy. Our results hold promise for advancing the design of THz phononic and phonon-integrated devices.

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

A heat and mass transfer model for evaluation of damaged cesium chloride radiological sources

Radiological sources are vital to many applications, and typically contain byproducts from waste streams which decay over time to form complex mixtures of elements. Daughter products resulting from these decay processes can introduce complicating factors to the integrity and safety of the vessel that contains the radiological material. Sources containing cesium 137 ( 137 Cs) are of particular interest, because the decay of this isotope produces barium metal which reacts readily and exothermically with oxygen. This work employs physics-based numerical tools to examine the thermal response of a radiological source containing a mixture of cesium and barium in the event that vessel walls are damaged and atmospheric gases contact source material. A parametric study was conducted to determine the sensitivity of the response to various factors, including vessel geometry, source material age, the degree of vessel damage, and other parameters. Here, it was found that the peak temperatures that occurred within the source material strongly depend on these parameters, particularly vessel geometry and age, which determine whether a breach would be a relatively minor accident or a catastrophic incident. Finally, the model’s sensitivity to uncertain thermophysical properties is discussed.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Advances in solar forecasting: Computer vision with deep learning

Renewable energy forecasting is crucial for integrating variable energy sources into the grid. It allows power systems to address the intermittency of the energy supply at different spatiotemporal scales. To anticipate the future impact of cloud displacements on the energy generated by solar facilities, conventional modeling methods rely on numerical weather prediction or physical models, which have difficulties in assimilating cloud information and learning systematic biases. Augmenting computer vision with machine learning overcomes some of these limitations by fusing real-time cloud cover observations with surface measurements acquired from multiple sources. This Review summarizes recent progress in solar forecasting from multisensor Earth observations with a focus on deep learning, which provides the necessary theoretical framework to develop architectures capable of extracting relevant information from data generated by ground-level sky cameras, satellites, weather stations, and sensor networks. Overall, machine learning has the potential to significantly improve the accuracy and robustness of solar energy meteorology; however, more research is necessary to realize this potential and address its limitations.

14 SOLAR ENERGY↗
DOE OSTI2021

A learning-augmented approach for AC optimal power flow

Because of the high nonlinearity of AC optimal power flow (OPF), numerous efforts have been made in recent decades to find efficient methods. Machine learning (ML) has proven to significantly reduce the computational costs in many real-world problems. Thus, this paper develops a learning-augmented method for solving AC OPF, which integrates both power network equations and ML to yield near-optimal solutions. More specifically, ML models are developed to first predict bus voltage magnitudes and angles. Then, physics-based network equations are employed to calculate the power injection at different buses. Three ML algorithms, i.e., random forest, multi-target decision tree, and extreme learning machine, are explored and compared. To evaluate the efficiency of the proposed learning-augmented AC OPF solver, the MATPOWER Interior Point Solver is adopted as a baseline. Case studies on both 500-bus and 4918-bus test networks show that the proposed learning-augmented method has reduced the computational time by 15–100 times depending on the network size with a minimal loss in optimality.

42 ENGINEERING↗
DOE OSTI2023

An implicit, conservative and asymptotic-preserving electrostatic particle-in-cell algorithm for arbitrarily magnetized plasmas in uniform magnetic fields

Here, we introduce a new electrostatic particle-in-cell algorithm capable of using large timesteps compared to particle gyro-period under a uniform external magnetic field. The algorithm extends earlier electrostatic fully implicit PIC implementations with a new asymptotic-preserving particle-push scheme that allows timesteps much larger than particle gyroperiods. In the large-timestep limit, the integrator preserves all particle drifts, while recovering the full orbit for small timesteps. The scheme allows for a seamless, efficient treatment of particles with coexisting magnetized and unmagnetized species, and conserves energy and charge exactly without spoiling implicit solver performance. We demonstrate by numerical experiment with several problems of variable species magnetization (diocotron instability, modified two-stream instability, and drift instability) that orders of magnitude wall-clock-time speedups vs. the standard fully implicit electrostatic PIC algorithm are possible without sacrificing solution accuracy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Multimode theory of electron hole transverse instability

We present Vlasov–Poisson three-dimensional linear stability analysis of an initially planar electron hole structure, solving for the distribution function by integration along unperturbed orbits. The non-sinusoidal potential perturbation shape (parallel to $B$ ) is expanded in eigenfunctions of the adiabatic Poisson operator, generalizing the prior assumption of a rigid shift of the equilibrium. We show that the shiftmode is then modified by a second discrete mode plus an integral over a continuum of wave-like modes. A rigorous treatment shows that the continuum can be approximated effectively by a single mode that satisfies the external wave dispersion relation, thus making the perturbation a weighted sum of three modes. We find numerically the solution for the complex instability frequency, and the corresponding three mode amplitudes determining the perturbation eigenmode. This multimode analysis refines the accuracy of the prior single-mode results, giving slightly higher growth rates at most parameters, as expected from the extra mode shape freedom. Oscillating modes near stability boundaries have larger mode distortions which help explain particle-in-cell simulations that observe instability up to ${\sim }20$ % beyond the prior shiftmode thresholds, and narrowing of the perturbation. At high magnetic field, the multimode analysis predicts a reduction of the already small growth rate.

Physics↗
DOE OSTI2024

Geologic Characterization of the South Georgia Rift Basin for Source Proximal CO2 Storage

The project Geologic Characterization of the South Georgia Rift Basin for Source Proximal CO2 Storage is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. The South Carolina Research Foundation and partners evaluated the feasibility of CCS in the Jurassic/ Triassic (J / TR) saline formations of the buried Mesozoic South Georgia Rift (SGR) Basin that extends from South Carolina into Georgia. The J / TR sequence, based on preliminary assessment of limited geologic and geophysical data, appears to have both the appropriate areal extent and multiple horizons to permanently and safely store CO2 The presence of several igneous rock layers within the sequence may potentially provide adequate seals to prevent upward CO2 migration into the Coastal Plain aquifer systems. Approximately 81 kilometers of 2-D seismic reflection data were collected by Bay Geophysical, Inc. to explore a portion of the SGR located in southern Georgia. The 81 kilometers were divided into two lines approximately 40.5 kilometers each, with Line 1 intersecting Georgia well GGS 3457. Line 2 intersects Line 1 at the southern portion of Line 1 to maximize the extent of coverage away from GGS-3457 (a deep well drilled in the 1980s for oil and gas exploration). This well had a set of usable logs, including gamma and neutron logs that provided promising results related to CO2 storage. Results showed sandstone with porosity values greater than 10 percent and a thickness of 120 meters. The design of the seismic shot was to extrapolate information away from the well and to better define the extent of the SGR and the necessary reservoir and caprock for a successful CO2 injection. A numerical simulation model of CO2 Injection and migration was developed based on the geology log for the GGS-3457 well. The simulation model was used to investigate the feasibility of injecting 30 million metric tons of CO2 into SGR J / TA sediments and integrity of the diabase layers as seals to prevent CO2 migration.

2-D seismic↗
DOE OSTI2023

From Deposition to Encapsulation: Roll-to-roll manufacturing of organic light emitting devices for lighting (Final Report)

Organic light emitting devices (OLEDs) are promising solid state light sources due to their high efficiency, high color quality and flexible form factors. The key to enable low cost OLED lighting, is to rapidly fabricate thin film organic layers on a continuous flexible roll, called R2R processing. This project aims to investigate the feasibility of mass production of OLED lighting using the R2R process, from deposition to encapsulation. A high efficiency white OLED is fabricated on 10 cm-wide substrate rolls in a pilot R2R tool comprising of two different organic deposition: vacuum thermal evaporation (VTE) and organic vapor phase deposition (OVPD). A high quality encapsulation process to package OLEDs is demonstrated using an atomic layer deposition tool integrated to the R2R system without air exposure. The method to achieve ultrahigh deposition rates required by R2R processing is demonstrated by OVPD. Uniform organic semiconductor thin films grown by OVPD at rates as high as 50 Å/s are achieved. A comprehensive numerical model that is capable of simulating complex, multilayer WOLED structures is developed to provide an alternative to experimental iterations of OLED design and tests. A cost estimate on the R2R production of WOLEDs for lighting is developed. Assuming a WOLED luminance of 10 klm/m 2 , the cost of a WOLED light engine is anticipated to be $\$ 12.5$ /klm. With incremental reduction in material and driver costs and improved luminance, the cost of WOLED lighting can be reduced to $\$ 6.3$ /klm in the near term. These findings suggest OLED lighting can be volume manufactured by R2R vapor deposition methods with much reduced costs compared to current batch processing methods, potentially positioning WOLEDs for use in numerous premium lighting applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Towards reverse mode automatic differentiation of Kokkos-based codes

Derivative computation is a key component of optimization, sensitivity analysis, uncertainty quantification, and the solving of nonlinear problems. Automatic differentiation (AD) is a powerful technique for evaluating such derivatives, and in recent years, has been integrated into programming environments such as Jax, PyTorch, and TensorFlow to support derivative computations needed for training of machine learning models, facilitating wide-spread use of these technologies. The C++ language has become the de facto standard for scientific computing due to numerous factors, yet language complexity has made the wide-spread adoption of AD technologies for C++ difficult, hampering the incorporation of powerful differentiable programming approaches into C++ scientific simulations. This is exacerbated by the increasing emergence of architectures, such as GPUs, with limited memory capabilities and requiring massive thread-level concurrency. C++ AD tools must effectively use these environments to bring novel scientific simulations to next-generation DOE experimental and observational facilities. In this project, we investigated source transformation-based automatic differentiation using LLVM compiler infrastructure to automatically generate portable and efficient gradient computations of Kokkos-based code. We have demonstrated that our proposed strategy is feasible by investigating the usage of a prototype LLVM-based source transformation tool to generate gradients of simple functions made of sequences of simple Kokkos parallel regions. Speedups of up to 500x compared to Sacado were observed on NVIDIA V100 GPU.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Potential applications of microbial genomics in nuclear non-proliferation

As nuclear technology evolves in response to increased demand for diversification and decarbonization of the energy sector, new and innovative approaches are needed to effectively identify and deter the proliferation of nuclear arms, while ensuring safe development of global nuclear energy resources. Preventing the use of nuclear material and technology for unsanctioned development of nuclear weapons has been a long-standing challenge for the International Atomic Energy Agency and signatories of the Treaty on the Non-Proliferation of Nuclear Weapons. Environmental swipe sampling has proven to be an effective technique for characterizing clandestine proliferation activities within and around known locations of nuclear facilities and sites. However, limited tools and techniques exist for detecting nuclear proliferation in unknown locations beyond the boundaries of declared nuclear fuel cycle facilities, representing a critical gap in non-proliferation safeguards. Microbiomes, defined as “characteristic communities of microorganisms” found in specific habitats with distinct physical and chemical properties, can provide valuable information about the conditions and activities occurring in the surrounding environment. Microorganisms are known to inhabit radionuclide-contaminated sites, spent nuclear fuel storage pools, and cooling systems of water-cooled nuclear reactors, where they can cause radionuclide migration and corrosion of critical structures. Microbial transformation of radionuclides is a well-established process that has been documented in numerous field and laboratory studies. These studies helped to identify key bacterial taxa and microbially-mediated processes that directly and indirectly control the transformation, mobility, and fate of radionuclides in the environment. Expanding on this work, other studies have used microbial genomics integrated with machine learning models to successfully monitor and predict the occurrence of heavy metals, radionuclides, and other process wastes in the environment, indicating the potential role of nuclear activities in shaping microbial community structure and function. Results of this previous body of work suggest fundamental geochemical-microbial interactions occurring at nuclear fuel cycle facilities could give rise to microbiomes that are characteristic of nuclear activities. These microbiomes could provide valuable information for monitoring nuclear fuel cycle facilities, planning environmental sampling campaigns, and developing biosensor technology for the detection of undisclosed fuel cycle activities and proliferation concerns.

59 BASIC BIOLOGICAL SCIENCES↗