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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Tilted lidar profiling: Development and testing of a novel scanning strategy for inhomogeneous flows

The most common profiling techniques for the atmospheric boundary layer based on a monostatic Doppler wind lidar rely on the assumption of horizontal homogeneity of the flow. This assumption breaks down in the presence of either natural or human-made obstructions that can generate significant flow distortions. The need to deploy ground-based lidars near operating wind turbines for the American WAKE experimeNt (AWAKEN) spurred a search for novel profiling techniques that could avoid the influence of the flow modifications caused by the wind farms. With this goal in mind, two well-established profiling scanning strategies have been retrofitted to scan in a tilted fashion and steer the beams away from the more severely inhomogeneous region of the flow. Results from a field test at the National Renewable Energy Laboratory's 135-m meteorological tower show that the accuracy of the horizontal mean flow reconstruction is insensitive to the tilt of the scan, although higher-order wind statistics are severely deteriorated at extreme tilts mainly due to geometrical error amplification. A numerical study of the AWAKEN domain based on the Weather Research and Forecasting Model and large-eddy simulation are also conducted to test the effectiveness of tilted profiling. It is shown that a threefold reduction of the error on inflow mean wind speed can be achieved for a lidar placed at the base of the turbine using tilted profiling.

17 WIND ENERGY↗
DOE OSTI2023

Online Adaptive Algorithm for Constraint Energy Minimizing Generalized Multiscale Discontinuous Galerkin Method

Here in this research, we propose an online basis enrichment strategy within the framework of a recently developed constraint energy minimizing generalized multiscale discontinuous Galerkin method. Combining the technique of oversampling, one makes use of the information of the current residuals to adaptively construct basis functions in the online stage to reduce the error of multiscale approximation. A complete analysis of the method is presented, which shows the proposed online enrichment leads to a fast convergence from multiscale approximation to the fine-scale solution. The error reduction can be made sufficiently large by suitably selecting oversampling regions and the number of oversampling layers. Further, the convergence rate of the enrichment algorithm depends on a factor of exponential decay regarding the number of oversampling layers and a user-defined parameter. Numerical results are provided to demonstrate the effectiveness and efficiency of the proposed online adaptive algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Quandary

Quandary numerically simulates and optimizes the time-evolution of open quantum systems. The underlying dynamics are modelled by Lindblad's master equation, a linear ordinary differential equation (ODE) describing quantum systems interacting with the environment. Quandary solves this ODE numerically by applying a time-stepping integration scheme, and utilizes a gradient-based optimization approach to determine optimal control pulses that drive the quantum system to a desired target state. Two optimization objectives are considered: (a) Unitary gate optimization that finds controls to realize a unitary gate transformation, and (b) optimal reset that aims to drive the quantum system to the ground states. Gradient-based optimization schemes utilizing Petsc's Tao optimization package are applied to generate control pulses that minimize the respective measure. To evaluate the gradient of the objective function, the discrete adjoint method is used while leveraging techniques from Algorithmic Differentiation to produce exact and consistent gradients. To mitigate excessive execution run times, the software can be build together with the XBraid software library which provides a parallelization strategy to distribute the time-evolution of the underlying dynamics onto multiple processor.

Petersson, NilsA.↗
DOE OSTI2021

Design of flexible proton beam imaging energy spectrometers (PROBIES)

A simple repeating step-filtered structure has been designed and deployed for the detection of MeV-energy proton beams produced from high-intensity short-pulse laser-driven experiments. Due to the diagnostic's ability to recover energy-resolved images of proton beams, it is called the proton beam imager and energy spectrometer (or PROBIES). The flexible design has two useful applications. First, it can be fielded in a configuration with a scintillator and imaging setup in order to collect proton beam information in high-repetition rate experiments. Here we show a configuration that collects the equivalent of a nine-film radiochromic film (RCF) stack, and a technique for using a neural network (NN) to rapidly analyze the images on millisecond time-scales. Second, when using RCF stacks as the detector, the energy resolution capability is greatly increased compared to traditional RCF stacks, while simultaneously reducing the total number of films required to make each measurement. We demonstrate a configuration that provides up to 90 spatially- and energy-resolved datapoints using just ten films. Both designs presented here can be easily customized for anticipated experimental outputs and manufactured at low cost. Such designs allow for rapid prototyping using modern high-accuracy 3D printers or on-demand computer numerical control machining to reduce the reliance on time-consuming film media, analysis, or bulky electro-magnetically dispersing spectrometers. The analysis of data is discussed and a methodology for constructing a NN that is capable of rapidly analyzing data on-the-fly with greater than 95% accuracy is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Generation of high-winding-number superfluid circulation in Bose-Einstein condensates

