Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical Methods”

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,189 records · Page 66

Diagnosing nuclear power plant pipe wall thinning due to flow accelerated corrosion using a passive, thermal non-destructive evaluation method: Feasibility assessment via numerical experiments

Flow accelerated corrosion (FAC) in nuclear power plant pipes is one of the leading causes of accidents, fatalities, damage and outages. Current FAC identification methods employ expensive sensing technology and are “active” methods, where the response of the piping system to an externally-generated thermal, mechanical or optical excitation must be measured. As a result, these techniques require a disruptive and time-consuming setup. Here we propose a method that utilizes pipe surface temperature measurements to passively monitor for FAC-induced pipe wall thinning without the need for expensive equipment or post-installation setup time. This diagnostic method utilizes a simulation data-driven diagnostic model to estimate the amount of thickness reduction in a pipe based on changes in measured steady-state pipe temperatures. In order to reduce the computational burden of generating large, simulation-based datasets, the behavior of the insulation of the pipe was modeled using a suitably calibrated heat transfer boundary parameter. Additionally, global sensitivity analysis was performed to determine system parameter(s), such as the temperature of water flowing inside the pipe, which significantly affect the steady state pipe wall temperature and could cause errors in diagnosis. Two diagnostic models, one using only the change in steady-state temperature as an indicator for FAC-induced pipe wall thinning and the other using water temperature as an additional diagnostic model input were evaluated for their ability to estimate thickness reductions in a pipe using simulated pipe wall temperature data. For the numerical experiments conducted in this work, both models estimated wall thickness with errors within 0.5 mm, indicating that the proposed technique can potentially be used as a low-cost, first-pass method for FAC monitoring.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Numerical Assessment of SNF Parameter Bandwidths and Muon Tomography Assay for Large Volume Casks - 18185

For the characterization of spent fuel and for the safe direct disposal of SNF we present: a) a method of numerical assessment of their radionuclides inventory that is derived from the initial fuel characteristics and individual burn-up / cooling time history of that particular fuel assembly, and b) how the numerical information can be experimentally assessed using muon tomography as a nondestructive assay of large and heterogeneous SNF casks. We combine both complementary methods, i.e. numerical assessment and muon tomography (MT) experimental assay, to provide a powerful tool for verifying the nuclear inventory without compromising the casks or radioactive material containers. Here, we shall present the numerical assessment technique and preliminary but very promising experimental results, and explain the scope of a versatile and mobile MT-detection system for the identification of nuclear material inside large casks containing spent fuel assemblies and/or other large volume heterogeneous material containers for which the European Commission has recently granted R and D funding within the CHANCE project to be developed. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Two-Loop Four-Gluon Amplitudes from Numerical Unitarity

We present the first numerical computation of two-loop amplitudes based on the unitarity method. As a proof of principle, we compute the four-gluon process in the leading-color approximation. We discuss the new method, analyze its numerical properties, and apply it to reconstruct the analytic form of the amplitudes. The numerical method is universal, and can be automated to provide multiscale two-loop computations for phenomenologically relevant signatures at hadron colliders.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Kolmogorov-Arnold wavefunctions

Here, this work investigates Kolmogorov-Arnold network-based (KAN) wave-function Ansätz as viable representations for quantum Monte Carlo simulations. Through systematic analysis of one-dimensional model systems, we evaluate their computational efficiency and representational power against established methods. Our numerical experiments suggest some efficient training methods and we explore how the computational cost scales with desired precision, particle number, and system parameters. Roughly speaking, KANs seem to be 10 times cheaper computationally than other neural-network-based Ansätz . We also introduce a novel approach for handling strong short-range potentials—a persistent challenge for many numerical techniques—which generalizes efficiently to higher-dimensional, physically relevant systems with short-ranged strong potentials common in atomic and nuclear physics.

1-dimensional systems↗

Planck 2015 results: VI. LFI mapmaking

