Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STRANGENESS”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Analysis of the shifted boundary method for the Stokes problem

The analysis of stability and accuracy of the shifted boundary method is developed for the Stokes flow equations. The key feature of the shifted boundary method, an embedded finite element method, is the shifting of the location where boundary conditions are applied from the true to a surrogate boundary, and an appropriate modification (shifting) of the value of the boundary conditions. An inf–sup condition is proved for the variational formulation associated to the shifted boundary method and we derive, by way of Strang’s second lemma, an optimal error estimate in the natural SBM norm. Here we also derive an L 2 -error estimate for the velocity field, by means of an extension of the Aubin–Nitsche approach to embedded, non-consistent, mixed finite element methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Time-discretization of a plasma-neutral MHD model with a semi-implicit leapfrog algorithm

The semi-implicit leapfrog time-discretization is a workhorse algorithm for initial-value MHD codes to bridge between vastly separated time scales. Inclusion of atomic interactions with neutrals breaks the functional structure of the MHD equations that exploited by the leapfrog. In this work, we address how to best integrate atomic physics into the semi-implicit leapfrog. Following the Crank-Nicolson method, one approach is to time-center the atomic interactions in the linear solver and use a Newton method to include the nonlinear contributions. Alternatively, another family of methods are based on operator-splitting the terms associated with the atomic interactions using a Strang-splitting technique. These methods naturally break equations into constituent ODE and PDE parts and preserve the structure exploited by the semi-implicit leapfrog. We study the accuracy and efficiency of these methods through a battery of 0D and 1D cases and show that a second-order-in-time Douglas-Rachford inspired coupling between the ODE and PDE advances is effective in reducing the time-discretization error to be comparable to that of Crank-Nicolson with Newton iteration of the nonlinear terms. Splitting ODE and PDE parts results in independent matrix solves for each field which reduces the computational cost considerably and provides parallelization over species relative to Crank-Nicolson.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗
DOE OSTI2020

Critique on “A Revisit of Navier-Stokes Equation”

In a recently published work (Sheng, 2020), Sheng identified inconsistencies in the concepts and principles leading to the Navier–Stokes equations and proposed a new asymmetric “friction tensor” to resolve these inconsistencies. We show that this friction tensor can be interpreted as the non-isotropic part of a vorticity-influenced viscous stress tensor following the analysis in Berdahl and Strang (1986). The asymmetric stress tensor however violates the principles of mechanical equilibrium in fluid flows that are free of external or body couples and is therefore non-physical in classical continuum hydrodynamics. We refute the arguments in favour of the asymmetric stress tensor in Sheng (2020) whilst also demonstrating that the friction tensor leads to the incorrect form of the Navier–Stokes equations in non-Cartesian co-ordinates and its ability to recover the correct governing equations in Cartesian co-ordinates is fortuitous.

asymmetric↗
DOE OSTI2025

NdPO 4 solubility and aqueous Neodymium speciation in supercritical fluids: An experimental study at 500–700 °C and 1.7 kbar

