Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mesh Adaptivity”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Magnetized interstellar molecular clouds – I. Comparison between simulations and Zeeman observations

The most accurate measurements of magnetic fields in star-forming gas are based on the Zeeman observations analysed by Crutcher et al. Here, we show that their finding that the 3D magnetic field scales approximately as density 0.65 can also be obtained from analysis of the observed line-of-sight fields. We present two large-scale adaptive-mesh-refinement magnetohydrodynamic simulations of several thousand M ⊙ of turbulent, isothermal, self-gravitating gas, one with a strong initial magnetic field (Alfven Mach number $\mathcal{M}_{A,0}=1$) and one with a weak initial field ($\mathcal{M}_{A,0}=10$). We construct samples of the 100 most massive clumps in each simulation and show that they exhibit a power-law relation between field strength and density ($\bar{n}_H$) in excellent agreement with the observed one. Our results imply that the average field in molecular clumps in the interstellar medium (ISM) is $\langle B_{tot}(\bar{n}_H)\rangle ≈ 43\bar{n}^{0.65}_{H,4}μG$. Furthermore, the median value of the ratio of the line-of-sight field to density 0.65 in the simulations is within a factor of about (1.3, 1.7) of the observed value for the strong- and weak-field cases, respectively. The median value of the mass-to-flux ratio, normalized to the critical value, is 70 per cent of the line-of-sight value. This is larger than the 50 per cent usually cited for spherical clouds because the actual mass-to-flux ratio depends on the volume-weighted field, whereas the observed one depends on the mass-weighted field. Our results indicate that the typical molecular clump in the ISM is significantly supercritical (! factor of 3). The results of our strong-field model are in very good quantitative agreement with the observations of Li et al., which show a strong correlation in field orientation between small and large scales. Because there is a negligible correlation in the weak-field model, we conclude that molecular clouds form from strongly magnetized (although magnetically supercritical) gas, in agreement with the conclusion of Li et al.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

No double detonations but core carbon ignitions in high-resolution, grid-based simulations of binary white dwarf mergers

We study the violent phase of the merger of massive binary white dwarf systems. Our aim is to characterize the conditions for explosive burning to occur, and identify a possible explosion mechanism of Type Ia supernovae. The primary components of our model systems are carbon–oxygen (C/O) white dwarfs, while the secondaries are made either of C/O or of pure helium. We account for tidal effects in the initial conditions in a self-consistent way, and consider initially well-separated systems with slow inspiral rates. We study the merger evolution using an adaptive mesh refinement, reactive, Eulerian code in three dimensions, assuming symmetry across the orbital plane. We use a corotating reference frame to minimize the effects of numerical diffusion, and solve for self-gravity using a multigrid approach. We find a novel detonation mechanism in C/O mergers with massive primaries. Here, the detonation occurs in the primary's core and relies on the combined action of tidal heating, accretion heating, and self-heating due to nuclear burning. The exploding structure is compositionally stratified, with a reverse shock formed at the surface of the dense ejecta. The existence of such a shock has not been reported elsewhere. The explosion energy (1.6 × 10 51 erg) and 56 Ni mass (0.86 M ⊙ ) are consistent with an SN Ia at the bright end of the luminosity distribution, with an approximated decline rate of Δm 15 (B) ≈ 0.99. Our study does not support double-detonation scenarios in the case of a system with a 0.6 M ⊙ helium secondary and a 0.9 M ⊙ primary. Although the accreted helium detonates, it fails to ignite carbon at the base of the boundary layer or in the primary's core.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Giant clumps in simulated high- z Galaxies: properties, evolution and dependence on feedback

Here, we study the evolution and properties of giant clumps in high- z disc galaxies using adaptive mesh refinement cosmological simulations at redshifts z ∼ 6–1. Our sample consists of 34 galaxies, of halo masses 10 11 –10 12 M⊙ at z = 2, run with and without radiation pressure (RP) feedback from young stars. While RP has little effect on the sizes and global stability of discs, it reduces the amount of star-forming gas by a factor of ∼2, leading to a similar decrease in stellar mass by z ∼ 2. Both samples undergo extended periods of violent disc instability continuously forming giant clumps of masses 10 7 –10 9 M ⊙ at a similar rate, though RP significantly reduces the number of long-lived clumps (LLCs). When RP is (not) included, clumps with circular velocity ≲ 40 (20) km s -1 , baryonic surface density ≲ 200 (100)M ⊙ pc -2 and baryonic mass ≲ 10 8.2 (10 7.3 ) M ⊙ are short-lived, disrupted in a few free-fall times. More massive and dense clumps survive and migrate towards the disc centre over a few disc orbital times. In the RP simulations, the distribution of clump masses and star formation rates (SFRs) normalized to their host disc is similar at all redshifts, exhibiting a truncated power law with a slope slightly shallower than -2. The specific SFR (sSFR) of the LLCs declines with age as they migrate towards the disc centre, producing gradients in mass, stellar age, gas fraction, sSFR and metallicity that distinguish them from the short-lived clumps which tend to populate the outer disc. Ex situ mergers comprise ~37 percent of the mass in clumps and ~29 percent of the SFR. They are more massive and with older stellar ages than the in situ clumps, especially near the disc edge. Lastly, roughly half the galaxies at redshifts z = 4–1 are clumpy, with ~3–30 percent of their SFR and ~0.1–3 percent of their stellar mass in clumps.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Bondi–Hoyle accretion in a turbulent, magnetized medium

