Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Unity”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Continuously differentiable PIC shape functions for triangular meshes

In this study, a new class of continuously-differentiable shape functions is developed and applied to two-dimensional electrostatic PIC simulation on an unstructured simplex (triangle) mesh. It is shown that troublesome aliasing instabilities are avoided for cold plasma simulation in which the Debye length is as small as 0.01 cell sizes. These new shape functions satisfy all requirements for PIC particle shape. They are non-negative, have compact support, and partition unity. They are given explicitly by cubic expressions in the usual triangle logical (areal) coordinates. The shape functions are not finite elements because their structure depends on the topology of the mesh, in particular, the number of triangles neighboring each mesh vertex. Nevertheless, they may be useful as approximations to solution of other problems in which continuity of derivatives is required or desired.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Taylor states in stellarators: A fast high-order boundary integral solver

We present a boundary integral equation solver for computing Taylor relaxed states in non-axisymmetric solid and shell-like toroidal geometries. The computation of Taylor states in these geometries is a key element for the calculation of stepped pressure stellarator equilibria. The integral representation of the magnetic field in this work is based on the generalized Debye source formulation, and results in a well-conditioned second-kind boundary integral equation. The integral equation solver is based on a spectral discretization of the geometry and unknowns, and the computation of the associated weakly-singular integrals is performed with high-order quadrature based on a partition of unity. The resulting scheme for applying the integral operator is then coupled with an iterative solver and suitable preconditioners. Several numerical examples are provided to demonstrate the accuracy and efficiency of our method, and a direct comparison with the leading code in the field is reported.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Improved C 1 shape functions for simplex meshes

Using continuously-differentiable (C 1 ) interpolation methods previously developed for computer visualization and other applications, a new class of shape functions is constructed and applied to electrostatic (ES) Particle-In-Cell (PIC) simulations on simplex (triangle or tetrahedral) meshes. The resulting shape functions are non-negative and partition unity,and lead to a reduction of shot noise and aliasing associated with particle interpolation to and from the mesh. For an energy-conserving implementation as chosen, C 1 interpolation of the potential leads to a continuous (C 0 ) electric field, with associated improved fidelity of particle motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A consistent and conservative model and its scheme for N -phase- M -component incompressible flows

Here, we propose a consistent and conservative model for multiphase and multicomponent incompressible flows, where there can be arbitrary numbers of phases and components. Each phase has a background fluid called the pure phase, each pair of phases is immiscible, and components are dissolvable in some specific phases. The model is developed based on the multiphase Phase-Field model including the contact angle boundary condition, the diffuse domain approach, and the analyses on the proposed consistency conditions for multiphase and multicomponent flows. The model conserves the mass of individual pure phases, the amount of each component in its dissolvable region, and thus the mass of the fluid mixture, and the momentum of the flow. It ensures that no fictitious phases or components can be generated and that the summation of the volume fractions from the Phase-Field model is unity everywhere so that there is no local void or overfilling. It satisfies a physical energy law and it is Galilean invariant. A corresponding numerical scheme is developed for the proposed model, whose formal accuracy is 2nd-order in both time and space. It is shown to be consistent and conservative and its solution is demonstrated to preserve the Galilean invariance and energy law. Numerical tests indicate that the proposed model and scheme are effective and robust to study various challenging multiphase and multicomponent flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Meshless discretization of the discrete-ordinates transport equation with integration based on Voronoi cells

The time-dependent, gray, linear radiation transport equation is discretized using the meshless local Petrov-Galerkin method with reproducing kernels. The integration is performed using a Voronoi tessellation, which creates a partition of unity that only depends on the position and extent of the kernels. The resolution of the integration automatically follows the particles and requires no manual adjustment. The discretization includes streamline-upwind Petrov-Galerkin stabilization to prevent oscillations and improve numerical conditioning. The angular quadrature is selectively refineable to increase angular resolution in chosen directions. The time discretization is done using backward Euler. The transport solve for each direction and the solve for the scattering source are both done using Krylov iterative methods. The results indicate first-order convergence in time and second-order convergence in space for linear reproducing kernels.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Data-driven acceleration of thermal radiation transfer calculations with the dynamic mode decomposition and a sequential singular value decomposition