A key aspect in the formation of rare earth elements (REE) deposits is the role of REE transport as aqueous REE complexes in supercritical hydrothermal solutions, where the nature of the aqueous complex is controlled by solution composition, temperature and pressure. Despite chloride being considered as one of the most abundant transporting ligands in magmatic-hydrothermal fluids, experimental investigations on the stability of aqueous REE chloride complexes are scarce above 300 °C. In this study, synthetic NdPO 4 crystals were reacted with non-saline and saline (0, 0.05 and 0.5 mNaCl), acidic (0.01 mHCl) aqueous solutions in a series of solubility experiments conducted at 500–700 °C and 1.7 kbar, where the solubilities were determined using a stable Nd isotope ( 145 Nd isotope spike) dilution technique. NdPO 4 solubility ranges between 28 ppm and 10,858 ppm, where solubility increases with both temperature and salinity. At 500 °C, log mNdPO 4 increases from –3.93 to –1.60 and there is a strong correlation between NdPO 4 solubility and NaCl concentrations (slope of 1.2 ± 0.3), indicating stabilization of the Nd chloride aqueous complexes with a stoichiometry corresponding to NdCl 2+ . At 600 °C, this correlation is weaker (slope of 0.4, log mNdPO 4 increases from –2.63 to –1.88) indicating the stabilization of both Nd chloride and hydroxyl species controlling solubility. At 700 °C, NdPO 4 solubility is largely independent of NaCl concentration indicating that solubility is controlled by Nd hydroxyl complexes, where stoichiometry suggests the neutral Nd(OH) 3 0 species is dominant. The solubility product (Ksp) of NdPO4 is derived from experimental data with the relation: log K sp = -41.81 – 0.057T – 20987/T, with T temperature in Kelvin. Comparison of the measured Nd phosphate solubility to thermodynamic predictions using the available Helgeson-Kirkham-Flowers equation of state parameters for aqueous Nd complexes indicate that predictions are up to three orders of magnitude lower compared to experimental observations. This discrepancy is most pronounced in saline solutions, suggesting that thermodynamic properties of the REE chloride species in supercritical fluids require revision. Numerical simulations of fluid-rock interaction between acidic, saline fluids and a Strange Lake felsic mineral assemblage demonstrates that NdPO 4 solubility predictions from models are four to six orders of magnitude lower than those calculated based on empirical fits from experiments, which suggests that acidic, saline fluids may play an important role in mobilizing large amounts of light REE from 450 to 700 °C.

58 GEOSCIENCES↗
DOE OSTI2011

A fast numerical scheme for causal relativistic hydrodynamics with dissipation

Highlights: {yields} We have developed a new multi-dimensional numerical scheme for causal relativistic hydrodynamics with dissipation. {yields} Our new scheme can calculate the evolution of dissipative relativistic hydrodynamics faster and more effectively than existing schemes. {yields} Since we use the Riemann solver for solving the advection steps, our method can capture shocks very accurately. - Abstract: In this paper, we develop a stable and fast numerical scheme for relativistic dissipative hydrodynamics based on Israel-Stewart theory. Israel-Stewart theory is a stable and causal description of dissipation in relativistic hydrodynamics although it includes relaxation process with the timescale for collision of constituent particles, which introduces stiff equations and makes practical numerical calculation difficult. In our new scheme, we use Strang's splitting method, and use the piecewise exact solutions for solving the extremely short timescale problem. In addition, since we split the calculations into inviscid step and dissipative step, Riemann solver can be used for obtaining numerical flux for the inviscid step. The use of Riemann solver enables us to capture shocks very accurately. Simple numerical examples are shown. The present scheme can be applied to various high energy phenomena of astrophysics and nuclear physics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Fourier-spectral element approximation of the ion–electron Braginskii system with application to tokamak edge plasma in divertor configuration

Due to the extreme conditions required to produce energy by nuclear fusion in tokamaks, simulating the plasma behavior is an important but challenging task. We focus on the edge part of the plasma, where fluid approaches are probably the best suited, and our approach relies on the Braginskii ion–electron model. Assuming that the electric field is electrostatic, this yields a set of 10 strongly coupled and non-linear conservation equations that exhibit multiscale and anisotropy features. The computational domain is a torus of complex geometrical section, that corresponds to the divertor configuration, i.e. with an “X-point” in the magnetic surfaces. To capture the complex physics that is involved, high order methods are used: The time-discretization is based on a Strang splitting, that combines implicit and explicit high order Runge–Kutta schemes, and the space discretization makes use of the spectral element method in the poloidal plane together with Fourier expansions in the toroidal direction. The paper thoroughly describes the algorithms that have been developed, provides some numerical validations of the key algorithms and exhibits the results of preliminary numerical experiments. In particular, we point out that the highest frequency of the system is intermediate between the ion and electron cyclotron frequencies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Dissipation-preserving spectral element method for damped seismic wave equations