We introduce simulations of accretion on to point masses embedded in an isothermal gas that is magnetized and supersonically turbulent, as occurs for protostars in molecular clouds. We use the ORION2 adaptive mesh refinement (AMR) code to carry out ideal magnetohydrodynamic simulations for an rms Mach number $\mathcal{M}$ rms =5 and a wide range of Alfvén Mach numbers. We find that the probability density functions for accretion rates in all models are very wide (±0.5 dex) and asymmetric, tilted to low accretion rates; the mean accretion rate is about twice the median. We find that the results of Lee et al. for magnetized Bondi-Hoyle accretion with the relative velocity parallel to the field describe our results to within a factor of 2, and we suggest that this should be valid at least for $\mathcal{M}$ rms ≲ 10. Our results demonstrate that turbulent magnetic fields of the strength observed in molecular clouds reduce the accretion rate relative to the classical Bondi-Hoyle rate by a factor of a few for Alfvén Mach numbers of order unity, but this is comparable to the reduction due to supersonic hydrodynamic turbulence alone. This reduction in accretion rates should be taken into account in analytic models of competitive accretion and analytic estimates of the accretion luminosities of young stellar objects in molecular clouds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Relativistic resistive magnetohydrodynamic reconnection and plasmoid formation in merging flux tubes

Herein we apply the general relativistic resistive magnetohydrodynamics code BHAC to perform a 2D study of the formation and evolution of a reconnection layer in between two merging magnetic flux tubes in Minkowski space–time. Small-scale effects in the regime of low resistivity most relevant for dilute astrophysical plasmas are resolved with very high accuracy due to the extreme resolutions obtained with adaptive mesh refinement. Numerical convergence in the highly non-linear plasmoid-dominated regime is confirmed for a sweep of resolutions. We employ both uniform resistivity and non-uniform resistivity based on the local, instantaneous current density. For uniform resistivity we find Sweet–Parker reconnection, from $η$ = 10-2 down to $η$ = 10 -4 , for a reference case of magnetization $σ$ = 3.33 and plasma-$β$ = 0.1. For uniform resistivity $η$ = 5 × 10 -5 the tearing mode is recovered, resulting in the formation of secondary plasmoids. The plasmoid instability enhances the reconnection rate to $v_{rec}$ ~ 0.03$c$ compared to $v_{rec}$ ~ 0.01$c$ for $η$ = 10 -4 . For non-uniform resistivity with a base level $η_0$ = 10 -4 and an enhanced current-dependent resistivity in the current sheet, we find an increased reconnection rate of $v_{rec}$ ~ 0.1$c$. The influence of the magnetization $σ$ and the plasma-$β$ is analysed for cases with uniform resistivity $η$ = 5 × 10 -5 and $η$ = 10 -4 in a range 0.5 ≤ $σ$ ≤ 10 and 0.01 ≤ $β$ ≤ 1 in regimes that are applicable for black hole accretion discs and jets. The plasmoid instability is triggered for Lundquist numbers larger than a critical value of $S_c$ ≈ 8000.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Magnetized interstellar molecular clouds – II. The large-scale structure and dynamics of filamentary molecular clouds

