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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Optimized Renewable Energy Integration for EV High-Power Dynamic Wireless Charging Systems

Dynamic wireless charging for electric vehicles is an emerging technology to reduce on-board battery size and extend driving range. Due to its unique characteristic of vehicle-speed-related pulse-like load profile, the high-power dynamic wireless charging system (DWCS) introduces high stress to the utility grid. In this paper, an optimization model for renewable energy integration in the DWCS is proposed to mitigate the grid impact and minimize the operation costs of the whole system. As the load profile of DWCS is related to the traffic volume and various approaching vehicle speeds, the annual average daily traffic data and a stochastic model are used to develop 24-hour load profile of DWCS. To find a tradeoff between grid impact mitigation and operation costs minimization, relationships among power demand from power grid, photovoltaic (PV) capacity, wind energy (WE) capacity and energy storage (ES) capacity are analyzed, and the optimization objective and constraints are developed. Numerical simulation results demonstrate that energy storage integration can greatly mitigate the grid impact of DWCS, and optimal ratio of PV and WE can significantly reduce the operation cost of DWCS.

Zeng, Rong↗
DOE OSTI2025

Computation of generalised magnetic coordinates asymptotically close to the separatrix

Integrals to calculate generalised magnetic coordinates from an input magnetic flux function asymptotically close to the separatrix are presented, and implemented in the GPEC/DCON code suite. These integrals allow characterisation of the magnetic equilibrium of a diverted tokamak, in magnetic coordinates, arbitrarily close to the last closed flux surface, avoiding the numerical issues associated with calculating diverging field-line integrals near a magnetic x-point. Finally, these methods may assist ongoing efforts to develop robust asymptotic equilibrium behaviour for spectral 3D MHD codes at the separatrix.

equilibrium edge truncation↗
DOE OSTI2022

Multi-element flow-driven spectral chaos (ME-FSC) method for uncertainty quantification of dynamical systems

The flow-driven spectral chaos (FSC) is a recently developed method for tracking and quantifying uncertainties in the long-time response of stochastic dynamical systems using the spectral approach. The method uses a novel concept called enriched stochastic flow maps as a means to construct an evolving finite-dimensional random function space that is both accurate and computationally efficient in time. In this paper, we present a multi-element version of the FSC method (the ME-FSC method for short) to tackle (mainly) those dynamical systems that are inherently discontinuous over the probability space. In ME-FSC, the random domain is partitioned into several elements, and then the problem is solved separately on each random element using the FSC method. Subsequently, results are aggregated to compute the probability moments of interest using the law of total probability. To demonstrate the effectiveness of the ME-FSC method in dealing with discontinuities and long-time integration of stochastic dynamical systems, four representative numerical examples are presented in this paper, including the Van-der-Pol oscillator problem and the Kraichnan-Orszag three-mode problem. Results show that the ME-FSC method is capable of solving problems that have strong nonlinear dependencies over the probability space, both reliably and at low computational cost.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Independent Review of Groundwater Remediation Strategy for Hexavalent Chromium and RDX Groundwater Plumes at Los Alamos National Laboratory (Rev. 1)