This article describes the extension of the conformal symplectic method to solve the damped acoustic wave equation and the elastic wave equations in the framework of the spectral element method. The conformal symplectic method is a variation of conventional symplectic methods to treat non-conservative time evolution problems, which has superior behaviors in long-time stability and dissipation preservation. To reveal the intrinsic dissipative properties of the model equations, we first reformulate the original systems in their equivalent conformal multi-symplectic structures and derive the corresponding conformal symplectic conservation laws. We thereafter separate each system into a conservative Hamiltonian system and a purely dissipative ordinary differential equation system. Based on the splitting methodology, we solve the two subsystems respectively. The dissipative one is cheaply solved by its analytic solution. While for the conservative system, we combine a fourth-order symplectic Nyström method in time and the spectral element method in space to cover the circumstances in realistic geological structures involving complex free-surface topography. The Strang composition method is adopted thereby to concatenate the corresponding two parts of solutions and generate the completed conformal symplectic method. A relative larger Courant number than that of the traditional Newmark scheme is found in the numerical experiments in conjunction with a spatial sampling of approximately 5 points per wavelength. A benchmark test for the damped acoustic wave equation validates the effectiveness of our proposed method in precisely capturing dissipation rate. The classical Lamb problem is used to demonstrate the ability of modeling Rayleigh wave in elastic wave propagation. More comprehensive numerical experiments are presented to investigate the long-time simulation, low dispersion and energy conservation properties of the conformal symplectic methods in both the attenuating homogeneous and heterogeneous media.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

A computational model for nanosecond pulse laser-plasma interactions

A multi-physics numerical model for laser-induced optical breakdown and laser-plasma interaction in a non-equilibrium gas is presented, accounting for: production of priming electrons via multi-photon ionization, energy absorption, cascade ionization, induced hydrodynamic response, and shock formation and propagation. The gas is governed by the Navier-Stokes equations, with non-equilibrium effects taken into account by means of a two-temperature model. The space-time dependence of the laser beam is modeled with a flux-tube formulation for the Radiative Transfer Equation. The flow governing equations are discretized in space using a second-order finite volume method. The semi-discrete equations are marched in time using an implicit-explicit (IMEX) dual time-stepping strategy, where diffusion and chemistry are solved implicitly, whereas convection is explicit. Application to a 20 ns long 50 mJ pulse laser-induced breakdown in quiescent O 2 shows the advantages of this temporal discretization during and just after the laser pulse, while a less-expensive symmetric Strang splitting (with implicit chemistry) is sufficient for the post-breakdown gas dynamics after ≃ 0.1 μs. The integrated model is shown to reproduce key features of corresponding experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Anomaly diffuse and dielectric relaxation in strontium doped lanthanum molybdate

Highlights: Black-Right-Pointing-Pointer The anomaly diffuse and dielectric relaxation behaviors are fitted by the Cole-Cole approach. Black-Right-Pointing-Pointer The peak in the LSMO is corresponding to different oxygen ion diffusion process. Black-Right-Pointing-Pointer We first give better explanation about the strange conductivity change caused by doping. Black-Right-Pointing-Pointer The oxygen ion diffusion is due to a combination of the dipolar relaxation and the motion of ions. -- Abstract: The dielectric properties of the La{sub 2-x}Sr{sub x}Mo{sub 2}O{sub 9-{delta}} (x = 0-0.2) ceramics were investigated in the temperature range of 300-800 K. Dielectric measurement reveals that two dielectric anomalies, associated with the oxygen ion diffusion, exist in frequency spectrum with x = 0.5. The broad dielectric peaks in tan {delta}({omega}) can be well fitted by a modified Cole-Cole approach. When x = 0.1, only one dielectric relaxation peak is observed, corresponding to different oxygen ion diffusion processes, as distinct from the only relaxation peak in the pure La{sub 2}Mo{sub 2}O{sub 9}. The relaxation parameters {tau}{sub 0}, the dielectric relaxation strength {Delta}, and the activation energy E{sub a} were obtained. The result of this work shows that, the conductivity change caused by doping between the two phases is due to the combination of the dipolar effects and motion of ions.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Resonances in QCD