This article describes the mapmaking procedure applied to Planck Low Frequency Instrument (LFI) data. The mapmaking step takes as input the calibrated timelines and pointing information. The main products are sky maps of I, Q, and U Stokes components. For the first time, we present polarization maps at LFI frequencies. The mapmaking algorithm is based on a destriping technique, which is enhanced with a noise prior. The Galactic region is masked to reduce errors arising from bandpass mismatch and high signal gradients. We apply horn-uniform radiometer weights to reduce the effects of beam-shape mismatch. The algorithm is the same as used for the 2013 release, apart from small changes in parameter settings. We validate the procedure through simulations. Special emphasis is put on the control of systematics, which is particularly important for accurate polarization analysis. We also produce low-resolution versions of the maps and corresponding noise covariance matrices. These serve as input in later analysis steps and parameter estimation. The noise covariance matrices are validated through noise Monte Carlo simulations. The residual noise in the map products is characterized through analysis of half-ring maps, noise covariance matrices, and simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Umbrella sampling: a powerful method to sample tails of distributions

We present the umbrella sampling (US) technique and show that it can be used to sample extremely low-probability areas of the posterior distribution that may be required in statistical analyses of data. In this approach sampling of the target likelihood is split into sampling of multiple biased likelihoods confined within individual umbrella windows. We show that the US algorithm is efficient and highly parallel and that it can be easily used with other existing Markov Chain Monte Carlo (MCMC) samplers. The method allows the user to capitalize on their intuition and define umbrella windows and increase sampling accuracy along specific directions in the parameter space. Alternatively, one can define umbrella windows using an approach similar to parallel tempering. We provide a public code that implements US as a standalone PYTHON package. We present a number of tests illustrating the power of the US method in sampling low-probability areas of the posterior and show that this ability allows a considerably more robust sampling of multimodal distributions compared to standard sampling methods. We also present an application of the method in a real world example of deriving cosmological constraints using the supernova type Ia data. We show that US can sample the posterior accurately down to the approximate to 15 sigma credible region in the Omega(m) - Omega(Lambda) plane, while for the same computational effort the affine-invariant MCMC sampling implemented in the emcee code samples the posterior reliably only to approximate to 3 sigma.

cosmology: cosmological parameters↗

Thermal Ratcheting Analysis of TEDS Packed-bed Thermocline Energy Storage Tank - Modeling Methodology and Data Validation

This report investigates numerical modeling methods for thermal ratcheting analysis of packed-bed thermal energy storage (TES) tank and discusses the validation results via comparison with experimental data. The experimental data obtained from various design characteristics of packed-bed thermocline tanks, including the Thermal Energy Distribution System (TEDS) TES tank at Idaho National Laboratory (INL), were used to validate thermal and mechanical models developed in this study to evaluate the thermal ratcheting potential. The thermal model was shown to predict the transient thermal propagation through the packed-bed thermocline tanks generally well. However, a larger discrepancy was observed during the comparison with the data from TEDS, presumably due to the uncertainty of boundary conditions given from the experiment. Based on the comparative study between the thermal model predictions and experimental data of various packed-bed thermocline tanks, potential improvements were suggested for the future TEDS experiments for more precise validation study. For mechanical (thermally induced stress) analysis, two different modeling approaches were tested to evaluate hoop stress applied to the packed-bed TES tank wall, which is a major cause of thermal ratcheting process: (i) infinite rigidity model and (ii) Drucker-Prager (DP) model. The ‘model (i)’ is a conservative method with infinite rigidity assumption of granular filler inside a TES tank, whereas the ‘model (ii)’ is a method that takes into account more realistic processes such as thermal expansion of filler and tank wall as well as inter-particle interactions during the cyclic operation of a packed-bed TES tank. The validity of each modeling method was examined by comparing the numerical simulation with the experimental data obtained from the packed-bed TES tank for Solar One pilot plant. Then, the effects of various model parameters were discussed to evaluate the thermal ratcheting potential of the TEDS TES tank. The preliminary thermal ratcheting analysis implies that the TEDS TES tank will hold its structural integrity during the normal operation cycles.

25 ENERGY STORAGE↗

Application of high-order numerical schemes and Newton-Krylov method to two-phase drift-flux model

