Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical approximation & analysis”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Exponential acceleration of macroscopic quantum tunneling in a Floquet Ising model

The exponential suppression of macroscopic quantum tunneling (MQT) in the number of elements to be reconfigured is an essential element of broken symmetry phases. This suppression is also a core bottleneck in quantum algorithms, such as traversing an energy landscape in optimization, and adiabatic state preparation more generally. In this work, we demonstrate exponential acceleration of MQT through Floquet engineering with the application of a uniform, high frequency transverse drive field. Using the ferromagnetic phase of the transverse field Ising model in one and two dimensions as a prototypical example, we identify three phenomenological regimes as a function of drive strength. For weak drives, the system exhibits exponentially decaying tunneling rates but robust magnetic order; in the crossover regime at intermediate drive strength, we find polynomial decay of tunnelling alongside vanishing magnetic order; and at very strong drive strengths both the Rabi frequency and time-averaged magnetic order are approximately constant with increasing system size. We support these claims with extensive full wavefunction and tensor network numerical simulations, and theoretical analysis. An experimental test of these results presents a technologically important and novel scientific question accessible on NISQ-era quantum computers.

Grattan, George↗
DOE OSTI2024

A Provably Accurate Randomized Sampling Algorithm for Logistic Regression

In statistics and machine learning, logistic regression is a widely-used supervised learning technique primarily employed for binary classification tasks. When the number of observations greatly exceeds the number of predictor variables, we present a simple, randomized sampling-based algorithm for logistic regression problem that guarantees high-quality approximations to both the estimated probabilities and the overall discrepancy of the model. Our analysis builds upon two simple structural conditions that boil down to randomized matrix multiplication, a fundamental and well-understood primitive of randomized numerical linear algebra. We analyze the properties of estimated probabilities of logistic regression when leverage scores are used to sample observations, and prove that accurate approximations can be achieved with a sample whose size is much smaller than the total number of observations. To further validate our theoretical findings, we conduct comprehensive empirical evaluations. Overall, our work sheds light on the potential of using randomized sampling approaches to efficiently approximate the estimated probabilities in logistic regression, offering a practical and computationally efficient solution for large-scale datasets.

Chowdhury, Agniva↗
DOE OSTI2022

Droplet formation simulation using mixed finite elements

Droplet formation happens in finite time due to the surface tension force. The linear stability analysis is useful to estimate the size of a droplet but fails to approximate the shape of the droplet. This is due to a highly nonlinear flow description near the point where the first pinch-off happens. A one-dimensional axisymmetric mathematical model was first developed by Eggers and Dupont [“Drop formation in a one-dimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] using asymptotic analysis. This asymptotic approach to the Navier–Stokes equations leads to a universal scaling explaining the self-similar nature of the solution. Numerical models for the one-dimensional model were developed using the finite difference [Eggers and Dupont, “Drop formation in a one-dimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] and finite element method [Ambravaneswaran et al., “Drop formation from a capillary tube: Comparison of one-dimensional and two-dimensional analyses and occurrence of satellite drops,” Phys. Fluids 14, 2606–2621 (2002)]. The focus of this study is to provide a robust computational model for one-dimensional axisymmetric droplet formation using the Portable, Extensible Toolkit for Scientific Computation. The code is verified using the Method of Manufactured Solutions and validated using previous experimental studies done by Zhang and Basaran [“An experimental study of dynamics of drop formation,” Phys. Fluids 7, 1184–1203 (1995)]. The present model is used for simulating pendant drops of water, glycerol, and paraffin wax, with an aspiration of extending the application to simulate more complex pinch-off phenomena.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

An asymptotic-preserving 2D-2P relativistic Drift-Kinetic-Equation solver for runaway electron simulations in axisymmetric tokamaks

We propose an asymptotic-preserving (AP), uniformly convergent numerical scheme for the relativistic collisional Drift-Kinetic Equation (rDKE) to simulate runaway electrons in axisymmetric toroidal magnetic field geometries typical of tokamak devices. The approach is derived from an exact Green's function solution with numerical approximations of quantifiable impact, and results in a simple, two-step operator-split algorithm, consisting of a collisional Eulerian step, and a Lagrangian orbit-integration step with analytically prescribed kernels. The AP character of the approach is demonstrated by analysis of the dominant numerical errors, as well as by numerical experiments. We demonstrate the ability of the algorithm to provide accurate answers regardless of plasma collisionality on a circular axisymmetric tokamak geometry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The eXtended virtual element method for elliptic problems with weakly singular solutions