We report on the EMMI Rapid Reaction Task Force meeting 'Resonances in QCD', which took place at GSI October 12-14, 2015 (Fig.~1). A group of 26 people met to discuss the physics of resonances in QCD. The aim of the meeting was defined by the following three key questions; what is needed to understand the physics of resonances in QCD?; where does QCD lead us to expect resonances with exotic quantum numbers?; and what experimental efforts are required to arrive at a coherent picture? For light mesons and baryons only those with up, down and strange quark content were considered. For heavy-light and heavy-heavy meson systems, those with charm quarks were the focus.This document summarizes the discussions by the participants, which in turn led to the coherent conclusions we present here.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Balance functions of (un)identified hadrons in Pb–Pb, p–Pb, and pp collisions at the LHC

In ultrarelativistic heavy-ion collisions, correlations of particles with opposite quantum numbers provide insight into quark production mechanisms and time scales, collective motion, and re-scattering in the hadronic phase. The longitudinal and azimuthal widths of balance functions for charged kaons and charged pions are used to examine the two-wave quark production model proposed to explain quark-antiquark production within the QGP, which predicts a large increase in up and down quark pairs relative to strange quark pairs around the time of hadronization. Balance functions are also analyzed in small collision systems such as p–Pb and pp to study fragmentation effects and possible collective effects in high-multiplicity events. A comprehensive set of balance functions has been measured using the ALICE detector, including results for unidentified hadrons in Pb–Pb and p–Pb collisions at √ s NN = 5.02 TeV, for charged pions in Pb–Pb collisions at √ s NN = 2.76 TeV, p–Pb collisions at √ s NN = 5.02 TeV and pp collisions at √ s NN = 7 TeV, and for charged kaons in Pb–Pb collisions at √ s NN = 2.76 TeV. The first balance function yield results are also presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Higher order cumulants of net baryon-number distributions at non-zero μ

Using recent results on higher order cumulants of conserved charge fluctuations from lattice QCD, we construct mean, variance, skewness, kurtosis, hyper-skewness and hyper-kurtosis of net-baryon number distributions for small baryon chemical potentials μ B . For the strangeness neutral case ( μ S = 0) at fixed ratio of electric charge to baryon number density ( n Q n B = 0.4 ), which is appropriate for a comparison with heavy ion collisions, we present results for κ B σ 2 B , S B σ 3 B / M B , κ H B σ 4 B and S H B σ 5 B / M B on the crossover line for the chiral transition, T pc ( μ B ). Furthermore, continuum extrapolations for this pseudo-critical transition line have recently been reported by HotQCD up to baryon chemical potentials μ B ≃ 300 MeV [A. Bazavov, et al., Chiral crossover in QCD at zero and non-zero chemical potentials, Phys. Lett. B 795 (2019) 15–21. arXiv:1812.08235 .]. These cumulant ratios are of direct relevance for comparisons with corresponding ratios measured by STAR in the BES-I and II runs at beam energies s NN ≥ 20 GeV . In particular, we point out that recent high statistics results on skewness and kurtosis of net-baryon number distributions obtained by STAR at s NN = 54.4 GeV put strong constraints on freeze-out parameters and are consistent with predictions from thermal QCD.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Chemical freeze-out parameters of net-kaons in heavy-ion collisions

In this work, we study chemical freeze-out parameters for heavy-ion collisions by performing two different thermal analyses. We analyze results from thermal fits for particle yields, as well as, net-charge fluctuations in order to characterize the chemical freeze-out. The Hadron Resonance Gas (HRG) model is employed for both methods. By separating the light hadrons from the strange hadrons in thermal fits, we study the proposed flavor hierarchy. For the net-charge fluctuations, we calculate the mean-over-variance ratio of the net-kaon fluctuations in the HRG model at the five highest energies of the RHIC Beam Energy Scan (BES) for different particle data lists. We compare these results with recent experimental data from the STAR collaboration in order to extract sets of chemical freeze-out parameters for each list. We focused on particle lists which differ largely in the number of resonant states. By doing so, our analysis determines the effect of the amount of resonances included in the HRG model on the freeze-out conditions. Our findings have potential impact on various other models in the field of relativistic heavy ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2013

