Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerics”

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

HPC-driven computational reproducibility in numerical relativity codes: a use case study with IllinoisGRMHD

Abstract Reproducibility of results is a cornerstone of the scientific method. Scientific computing encounters two challenges when aiming for this goal. Firstly, reproducibility should not depend on details of the runtime environment, such as the compiler version or computing environment, so results are verifiable by third-parties. Secondly, different versions of software code executed in the same runtime environment should produceconsistent numerical results for physical quantities. In this manuscript, we test the feasibility of reproducing scientific results obtained using theIllinoisGRMHDcode that is part of an open-source community software for simulation in relativistic astrophysics, theEinstein Toolkit. We verify that numerical results of simulating a single isolated neutron star withIllinoisGRMHDcan be reproduced, and compare them to results reported by the code authors in 2015. We use two different supercomputers: Expanse at SDSC, and Stampede2 at TACC. By compiling the source code archived along with the paper on both Expanse and Stampede2, we find thatIllinoisGRMHDreproduces results published in its announcement paper up to errors comparable to round-off level changes in initial data parameters. We also verify that a current version ofIllinoisGRMHDreproduces these results once we account for bug fixes which have occurred since the original publication.

Astronomy & Astrophysics↗
DOE OSTI2019

Comment on ‘Numerical modeling of tokamak breakdown phase driven by pure Ohmic heating under ideal conditions’

In this comment, we point out possible critical numerical flaws of recent particle simulation studies (Jiang et al 2016 Nucl. Fusion 56 126017, Peng et al 2018 Nucl. Fusion 58 026007) on the electrical gas breakdown in a simple one-dimensional periodic slab geometry. We show that their observations on the effects of the ambipolar electric fields during the breakdown, such as the sudden reversal of the ion flow direction, could not be real physical phenomena but resulting from numerical artifacts violating the momentum conservation law. We show that an incomplete implementation of the direct-implicit scheme can cause the artificial electric fields and plasma transports resulting in fallacies in simulation results. We also discuss that their simple plasma model without considering poloidal magnetic fields seriously mislead the physical mechanism of the electrical gas breakdown because it cannot reflect important dominant plasma dynamics in the poloidal plane (Yoo et al 2018 Nat. Commun. 9 3523).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Numerical quantum clock simulations for measuring tunneling times

Here, we numerically study two methods of measuring tunneling times using a quantum clock. In the conventional method using the Larmor clock, we show that the Larmor tunneling time can be shorter for higher tunneling barriers. In the second method, we study the probability of a spin-flip of a particle when it is transmitted through a potential barrier including a spatially rotating field interacting with its spin. According to the adiabatic theorem, the probability depends on the velocity of the particle inside the barrier. It is numerically observed that the probability increases for higher barriers, which is consistent with the result obtained by the Larmor clock. By comparing outcomes for different initial spin states, we suggest that one of the main causes of the apparent decrease in the tunneling time can be the filtering effect occurring at the end of the barrier.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Numerical simulation of RHIC polarized proton run 17 spin flipper experiments

RHIC nine-magnet spin flipper has been operated successfully during RHIC polarized proton Run 17, with 97% spin flip efficiency achieved. The results show the importance of mirror resonance removal, small spin tune spread, and proper spin flipper driving tune sweep speed. Detailed spin tracking simulations, based on a Lorentz force and Thomas-BMT differential equation numerical solver code for accuracy, have been carried out to understand the experimental results. Agreement within measurement accuracy is obtained at injection energy, 23.8 GeV. It is not as tight at 255 GeV, reasons for that are exposed. These measurements and numerical studies allow to determine the sensitivity of spin-flip efficiency to the dispersion slopes at the two Siberian snakes and to the ac dipole frequency sweep speed. They also provide guidance for future developments at BNL’s electron-ion collider. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Numerically exact generalized Green's function cluster expansions for electron-phonon problems

