Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical”

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,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

NUMERICAL SIMULATIONS OF CORONAL HEATING THROUGH FOOTPOINT BRAIDING

Advanced three-dimensional (3D) radiative MHD simulations now reproduce many properties of the outer solar atmosphere. When including a domain from the convection zone into the corona, a hot chromosphere and corona are self-consistently maintained. Here we study two realistic models, with different simulated areas, magnetic field strength and topology, and numerical resolution. These are compared in order to characterize the heating in the 3D-MHD simulations which self-consistently maintains the structure of the atmosphere. We analyze the heating at both large and small scales and find that heating is episodic and highly structured in space, but occurs along loop-shaped structures, and moves along with the magnetic field. On large scales we find that the heating per particle is maximal near the transition region and that widely distributed opposite-polarity field in the photosphere leads to a greater heating scale height in the corona. On smaller scales, heating is concentrated in current sheets, the thicknesses of which are set by the numerical resolution. Some current sheets fragment in time, this process occurring more readily in the higher-resolution model leading to spatially highly intermittent heating. The large-scale heating structures are found to fade in less than about five minutes, while the smaller, local, heating shows timescales of the order of two minutes in one model and one minutes in the other, higher-resolution, model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

LIGHT BRIDGE IN A DEVELOPING ACTIVE REGION. II. NUMERICAL SIMULATION OF FLUX EMERGENCE AND LIGHT BRIDGE FORMATION

Light bridges, the bright structure dividing umbrae in sunspot regions, show various activity events. In Paper I, we reported on an analysis of multi-wavelength observations of a light bridge in a developing active region (AR) and concluded that the activity events are caused by magnetic reconnection driven by magnetconvective evolution. The aim of this second paper is to investigate the detailed magnetic and velocity structures and the formation mechanism of light bridges. For this purpose, we analyze numerical simulation data from a radiative magnetohydrodynamics model of an emerging AR. We find that a weakly magnetized plasma upflow in the near-surface layers of the convection zone is entrained between the emerging magnetic bundles that appear as pores at the solar surface. This convective upflow continuously transports horizontal fields to the surface layer and creates a light bridge structure. Due to the magnetic shear between the horizontal fields of the bridge and the vertical fields of the ambient pores, an elongated cusp-shaped current layer is formed above the bridge, which may be favorable for magnetic reconnection. The striking correspondence between the observational results of Paper I and the numerical results of this paper provides a consistent physical picture of light bridges. The dynamic activity phenomena occur as a natural result of the bridge formation and its convective nature, which has much in common with those of umbral dots and penumbral filaments.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2013

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↗
DOE OSTI2013

Numerical evaluation of the three-point scalar-tensor cross-correlations and the tensor bi-spectrum

Utilizing the Maldacena formalism and extending the earlier efforts to compute the scalar bi-spectrum, we construct a numerical procedure to evaluate the three-point scalar-tensor cross-correlations as well as the tensor bi-spectrum in single field inflationary models involving the canonical scalar field. We illustrate the accuracy of the adopted procedure by comparing the numerical results with the analytical results that can be obtained in the simpler cases of power law and slow roll inflation. We also carry out such a comparison in the case of the Starobinsky model described by a linear potential with a sudden change in the slope, which provides a non-trivial and interesting (but, nevertheless, analytically tractable) scenario involving a brief period of deviation from slow roll. We then utilize the code we have developed to evaluate the three-point correlation functions of interest (and the corresponding non-Gaussianity parameters that we introduce) for an arbitrary triangular configuration of the wavenumbers in three different classes of inflationary models which lead to features in the scalar power spectrum, as have been recently considered by the Planck team. We also discuss the contributions to the three-point functions during preheating in inflationary models with a quadratic minimum. We conclude with a summary of the main results we have obtained.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Numerical evaluation of the tensor bispectrum in two field inflation