A quenched study of the Schroedinger functional with chirally rotated boundary conditions: applications

In a previous paper (González López, et al., 2013) [1], we have discussed the non-perturbative tuning of the chirally rotated Schrödinger functional (χSFχSF). This tuning is required to eliminate bulk O(a) cutoff effects in physical correlation functions. Using our tuning results obtained in González López et al. (2013) [1] we perform scaling and universality tests analyzing the residual O(a) cutoff effects of several step-scaling functions and we compute renormalization factors at the matching scale. As an example of possible application of the χSFχSF we compute the renormalized strange quark mass using large volume data obtained from Wilson twisted mass fermions at maximal twist.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Improving Dark Matter Searches by Measuring the Nucleon Axial Form Factor: Perspectives from MicroBooNE

The MicroBooNE neutrino experiment at Fermilab is constructing a liquid-argon time-projection chamber for the Booster Neutrino Beam to study neutrino oscillations and interactions with nucleons and nuclei, starting in 2014. We describe the experiment and focus on its unique abilities to measure cross sections at low values of $Q^2$. In particular, the neutral-current elastic scattering cross section is especially interesting, as it is sensitive to the contribution of the strange sea quark spin to the angular-momentum of the nucleon, $\Delta s$. Implications for dark-matter searches are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Two-nucleon higher partial-wave scattering from lattice QCD

We present a determination of nucleon-nucleon scattering phase shifts for ℓ ≥ 0. The S, P, D and F phase shifts for both the spin-triplet and spin-singlet channels are computed with lattice Quantum ChromoDynamics. For ℓ > 0, this is the first lattice QCD calculation using the Lüscher finite-volume formalism. This required the design and implementation of novel lattice methods involving displaced sources and momentum-space cubic sinks. To demonstrate the utility of our approach, the calculations were performed in the SU(3)-flavor limit where the light quark masses have been tuned to the physical strange quark mass, corresponding to m π = m K ≈ 800MeV. In this work, we have assumed that only the lowest partial waves contribute to each channel, ignoring the unphysical partial wave mixing that arises within the finite-volume formalism. This assumption is only valid for sufficiently low energies; we present evidence that it holds for our study using two different channels. Two spatial volumes of V ≈ (3.5fm) 3 and V ≈ (4.6fm) 3 were used. The finite-volume spectrum is extracted from the exponential falloff of the correlation functions. Said spectrum is mapped onto the infinite volume phase shifts using the generalization of the Lüscher formalism for two-nucleon systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Measurement of B s 0 meson production in pp and PbPb collisions at s NN = 5.02 TeV

The production cross sections of B$_s^0$ mesons and charge conjugates are measured in proton-proton (pp) and PbPb collisions via the exclusive decay channel B$_s^0$→J/ψ&straightphi;→μ + μ - K + K - at a center-of-mass energy of 5.02 TeV per nucleon pair and within the rapidity range |y|<2.4 using the CMS detector at the LHC. The pp measurement is performed as a function of transverse momentum ( p T ) of the B$_s^0$ mesons in the range of 7 to 50 GeV/ c and is compared to the predictions of perturbative QCD calculations. The B$_s^0$ production yield in PbPb collisions is measured in two p T intervals, 7 to 15 and 15 to 50 GeV/c, and compared to the yield in pp collisions in the same kinematic region. The nuclear modification factor ( R AA ) is found to be 1.5±0.6(stat)±0.5(syst) for 7–15 GeV/c, and 0.87±0.30(stat)±0.17(syst) for 15–50 GeV/c, respectively. Within current uncertainties, the B$_s^0$ results are consistent with models of strangeness enhancement, and suppression by parton energy loss, as observed for the B + mesons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