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

On the convergence of statistics in simulations of stationary incompressible turbulent flows

When reporting statistics from simulations of statistically stationary chaotic phenomenon, it is important to verify that the simulations are time-converged. This condition is connected with the statistical error or number of digits with which statistics can be reliably reported. In this work we consider numerical experiments of low Reynolds number incompressible homogeneous and isotropic turbulence as a model problem to investigate statistical convergence over finite simulation times. Specifically, we investigate the time integration requirements that allow meaningful reporting of the statistical error associated with finiteness of the temporal domain. We address two key questions: (1) How long should a simulation be performed in terms of large eddy time, and (2) How should the simulation time be divided among temporal windows over which a quantity of interest is estimated so that its statistical error could be reliably reported? We find that reliable reporting of statistical errors requires simulations on the order of 10 4 large eddy times, which is orders of magnitude longer than typically performed. Additionally, data post-processing should employ windows of at least ten times the large eddy time scale, with the most robust computation of statistical error of the mean requiring window sizes of an additional factor of ten. For practical simulations, we demonstrate that it is possible to estimate the statistical error within a factor of two under a less stringent condition in which a minimum of four windows with size at least ten large eddy times are used. In conclusion, our observations for homogeneous isotropic turbulence are also shown to hold in turbulent channel flow.

42 ENGINEERING↗
DOE OSTI2021

Coupling Surface Flow with High-performance Subsurface Reactive Flow and Transport Code PFLOTRAN

Water exchange between the surface and subsurface is important for both water resource management and environmental protection. In this paper, we develop coupled surface and subsurface flow simulation capability in a parallel subsurface flow and reactive transport code PFLOTRAN. We sequentially couple the diffusion wave-based surface flow with the subsurface flow governedby the Richards equation in PFLOTRAN. These two flow domains are linked with a boundary condition switching method that ensures continuity of pressure and flux at the surface-subsurface interface. We verify the coupled code against other existing hydrologic models and observation data using a number of numerical experiments. The coupled hydrological model exhibits good performance in strong parallel scaling tests. The new coupled surface and subsurface simulator significantly advance community simulation capability towards improving integrated hydrologic and biogeochemical understanding of complex systems such as watersheds and river corridors. Keywords: Surface flow, Integrated hydrological modeling, Boundary condition switching, Parallel computing

Wu, Runjian↗
DOE OSTI2024

A high-order, localized-artificial-diffusivity method for Eulerian simulation of multi-material elastic-plastic deformation with strain hardening

A high-order method for Eulerian simulation of material undergoing large elastic–plastic deformation is developed. Thermodynamically consistent hyperelastic constitutive relations are assumed, facilitating the treatment of solids, liquids, and gases in a unified manner. Here, the method enables the simulation of multi-material interactions using a diffuse interface approach. Numerical capturing of material interfaces, shock waves, contact surfaces, and elastic-plastic strain discontinuities using high-order compact-difference schemes is assisted by Localized Artificial Diffusivity (LAD). In the new setting involving elastic–plastic deformation, the previously established terms for the artificial properties are verified to effectively regularize normal shocks. Additional LAD terms are introduced to the elastic and plastic kinematic equations to regularize shear shocks and other strain discontinuities, improving solution stability. Other important features of the method that improve robustness include the numerical treatment of compatibility terms in the kinematic equations, and the treatment of rotation. Particular emphasis is focused toward new advancements of the methods for plastic-deformation integration and the associated strain hardening of the material, including rate-dependent plasticity. The method is demonstrated on a variety of test problems, including 1-D impacts, a variant of the Shu-Osher problem, a Taylor impact, and a Richtmyer-Meshkov instability between two elastic–plastic solids with strain hardening.

42 ENGINEERING↗
DOE OSTI2024

Tensor Hypercontraction of Cluster Perturbation Theory: Quartic Scaling Perturbation Series for the Coupled Cluster Singles and Doubles Ground-State Energies