This study concerns the application and solver robustness of the Newton-Krylov method in solving two-phase flow drift-flux model problems using high-order numerical schemes. In our previous studies, the Newton-Krylov method has been proven as a promising solver for two-phase flow drift-flux model problems. However, these studies were limited to use first-order numerical schemes only. Moreover, the previous approach to treating the drift-flux closure correlations was later revealed to cause deteriorated solver convergence performance, when the mesh was highly refined, and also when higher-order numerical schemes were employed. In this study, a second-order spatial discretization scheme that has been tested with two-fluid two-phase flow model was extended to solve drift-flux model problems. In order to improve solver robustness, and therefore efficiency, a new approach was proposed to treating the mean drift velocity of the gas phase as a primary nonlinear variable to the equation system. With this new approach, significant improvement in solver robustness was achieved. With highly refined mesh, the proposed treatment along with the Newton-Krylov solver were extensively tested with two-phase flow problems that cover a wide range of thermal-hydraulics conditions. Satisfactory convergence performances were observed for all test cases. Numerical verification was then performed in the form of mesh convergence studies, from which expected orders of accuracy were obtained for both the first-order and the second-order spatial discretization schemes. Finally, the drift-flux model, along with numerical methods presented, were validated with three sets of flow boiling experiments that cover different flow channel geometries (round tube, rectangular tube, and rod bundle), and a wide range of test conditions (pressure, mass flux, wall heat flux, inlet subcooling and outlet void fraction).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Edge Equilibrium Code (EEC) For Tokamaks

The edge equilibrium code (EEC) described in this paper is developed for simulations of the near edge plasma using the finite element method. It solves the Grad-Shafranov equation in toroidal coordinate and uses adaptive grids aligned with magnetic field lines. Hermite finite elements are chosen for the numerical scheme. A fast Newton scheme which is the same as implemented in the equilibrium and stability code (ESC) is applied here to adjust the grids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: exploring the halo occupation distribution model for emission line galaxies

We study the modelling of the halo occupation distribution (HOD) for the eBOSS DR16 emission line galaxies (ELGs). Motivated by previous theoretical and observational studies, we consider different physical effects that can change how ELGs populate haloes. We explore the shape of the average HOD, the fraction of satellite galaxies, their probability distribution function (PDF), and their density and velocity profiles. Our baseline HOD shape was fitted to a semi-analytical model of galaxy formation and evolution, with a decaying occupation of central ELGs at high halo masses. We consider Poisson and sub/super-Poissonian PDFs for satellite assignment. We model both Navarro–Frenk–White and particle profiles for satellite positions, also allowing for decreased concentrations. We model velocities with the virial theorem and particle velocity distributions. Additionally, we introduce a velocity bias and a net infall velocity. We study how these choices impact the clustering statistics while keeping the number density and bias fixed to that from eBOSS ELGs. The projected correlation function, w p , captures most of the effects from the PDF and satellites profile. The quadrupole, ξ 2 , captures most of the effects coming from the velocity profile. We find that the impact of the mean HOD shape is subdominant relative to the rest of choices. Here, we fit the clustering of the eBOSS DR16 ELG data under different combinations of the above assumptions. The catalogues presented here have been analysed in companion papers, showing that eBOSS RSD+BAO measurements are insensitive to the details of galaxy physics considered here. These catalogues are made publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of a Wave-Powered, Reverse-Osmosis System and Its Economic Availability in the United States: Preprint

A wave energy converter (WEC) system has the potential to convert the wave energy resource directly into the high-pressure flow that is needed by the desalination system to permeate saltwater through the reverse-osmosis membrane to generate clean water. In this study, a wave-to-water numerical model was developed to investigate the potential use of a wave-powered desalination system (WPDS) for water production in the United States. The model was developed by coupling a time-domain radiation-and-diffraction-method-based numerical tool (WEC-Sim) for predicting the hydrodynamic performance of WECs with a solution-diffusion model that was used to simulate the reverse-osmosis process. To evaluate the feasibility of the WPDS, the wave-to-water numerical model was applied to simulate a desalination system that used an oscillating surge WEC device to pump seawater through the system. The annual water production was estimated based on the wave resource at a reference site on the coast of northern California to investigate the potential cost of water in that area, where the cost of water and electricity is high compared to other regions. In the scenario evaluated, for a 100-unit utility-scale electricity-producing array, the estimated levelized cost of energy for these WECs is about 3-6 times the U.S.'s current, unsubsidized electricity rates. However, with clean water as an end product and by directly producing pressurized water with WECs, rather than electricity as an intermediary, it is presently only 12 percent greater than typical water cost in California. This study suggests that a WEC array that produces water may be a viable, near-term solution to the nation's water supply, and the niche application of the WPDS may also provide developers with new opportunities to further develop technologies that benefit both the electric and drinking water markets.