Site operations at the Los Alamos National Laboratory (LANL) resulted in the release of oxidized chromium, Cr(VI), into Sandia Canyon from cooling tower effluent from 1956 until 1972. The chromium traveled with the surface water approximately 3 miles downstream before migrating below ground surface. Chromium concentrations exceed 50 μg/L in the upper portion of the aquifer. Another LANL groundwater plume of concern is associated with RDX (Royal Demolition Explosives, 1,3,5-trinitro-1,3,5-triazine). Between 1951 and 1996, RDX was released to the mesa-top facilities' process water outfall, adjacent and underlying soils, and alluvial sediments, along with surface water in Cañon de Valle. Between 2000 and 2010, two remedial actions were deployed, removing much of the near-surface RDX, however, recharge due to precipitation has transported RDX into the perched-intermediate zone and into the regional aquifer. The report documents an independent technical review by scientists from the Department of Energy (DOE) Network of National Laboratories for Environmental Management and Stewardship (NNLEMS) to provide recommendations for potential near term actions to address and optimize remediation for both the Cr(VI) and RDX plumes. The proposed near-term remedial actions include design of pump and treat systems for Cr(VI) and monitoring and study for natural attenuation for RDX. The review assesses existing data, conceptual and numerical modeling, and it recommends a technical integration process to support identifying and implementing strategic, effective and efficient remedies. The DOE Environmental Management Los Alamos Field Office (EM-LA) and their cleanup contractor Newport News Nuclear-BWTX, LLC Los Alamos (N3B) provided the information required for the review. Interviews were also conducted with regulators to obtain the full spectrum of technical, regulatory and scientific perspectives. The independent review team was impressed by the capabilities, experiences, innovativeness, and insightfulness of the technical representatives from both the regulator, the New Mexico Environment Department (NMED) and N3B. Incorporation of vadose zone flow pathways in the conceptual site model (CSM) and configuring the numerical modeling for the site was generally state-of-the-practice (or better). This could be considered state-of-the-art by addressing uncertainties related to spatial extent of hydraulic windows. The reviews from the regulators were thorough and often provided useful concepts for consideration and future/study resolutions. The overarching consensus recommendation of independent review team is that the LANL groundwater plumes should be addressed in context of the emerging "management of complex sites" paradigm.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Superdiffusion from Nonabelian Symmetries in Nearly Integrable Systems

The Heisenberg spin chain is a canonical integrable model. As such, it features stable ballistically propagating quasiparticles, but spin transport is subballistic at any nonzero temperature: An initially localized spin fluctuation spreads in time t to a width t 2/3 . This exponent as well as the functional form of the dynamical spin correlation function suggest that spin transport is in the Kardar–Parisi–Zhang (KPZ) universality class. However, the full counting statistics of magnetization is manifestly incompatible with KPZ scaling. A simple two-mode hydrodynamic description, derivable from microscopic principles, captures both the KPZ scaling of the correlation function and the coarse features of the full counting statistics, but remains to be numerically validated. These results generalize to any integrable spin chain invariant under a continuous nonabelian symmetry and are surprisingly robust against moderately strong integrability-breaking perturbations that respect the nonabelian symmetry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Microbial carbon limitation: The need for integrating microorganisms into our understanding of ecosystem carbon cycling

Abstract Numerous studies have demonstrated that fertilization with nutrients such as nitrogen, phosphorus, and potassium increases plant productivity in both natural and managed ecosystems, demonstrating that primary productivity is nutrient limited in most terrestrial ecosystems. In contrast, it has been demonstrated that heterotrophic microbial communities in soil are primarily limited by organic carbon or energy. While this concept of contrasting limitations, that is, microbial carbon and plant nutrient limitation, is based on strong evidence that we review in this paper, it is often ignored in discussions of ecosystem response to global environment changes. The plant‐centric perspective has equated plant nutrient limitations with those of whole ecosystems, thereby ignoring the important role of the heterotrophs responsible for soil decomposition in driving ecosystem carbon storage. To truly integrate carbon and nutrient cycles in ecosystem science, we must account for the fact that while plant productivity may be nutrient limited, the secondary productivity by heterotrophic communities is inherently carbon limited. Ecosystem carbon cycling integrates the independent physiological responses of its individual components, as well as tightly coupled exchanges between autotrophs and heterotrophs. To the extent that the interacting autotrophic and heterotrophic processes are controlled by organisms that are limited by nutrient versus carbon accessibility, respectively, we propose that ecosystems by definition cannot be ‘limited’ by nutrients or carbon alone. Here, we outline how models aimed at predicting non‐steady state ecosystem responses over time can benefit from dissecting ecosystems into the organismal components and their inherent limitations to better represent plant–microbe interactions in coupled carbon and nutrient models.

Soong, Jennifer L.↗
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↗
DOE OSTI2022

Code verification for practically singular equations