We generalize the family of approximate momentum average methods to formulate a numerically exact, convergent hierarchy of equations whose solution provides an efficient algorithm to compute the Green's function of a particle dressed by bosons suitable in the entire parameter regime. We use this approach to extract ground-state properties and spectral functions. Our approximation-free framework, dubbed the generalized Green's function cluster expansion (GGCE), allows access to exact numerical results in the extreme adiabatic limit, where many standard methods struggle or completely fail. We showcase the performance of the method, specializing three important models of charge-boson coupling in solids and molecular complexes: the molecular Holstein model, which describes coupling between charge density and local distortions, the Peierls model, which describes modulation of charge hopping due to intersite distortions, and a more complex Holstein + Peierls system with couplings to two different phonon modes, paradigmatic of charge-lattice interactions in organic crystals. Furthermore, the GGCE serves as an efficient approach that can be systematically extended to different physical scenarios, thus providing a tool to model the frequency dependence of dressed particles in realistic settings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Numerically modeling stochastic inflation in slow-roll and beyond

In this work, we present a complete numerical treatment of inflationary dynamics under the influence of stochastic corrections from sub-Hubble modes. We discuss how to exactly model the stochastic noise terms arising from the sub-Hubble quantum modes that give rise to the coarse-grained inflaton dynamics in the form of stochastic differential equations. The stochastic differential equations are solved event by event on a discrete time grid. We then compute the power spectrum of curvature perturbations that can be compared with the power spectrum computed in the traditional fashion using the Mukhanov-Sasaki equation by canonically quantizing the inflaton fluctuations. Our numerical procedure helps us to easily extend the formalism to ultra slow-roll inflation and study the possibility of primordial black hole formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Numerical relativity higher order gravitational waveforms of eccentric, spinning, nonprecessing binary black hole mergers

We use the open source, community-driven, numerical relativity software, the EINSTEIN TOOLKIT to study the physics of eccentric, spinning, nonprecessing binary black hole mergers with mass-ratios q = {2; 4; 6}, individual dimensionless spin parameters chi 1z = +/- 0.6, chi 2z = +/- 0.3, that include higher order gravitational wave modes l <= 4, except for memory modes. Assuming stellar mass binary black hole mergers that may be detectable by the advanced LIGO detectors, we find that including modes up to l = 4 increases the signal-to-noise of compact binaries between 3.5% to 35%, compared to signals that only include the l = ImI = 2 mode. We use two waveform models, TEOBResumS and SEOBNRE, which incorporate spin and eccentricity corrections in the waveform dynamics, to quantify the orbital eccentricity of our numerical relativity catalog in a gauge-invariant manner through fitting factor calculations. Our findings indicate that the inclusion of higher order wave modes has a measurable effect in the recovery of moderately and highly eccentric black hole mergers, and thus it is essential to develop waveform models and signal processing tools that accurately describe the physics of these astrophysical sources.

Joshi, Abhishek V.↗
DOE OSTI2024

Simulating a numerical UV completion of quartic Galileons

The Galileon theory is a prototypical effective field theory that incorporates the Vainshtein screening mechanism—a feature that arises in some extensions of general relativity, such as massive gravity. The Vainshtein effect requires that the theory contain higher order derivative interactions, which results in Galileons, and theories like them, failing to be technically well posed. While this is not a fundamental issue when the theory is correctly treated as an effective field theory, it nevertheless poses significant practical problems when numerically simulating this model. These problems can be tamed using a number of different approaches: introducing an active low-pass filter and/or constructing a UV completion at the level of the equations of motion, which controls the high momentum modes. These methods have been tested on cubic Galileon interactions, and have been shown to reproduce the correct low-energy behavior. Here we show how the numerical UV-completion method can be applied to quartic Galileon interactions, and present the first simulations of the quartic Galileon model using this technique. We demonstrate that our approach can probe physics in the regime of the effective field theory in which the quartic term dominates, while successfully reproducing the known results for cubic interactions. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2020

Fluid-driven fractures in granular media: Insights from numerical investigations