In this work, we present a method for accelerating discrete ordinates radiative transfer calculations for radiative transfer. Our method works with nonlinear positivity fixes, in contrast to most acceleration schemes. The method is based on the dynamic mode decomposition (DMD) and using a sequence of rank-one updates to compute the singular value decomposition needed for DMD. Using a sequential method allows us to automatically determine the number of solution vectors to include in the DMD acceleration. We present results for slab geometry discrete ordinates calculations with the standard temperature linearization. Compared with positive source iteration, our results demonstrate that our acceleration method reduces the number of transport sweeps required to solve the problem by a factor of about 3 on a standard diffusive Marshak wave problem, a factor of several thousand on a cooling problem where the effective scattering ratio approaches unity, and a factor of 20 improvement in a realistic, multimaterial radiating shock problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

FSISPH: An SPH formulation for impacts between dissimilar materials

Here, we present an SPH formulation with several new features designed to better model the fully-compressible interaction of dissimilar materials. We developed the new method to simulate the atmospheric entry and break-up of small celestial bodies in planetary atmospheres. The formulation uses a unity-based, density-energy discretization of the hydrodynamic conservation laws with linear-corrected kernel gradients. To account for variations in compressibility, we use an HLLC approximate Riemann solver to adjust the velocity gradient at material interfaces. To handle large transverse velocity discontinuities, we introduce a simple slip interface model that limits the artificial viscosity at material interfaces. Diffusion is optionally applied through the velocity gradient and this allows the density and specific thermal energy to evolve in a manner more consistent with the first law of thermodynamics in comparison to other more direct diffusion schemes. We also introduce a material-local second-order artificial conduction scheme used to smooth the specific thermal energy field. Material damage fits neatly under this framework by treating the damage front as a material interface. The method has been implemented as a solver, FSISPH, within the code, Spheral++, and is publicly available on github. We test our new solver on a number of classic shock, mixing, and multi-material problem. The components we outline can significantly improve accuracy of SPH for problems with sharp contact discontinuities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Analysis of the High-Pressure High-Temperature (HPHT) growth of single crystal diamond

A multi-scale, computational model is developed to describe the growth characteristics of single-crystal diamond in the High-Pressure, High-Temperature (HPHT) process. This model is the first to connect phase-change kinetics governing crystal growth to the continuum transport of carbon through the growth cell. Results show the importance of convective transport driven by buoyant flow in the metallic solvent, which increases the growth rate by nearly an order of magnitude over that obtained under diffusion alone. Parametric studies show how crystal growth may be kinetically-limited or transport-limited, depending on the value of the macroscopic kinetic coefficient. Estimating this kinetic coefficient from growth experiments yields a phase-change Damköhler number of unity, indicating a mixed regime where phase-change kinetics and transport are comparable and strongly coupled in this system. Mechanisms responsible for slowing growth as the crystal size increases are explained. Finally, supersaturation inhomogeneities along the facets of larger crystals are predicted, which may be relevant to solvent inclusion formation during growth.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Comparison of fiber-optic distributed temperature sensing and high-sensitivity sensor spatial surveying of stream temperature