Even though Cluster Perturbation has been shown to be a robust non-iterative alternative to Coupled Cluster Theory, it is still plagued by high order polynomial computational scaling and the storage of higher order tensors. Here we present a proof-of-concept strategy for implementing Cluster Perturbation Theory ground state energy series for the coupled cluster singles and doubles energy with N 4 computational scaling using Tensor Hypercontraction (THC). The reduction in computational scaling by two orders is achieved by decomposing two electron repulsion integrals, doubles amplitudes and multipliers, as well as selected doubles intermediates to Tensor Hypercontraction format. Using the outlined strategy, we showcase that the Tensor Hypercontraction pilot implementations retain numerical accuracy to within 1 kcal/mol relative to corresponding conventional and density fitting implementations and we empirically verify the N 4 scaling.

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

Super-strong magnetic field-dominated ion beam dynamics in focusing plasma devices

High energy density physics is the field of physics dedicated to the study of matter and plasmas in extreme conditions of temperature, densities and pressures. It encompasses multiple disciplines such as material science, planetary science, laboratory and astrophysical plasma science. For the latter, high energy density states can be accompanied by extreme radiation environments and super-strong magnetic fields. The creation of high energy density states in the laboratory consists in concentrating/depositing large amounts of energy in a reduced mass, typically solid material sample or dense plasma, over a time shorter than the typical timescales of heat conduction and hydrodynamic expansion. Laser-generated, high current–density ion beams constitute an important tool for the creation of high energy density states in the laboratory. Focusing plasma devices, such as cone-targets are necessary in order to focus and direct these intense beams towards the heating sample or dense plasma, while protecting the proton generation foil from the harsh environments typical of an integrated high-power laser experiment. A full understanding of the ion beam dynamics in focusing devices is therefore necessary in order to properly design and interpret the numerous experiments in the field. In this work, we report a detailed investigation of large-scale, kilojoule-class laser-generated ion beam dynamics in focusing devices and we demonstrate that high-brilliance ion beams compress magnetic fields to amplitudes exceeding tens of kilo-Tesla, which in turn play a dominant role in the focusing process, resulting either in a worsening or enhancement of focusing capabilities depending on the target geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Linearised Fokker–Planck collision model for gyrokinetic simulations

We introduce a gyrokinetic, linearised Fokker–Planck collision model that satisfies conservation laws and is accurate at arbitrary collisionalities. The differential test-particle component of the operator is exact; the integral field-particle component is approximated using a spherical harmonic and a modified Laguerre polynomial expansion developed by Hirshman and Sigmar (1976 Phys. Fluids 19 1532). The numerical methods of the implementation in the δf-gyrokinetic code stella (Barnes et al 2019 J. Comput. Phys. 391 365–80) are discussed, and conservation properties of the operator are demonstrated. The collision model is then benchmarked against the collision model of the gyrokinetic solver GS2 in the limiting cases of a reduced test-particle collision operator and energy- and momentum-conserving operator. The accuracy of the full collision model is investigated by solving the parallel Spitzer-Härm problem for the transport coefficients. It is shown that retaining collisional energy flux and higher-order terms in the field-particle operator reduces errors in the transport coefficients from 10%–25% for a simple momentum- and energy-conserving model to under 1%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Inverse-designed metasurfaces with facile fabrication parameters

Abstract Optical metasurfaces are planar nanostructured devices that are industrially attractive in part because they utilize high-throughput microelectronic fabrication techniques for implementation. It is therefore critical to develop design paradigms that can balance the realization of highly efficient wavefront responses together with device manufacturability. We introduce a gradient-based design framework for freeform metasurfaces in which nanoscale elements are explicitly constrained to feature basic shapes, nearly uniform feature sizes, and exceptionally low aspect ratios. In spite of the apparent uniformity of the metasurface geometric features, the devices are able to utilize nonlocal near-field optical coupling to achieve highly efficient and extreme wavefront scattering beyond conventional design methodologies. Utilizing this approach, we design facile high-numerical-aperture devices such as beam deflectors and large-area metalenses capable of diffraction-limited focusing. We anticipate that these concepts can facilitate the design and integration of metasurfaces into monolithic optical systems.

Optics↗
DOE OSTI2024

Effective many-body interactions in reduced-dimensionality spaces through neural network models