We report the method-of-moments implementation of the electric-field integral equation (EFIE) yields many code-verification challenges due to the various sources of numerical error and their possible interactions. Matters are further complicated by singular integrals, which arise from the presence of a Green's function. To address these singular integrals, an approach is presented in wherein both the solution and Green's function are manufactured. Because the arising equations are poorly conditioned, they are reformulated as a set of constraints for an optimization problem that selects the solution closest to the manufactured solution. In this paper, we demonstrate how, for such practically singular systems of equations, computing the truncation error by inserting the exact solution into the discretized equations cannot detect certain orders of coding errors. On the other hand, the discretization error from the optimal solution is a more sensitive metric that can detect orders less than those of the expected convergence rate.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Verification of Bison fission product species conservation under TRISO reactor conditions

When assessing the reliability and predictive capabilities of a simulation tool, code verification is used to ensure that the implemented numerical algorithm is a faithful representation of its underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. During this process, numerical results in a discrete solution are compared to the analytical solution of the mathematical model. Here, in this paper, the code verification process is applied to one-dimensional spatiotemporal problems that exercise partial differential equation governing the conservation of fission product species (or mass diffusion). Numerical experiments were performed in the Bison fuel performance code to evaluate its predictive capability under various TRISO reactor conditions such as base irradiation and safety heating test conditions for either short- or long-lived fission product species, as well as a case concerning evaporation from the outer surface of a particle. The code predictions were compared with the expected exact results obtained from the analytical expressions, and the fact that they demonstrate the correct analytical behavior provides strong evidence of proper numerical algorithm implementation.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2021

Impact of Integration Scheme on Performance of Anisotropic Plasticity Models

Given the prevalent role of metals in a variety of industries, schemes to integrate corresponding constitutive models in finite element applications have long been studied. A number of formulations have been developed to accomplish this task; each with their own advantages and costs. Often the focus has been on ensuring the accuracy and numerical stability of these algorithms to enable robust integration. While important, emphasis on these performance metrics may often come at the cost of computational expense potentially neglecting the needs of individual problems. In the current work, the performance of two of the most common integration methods for anisotropic plasticity -- the convex cutting plane (CCP) and closest point projection (CPP) -- across a variety of metrics is assessed; including accuracy and cost. A variety of problems are considered ranging from single elements to large representative simulations including both implicit quasistatic and explicit transient dynamic type responses. The relative performance of each scheme in the different instances is presented with an eye towards guidance on when the different algorithms may be beneficial.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Report on FY 2023 Research and Development on Specially Designed Creep-fatigue Experiments on Alloy 617 in Support of Improving Creep-fatigue Evaluation Approaches

Experimental and numerical studies in support of developing the integrated Elastic–Perfectly Plastic (EPP) plus Simplified Model Test (SMT) design methodology, referred to as the EPP+SMT method, continued in FY 2023. This report focuses on the methods for extrapolating the EPP+SMT creep-fatigue life curves at long hold times and low strain ranges at elevated temperatures. In this work, the available uniaxial creep-fatigue failure data on Alloy 617 at temperatures of 950°C and 850°C were analyzed to provide a guidance on the development of the extrapolation method and the creep-fatigue failure criteria. A viscoplastic constitutive model for Alloy 617 was adopted to extrapolate the mechanical responses to low strain ranges and long hold times. A set of design curves of Alloy 617 at temperatures of 950°C, 850°C, and 800°C with tensile hold times of 1 hr, 100 hr, and 1,000 hr are developed. Furthermore, creep-fatigue testing on two notch specimen geometries, shallow-notch and sharp V-notch, on Alloy 617 was performed 950°C to understand the multiaxial stress relaxation behavior. The experimental and numerical results on the notch specimens were compared with those on the standard uniaxial smooth bar specimens. The effect of multi-axial stress state combined with elastic follow-up on the stress relaxation behavior was investigated in this report.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A large deformation multiphase continuum mechanics model for shock loading of soft porous materials