We investigate the mechanisms of opening-mode fracture initiation in granular media. The study is based on a simulation of grain-scale fluid-grain interactions through a coupled numerical approach in which the discrete element method is used to solve for the mechanics of a solid granular medium, and computational fluid dynamics is used to model fluid flow and drag forces. In this work, we present benchmark problems with analytical solutions and validate this numerical model against experiments on a viscous-drag-driven cavity in the literature. Additional simulation results show fracture initiation mechanisms in a random granular packing subjected to constant boundary stresses and to fluid injection with a localized source. The dimensionless variable $F_s/F_{sk}$ (ratio of seepage force $F_s$ and skeletal force $F_{sk}$) incorporates the impacts of physical properties and injection parameters including fluid viscosity, injection velocity, grain size, and effective stresses, and it has been used as a criterion separating regimes of fluid invasion and drag-driven fracture opening. Our simulation results show that $F_s/F_{sk}$ in combination with $τ_1$ (ratio of diffusion time from hydromechanical coupling and injection time) serves as a prediction of fracture opening within granular packing. We suggest a simple criterion ($F_s/F_{sk}$ > 1 or $τ_1$ > 0.17) that is valid for various types of granular media and injection conditions to determine if fracture opening will occur. Among other applications, this study is useful to predict the initiation and propagation of fractures in natural sediments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Direct numerical simulation and Reynolds-averaged Navier-Stokes modeling of the sudden viscous dissipation for multicomponent turbulence

Simulations of a turbulent multicomponent fluid mixture undergoing isotropic deformations are carried out to investigate the sudden viscous dissipation. This dissipative mechanism was originally demonstrated using simulations of an incompressible single-component fluid [S. Davidovits and N. J. Fisch, Phys. Rev. Lett. 116, 105004 (2016)]. By accounting for the convective and diffusive transfer of various species, the current work aims to increase the physical fidelity of previous simulations and their relevance to inertial confinement fusion applications. Direct numerical simulations of the compressed fluid show that the sudden viscous dissipation of turbulent kinetic energy is unchanged from the single-component scenario. More importantly, the simulations demonstrate that the mass fraction variance and covariance for the various species also exhibit a sudden viscous decay. Reynolds-averaged Navier-Stokes simulations were carried out using the $k - l$ model to assess its ability to reproduce the sudden viscous dissipation. Results show that the standard $k - l$ formulation does not capture the sudden decay of turbulent kinetic energy, mass-fraction variance, and mass-fraction covariance for simulations with various compression and expansion rates, or different exponents for the power-law model of viscosity. A new formulation of the $k - l$ model that is based on previous improvements to the $k - ε$ family of models is proposed, which leads to consistently good agreement with the direct numerical simulations for all the isotropic deformations under consideration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Grid resolution requirement for resolving rare and high intensity wall-shear stress events in direct numerical simulations

Turbulent signals are intermittent with large instantaneous fluctuations. Such large fluctuations lead to small Kolmogorov scales that are hard to resolve in numerical simulations [P. K. Yeung, K. R. Sreenivasan, and S. B. Pope, Effects of finite spatial and temporal resolution in direct numerical simulations of incompressible isotropic turbulence, Phys. Rev. Fluids 3, 064603 (2018)]. Here, the present paper follows the above basic logic, but instead of dissipation events in isotropic turbulence, we study wall-shear stress events in plane channel flow. Wall-shear stress fluctuations are increasingly more intermittent as the Reynolds number increases. Hence, one has to employ higher grid resolutions as the Reynolds number increases in order to resolve a given percentage of wall-shear stress events. The objective of this paper is to quantify effects of the grid resolutions on the rare and high intensity wall-shear stress events. We find that the standard grid resolution resolves about 99% of the wall-shear stress events at Re τ = 180. A slightly higher grid resolution has to be employed in order to resolve 99% of the wall-shear stress events at higher Reynolds numbers, and if the standard grid resolution is used for, e.g., a Re τ = 10000 channel flow, one resolves about 90%–95% wall-shear stress events.

42 ENGINEERING↗
DOE OSTI2021

Numerical calculation of the particle–fluid–particle stress in random arrays of fixed particles