Abstract We perform ideal magnetohydrodynamics high-resolution adaptive mesh refinement simulations with driven turbulence and self-gravity and find that long filamentary molecular clouds are formed at the converging locations of large-scale turbulence flows and the filaments are bounded by gravity. The magnetic field helps shape and reinforce the long filamentary structures. The main filamentary cloud has a length of ∼4.4 pc. Instead of a monolithic cylindrical structure, the main cloud is shown to be a collection of fibre/web-like substructures similar to filamentary clouds such as L1495. Unless the line-of-sight is close to the mean field direction, the large-scale magnetic field and striations in the simulation are found roughly perpendicular to the long axis of the main cloud, similar to L1495. This provides strong support for a large-scale moderately strong magnetic field surrounding L1495. We find that the projection effect from observations can lead to incorrect interpretations of the true three-dimensional physical shape, size, and velocity structure of the clouds. Helical magnetic field structures found around filamentary clouds that are interpreted from Zeeman observations can be explained by a simple bending of the magnetic field that pierces through the cloud. We demonstrate that two dark clouds form a T-shaped configuration that is strikingly similar to the infrared dark cloud SDC13, leading to the interpretation that SDC13 results from a collision of two long filamentary clouds. We show that a moderately strong magnetic field (${{\cal M}_{\rm A}}\sim 1$) is crucial for maintaining a long and slender filamentary cloud for a long period of time ∼0.5 Myr.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Langevin and Fokker-Planck analyses for diffusion-mediated passing of circular and discorectangular species in two-dimensional channels

The propensity for pairs of diffusing species to pass each other within narrow channels or pores is of basic interest as a first-passage-type problem. It is also of relevance for solution-phase transport in nanoporous materials, and in particular for catalytic conversion reactions where high yield requires that product species can efficiently pass reactant species to exit the pores. Here, we analyze a two-dimensional model with nonoverlapping circular and discorectangular species confined to a rectangular channel, and where passing is mediated by Brownian dynamics in an implicit solvent. For narrower channels where passing is still possible, the discorectangle must align with the channel to pass the circular species. Behavior of the passing propensity, 𝑃, can be assessed by strongly damped Langevin simulations, or within an equivalent Fokker-Planck equation (FPE) formalism. The latter corresponds to a diffusion problem in a “higher-dimensional channel” with a constriction. We assess the variation of the passing propensity, 𝑃, for a broad range of channel width including its scaling just above the threshold where passing is sterically blocked. Analysis of 𝑃 versus the rotational diffusion coefficient 𝐷 𝑟 of the discorectangle reveals a significant decrease in 𝑃 for lower 𝐷 𝑟 for moderate channel width. This prompts a direct analysis of the regime where 𝐷 𝑟 → 0 , for which the FPE can be reduced to a three-dimensional diffusion problem, precise analysis of which is facilitated by adaptive-mesh finite element methods. The dependence of 𝑃 on the aspect ratio of the discorectangle is also assessed.