We experimentally and numerically demonstrate a method to generate multiply quantized superfluid circulation about an obstacle in highly oblate Bose-Einstein condensates (BECs). We experimentally achieve pinned superflow with winding numbers as high as 11, which persists for at least 4 s. Our method conceptually involves spiraling a blue-detuned laser beam, around and towards the center of the BEC, and is experimentally implemented by moving the BEC in a spiral trajectory around a stationary laser beam. This optical potential serves first as a repulsive stirrer to initiate superflow, and then as a pinning potential to transport the superfluid circulation within the BEC. The spiral technique can be used either to generate a high-winding-number persistent current, or for controlled placement of a cluster of singly quantized vortices of the same circulation. Thus, the technique may serve as a building block in experimental architectures to create on-demand vortex distributions in BECs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Collisional simulations of the modulator section in coherent electron cooling

The first section of any coherent electron cooling (CeC) system is the modulator, where the density of the electron beam is modulated by the copropagating ion beam. This density modulation is a result of Coulomb collisions between the individual particles of the two beams. The pairwise, stochastic part of the interactions impacts the overall performance of the CeC process. We present the first simulations of the density modulations of the electron beams from a collisional picture of the dynamics, considering the proof-of-principle CeC experiments at Brookhaven National Laboratory. These simulations were performed using PHAD, which is the first efficient, large-scale collisional numerical method in beam physics that we have previously developed and benchmarked. Realistic beam distributions and external fields have been optimized to provide strong modulation signals necessary for variations of coherent electron cooling systems. Cooling performance limits and potential collisionless simulation pitfalls are pointed out. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Second-order magnetic responses in quantum magnets: Magnetization under ac magnetic fields

Here we investigate second-order magnetic responses of quantum magnets against ac magnetic fields. We focus on the case where the z component of the spin is conserved in the unperturbed Hamiltonian and the driving field is applied in the xy plane. We find that linearly polarized driving fields induce a second-harmonic response, while circularly polarized fields generate only a zero-frequency response, leading to a magnetization with a direction determined by the helicity. Employing an unbiased numerical method, we demonstrate the nonlinear magnetic effect driven by the circularly polarized field in the XXZ model and show that the magnitude of the magnetization can be predicted by the dynamical spin structure factor in the linear response regime.

36 MATERIALS SCIENCE↗
DOE OSTI2024

PETSc/TAO Users Manual V.3.21

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Toolkit for Advanced Optimization (TAO) for the numerical solution of partial differential equations (PDEs) and related problems on high-performance computers. PETSc/TAO is a suite of data structures and routines that provide the building blocks for implementing large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication. PETSc/TAO includes a large suite of parallel linear solvers, nonlinear solvers, time integrators, and optimizers that may be used in application codes written in Fortran, C, C++, and Python (via petsc4py; see Getting Started ). The library is organized hierarchically, enabling users to employ the abstraction level most appropriate for a particular problem. By using techniques of object-oriented programming, PETSc provides enormous flexibility for users.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Simulation-Based Inference for Neutrino Interaction Model Parameter Tuning

High-energy physics experiments studying neutrinos rely heavily on simulations of their interactions with atomic nuclei. Limitations in the theoretical understanding of these interactions typically necessitate ad hoc tuning of simulation model parameters to data. Traditional tuning methods for neutrino experiments have largely relied on simple algorithms for numerical optimization. While adequate for the modest goals of initial efforts, the complexity of future neutrino tuning campaigns is expected to increase substantially, and new approaches will be needed to make progress. In this paper, we examine the application of simulation-based inference (SBI) to the neutrino interaction model tuning for the first time. Using a previous tuning study performed by the MicroBooNE experiment as a test case, we find that our SBI algorithm can correctly infer the tuned parameter values when confronted with a mock data set generated according to the MicroBooNE procedure. This initial proof-of-principle illustrates a promising new technique for next-generation simulation tuning campaigns for the neutrino experimental community.

Tame-Narvaez, Karla Maria [Fermilab]↗
DOE OSTI2025

Simulation-based inference for neutrino interaction model parameter tuning

High-energy physics experiments studying neutrinos rely heavily on simulations of their interactions with atomic nuclei. Limitations in the theoretical understanding of these interactions typically necessitate ad hoc tuning of simulation model parameters to data. Traditional tuning methods for neutrino experiments have largely relied on simple algorithms for numerical optimization. While adequate for the modest goals of initial efforts, the complexity of future neutrino tuning campaigns is expected to increase substantially, and new approaches will be needed to make progress. In this paper, we examine the application of simulation-based inference (SBI) to the neutrino interaction model tuning for the first time. Using a previous tuning study performed by the MicroBooNE experiment as a test case, we find that our SBI algorithm can correctly infer the tuned parameter values when confronted with a mock data set generated according to the MicroBooNE procedure. This initial proof-of-principle illustrates a promising new technique for next-generation simulation tuning campaigns for the neutrino experimental community.