A large deformation, coupled finite-element (FE) model is developed to simulate the multiphase response of soft porous materials subjected to high strain-rate loading. The approach is based on the theory of porous media (TPM) at large deformations. Simplifications to the one-dimensional regime studied in the numerical simulations follow. An overview of several different time integration schemes is presented for the purpose of solving the nonlinear dynamic coupled balance of momenta (mixture and fluid) and balance of mass of the mixture equations. Numerical examples are presented for (i) verification against closed-form analytical solutions assuming small loads, (ii) demonstrating large deformation effects at high strain-rate, and (iii) showing differences in deformations between a single-phase elastodynamics model with occluded compressible pore fluid and a multiphase poroelastodynamics model at high strain-rate. The multiphase model shows that the relative motion of the pore fluid significantly dampens the deformation response of the solid skeleton as compared to the single-phase model, and makes it possible to extract quantitative values for the stresses of the different constituents, thereby allowing one to form preliminary conclusions about the onset of damage in the solid skeleton. The novelty of the current work is developing a multiphase, large deformation, mixture theory numerical model for high strain-rate loading of soft porous materials. It was discovered that explicit, adaptive time-stepping Runge–Kutta schemes offer high accuracy at relatively low cost when compared to traditional implicit or explicit central difference time-stepping schemes for shock-like loadings. Here, shock viscosity is added to the mixture momentum balance equation to regularize the shock front, and a stabilization term is added to the mixture mass balance equation to stabilize equal order interpolation finite elements for the coupled finite element solution of multiphase materials.

Engineering↗
DOE OSTI2023

Combining machine-learned and empirical force fields with the parareal algorithm: application to the diffusion of atomistic defects

We numerically investigate an adaptive version of the parareal algorithm in the context of molecular dynamics. This adaptive variant has been originally introduced in [1]. We focus here on test cases of physical interest where the dynamics of the system is modelled by the Langevin equation and is simulated using the molecular dynamics software LAMMPS. In this work, the parareal algorithm uses a family of machine-learning spectral neighbor analysis potentials (SNAP) as fine, reference, potentials and embedded-atom method potentials (EAM) as coarse potentials. We consider a self-interstitial atom in a tungsten lattice and compute the average residence time of the system in metastable states. Our numerical results demonstrate significant computational gains using the adaptive parareal algorithm in comparison to a sequential integration of the Langevin dynamics. We also identify a large regime of numerical parameters for which statistical accuracy is reached without being a consequence of trajectorial accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Jacobian-free Newton–Krylov method for the simulation of non-thermal plasma discharges with high-order time integration and physics-based preconditioning

A preconditioning framework for the numerical simulation of non-thermal streamer discharges is developed using the Jacobian-free Newton-Krylov (JFNK) method. A reduced plasma fluid model is considered, consisting of electrons, one positive ion, one negative ion, and the electrostatic potential. Here, the plasma kinetics model includes ionization, electron-ion recombination, electron attachment, electron detachment, and ion-ion recombination. The governing equations are made dimensionless, discretized in space with finite differences, and integrated in time with a fully implicit method based on high-order backward differentiation formulas. The preconditioning framework is based on a linearized form of the governing equations and physics-based operator splitting. The efficiency of the preconditioning strategy is assessed through two test cases: streamer propagation between parallel plates and an axisymmetric pin-to-pin discharge. The fully implicit approach overcomes traditional restrictions in the time step size due to processes such as electron drift, electron diffusion, and dielectric relaxation. Excellent performance is observed through relevant statistics of the JFNK solver, although the number of linear iterations increases for the pin-to-pin discharge when nonlinear numerical boundary conditions are imposed at the electrodes. Performance studies show scalability with O(100-1000) processors for O(10M) unknowns with ample room for optimization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Beyond Limber: efficient computation of angular power spectra for galaxy clustering and weak lensing

Angular two-point statistics of large-scale structure observables are important cosmological probes. To reach the high accuracy required by the statistical precision of future surveys, some of these statistics may need to be computed without the commonly employed Limber approximation; the exact computation however requires integration over Bessel functions, and a brute-force evaluation is slow to converge. Here, we present a new method based on our generalized FFTLog algorithm for the efficient computation of angular power spectra beyond the Limber approximation. The new method significantly simplifies the calculation and improves the numerical speed and stability. It is easily extended to handle integrals involving derivatives of Bessel functions, making it equally applicable to numerically more challenging cases such as contributions from redshift-space distortions and Doppler effects. We implement our method for galaxy clustering and galaxy-galaxy lensing power spectra. We find that using the Limber approximation for galaxy clustering in future analyses like LSST Year 1 and DES Year 6 may cause significant biases in cosmological parameters, indicating that going beyond the Limber approximation is necessary for these analyses.