Rahman, Md Khaledur [Ames Laboratory (AMES), Ames,↗
DOE OSTI2026

Seamlessly joining length scales: From atomistic thermal graphs to anisotropic continuum conductivity

Thermal transport in complex solids is governed by local structure, defects, and anisotropy, yet most continuum models still rely on oversimplified and homogenized conductivities. Here, we bridge atomistic and continuum descriptions by building finite element (FE) models directly from the site-projected thermal conductivity (SPTC), an atomic-level decomposition of the Green–Kubo thermal conductivity. We introduce a toolkit, the “Simulator Collection for Atomic-to-Continuum Scales (SCACS)”, which uses a graph neural network to predict SPTC on large atomic structures, coarse-grains these fields into anisotropic conductivity tensors, and embeds them into the heat-flow FE equation with a customized, anisotropy-aware adaptive mesh refinement scheme. Applied to silicon nanostructures, the resulting FE models act as representative volume elements, reproduce bulk conductivities, and capture interfacial and defect-driven anisotropy while maintaining thermodynamic consistency. Additionally, SCACS predicts experimental conductance trends and fields. This work demonstrates a general route for transferring atomistic transport information into device-scale thermal simulations with physics-based approximations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Effective Demagnetizing Factors of Diamagnetic Samples of Various Shapes

Effective demagnetizing factors that connect the sample magnetic moment with the applied magnetic field are calculated numerically for perfectly diamagnetic samples of various nonellipsoidal shapes. The procedure is based on calculating the total magnetic moment by integrating the magnetic induction obtained from a full three-dimensional (3D) solution of the Maxwell equations using an adaptive mesh. The results are relevant for superconductors (and conductors in ac fields) when the London penetration depth (or the skin depth) is much smaller than the sample size. In conclusion, simple but reasonably accurate approximate formulas are given for practical shapes including rectangular cuboids, finite cylinders in axial and transverse fields, as well as infinite rectangular and elliptical cross-section strips.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Real time time-dependent density functional theory using higher order finite-element methods

We present a computationally efficient approach to solve the time-dependent Kohn-Sham equations in real time using higher order finite-element spatial discretization, applicable to both pseudopotential and all-electron calculations. To this end, we develop an a priori mesh adaption technique, based on the semidiscrete (discrete in space but continuous in time) error estimate on the time-dependent Kohn-Sham orbitals, to construct an efficient finite-element discretization. Subsequently, we obtain the full-discrete error estimate to guide our choice of the time step. We employ spectral finite elements along with special reduced order quadrature to render the overlap matrix diagonal, thereby simplifying the inversion of the overlap matrix that features in the evaluation of the discrete time-evolution operator. We use the second-order Magnus operator as the time-evolution operator, wherein the action of the discrete Magnus operator, expressed as exponential of a matrix, on the Kohn-Sham orbitals is obtained efficiently through an adaptive Lanczos iteration. We observe close to optimal rates of convergence of the dipole moment with respect to spatial and temporal discretization, for both pseudopotential and all-electron calculations. We demonstrate a staggering 100-fold reduction in the computational time afforded by higher order finite elements over linear finite elements, for both pseudopotential and all-electron calculations. Further, for similar level of accuracy, we obtain significant computational savings by our approach as compared to state-of-the-art finite-difference methods. Here, we also demonstrate the competence of higher order finite elements for all-electron benchmark systems. Lastly, we observe good parallel scalability of the proposed method on many hundreds of processors.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Handing off the outcome of binary neutron star mergers for accurate and long-term postmerger simulations

In this work, we perform binary neutron star (BNS) merger simulations in full dynamical general relativity with IllinoisGRMHD, on a Cartesian grid with adaptive-mesh refinement. After the remnant black hole has become nearly stationary, the evolution of the surrounding accretion disk on Cartesian grids over long timescales (~1s) is suboptimal, as Cartesian coordinates over-resolve the angular coordinates at large distances, and the accreting plasma flows obliquely across coordinate lines dissipating angular momentum artificially from the disk. To address this, we present the HandOff, a set of computational tools that enables the transfer of general relativistic magnetohydrodynamic (GRMHD) and spacetime data from IllinoisGRMHD to harm3d, a GRMHD code that specializes in modeling black hole accretion disks in static spacetimes over long timescales, making use of general coordinate systems with spherical topology. We demonstrate that the handoff allows for a smooth and reliable transition of GRMHD fields and spacetime data, enabling us to efficiently and reliably evolve BNS dynamics well beyond merger. We also discuss future plans, which involve incorporating advanced equations of state and neutrino physics into BNS simulations using the HandOff approach.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Addition of tabulated equation of state and neutrino leakage support to illinoisgrmhd

Here we have added support for realistic, microphysical, finite-temperature equations of state (EOS) and neutrino physics via a leakage scheme to illinoisgrmhd, an open-source GRMHD code for dynamical spacetimes in the einstein toolkit. These new features are provided by two new, nrpy+-based codes: nrpyeos, which performs highly efficient EOS table lookups and interpolations, and nrpyleakage, which implements a new, adaptive mesh refinement (AMR)-capable neutrino leakage scheme in the einstein toolkit. We have performed a series of strenuous validation tests that demonstrate the robustness of these new codes, particularly on the Cartesian AMR grids provided by carpet. Furthermore, we show results from fully dynamical GRMHD simulations of single unmagnetized neutron stars, and magnetized binary neutron star mergers. This new version of illinoisgrmhd, as well as nrpyeos and nrpyleakage, is pedagogically documented in jupyter notebooks and fully open source. The codes will be proposed for inclusion in an upcoming version of the einstein toolkit.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Ameliorating the Courant-Friedrichs-Lewy condition in spherical coordinates: A double FFT filter method for general relativistic MHD in dynamical spacetimes

Numerical simulations of merging compact objects and their remnants form the theoretical foundation for gravitational wave and multimessenger astronomy. While Cartesian-coordinate-based adaptive mesh refinement is commonly used for simulations, spherical-like coordinates are more suitable for nearly spherical remnants and azimuthal flows due to lower numerical dissipation in the evolution of fluid angular momentum, as well as requiring fewer numbers of computational cells. However, the use of spherical coordinates to numerically solve hyperbolic partial differential equations can result in severe Courant-Friedrichs-Lewy (CFL) stability condition time step limitations, which can make simulations prohibitively expensive. This paper addresses this issue for the numerical solution of coupled spacetime and general relativistic magnetohydrodynamics evolutions by introducing a double fast Fourier transform (FFT) filter and implementing it within the fully message passing interface (mpi)-parallelized sphericalnr framework in the einstein toolkit. In conclusion, we demonstrate the effectiveness and robustness of the filtering algorithm by applying it to a number of challenging code tests, and show that it passes these tests effectively, demonstrating convergence while also increasing the time step significantly compared to unfiltered simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Signatures of primordial energy injection from axion strings

Axion strings are horizon-size topological defects that may be produced in the early Universe. Ultralight axion-like particles may form strings that persist to temperatures below that of big bang nucleosynthesis. Such strings have been considered previously as sources of gravitational waves and cosmic microwave background (CMB) polarization rotation. In this work we show, through analytic arguments and dedicated adaptive mesh refinement cosmological simulations, that axion strings deposit a subdominant fraction of their energy into high-energy Standard Model (SM) final states, for example, by the direct production of heavy radial modes that subsequently decay to SM particles. This high-energy SM radiation is absorbed by the primordial plasma, leading to novel signatures in precision big bang nucleosynthesis, the CMB power spectrum, and gamma-ray surveys. In particular, we show that CMB power spectrum data constrains axion strings with decay constants f a ≲ 10 12 GeV , up to model dependence on the ultraviolet completion, for axion masses m a ≲ 10 − 29 eV ; future CMB surveys could find striking evidence of axion strings with lower decay constants. Published by the American Physical Society 2024

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2010

Hydrodynamic study of plasma amplifiers for soft-x-ray lasers: A transition in hydrodynamic behavior for plasma columns with widths ranging from 20 {mu}m to 2 mm

Plasma-based seeded soft-x-ray lasers have the potential to generate high energy and highly coherent short pulse beams. Due to their high density, plasmas created by the interaction of an intense laser with a solid target should store the highest amount of energy density among all plasma amplifiers. Our previous numerical work with a two-dimensional (2D) adaptive mesh refinement hydrodynamic code demonstrated that careful tailoring of plasma shapes leads to a dramatic enhancement of both soft-x-ray laser output energy and pumping efficiency. Benchmarking of our 2D hydrodynamic code in previous experiments demonstrated a high level of confidence, allowing us to perform a full study with the aim of the way for 10-100 {mu}J seeded soft-x-ray lasers. In this paper, we describe in detail the mechanisms that drive the hydrodynamics of plasma columns. We observed transitions between narrow plasmas, where very strong bidimensional flow prevents them from storing energy, to large plasmas that store a high amount of energy. Millimeter-sized plasmas are outstanding amplifiers, but they have the limitation of transverse lasing. In this paper, we provide a preliminary solution to this problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

High-resolution Navier-Stokes simulations of Richtmyer-Meshkov instability with reshock

The interaction of a Mach-1.45 shock wave with a perturbed planar interface between sulfur hexafluoride and air is studied through high-resolution two-dimensional (2D) and three-dimensional (3D) shock-capturing adaptive mesh refinement simulations of multispecies Navier-Stokes equations. The sensitivities of time-dependent statistics on grid resolution for 2D and 3D simulations are evaluated to ensure the accuracy of the results. Here, the numerical results are used to examine the differences between the development of 2D and 3D Richtmyer-Meshkov instability during two different stages: (1) initial growth of hydrodynamic instability from multimode perturbations after first shock and (2) transition to chaotic or turbulent state after reshock. The effects of the Reynolds number on the mixing in 3D simulations are also studied through varying the transport coefficients.

42 ENGINEERING↗
DOE OSTI2022

On Using Linux Kernel Huge Pages with FLASH, an Astrophysical Simulation Code

We present efforts at improving the performance of FLASH, a multi-scale, multi-physics simulation code principally for astrophysical applications, by using huge pages on Ookami, an HPE Apollo 80 A64FX platform. FLASH is written principally in modern Fortran and makes use of the PARAMESH library to manage a block-structured adaptive mesh. We explored options for enabling the use of huge pages with several compilers, but we were only able to successfully use huge pages when compiling with the Fujitsu compiler. As a result, the use of huge pages substantially reduced the number of translation lookaside buffer misses, but overall performance gains were marginal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

An Evaluation of An Asynchronous Task Based Dataflow Approach For Uintah

The challenge of running complex physics code on the largest computers available has led to dataflow paradigms being explored. While such approaches are often applied at smaller scales, the challenge of extreme-scale data flow computing remains. The Uintah dataflow framework has consistently used dataflow computing at the largest scales on complex physics applications. At present Uintah contains two main dataflow models. Both are based upon asynchronous communication. One uses a static graph-based approach with asynchronous communication and the other uses a more dynamic approach that was introduced almost a decade ago. Subsequent changes within the Uintah runtime system combined with many more large scale experiments, has necessitated a reevaluation of these two approaches, comparing them in the context of large scale problems. While the static approach has worked well for some large-scale simulations, the dynamic approach is seen to offer performance improvements over the static case for a challenging fluid-structure interaction problem at large scale that involves fluid flow and a moving solid represented using particle method on an adaptive mesh.

Humphrey, Alan↗