This paper introduces a novel eXtended virtual element method, an extension of the conforming virtual element method. The X-VEM is formulated by incorporating appropriate enrichment functions in the local spaces. The method is designed to handle highly generic enrichment functions, including singularities arising from fractured domains. By achieving consistency on the enrichment space, the method is proven to achieve arbitrary approximation orders even in the presence of singular solutions. The paper includes a complete convergence analysis under general assumptions on mesh regularity, and numerical experiments validating the method’s accuracy on various mesh families, demonstrating optimal convergence rates in the L 2 - and H 1 - norms on fractured or L-shaped domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Forward and inverse modeling of fault transmissibility in subsurface flows

Characterizing physical properties of faults, such as their transmissibility, is crucial for performing predictive numerical simulation of subsurface flows, such as those encountered in petroleum engineering and remediation of subsurface contamination. Here, this paper provides a complete investigation of the inverse problem for fault transmissibility in subsurface flow models, under appropriate assumptions on fault structure. In particular, the following aspects are considered: 1) fault modeling and well-posedness of the forward problem; 2) finite element (FEM) discretizations of the forward problem and their rigorous a priori convergence analysis; 3) Well-posedness of the Bayesian inverse problem, FEM discretization of the infinite dimensional Bayesian inverse formulation, and its rigorous a priori analysis. Moreover, computation of the maximum a posteriori (MAP) point via fast inexact Newton-conjugate gradient optimization and a Laplace approximation of the Bayesian posterior are also presented. Numerical results illustrate the use of the proposed fault model in forward and inverse problems for subsurface flows in two dimensional domains with multiple faults.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Numerical methods for nonlocal and fractional models

Partial differential equations (PDEs) are used with huge success to model phenomena across all scientific and engineering disciplines. However, across an equally wide swath, there exist situations in which PDEs fail to adequately model observed phenomena, or are not the best available model for that purpose. On the other hand, in many situations, nonlocal models that account for interaction occurring at a distance have been shown to more faithfully and effectively model observed phenomena that involve possible singularities and other anomalies. Here, we consider a generic nonlocal model, beginning with a short review of its definition, the properties of its solution, its mathematical analysis and of specific concrete examples. We then provide extensive discussions about numerical methods, including finite element, finite difference and spectral methods, for determining approximate solutions of the nonlocal models considered. In that discussion, we pay particular attention to a special class of nonlocal models that are the most widely studied in the literature, namely those involving fractional derivatives. The article ends with brief considerations of several modelling and algorithmic extensions, which serve to show the wide applicability of nonlocal modelling.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Radiative heat transfer in FLiBe molten salt participating medium in a vertical heated tube under forced and mixed convection laminar flows

The contribution of radiative heat transfer (RHT) to convective heat transfer in a heated vertical pipe with laminar flow and constant wall temperature, for FLiBe molten salt, is investigated computationally. The combined effects of conduction in the fluid, forced convection, buoyancy, temperature-dependent physical properties, and thermal radiation are investigated. The P1 approximation is employed in the discretization of the Radiative Transfer Equation (RTE). The COMSOL Multiphysics software is used to generate the numerical solutions. The theoretical analysis developed here demonstrates that only for intermediate values of the optical thickness, i.e. τ D ~ O 1 , the participating media effects are expected to be important. This analysis is confirmed computationally and a value of τ D ≈ 4 is shown to lead to the highest increase in overall heat transfer behavior. Under forced convection, the Nusselt ratio N u total / N u no - rad reaches a peak value of 1.76 at τ D = 4 and z/D = 200 , and is less than 1.1 for τ D < 0 . 1 and τ D > 60 . Under aiding-flow mixed convection, N u total / N u no - rad reaches a peak of 1.34 at τ D = 4 . 2 and is less than 1.1 for τ D < 0 . 4 and τ D > 36 . Under opposing-flow mixed convection, Nu total / Nu no - rad reaches a peak of 1.81 at τ D = 4 . 2 and is less than 1.1 for τ D < 0 . 2 and τ D > 50 . RHT effects on the overall heat transfer are most pronounced in opposing mixed convection, where N u total / N u no - rad is 2.26 observed at z/D = 50 . The sensitivity to wall emissivity is evaluated; in forced convection, the peak N u total / N u no - rad is 1.66 at ε = 0 (reflective wall) and 2.04 at ε = 1 (absorptive wall). Finally, the sensitivity to pipe diameter is also discussed. Overall, for a vertical heated tube, radiative heat transfer effects lead to enhancement of heat transfer by as high as a factor of two, and they depend on the optical thickness of the flow, mixed convection environment ( Gr / R e 2 and direction of flow relative to the gravitational force), surface emissivity, and entrance effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