Measuring surface water temperature spatial variability is needed to estimate the interaction between surface water and groundwater, evaluate fish habitat and thermal inertia, and to estimate streamflow frequency and duration. Fiber optic distributed temperature sensing (FO-DTS) has been used in rivers and lakes, providing high-resolution and sensitive temperature monitoring over large temporal and spatial scales. However, in streams with cobbly or bedrock-lined streambeds and variable bathymetry, use of FO-DTS to measure temperature close to the surface water and groundwater interface can be challenging if even feasible. FO-DTS can also be costly, involve difficult installations, and require an advanced understanding of the technology, calibration, and data processing. In this study, we compared FO-DTS stream temperature survey results to an alternative temperature survey method employing a towed transect of high-resolution temperature loggers spaced at 1-m and transported in the stream along the study reach, to measure the spatial distribution of stream-water temperature in East Fork Poplar Creek near Oak Ridge, Tennessee, USA. We assessed the applicability and limitations of the two methods, and quantitatively compared in-situ temperature survey results measured simultaneously with each method. Regression results showed strong temporal and spatial correlation between the two methods. Differences were only elevated near the stream banks in areas that were coincident with correlation slope deviations from unity, which was attributed to shallower water and lower data density. Kriging standard errors were also low at channel center with minor increases near the stream banks. Furthermore, the results suggested that the array of the individual temperature sensors can provide a practical alternative to FO-DTS for thermal characterization of surface water, providing slightly lower spatial and temporal resolution, but with higher accuracy of temperature measurement, with greater simplicity, and with a broader range of conditions where it may be applied.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Multi-physics melt pool modeling and process optimization for laser direct energy deposition of Nb-based refractory C103: Defect formation, geometric precision, and process mapping

Recent developments in additive manufacturing (AM) technology have reignited interest in the fabrication of the Nb-based refractory C103 alloy offering solutions to the challenges posed by traditional manufacturing methods. However, the limited numerical and experimental studies on laser direct energy deposition (DED) of C103 have hindered the understanding of the relationships between process parameters and build quality. This has made it challenging to consistently produce parts with the desired quality and microstructure suitable for critical applications. In this study, we focus on optimizing the laser DED process for C103 by employing a hybrid approach that combines experimental techniques and computational fluid dynamics (CFD). This approach facilitates the development of process maps for defect detection and geometric precision. To achieve this, multi-layer C103 samples were fabricated using laser DED under various process parameters, enabling the creation of a process map for defect detection. Additionally, a multi-physics, multiphase simulation framework was developed within a high-performance computing (HPC) environment to establish process maps for geometric precision. Using these process maps, printability windows were identified for achieving both the desired geometric accuracy and defect-free prints. It was observed that prints with a power-to-velocity (P/V) ratio close to unity resulted in defect-free outcomes. This study provides a foundation for reducing design lead time and rejected parts, ultimately optimizing the laser DED process for C103.

Defect formation and geometric precision↗
DOE OSTI2024

A dynamically bare metal interface enables reversible magnesium electrodeposition at 50 mAh cm −2

Understanding and facilitating pure magnesium nucleation/growth electrodeposition behavior with ultrahigh Coulombic efficiency is complicated by the phenomenon of solid electrolyte interphase (SEI) formation in state-of-the-art, halogen-free magnesium electrolytes. Defining the electrolyte properties necessary to achieve ideal electrodeposition/stripping (E/S) thus remains elusive. Here, in this work, we reveal for the first time, rapid magnesium electrodeposition behavior that forms densely aligned, micron-sized thin platelets by establishing a dynamic bare magnesium/electrolyte interface during high-rate net plating. This effectively “SEI-free” interface allows facile magnesium diffusion and migration in stripping with near-unity E/S efficiency under demanding conditions over long-term cycling. The intrinsic electrolyte stability of the salt/solvent at the molecular level is the key to forming such an interface. The efficacy of the dynamic bare interface and an electrodeposited, free-standing magnesium anode is demonstrated in a high-areal-capacity full cell. These findings provide new design principles and fundamental understanding of interfacial chemistry in multivalent metal batteries.

Mg full cell↗
DOE OSTI2023

Dense Ru single-atoms integrated with sulfoacids for cellulose valorization to isosorbide