Tame-Narvaez, Karla [Fermilab] (ORCID:000000022249↗
DOE OSTI2022

An Online Dynamic Amplitude-Correcting Gradient Estimation Technique to Align X-ray Focusing Optics

High-brightness X-ray pulses, as generated at synchrotrons and X-ray free electron lasers (XFELs), are used in a variety of scientific experiments. At these facilities, measurements often require optical equipment, e.g Compound Refractive Lenses (CRLs) to be precisely aligned and focused. The lateral alignment of CRLs to a beamline requires precise positioning along four axes: two translational, and the two rotational. At a synchrotron, alignment is often accomplished manually. However, XFEL beamlines present a beam brightness that fluctuates stochastically, making manual alignment a time-consuming endeavor. Automation using simplex or classic stochastic descent often fails, given the errant gradient estimates. Herein we present a dynamic-amplitude correction to the usual gradient based on the combination of a generalized finite difference stencil and a time-dependent sampling pattern. Intensity is recorded periodically, then used to normalize numerical derivatives against fluctuations. Error expectation is analyzed, and efficacy is demonstrated on classic benchmarks. We provide a proof of concept by laterally aligning optics on a simulated XFEL beamline using data recorded at both synchrotron and XFEL facilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Efficient kinetic thermal inverse modeling for organic material decomposition

The prevalent use of organic materials in manufacturing is a fire safety concern, and motivates the need for predictive thermal decomposition models. A critical component of predictive modeling is numerical inference of kinetic parameters from bench scale data. Currently, an active area of computational pyrolysis research focuses on identifying efficient, robust methods for optimization. This paper demonstrates that kinetic parameter calibration problems can successfully be solved using classical gradient-based optimization. We explore calibration examples that exhibit characteristics of concern: high nonlinearity, high dimensionality, complicated schemes, overlapping reactions, noisy data, and poor initial guesses. The examples demonstrate that a simple, non-invasive change to the problem formulation can simultaneously avoid local minima, avoid computation of derivative matrices, achieve a computational efficiency speedup of 10x, and make optimization robust to perturbations of parameter components. Techniques from the mathematical optimization and inverse problem communities are employed. By re-examining gradient-based algorithms, we highlight opportunities to develop kinetic parameter calibration methods that should outperform current methods.

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

Approximating accelerator impedances with resonator networks

It is common in the accelerator community to use the impedance of accelerator components to describe wake interactions in the frequency domain. However, it is often desirable to understand such wake interactions in the time domain in a general manner for excitations that are not necessarily Gaussian in nature. Here, the conventional method for doing this involves taking the inverse Fourier Transform of the component impedance, obtaining the Green's Function, and then convolving it with the desired excitation distribution. This method can prove numerically cumbersome, for a convolution integral must be evaluated for each individual point in time when the wake function is desired. An alternative to this method would be to compute the wake function analytically, which would sidestep the need for repetitive integration. Only a handful of cases, however, are simple enough for this method to be tenable. One of these cases is the case where the component in question is an RLC resonator, which has a closed-form analytical wake function solution. This means that a component which can be represented in terms of resonators can leverage this solution. As it happens, common network synthesis techniques may be used to map arbitrary impedance profiles to RLC resonator networks in a manner the accelerator community has yet to take advantage of. In this work, we will use Foster Canonical Resonator Networks and partial derivative descent optimization to develop a technique for synthesizing resonator networks that well approximate the impedances of real-world accelerator components. We will link this synthesis to the closed-form resonator wake function solution, giving rise to a powerful workflow that may be used to streamline beam dynamics simulations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Multiscale Ecosystem for solving Maxwell-Schrodinger equations of open quantum systems (OpenMS)

Light-matter interactions play an important role in many branches of physics, chemistry, energy, and materials science. In the strong coupling regime, light-matter interactions are able to tune the materials properties via the formation of new quasiparticles, such as plasmons and polaritons. However, theoretical and numerical modeling of the light-matter interaction-mediated processes remain a big challenge because light-matter interactions are fundamentally multiscale and multiphysics problems involving multiple interactions between electrons, nuclei, and photons at different time/length scales. In the current community, light-matter interactions were treated at different levels of theoretical complexity in quantum chemistry and quantum optics. In quantum optics or quantum photonics, the matter is usually simplified as a few-level system, and the light is treated quantum-mechanically. On the other hand, quantum chemistry explores first-principles methods, including both single-particle and many-body-based techniques, to describe the electronic properties of matter in detail. However, light is usually prescribed as a classical electromagnetic field, and the light-matter interaction is taken into account as an external potential via classical approximations. This software is designed to fill current modeling shortcomings by delivering a first-ever scalable multiscale platform for simulating light-matter interactions in realistic electromagnetic environments. The software solves Maxwell and Schrodinger equations self-consistent on the heterogeneous platforms. It implements HF/DFT, TDDFT, and coupled-cluster counterparts for light-matter interactions and adopts modular programming to offload massively parallel algorithms on a large number of CPU/GPUs.

Zhang, Yu↗
DOE OSTI2025

Advancements in Heat Transfer and Fluid Mechanics (Fundamentals and Applications)

Thermo-fluid science is a foundational discipline for numerous mechanical systems, particularly in energy production and building equipment, where thermal and mechanical energy transfer play critical roles. Advancements in heat transfer and fluid mechanics have significantly enhanced these systems, driving progress in associated market sectors. For instance, evaporator and condenser coils are essential for optimizing vapor compression cycles in building equipment. Heat pumps used for space and water heating constitute a major share of building systems, with microchannel heat exchanger technology at the forefront of these innovations. It is worth mentioning that advancements in heat transfer and fluid mechanics significantly alleviate the challenge of designing energy-efficient building equipment. The development of research techniques has contributed significantly to the advancement of science; for example, in the experimental field, non-intrusive measurement techniques such as Particle Image Velocimetry (PIV) are now capable of resolving flow behavior in a 3D format for different length and time scales.

42 ENGINEERING↗
DOE OSTI2021

Synthetic Absorption Lines from Simulations of Multiphase Gas in Galactic Winds

Supernova-driven galactic winds are multiphase streams of gas that are often observed flowing at a range of velocities out of star-forming regions in galaxies. In this study, we use high-resolution 3D simulations of multiphase galactic winds modeled with the hydrodynamics code Cholla to investigate the connection between numerical studies and observations. Using a simulated interaction between a hot $T$ ~ 10 4 K supernova-driven wind and a cool $T$ ~ 10 7 K cloud of interstellar material, we create mock observables, including the optical depth ($\tau$) and covering fraction ($C$ $f$ ) of six commonly observed ions (Si II , C II , Si IV , C IV , N V , and O VI ) as a function of gas velocity. We compare our mock observables to surveys of galactic winds in the literature, finding good agreement with velocities and profiles of the low ions. We then compute "empirical" values for the optical depth and covering fraction following observational techniques and compare them to the values calculated directly from the simulation data. We find that the empirically computed values tend to underestimate the "true" value of $\tau$ for ions with high optical depth and overestimate the "true" value of $\tau$ for ions with low optical depth relative to the simulated data. In conclusion, the empirically estimated covering fractions match our direct calculations very well for the low-ionization ions; for the high-ionization ions, the empirical covering fractions underestimate the directly computed values by up to ~40%.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

On differentiable local bounds preserving stabilization for Euler equations

This work presents the design of nonlinear stabilization techniques for the finite element discretization of Euler equations in both steady and transient form. Implicit time integration is used in the case of the transient form. A differentiable local bounds preserving method has been developed, which combines a Rusanov artificial diffusion operator and a differentiable shock detector. Nonlinear stabilization schemes are usually stiff and highly nonlinear. This issue is mitigated by the differentiability properties of the proposed method. Moreover, in order to further improve the nonlinear convergence, we also propose a continuation method for a subset of the stabilization parameters. The resulting method has been successfully applied to steady and transient problems with complex shock patterns. Numerical experiments show that it is able to provide sharp and well resolved shocks. Furthermore, the importance of the differentiability is assessed by comparing the new scheme with its non-differentiable counterpart. Numerical experiments suggest that, for up to moderate nonlinear tolerances, the method exhibits improved robustness and nonlinear convergence behavior for steady problems. Additionally, in the case of transient problem, we also observe a reduction in the computational cost.

42 ENGINEERING↗
DOE OSTI2025

A hereditary integral, transient network approach to modeling permanent set and viscoelastic response in polymers

An efficient numerical framework is presented for modeling viscoelasticity and permanent set of polymers. It is based on the hereditary integral form of transient network theory, in which polymer chains belong to distinct networks each with different natural equilibrium states. Chains continually detach from previously formed networks and reattach to new networks in a state of zero stress. The free energy of these networks is given in terms of the deformation gradient relative to the configuration at which the network was born . A decomposition of the kernel for various free energies allows for a recurrence relationship to be established, bypassing the need to integrate over all time history. The technique is established for both highly compressible and nearly incompressible materials through the use of neo-Hookean, Blatz–Ko, Yeoh, and Ogden-Hill material models. Multiple examples are presented showing the ability to handle rate-dependent response and residual strains under complex loading histories.

Finite element method↗