X-ray Computed Tomography of Irradiated and Unirradiated AGR-3/4 Compacts

X-ray Computed Tomography (XCT) has been utilized to image and characterize compacts from the combined third and fourth irradiation of the Advanced Gas Reactor (AGR) Program, AGR-3/4, fuel. The experiment contained tristructural isotropic (TRISO)-coated fuel particles as well as designed-to-fail (DTF) fuel particles. Two irradiated compacts, representing the lower and higher range of AGR-3/4 burnup (4.85% and 14.92% fissions per initial heavy metal atom FIMA) were examined. These represent the first known highly irradiated TRISO fuel compacts to be examined via X-ray CT. Additionally, two unirradiated compacts from the same production batch as the examined irradiation compacts were also imaged for a baseline comparison. As XCT of irradiated TRISO compacts is not a commonly implemented characterization technique, a significant portion of the report focuses on developed methodology and imaging conditions. A specialized sample shielding device was developed and fabricated specifically to limit received dose to staff during sample preparation for XCT and to minimize excess gamma radiation dose to sensitive electronic components with the utilized X-ray system. Significant penetration through the uranium oxycarbide fuel kernels by significantly hardening the X-ray beam with specialized proprietary filters acquired from Carl Zeiss NTS Ltd. The filter utilized resulted in an average X ray photon energy of ~110 keV which approaches uranium’s K-edge (~115 keV), maximizing penetration for a microfocus X-ray source. The gamma-radiation emitted from the irradiated AGR-3/4 TRISO compacts, has the same properties and mechanisms for interaction with matter as X-rays, thus the detection of gamma-radiation by the utilized X-ray detectors was initially a concern. However, although ?-rays did produce an observable signal on the X-ray detector, its contribution to the overall imaging results appeared negligible upon 3D reconstruction. The neglibile impact on the resulting 3D reconstructed volumes were likely the result of: (1) a significantly lower detection efficiency for ?-rays relative to X-rays; (2) An X-ray flux at the detector several orders of magnitude higher than that of the impinging ?-rays from the irradiated compacts. These results suggest that irradiated compacts with significantly higher radiation fields can be examined in the future if an acceptable route for sample handling and preparation can be determined. Additionally, the 3D imaging results of XCT can provide a valuable means of assessing compacts. While in many ways complimentary to traditional post irradiation examination techniques such as optical ceramography, XCT can provide additional insight into compact features traditionally difficult to discern directly from cross-sectional imaging alone. Preliminary analyses on kernel size, morphology (aspect ratio and sphericity), and kernel orientation were presented. Sphericity, a simple morphological shape descriptor, was utilized to screen for kernel extrusions within the high burnup compact. The number of kernel extrusions identified via XCT represented an approximate two-fold increase from the quantity of extruded particles observed (via optical ceramography) in adjacent compacts from the same irradiation capsule. While numerical analysis of the compact datasets was highly preliminary, initial results show promise for providing complimentary metrics to current AGR-3/4 PIE and potentially additional insight into the processes driving TRISO fuel degradation during reactor operation. Additional analyses to be performed at a later date include a more detailed examination of kernel size, kernel sphericity (and observed kernel extrusions), and sphericity. Given all particles can be observed in a single data volume possible correlation of spatial position with observed kernel features will also be made at a later date.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Dynamic Strength and Equation of State of Epon 828 and Diethanolamine (DEA) Polymer Epoxy Under Shock Loading