16 TIDAL AND WAVE POWER↗

A modal approach to the numerical calculation of primordial non-Gaussianities

We propose a new method to numerically calculate higher-order correlation functions of primordial fluctuations generated from any early-universe scenario. Our key-starting point is the realization that the tree-level In-In formalism is intrinsically separable. This enables us to use modal techniques to efficiently calculate and represent non-Gaussian shapes in a separable form well suited to data analysis. We prove the feasibility and the accuracy of our method by applying it to simple single-field inflationary models in which analytical results are available, and we perform non-trivial consistency checks like the verification of the single field consistency relation. We also point out that the iε prescription is automatically taken into account in our method, preventing the need for ad-hoc tricks to implement it numerically.

79 ASTRONOMY AND ASTROPHYSICS↗

PMU-Aided Voltage Security Assessment for a Wind Power Plant

Because wind power penetration levels in electric power systems are continuously increasing, voltage stability is a critical issue for maintaining power system security and operation. The traditional methods to analyze voltage stability can be classified into two categories: dynamic and steady-state. Dynamic analysis relies on time-domain simulations of faults at different locations; however, this method needs to exhaust faults at all locations to find the security region for voltage at a single bus. With the widely located phasor measurement units (PMUs), the Thevenin equivalent matrix can be calculated by the voltage and current information collected by the PMUs. This paper proposes a method based on a Thevenin equivalent matrix to identify system locations that will have the greatest impact on the voltage at the wind power plant's point of interconnection. The number of dynamic voltage stability analysis runs is greatly reduced by using the proposed method. The numerical results demonstrate the feasibility, effectiveness, and robustness of the proposed approach for voltage security assessment for a wind power plant.

17 WIND ENERGY↗

Discrete Hahn polynomials for numerical solution of two-dimensional variable-order fractional Rayleigh–Stokes problem

The main aim of this paper is to find the numerical solutions of 2D Rayleigh–Stokes problem with the variable-order fractional derivatives in the Riemann–Liouville sense. The presented method is based on collocation procedure in combination with the new operational matrix of the variable-order fractional derivatives, in the Caputo sense, for the discrete Hahn polynomials. The main advantage of the proposed method is obtaining a global approximation for spatial and temporal discretizations, and it reduced the problem to an algebraic system, which is easier to solve. Also, the profit of approximating a continuous function by Hahn polynomials is that for computing the coefficients of the expansion, we only have to compute a summation and the calculation of coefficients is exact. The error bound for the approximate solution is estimated. Finally, we evaluate results of the presented method with other numerical methods.

97 MATHEMATICS AND COMPUTING↗

Towards modelling AR Sco: generalized particle dynamics and strong radiation-reaction regimes

Numerical simulations of relativistic plasmas have become more feasible, popular, and crucial for various astrophysical sources with the availability of computational resources. The necessity for high-accuracy particle dynamics is especially highlighted in pulsar modelling due to the extreme associated electromagnetic fields and particle Lorentz factors. Including the radiation-reaction force in the particle dynamics adds even more complexity to the problem, but is crucial for such extreme astrophysical sources. We have also realized the need for such modelling concerning magnetic mirroring and particle injection models proposed for AR Sco, the first white dwarf pulsar. This paper demonstrates the benefits of using higher-order explicit numerical integrators with adaptive time-step methods to solve the full particle dynamics with radiation-reaction forces included. We show that for standard test scenarios, namely various combinations of uniform E- and B-fields and a static dipole B-field, the schemes we use are equivalent to and in extreme field cases outperform standard symplectic integrators in accuracy. We show that the higher-order schemes have massive computational time improvements due to the adaptive time-steps we implement, especially in non-uniform field scenarios and included radiation reaction where the particle gyro-radius rapidly changes. When balancing accuracy and computational time, we identified the adaptive Dormand–Prince eighth-order scheme to be ideal for our use cases. The schemes we use maintain accuracy and stability in describing the particle dynamics and we indicate how a charged particle enters radiation-reaction equilibrium and conforms to the analytical Aristotelian Electrodynamics expectations.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimization performance, fidelity, and cost: SIAM VQE