Accurately describing properties of challenging problems in physical sciences often requires complex mathematical models that are unmanageable to tackle head on. Therefore, developing reduced-dimensionality representations that encapsulate complex correlation effects in many-body systems is crucial to advance the understanding of these complicated problems. However, a numerical evaluation of these predictive models can still be associated with a significant computational overhead. To address this challenge, in this paper we discuss a combined framework that integrates recent advances in the development of active-space representations of coupled cluster (CC) downfolded Hamiltonians with neural network approaches. The primary objective of this effort is to train neural networks to eliminate the computationally expensive steps required for evaluating hundreds or thousands of Hugenholtz diagrams, which correspond to multidimensional tensor contractions necessary for evaluating a many-body form of downfolded effective Hamiltonians. Using small molecular systems (the H 2 O and HF molecules) as examples, we demonstrate that training neural networks employing effective Hamiltonians for a few nuclear geometries of molecules can accurately interpolate or extrapolate their forms to other geometrical configurations characterized by different intensities of correlation effects. We also discuss differences between effective interactions that define CC downfolded Hamiltonians with those of bare Hamiltonians defined by Coulomb interactions in the active spaces. Published by the American Physical Society 2024

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Coordinated Inverter Control to Increase Dynamic PV Hosting Capacity: A Real-Time Optimal Power Flow Approach

High penetrations of distributed photovoltaics (PV) could cause adverse grid impacts, such as voltage violations. The recent development in inverter technologies provides the opportunity to develop control systems to realize effective PV governance and thus to improve dynamic PV hosting capacity for distribution grids. In this paper, a novel distributed energy resource management system (DERMS) solution is proposed by adopting the real-time optimal power flow (OPF) approach for coordinated control of the distributed PV inverters. The proposed approach eliminates the dependence on load knowledge via measurement feedback correction, and it can be implemented in real time. One challenge is that the technique is sensitive to the data availability and integrity of voltage measurements. Therefore, a decentralized DERMS approach is developed by leveraging the concepts of state estimation. The framework and the effectiveness of the solution approach are numerically demonstrated on a real distribution feeder in Southern California.

14 SOLAR ENERGY↗
DOE OSTI2024

On a High-Mirror Stellarator Reactor Exploratory Concept With Neutrons Concentrated on Centrifuge Liquids

In the framework of fusion energy research based on magnetic confinement, stellarators allow numerous degrees of freedom for the design of the magnetic trap and plasma shape. Taking advantage of these features, some plasma shapes might benefit several of the many integrated elements involved in commercial fusion reactors, like, e.g., decreasing the number, mass and complexity of replaced activated in-vessel components (IVC) (i.e., by using liquids), extraction of large power, tritium generation, and remote maintenance. Certainly, free-surface liquid materials were proposed for tokamaks and field-reversed configuration (FRC) to try to improve some of such elements, i.e., in advanced power extraction (APEX) studies. Some reactor-relevant quasi-isodynamic (QI) magnetic configurations exhibit a relatively straight sector of plasma and high magnetic mirror. The combination of those elements and possibilities in a single stellarator reactor concept might have some advantages, in spite of the uncertainties due to the current low technological readiness level (TRL). The proposed and studied reactor concept is based on a vacuum vessel having short curved sectors and longer wide cylindrical sectors, which encloses a high-mirror low-vertical-excursion magnetic configuration, and swirling liquids or rotating cylinders, which centrifuge molten Li salts located at the low field region. Thus, the molten salts (if possible covered by a thin layer of liquid lithium) would be located on the internal perimeter of the cylinder, to act as particle exhaust (except for helium), neutron power extractor, and tritium breeder. The high-mirror feature tries to concentrate the neutron power at the cylindrical sectors, which might avoid using breeding materials at the curved sectors. The different elements of the concept are exploratorily studied and defined, and the difficulties assessed.

centrifuge↗
DOE OSTI2021

Sensitivity Analysis of Tunable Equation of State Material Model In Pulsed Mercury Target Simulation