Metal-acid bifunctional catalysts are the unity of two opposites (metal-acid repelling) for cellulosic biomass valorization to isosorbide. It is challenging to improve the selectivity of metal hydrogenation catalysts and their synergism with acids for catalytic hydrolysis and dehydration to achieve substantial isosorbide. Herein, dense Ru single-atoms (10.1 wt% of Ru SAs) are anchored on sulfoacid-functionalized hollow mesoporous carbon shells, designed by assembling silica and 8-hydroxyquinoline-modified chitosan (HQ-CTS) through in situ Stöber templating strategy before pyrolysis and acid treatment. Based on X-ray absorption fine structure and computational modeling results, the structure of Ru SAs is determined as RuN 4 , which is more selective for a transitional glucose hydrogenation to sorbitol than Ru001 of Ru clusters. A lower-energy barrier of 1.21 (0.72) eV is delivered over RuN 4 (Ru001) for glucose hydrogenation (isomerization). These Ru SAs are integrated with sulfoacids (SO 3 H) but resistent against acids, rendering enhanced isosorbide yield in water as compared to Ru clusters, via a one-pot cascade reaction under harsh conditions (220 °C, 6 MPa H 2 ). The elaborately fabricated dense Ru SAs and sulfoacids, achieved by varying the addition time of HQ-CTS during the in situ Stöber templating process, improve the synergism of glucose hydrogenation with cellulose hydrolysis and sorbitol dehydration. In conclusion, this study provides a new idea for rational design of high-performance metal-acid bifunctional catalysts toward one-pot conversion of cellulose to isosorbide.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Superconducting nanowires as high-rate photon detectors in strong magnetic fields

Superconducting nanowire single photon detectors are capable of single-photon detection across a large spectral range, with near unity detection efficiency, picosecond timing jitter, and sub-10 µm position resolution at rates as high as 10 9 counts/s. In an effort to bring this technology into nuclear physics experiments, we fabricate Niobium Nitride nanowire detectors using novel material preparation methods and test their performance in strong magnetic fields. Furthermore, we demonstrate that these devices are capable of detection of 400 nm wavelength photons with saturated internal quantum efficiency at liquid helium temperatures and in magnetic fields potentially up to 8 T at high rates and with nearly zero dark counts.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Transverse-momentum and event-shape dependence of D-meson flow harmonics in Pb–Pb collisions at s NN = 5.02 TeV

The elliptic and triangular flow coefficients v 2 and v 3 of prompt D 0 , D + , and D * + mesons were measured at midrapidity ( | y | < 0.8 ) in Pb–Pb collisions at the centre-of-mass energy per nucleon pair of s NN = 5.02 TeV with the ALICE detector at the LHC. The D mesons were reconstructed via their hadronic decays in the transverse momentum interval 1 < p T < 36 GeV / c in central (0–10%) and semi-central (30–50%) collisions. Compared to pions, protons, and J / ψ mesons, the average D-meson v n harmonics are compatible within uncertainties with a mass hierarchy for p T ≲ 3 GeV / c , and are similar to those of charged pions for higher p T . The coupling of the charm quark to the light quarks in the underlying medium is further investigated with the application of the event-shape engineering (ESE) technique to the D-meson v 2 and p T -differential yields. The D-meson v 2 is correlated with average bulk elliptic flow in both central and semi-central collisions. Within the current precision, the ratios of per-event D-meson yields in the ESE-selected and unbiased samples are found to be compatible with unity. All the measurements are found to be reasonably well described by theoretical calculations including the effects of charm-quark transport and the recombination of charm quarks with light quarks in a hydrodynamically expanding medium.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Measurement of cold nuclear matter effects for inclusive J / ψ in p +Au collisions at s NN = 200 GeV

Measurement by the STAR experiment at RHIC of the cold nuclear matter (CNM) effects experienced by inclusive J / ψ at mid-rapidity in 0-100% p +Au collisions at s NN = 200 GeV is presented. Such effects are quantified utilizing the nuclear modification factor, R p Au , obtained by taking a ratio of J / ψ yield in p +Au collisions to that in p + p collisions scaled by the number of binary nucleon-nucleon collisions. The differential J / ψ yield in both p + p and p +Au collisions is measured through the dimuon decay channel, taking advantage of the trigger capability provided by the Muon Telescope Detector in the RHIC 2015 run. Consequently, the J / ψ R p Au is derived within the transverse momentum ( p T ) range of 0 to 10 GeV/ c . A suppression of approximately 30% is observed for p T < 2 GeV/ c , while J / ψ R p Au becomes compatible with unity for p T greater than 3 GeV/ c , indicating the J / ψ yield is minimally affected by the CNM effects at high p T . Comparison to a similar measurement from 0-20% central Au+Au collisions reveals that the observed strong J / ψ suppression above 3 GeV/ c is mostly due to the hot medium effects, providing strong evidence for the formation of the quark-gluon plasma in these collisions. Several model calculations show qualitative agreement with the measured J / ψ R p Au , while their agreement with the J / ψ yields in p + p and p +Au collisions is worse.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Elliptic flow of charged particles at midrapidity relative to the spectator plane in Pb–Pb and Xe–Xe collisions