Based on the nearest particle statistics [Zhang, J. Fluid Mech. 910, A16 (2021)], the phase interaction force in a multiphase flow is decomposed into a particle–mean-field force and the divergence of the particle–fluid–particle (PFP) stress. The PFP stress is proportional to the correlation product of the distance from a particle to its nearest neighbor and the force on the particle conditionally averaged on the nearest-neighbor location. In this work, a functional form of the stress is obtained corrected to the first order of the ratio between the interparticle distance and the macroscopic length scale of the flow. Particle-resolved numerical simulations are used to calculate the PFP stress in random arrays of fixed particles and to explore the physics represented by the stress. The numerical results show that the PFP stress is attractive along the direction of the flow and is repulsive in the directions perpendicular to the flow. In the flow regime simulated, this PFP stress can be considered as a macroscopic representation of the drafting–kissing–tumbling mechanism. Finally, the Reynolds stress for the fluid phase is also calculated and compared with the PFP stress.

42 ENGINEERING↗
DOE OSTI2020

A Digital Twin Approach to Study Additive Manufacturing Processing Using Embedded Optical Fiber Sensors and Numerical Modeling

One of the major challenges for metal-powder-based additive manufacturing is measuring and mitigating residual strain induced during the manufacturing processes. This article reports distributed fiber optic sensors embedded in Inconel alloy components as experimental means to validate numerical models of additive manufacturing process. Electroplating was used to deposit a metal protective jacket onto standard telecom single-mode fibers for strain measurements, Fiber sensors were embedded in an Inconel alloy substrate using the laser engineered net shaping (LENS) process. Here by using a Rayleigh-scattering optical frequency domain reflectometer (OFDR), temperature changes, and residual strain in the metal substrate were monitored with 5 mm spatial resolution during the LENS process. Using finite element analysis, temperature and strain profiles induced by the LENS deposition processes were also numerically studied. Discrepancies between the simulated temperature and strain profiles and those measured directly were less than 10%. Results presented in this article demonstrates a digital twin approach to fuse modeling results with distributed fiber sensor measurement data to study additive manufacturing process toward design and fabrication process optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Numerical Oscillation Prevention for PV Inverter Controllers in Quasi-Steady-State Simulators

Droop-based PV inverter control methods are popular for preventing overvoltage in low voltage networks with high installations of PV solar. Quasi-steady-state simulators (QSTS) can be used to reduce the computational burden of conducting long-term techno-economic analyses of such PV inverter controllers. If dynamic models for droop-based controllers are directly implemented in QSTS, numerical oscillation occurs and the controllers fail to converge to their steady-state operation. In this paper, a power network sensitivity-based algorithm is proposed for preventing numerical oscillation while implementing active and reactive droop-based PV inverter controllers in quasisteady state simulation tool (e.g. in GridLAB-D). Convergence of the sensitivity-based algorithm is verified using a typical radial distribution feeder. It is found that all three controllers reached converged state in less than 15 iterations.

GridLAB-D, PV inverter control, quasi-steady-state↗
DOE OSTI2020

The Numerical Computation of Clear-Sky Transmittance and Reflectance of Solar Radiation for the Modeling of PV System Performance

Modeling of photovoltaics (PV) system performance requires numerical computation of plane-of-array (POA) irradiances that are influenced by atmospheric absorption and scattering. Clear-sky transmittance and reflectance of diffuse irradiance are essential for the understanding of solar radiation related to multiple reflection between land surface and clouds. This study analytically solves the radiative transfer equation for clear-sky conditions and numerically computes the transmittance and reflectance of diffuse irradiance as functions of the optical thickness of the atmosphere. The computation will be implemented in all-sky radiative transfer models to improve their computational efficiency.

14 SOLAR ENERGY↗
DOE OSTI2020

Numerical Investigation of a Central Fuel Property Hypothesis Under Boosted Spark-Ignition Conditions