Polymers are increasingly utilized in engineering applications that can experience high loading rates, necessitating increased understanding of their response under such conditions. The tamped Richtmyer-Meshkov Instability (RMI) method was used to characterize the equation of state and dynamic strength of the polymer Epon 828 cured with Diethanolamine (DEA). Plate impact experiments that drove a uniaxial shock compression wave across a sinusoidally corrugated metal-polymer interface were performed to generate shock stresses from 4-12 GPa and strain rates of approximately 1/s in the polymer. X-ray phase contrast imaging recorded the shock motion in the polymer and subsequent interface evolution. Analysis of this data yielded the polymer equation of state and, in conjunction with numerical modeling, the dynamic strength. The equation of state was validated against one-dimensional plate impact experiments from existing literature. The dynamic strength was compared to prior data for Epon 828 and related polymers at lower strain rates and found to exhibit significant strain rate and pressure-hardening effects. The strength of the Epon 828 polymer at 10 6 1/s was found to be approximately 1.5 GPa, suggesting that it is comparable to the strength of high strength metals at these dynamic conditions.

Equation of state↗
DOE OSTI2021

Scaling patch analysis of turbulent planar plume

Proper scaling in turbulent planar plumes is investigated in this study using a scaling patch approach. Based on the scaled boundary conditions, a proper velocity scale for the mean axial flow is the plume centerline velocity U ref = U ctr , and a proper temperature scale for the temperature excess is Θ ref = T ctr – T ∞ , where T ctr is the plume centerline temperature and T ∞ is the ambient fluid temperature. By seeking an admissible scaling, a key concept in the scaling patch approach, for the mean continuity, mean momentum, and mean energy equations, respectively, the following is found: (1) a proper scale for the mean transverse flow is V ref = ( d δ / d x ) U ctr , where d δ / d x is the growth rate of the plume width. (2) A proper scale for the Reynolds shear stress is R vu , ref = U ctr V ref = ( d δ / d x ) U ctr 2 , a mix of the scales for the mean axial and transverse flows. (3) A proper scale for the turbulent heat flux is R v θ , ref = V ref Θ ctr , a mix of the scales for the mean transverse flow and mean temperature excess. The mean transverse flow thus plays a critical role in the scaling of turbulent planar plumes. Approximate functions are developed for the scaled mean transverse flow, Reynolds shear stress, and turbulent temperature flux, and are found to agree favorably with experimental and numerical simulation data. The integral analysis of the mean momentum equation yields a Richardson number Ri, which remains invariant in the axial direction. The Richardson number is defined as Ri = def g β Θ ctr δ t / ( U ctr V ref ) ≈ 1 / 2 , where g is the gravitational acceleration, β is the thermal expansion coefficient, and δ t is the plume half-width based on the mean temperature profile. This Richardson number arises directly from the scaling patch analysis of the mean momentum equation, including both the streamwise and transverse velocity scales.

42 ENGINEERING↗
DOE OSTI2025

Massively parallel and universal approximation of nonlinear functions using diffractive processors

Nonlinear computation is essential for a wide range of information processing tasks, yet implementing nonlinear functions using optical systems remains a challenge due to the weak and power-intensive nature of optical nonlinearities. Overcoming this limitation without relying on nonlinear optical materials could unlock unprecedented opportunities for ultrafast and parallel optical computing systems. Here, we demonstrate that large-scale nonlinear computation can be performed using linear optics through optimized diffractive processors composed of passive phase-only surfaces. In this framework, the input variables of nonlinear functions are encoded into the phase of an optical wavefront—e.g., via a spatial light modulator (SLM)—and transformed by an optimized diffractive structure with spatially varying point-spread functions to yield output intensities that approximate a large set of unique nonlinear functions–all in parallel. We provide proof establishing that this architecture serves as a universal function approximator for an arbitrary set of bandlimited nonlinear functions, also covering wavelength-multiplexed nonlinear functions as well as multi-variate and complex-valued functions that are all-optically cascadable. Our analysis also indicates the successful approximation of typical nonlinear activation functions commonly used in neural networks, including the sigmoid, tanh, ReLU (rectified linear unit), and softplus. We numerically demonstrate the parallel computation of one million distinct nonlinear functions, accurately executed at wavelength-scale spatial density at the output of a diffractive optical processor. Furthermore, we experimentally validated this framework using in situ optical learning and approximated 35 unique nonlinear functions in a single shot using a compact setup consisting of an SLM and an image sensor. These results establish diffractive optical processors as a scalable platform for massively parallel universal nonlinear function approximation, paving the way for new capabilities in analog optical computing based on linear materials.