Measurements of the elliptic flow coefficient relative to the collision plane defined by the spectator neutrons v 2 {$Ψ_{SP}$} in collisions of Pb ions at center-of-mass energy per nucleon–nucleon pair $\sqrt{s_{NN}}=$ 2.76 TeV and Xe ions at $\sqrt{s_{NN}}=$ 5.44 TeV are reported. The results are presented for charged particles produced at midrapidity as a function of centrality and transverse momentum for the 5–70% and 0.2–6 GeV/c ranges, respectively. The ratio between v 2 {$Ψ_{SP}$} and the elliptic flow coefficient relative to the participant plane v 2 {4}, estimated using four-particle correlations, deviates by up to 20% from unity depending on centrality. This observation differs strongly from the magnitude of the corresponding eccentricity ratios predicted by the T R ENTo and the elliptic power models of initial state fluctuations that are tuned to describe the participant plane anisotropies. The differences can be interpreted as a decorrelation of the neutron spectator plane and the reaction plane because of fragmentation of the remnants from the colliding nuclei, which points to an incompleteness of current models describing the initial state fluctuations. A significant transverse momentum dependence of the ratio v 2 {$Ψ_{SP}$} /v 2 {4} is observed in all but the most central collisions, which may help to understand whether momentum anisotropies at low and intermediate transverse momentum have a common origin in initial state fluctuations. The ratios of v 2 {$Ψ_{SP}$} and v 2 {4} to the corresponding initial state eccentricities for Xe–Xe and Pb–Pb collisions at similar initial entropy density show a difference of (7.0 ± 0.9)% with an additional variation of +1.8% when including RHIC data in the T R ENTo parameter extraction. These observations provide new experimental constraints for viscous effects in the hydrodynamic modeling of the expanding quark–gluon plasma produced in heavy-ion collisions at the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Measurement of beauty-strange meson production in Pb–Pb collisions at s NN = 5.02 TeV via non-prompt D s + mesons

The production yields of non-prompt ${D}^{+}_{s}$ mesons, namely ${D}^{+}_{s}$ mesons from beauty-hadron decays, were measured for the first time as a function of the transverse momentum ($p_T$) at midrapidity (|$y$| < 0.5) in central and semi-central Pb–Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s{NN}}$ = 5.02 TeV with the ALICE experiment at the LHC. The ${D}^{+}_{s}$ mesons and their charge conjugates were reconstructed from the hadronic decay channel ${D}^{+}_{s}$ → $φπ^+$, with $φ$ → K - K + , in the 4 < $p_T$ < 36 GeV/c and 2 < $p_T$ < 24 GeV/ $c$ intervals for the 0–10% and 30–50% centrality classes, respectively. The measured yields of non-prompt ${D}^{+}_{s}$ mesons are compared to those of prompt ${D}^{+}_{s}$ and non-prompt $D^0$ mesons by calculating the ratios of the production yields in Pb–Pb collisions and the nuclear modification factor $R_{AA}$. The ratio between the $R_{AA}$ of non-prompt ${D}^{+}_{s}$ prompt ${D}^{+}_{s}$ mesons, and that between the $R_{AA}$ of non-prompt ${D}^{+}_{s}$ and non-prompt $D^0$ mesons in central Pb–Pb collisions are found to be on average higher than unity in the 4 < $p_T$ < 12 GeV/ $c$ interval with a statistical significance of about 1.6σ and 1.7σ, respectively. The measured $R_{AA}$ ratios are compared with the predictions of theoretical models of heavy-quark transport in a hydrodynamically expanding QGP that incorporate hadronisation via quark recombination.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