In the present work, a central fuel property hypothesis (CFPH), which states that fuel properties are sufficient to provide an indication of a fuel’s performance irrespective of its chemical composition, was numerically investigated. In particular, the objective of the study was to determine whether Research Octane Number (RON) and Motor Octane Number (MON), as fuel properties, are sufficient to describe a fuel’s knock-limited performance under boosted spark-ignition (SI) conditions within the framework of CFPH. To this end, four TPRF-bioblendstock surrogates having different compositions but matched RON (=98) and MON (=90), were first generated using a non-linear regression model based on artificial neural network (ANN). Additionally, three unconventional bioblendstocks were included in the analysis: di-isobutylene (DIB), isobutanol, and Anisole. Skeletal reaction mechanisms were generated for the TPRF-DIB, TPRF-isobutanol, and TPRF-anisole blends from a detailed kinetic mechanism. Thereafter, numerical simulations were performed for the fuel surrogates using the skeletal mechanisms and a virtual cooperative fuel research (CFR) engine model, under a representative boosted operating condition. In the computational fluid dynamics (CFD) model, the G-equation approach was employed to track the turbulent flame front and the well-stirred reactor model combined with the multi-zone binning strategy was used to capture auto-ignition in the end-gas. In addition, laminar flame speed (LFS) was tabulated for each blend as a function of pressure, temperature, and equivalence ratio a priori, and the lookup tables were used to prescribe laminar flame speed as an input to the G-equation model. Parametric spark timing sweeps were performed for each fuel blend to determine the corresponding knock-limited spark advance (KLSA) and 50% burn point (CA50) at the respective KLSA timing. It was observed that despite same RON, MON, and engine operating conditions, the TPRF-anisole blend exhibited markedly different knock-limited performance from the other three blends. This deviation from the octane index (OI) expectation was shown to be caused by differences in laminar flame speed. However, it was found that relatively large fuel-specific differences in LFS (>20%) would have to be present to cause any appreciable deviation from the OI framework. Otherwise, RON and MON would still be robust enough to predict a fuel’s knock-limited performance.

42 ENGINEERING↗
DOE OSTI2022

Experimental and Numerical Study on the Effect of Annular Combustor Design on Thermal Uniformity Jet in Crossflow

Abstract Non-uniformity of the exit flow temperature represents one of the significant damages to gas-turbine components, particularly turbine blades. This may occur in the course of gas-turbine operation. This paper aims to provide passive techniques by modifying the combustor design rather than changing the flow parameters to improve the thermal uniformity and turbine blades to reduce thermal stresses and increase turbine blades’ life span. An acceptable agreement between the numerical and experimental results has been achieved, and the agreement includes the velocity and temperature profile. Four different angles have been tested numerically and experimentally with a maximum error of 5% at two different Reynolds numbers. Designing the outer combustor surface with a 45-deg angle bend can give a more uniform temperature distribution of 37% higher than the basic design with only a 0.5% higher pressure drop.

Energy & Fuels↗
DOE OSTI2023

Direct Numerical Simulation of Transitional and Turbulent Flows Over Multi-Scale Surface Roughness—Part I: Methodology and Challenges

Abstract High-fidelity simulation of transitional and turbulent flows over multi-scale surface roughness presents several challenges. For instance, the complex and irregular geometrical nature of surface roughness makes it impractical to employ conforming structured grids, commonly adopted in large-scale numerical simulations due to their high computational efficiency. One possible solution to overcome this problem is offered by immersed boundary methods, which allow wall boundary conditions to be enforced on grids that do not conform to the geometry of the solid boundary. To this end, a three-dimensional, second-order accurate boundary data immersion method (BDIM) is adopted. A novel mapping algorithm that can be applied to general three-dimensional surfaces is presented, together with a newly developed data-capturing methodology to extract and analyze on-surface flow quantities of interest. A rigorous procedure to compute gradient quantities such as the wall shear stress and the heat flux on complex non-conforming geometries is also introduced. The new framework is validated by performing a direct numerical simulation (DNS) of fully developed turbulent channel flow over sinusoidal egg-carton roughness in a minimal-span domain. For this canonical case, the averaged streamwise velocity profiles are compared against results from the literature obtained with a body-fitted grid. General guidelines on the BDIM resolution requirements for multi-scale roughness simulation are given. Momentum and energy balance methods are used to validate the calculation of the overall skin friction and heat transfer at the wall. The BDIM is then employed to investigate the effect of irregular homogeneous surface roughness on the performance of an LS-89 high-pressure turbine blade at engine-relevant conditions using DNS. This is the first application of the BDIM to realize multi-scale roughness for transitional flow in transonic conditions in the context of high-pressure turbines. The methodology adopted to generate the desired roughness distribution and to apply it to the reference blade geometry is introduced. The results are compared to the case of an equivalent smooth blade.

Engineering↗