weak gravitational lensing↗
DOE OSTI2022

Novel Co‐Catalytic Activities of Solid and Liquid Phase Catalysts in High‐Rate Li‐Air Batteries

Abstract Li‐air batteries are considered strong candidates for the next‐generation energy storage systems designed for electrical transportation. However, low cyclability and current rates are two major drawbacks that hinder them from further realization. These issues necessitate the discovery of novel materials to significantly enhance the redox process of discharge products. In this study, a novel catalytic system comprised of tin sulfide (SnS) nanoflakes as a solid catalyst and tin iodide (SnI 2 ) as a dual‐functional electrolyte additive is discovered. This system enables operating the battery at high current rates up to 10 000 mA g −1 (corresponding to 1 mA cm −2 ). The SnS catalyst shows outstanding catalytic activity for both oxygen reduction and evolution reactions compared to carbon, noble metals, and other transition metal dichalcogenides. It also exhibits good structural integrity at high rates. The computations indicate numerous possible oxygen reduction sites without oxygen dissociations on the SnS surface through solution mechanism that is likely responsible for the formation of Li 2 O 2 . The calculations also indicate that the role of the SnI 2 is not only reacting with the lithium anode to provide protection but reducing the charge potential by promoting catalytic decomposition of the Li 2 O 2 . This work provides new novel additives for designing high‐rate Li‐air batteries.

25 ENERGY STORAGE↗
DOE OSTI2022

Toward control co-design of utility-scale wind turbines: Collective vs. individual blade pitch control

A large-eddy simulation framework has been coupled with controller modules to systematically investigate the impacts of collective (CPC) and individual (IPC) pitch control strategies on utility-scale wind turbine energy production and fatigue loads. Wind turbine components were parameterized using an actuator surface model to simulate the rotor blades and the turbine nacelle. The baseline CPC and IPC algorithms, consisting of single-input single-output proportional–integral controllers and two integral controllers, respectively, were incorporated into the numerical framework. A series of simulations were carried out to investigate the relative performance of the two controllers under various turbulent inflow conditions, spanning hub-height velocities of 7 to 14 m/s. The numerical simulation results of this study showed that, in comparison to the CPC, the IPC controller could successfully reduce the damage equivalent loads of utility-scale turbines at regions 2 and 3 of turbine operation by about 3% and 40%, respectively, without any penalty on the power production of the turbine. It was also shown that, despite its minor impact on the turbulence kinetic energy of the wake, the IPC controller did not influence the recovery of the turbine wake.

17 WIND ENERGY↗
DOE OSTI2021

On the computational solution of vector-density based continuum dislocation dynamics models: A comparison of two plastic distortion and stress update algorithms

Continuum dislocation dynamics models of mesoscale plasticity consist of dislocation transport-reaction equations coupled with crystal mechanics equations. The coupling between these two sets of equations is such that dislocation transport gives rise to the evolution of plastic distortion (strain), while the evolution of the latter fixes the stress from which the dislocation velocity field is found via a mobility law. Earlier solutions of these equations employed a staggered solution scheme for the two sets of equations in which the plastic distortion was updated via time integration of its rate, as found from Orowan's law. In this work, we show that such a direct time integration scheme can suffer from accumulation of numerical errors. We introduce an alternative scheme based on field dislocation mechanics that ensures consistency between the plastic distortion and the dislocation content in the crystal. The new scheme is based on calculating the compatible and incompatible parts of the plastic distortion separately, and the incompatible part is calculated from the current dislocation density field. Stress field and dislocation transport calculations were implemented within a finite element based discretization of the governing equations, with the crystal mechanics part solved by a conventional Galerkin method and the dislocation transport equations by the least squares method. A simple test was first performed to show the accuracy of the two schemes for updating the plastic distortion, which shows that the solution method based on field dislocation mechanics is more accurate. This method then was used to simulate an austenitic steel crystal under uniaxial loading and multiple slip conditions. By considering dislocation interactions caused by junctions, a hardening rate similar to discrete dislocation dynamics simulation results was obtained. Finally, the simulations show that dislocations exhibit some self-organized structures as the strain is increased.

36 MATERIALS SCIENCE↗