A pulsed neutron spallation target is subjected to very short but intense loads from repeated proton pulses. Approximately 60% of the energy from each proton pulse is deposited into the mercury target material and the stainless-steel target structure, leading to a high-pressure region in both the stationary target structure and the flowing mercury. The high-pressure region propagates and leads to fluid-structure interaction. The resultant loading on the target structure containing liquid mercury is difficult to predict, although various simulation approaches and material models for the mercury have been tried. To date, the best match of simulation to experimental data is obtained by using an equation of state (EOS) material model with a specified tensile cutoff pressure, which simulates the cavitation threshold. The inclusion of a threshold to represent cavitation is key to the successful predictions of stress waves triggered by the high-energy pulse striking the mercury and vessel. However, recent measurements of target structure strain show that significant discrepancies remain between the measured and simulated strain values in the EOS mercury model. These differences grow when noncondensable helium gas is intentionally injected into the flowing mercury to reduce the loading on the structure. An EOS-based proportional–integral–derivative (PID) mercury model has been proposed to reduce the gap between the measured and simulated vessel strain responses for targets with gas injection. The conceptual and numerical description and initial investigation of the PID model are presented in previous work. Further studies of this PID model — including the sensitivity of the structure’s strain response to model parameters (the tensile cutoff, PID parameters Kp, Ki, and Kd) — are reported in this article. Results show the strain response is more sensitive to changes in the tensile cutoff value than to changes in the model parameters Kp, Ki, and Kd. These results will aid in future work where the model parameters will be optimized to match simulation data to strain measurements.

Lin, Lianshan↗
DOE OSTI2024

Implementation of the Glued Sphere Discrete Element Model for Non-Spherical Particles in MFiX Software

To enhance solver capabilities, simulation flexibility and model validation within the MFiX software, the U.S. Department of Energy (DOE) is funding efforts to develop and integrate the glued-sphere discrete element method into the latest version of MFiX as a dedicated computational module. The glued-sphere discrete element method is a numerical technique to depict the behavior of non-spherical particles in granular flows or particulate systems by representing them as a collection of component spheres. These spheres are bonded together to approximate the shape and mechanical/chemical properties of a more complex particle. The method effectively reuses the existing sphere-sphere collision algorithm, interphase momentum and heat transfer calculations utilized in the traditional discrete element method, extending these capabilities to non-spherical particles. Additionally, this method explicitly resolves intra-particle temperature and species distributions. The MFiX glued-sphere computational module includes tools for generating glued sphere configurations, a dedicated solver, and visualization capabilities in post-processing. More specifically within the computational module, collision detection and calculations were first performed on component spheres and then mapped onto non-spherical particles. The linear spring-dashpot model was utilized to simulate the sphere-sphere interactions.

Ke, Renjie↗
DOE OSTI2025

Ab initio equation of state for vanadium

Accurate materials’ equations of state (EOS) are essential for understanding materials properties as well as for use in multiphysics simulations. In particular, hydrodynamics simulations are based on three fundamental conservation laws (mass, momentum and energy) that form an under-determined system of equations. The equation of state serves as an additional closure relation between thermodynamic variables for a given material that enables numerical hydrodynamics simulation. In this report, we focus on the development of an ab initio EOS for the body centered cubic (BCC) phase of Vanadium (V) for eventual integration into a multiphase EOS in the OpenSesame EOS database.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The Role of the Digital Thread for Security, Resilience, and Adaptability in Manufacturing

The manufacturing sector is dramatically evolving with recent technical and digital advances. Among these, the digital thread is revolutionizing manufacturing operations well beyond the historic downloading of programs to computer numerically controlled (CNC) machine tools and the uploading of edited programs from CNC controllers. The digital thread is “the communication framework that enables a connected data flow and integrated view of the asset’s data throughout its lifecycle across traditionally siloed functional perspectives” (Leiva 2016). It is pervasive around the world and has changed the way society operates; for example, map apps are used to provide time-optimal directions to destinations using real-time traffic feedback. The digital thread makes it possible to digitally verify products, ensure that the latest technologies are deployed across the entire manufacturing ecosystem, and strengthen the workforce by making each individual more efficient and effective. These three foundational advanced manufacturing concepts will ensure a next-generation secure, resilient, and adaptable manufacturing ecosystem and sustainably support current and future needs of society (Lynn et al. 2020).

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Jamming Detection for Low-Resolution SC-FDE Systems: A Machine Learning Approach