Rahman, Md Sadman Sakib [University of California,↗
DOE OSTI2024

Efficient sensitivity analysis of the thermal profile in powder bed fusion of metals using hypercomplex automatic differentiation finite element method

Rapid cyclic temperature fluctuation occurring in powder bed fusion of metals using a laser beam (PBF-LB/M) influences the formation of flaws in printed parts. Consequently, there is a pressing need to enhance the quality of printed parts by developing innovative methodologies that can predict thermal histories and help uncover the intricate relationships between process parameters and thermal profiles. Sensitivity Analysis (SA) emerges as an essential tool for this, offering the potential for process optimization and enhanced quality control. Nonetheless, conventional SA methodologies often incur in excessive computational costs and potential numerical approximation errors. Here, to address this technical challenge, we present a novel method for SA that integrates the HYPercomplex-based Automatic Differentiation (HYPAD) technique with transient thermal simulations conducted via the finite element method (FEM). Leveraging this methodology, we efficiently and accurately perform SA for PBF-LB/M processes in a post-processing step. Compared to traditional methods like Finite Differences (FD), HYPAD-FEM required 96 % less computational time for obtaining sensitivities for 22 process parameters, under a comparative study conducted within the context of the 2018–02 AM benchmark of the National Institute of Standards and Technology. In summary, HYPAD-FEM offers superior efficiency and accuracy in SA over conventional methods, delivering the best sensitivity of a model without the need for step-size selection and problem or parameter-based implementations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multiscale analysis in solids with unseparated scales: fine-scale recovery, error estimation, and coarse-scale adaptivity

There are several engineering applications in which the assumptions of homogenization and scale separation may be violated, in particular, for metallic structures constructed through additive manufacturing. Instead of resorting to direct numerical simulation of the macroscale system with an embedded fine scale, an alternative approach is to use an approximate macroscale constitutive model, but then estimate the model-form error using a posteriori error estimation techniques and subsequently adapt the macroscale model to reduce the error for a given boundary value problem and quantity of interest. Here, we investigate this approach to multiscale analysis in solids with unseparated scales using the example of an additively manufactured metallic structure consisting of a polycrystalline microstructure that is neither periodic nor statistically homogeneous. As a first step to the general nonlinear case, we focus here on linear elasticity in which each grain within the polycrystal is linear elastic but anisotropic.

42 ENGINEERING↗
DOE OSTI2024

Quasiparton distributions in massive QED2: Toward quantum computation

We analyze the quasiparton distributions of the lightest 𝜂′ meson in massive two-dimensional quantum electrodynamics (QED2) by exact diagonalization. The Hamiltonian and boost operators are mapped onto spin qubits in a spatial lattice with open boundary conditions. The lowest excited state in the exact diagonalization is shown to interpolate continuously between an anomalous 𝜂′ state at strong coupling, and a nonanomalous heavy meson at weak coupling, with a cusp at the critical point. The boosted 𝜂′ state follows relativistic kinematics but with large deviations in the luminal limit. The spatial quasiparton distribution function and amplitude for the 𝜂′ state are computed numerically for increasing rapidity both at strong and weak coupling, and compared to the exact light front results. The numerical results from the boosted form of the spatial parton distributions, compare fairly with the inverse Fourier transformation of the luminal parton distributions, derived in the lowest Fock space approximation. Our analysis points out some of the limitations facing the current lattice program for the parton distributions.

Lattice field theory↗
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 OSTI2020

Quantum approximate optimization of the long-range Ising model with a trapped-ion quantum simulator

Quantum computers and simulators may offer significant advantages over their classical counterparts, providing insights into quantum many-body systems and possibly improving performance for solving exponentially hard problems, such as optimization and satisfiability. Here, we report the implementation of a low-depth Quantum Approximate Optimization Algorithm (QAOA) using an analog quantum simulator. We estimate the ground-state energy of the Transverse Field Ising Model with long-range interactions with tunable range, and we optimize the corresponding combinatorial classical problem by sampling the QAOA output with high-fidelity, single-shot, individual qubit measurements. We execute the algorithm with both an exhaustive search and closed-loop optimization of the variational parameters, approximating the ground-state energy with up to 40 trapped-ion qubits. We benchmark the experiment with bootstrapping heuristic methods scaling polynomially with the system size. We observe, in agreement with numerics, that the QAOA performance does not degrade significantly as we scale up the system size and that the runtime is approximately independent from the number of qubits. We finally give a comprehensive analysis of the errors occurring in our system, a crucial step in the path forward toward the application of the QAOA to more general problem instances.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