This dataset contains files storing results from classically-simulated quantum subroutines within a dynamical mean-field theory workflow, and jupyter notebooks processing the data in these files to generate plots. The files store: (1) Results from variational quantum eigensolver (VQE) simulations searching for optimal parameters allowing parametrized quantum circuits to prepare approximations to ground states of different Anderson impurity models (AIMs) (2) Results from simulations of a quantum Lanczos algorithm (QLA) estimating the Lanczos coefficients defining the continued-fraction representation of an (AIM) Green’s function Description: Any file named vqe_gs_results* stores approximations to the ground state and energy of a given AIM estimated using three different methods: (1) Numerical diagonalization (2) Ideal VQE simulation (3) VQE simulation with sampling noise For each VQE simulations metadata about the optimization (optimization results plus number of quantum circuits that would have been executed on real hardware) is also stored. Any file named qla_dos_results* estimations for the Lanczos coefficients defining the Green’s function of an AIM. The stored estimations are achieved using different methods: (1) Numerical Lanczos algorithm from initial states obtained from numerical diagonalization (2) Simulated quantum Lanczos algorithm from initial states prepared from parametrized quantum circuits yielded by corresponding ideal and noisy VQE subroutines. The dataset is used and described in M. Karabin et al., "Quantum solver for single-impurity Anderson models with particle-hole symmetry", Phys. Rev. Research 8, 033066 (2026). DOI: https://doi.org/10.1103/7ys3-tl4l

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A semi-analytical method for simulating matrix diffusion in numerical transport models

We report a semi-analytical approximation for transient matrix diffusion is developed for use in numerical contaminant transport simulators. This method is an adaptation and extension of the heat conduction method of Vinsome and Westerveld (1980) used to simulate heat losses during thermally enhanced oil recovery. The semi-analytical method is used in place of discretization of the low permeability materials, and it represents the concentration profile in the low permeability materials with a fitting function that is adjusted in each element at each time-step. The resulting matrix diffusion fluxes are added to the numerical model as linear concentration-dependent source/sink terms. Since only the high permeability zones need to be discretized, the numerical formulation is extremely efficient compared to traditional approaches that require discretization of both the high and low permeability zones. The semi-analytical method compares favorably with the analytical solution for transient one-dimensional diffusion with first order decay, with a two-layer aquifer/aquitard solution, with the solution for transport in a fracture with matrix diffusion and decay, and with a fully numerical solution for transport in a thin sand zone bounded by clay with variable decay rates.

54 ENVIRONMENTAL SCIENCES↗

Numerical and experimental evaluation of mechanical performance of the multifunctional energy storage composites

This work presents numerical simulation methods to model the mechanical behavior of the multifunctional energy storage composites (MESCs), which consist of a stack of multiple thin battery layers reinforced with through-the-hole polymer rivets and embedded inside carbon fiber composite laminates. MESC has been demonstrated through earlier experiments on its exceptional behavior as a structural element as well as a battery. However, the inherent complex infrastructure of the MESC design has created significant challenges in simulation and modeling. A novel homogenization technique was adopted to characterize the multi-layer properties of battery material using physics-based constitutive equations combined with nonlinear deformation theories to handle the interface between the battery layers. Second, mechanical damage and failure modes among battery materials, polymer reinforcements, and carbon fiber-polymer interfaces were characterized through appropriate models and experiments. The model of MESCs has been implemented in a commercial finite element code in ABAQUS. A comparison of structural response and failure modes from numerical simulations and experimental tests are presented. The results of the study showed that the predictions of elastic and damage responses of MESCs at various loading conditions agreed well with the experimental data. © 2021

Materials Science↗