We evaluate the dimensionless non-Gaussianity parameter h {sub NL}, that characterizes the amplitude of the tensor bispectrum, numerically for a class of two field inflationary models such as double inflation, hybrid inflation and aligned natural inflation. We compare the numerical results with the slow roll results which can be obtained analytically. In the context of double inflation, we also investigate the effects on h {sub NL} due to curved trajectories in the field space. We explicitly examine the validity of the consistency relation governing the tensor bispectrum in the squeezed limit. Lastly, we discuss the contribution to h {sub NL} due to the epoch of preheating in two field models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Numerical analysis of ELM stability with rotation and ion diamagnetic drift effects in JET

Stability to the type-I edge localized mode (ELM) in JET plasmas was investigated numerically by analyzing the stability to a peeling-ballooning mode with the effects of plasma rotation and ion diamagnetic drift. The numerical analysis was performed by solving the extended Frieman–Rotenberg equation with the MINERVA-DI code. To take into account these effects in the stability analysis self-consistently, the procedure of JET equilibrium reconstruction was updated to include the profiles of ion temperature and toroidal rotation, which are determined based on the measurement data in experiments. With the new procedure and MINERVA-DI, it was identified that the stability analysis including the rotation effect can explain the ELM trigger condition in JET with ITER like wall (JET-ILW), though the stability in JET with carbon wall (JET-C) is hardly affected by rotation. The key difference is that the rotation shear in JET-ILW plasmas analyzed in this study is larger than that in JET-C ones, the shear which enhances the dynamic pressure destabilizing a peeling-ballooning mode. In addition, the increase of the toroidal mode number of the unstable MHD mode determining the ELM trigger condition is also important when the plasma density is high in JET-ILW. Though such modes with high toroidal mode number are strongly stabilized by the ion diamagnetic drift effect, it was found that plasma rotation can sometimes overcome this stabilizing effect and destabilizes the peeling-ballooning modes in JET-ILW.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2018

Numerical discretization of completely integrable nonlinear Hamiltonian systems

Completely integrable nonlinear Hamiltonian systems offer attractive models for novel accelerator systems that are being developed to push the frontiers of high intensity and high power beam applications. However, all known numerical methods used to discretize and simulate these systems spoil the integrability property of the original, time-continuous, theoretically integrable system. Since the hallmark of integrability is the existence of sufficiently many functionally independent first integrals in involution, we investigated several methods’ performance in terms of invariant preservation over longer timescales. We conclude that, while novel structure preserving integrators would be ideal, symplectic methods under mild restrictions offer adequate invariant preservation. We especially highlight the very good performance of the Stormer-Verlet method at low accuracy and a modified Lobatto method at high accuracy for the Integrable Optics Test Accelerator (IOTA). However, we show that symplectic integrators are limited by numerically induced resonances, which make interpretation of some simulation results challenging. These limitations and workarounds are discussed.

43 PARTICLE ACCELERATORS↗
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 OSTI2014

Numerical exploration of spontaneous broken symmetries in multiorbital Hubbard models

Here, we study three proposals for broken symmetry in the cuprate pseudogap-oxygen antiferromagnetism, Theta(II) orbital loop currents, and circulating currents involving apex oxygens-through numerical exploration of multiorbital Hubbard models. Our numerically exact results show no evidence for the existence of oxygen antiferromagnetic order or the Theta(II) phase in the three-orbital Hubbard model. The model also fails to sustain an ordered current pattern even with the presence of additional apex oxygen orbitals. Thus, we conclude that it is difficult to stabilize the aforementioned phases in the multiorbital Hubbard models for parameters relevant to cuprate superconductors. However, the Theta(II) phase might be stabilized through explicit flux terms. We also found an enhanced propensity for circulating currents with such terms in calculations simulating applied stress or strain, which skew the copper-oxygen plane to resemble a kagome lattice. We propose an experimental viewpoint to shed additional light on this problem.

36 MATERIALS SCIENCE↗
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↗