Jammers interfere with communication between base stations (BSs) and legitimate users, leading to degradation of wireless system performance. Our study focuses on jamming detection for wideband single-carrier frequency domain equalization (SC-FDE) systems with low-resolution analog-to digital converters (ADCs). In such systems, jamming detection is challenging because traditional analytical approaches cannot be directly applied due to the delay dispersion in wideband channels and the non-linearity induced by low-resolution ADCs. We propose a machine learning (ML)-based jamming detection method that directly uses the quantized receive signals. Significantly, our ML-based detector can be integrated into existing standard frameworks, such as unique word (UW)-based SC-FDE systems, as it uses existing pilots without requiring additional pilots for jamming detection. Through numerical simulations, we show that two or more bits provide satisfactory performance compared to unquantized scenarios. Additionally, we demonstrate that using more and well-separated pilot symbols improves performance.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2026

Restricted Partition Function Semiclassical Transition State Theory (RPF-SCTST): Applications to Reactions Involving Hydrogen Cyanide, Hydrogen Peroxide, and Formaldehyde Oxide

We demonstrate the efficiency of the numerical calculation of thermal semiclassical transition state theory (SCTST) rates across several representative chemical reactions using the restricted partition function (RPF)─viz a function that depends only on the imaginary action associated with a reaction. Here, we treat the potential energy surfaces (PESs) through fourth order expansions around the saddle point and integrate the Hamiltonian up to second order in employing vibrational perturbation theory. We apply this formalism to uni- and bimolecular reactions with different─though relatively high─barrier heights to uncover the influence of the barrier properties, minimal energies, and separability of the rotational motion on rate constants and tunneling corrections. Although all modes are coupled within the RPF, we found in our numerical examples that the rotational component in the absence of significant rotational distortions can be separated from the vibrational contributions. Moreover, a classical treatment of the rotational contribution is adequate over the temperature range from approximately 100 to 1000 K. We also found that the choice of DFT basis set can lead to variations in rate constants of up to an order of magnitude. Lastly, different schemes for counting eligible energy levels result in rate constants that differ by no more than a factor of 2, with the remaining discrepancies attributed to the treatment of near-convergent levels in perturbation theory.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Godiva IV Shielding Benchmark Overview [Slides]

This presentation shows that Integral experiments increase the level of realism from that of data measurements. Further, Godiva IV is capable of intensities to activate foils. IER 498 uses a relative measurement to provide an intermediate level of realism. With reproducibility and room return management are key components. The test matrix put forward in CED-2 enables numerous cross-comparisons. IER 498 facilitates various additional works using a room return shield to isolate the source. IER 557 demonstrated burst reproducibility sufficient for IER 498.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Experimentation on Finned-Tube Microchannel Heat Exchanger Incorporating Phase Change Material and R-410A

Integrating latent heat thermal energy storage, particularly phase change materials (PCM), in building heating, ventilation, air conditioning, and refrigeration (HVAC&R) systems, offers substantial benefits for peak load shifting and energy efficiency. While numerous studies explore PCM-embedded heat exchangers (PCM-HX) using single-phase heat transfer fluids, there needs to be more investigation regarding PCM-HX utilizing two-phase refrigerant flow. This gap is even more pronounced for microchannel heat exchanger (MCHX) configurations. In this study, an MCHX-based PCM-HX was built and tested to analyze and evaluate the performance of a PCM-HX operating with a two-phase refrigerant as the working fluid. A commercially available MCHX condenser unit was embedded with 2.6 kg of RT35 PCM and tested using R-410A as the working fluid, using an in-house test loop to control the refrigerant inlet conditions precisely. Thin and minimally invasive T-type thermocouples traced the temperature evolution along MCHX channels and tube banks. Additional larger T-type thermocouple probes at the refrigerant inlet and outlet measured the average power of the PCM-HX during condensing tests, ranging from 0.7 kW to 1.6 kW based on mass flow rates varying from 3.5 g/s to 8 g/s. This experimental study on MCHX-based PCM-HX with a two-phase refrigerant has the potential to offer practical design guidelines for the direct integration of PCM-TES in HVAC&R systems without the use of a secondary loop, while caution is advised to avoid a substantial increase in the refrigerant charge amount.

42 ENGINEERING↗