Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “favorability”

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,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Impacts of point defects on shallow doping in cubic boron arsenide: A first principles study

Cubic boron arsenide (BAs) stands out as a promising material for advanced electronics, thanks to its exceptional thermal conductivity and ambipolar mobility. However, effective control of p- and n-type doping in BAs poses a significant challenge, mostly as a result of the influence of defects. In the present study, we employed density functional theory (DFT) to explore the impacts of the common point defects and impurities on p-type doping of Be B and Si As , and on n-type doping of Si B and Se As . We found that the most favorable point defects formed by C, O, and Si are C As , O B O As , Si As , C As Si B , and O B Si As , which have formation energies of less than 1.5 eV. While the O impurity detrimentally affects both p- and n-type dopings, C and Si impurities are harmful for n-type dopings, making n-type doping a potential challenge. Interestingly, the antisite defect pair A s B B As benefits both p- and n-type doping. Finally, the doping limitation analysis presented in this study can potentially pave the way for strategic development in the area of BAs-based electronics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Transferable predictions of energetic and structural properties for refractory solid solution alloys across chemical compositions

We present a data-efficient approach to train graph neural networks (GNNs) on density functional theory (DFT) data for accurate and transferable predictions of energetic and structural properties of refractory solid solution alloys in the niobium-tantalum-vanadium (Nb-Ta-V) chemical space. We start by training the GNN model only on DFT data that describes refractory binary alloys niobium-tantalum (Nb-Ta), niobium-vanadium (Nb-V), and tantalum-vanadium (Ta-V) to predict formation enthalpy and root mean squared displacement. Once trained, the GNN predictions are tested on DFT data describing refractory ternary alloys Nb-Ta-V. While, unsurprisingly, direct transferability from binary to ternary is not sufficiently accurate, augmenting the training with only 1% of the available ternary data (uniformly distributed across the entire range of chemical compositions) improves significantly the quality of the GNN predictions. For comparison, we assess the transferability in the opposite direction by training GNN models on ternary Nb-Ta-V data and making predictions on binaries Nb-Ta, Nb-V, and Ta-V, which exhibits notably higher predictive errors. The proposed methodology, which favors transferability from lower-component to higher-component alloys, offers an efficient path towards avoiding the curse of dimensionality incurred when collecting DFT data for discovery and design of multi-component disordered alloys.

Density functional theory calculations↗
DOE OSTI2025

Impact of amorphous pockets on displacement damage evolution in silicon

Silicon has long been known to exhibit amorphization in response to heavy particle bombardment. For doses below the total amorphization threshold, partial amorphization is observed in the form of scattered amorphous pockets. While extensive research has gone into modeling the formation and evolution of amorphous pockets in response to irradiation, no studies yet investigate their impact on the evolution of other damage such as interstitial supersaturation and clustering. In this study, we survey the impact of amorphous pockets on defect evolution in silicon when treated as static sinks. MD is first used to show that amorphous pockets provide energetically favorable sites for point defects relative to the crystalline bulk, supporting the hypothesis that they act as sinks. A 0-D cluster dynamics model is then constructed, taking an interstitial clustering model from the literature and including amorphous pockets as a sink species. We conduct our survey for temperatures between 30 and 400 °C and sink strengths between 1 to 6 x 10 10 cm −2 . Both implantation- and radiation-induced damage states are investigated using interstitial and vacancy concentrations as initial condition variables. We find that, due to the differing migration rates of the interstitial and the vacancy, amorphous pockets have a non-monotonic impact on the final damage state depending on the effective sink strength of the amorphous pockets, resulting in increased damage formation in regimes of intermediate amorphization. In conclusion, this result emphasizes the important role of amorphous pockets in governing the evolution of damage in partially amorphized crystalline materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A new interatomic potential for mixed Mg-Al-Ga-In spinels

While density functional theory (DFT) has become the de facto approach for accurate simulation of materials at the atomic scale, there are many aspects of materials that are simply out of reach of DFT methods. In particular, finite temperature properties such as diffusivities, the structure and properties of grain boundaries and interfaces, and the study of defect properties in complex alloys are computationally challenging for DFT methods. Recently, a new class of spinels in which three cations order over two sublattices was discovered. In order to predict the properties of these types of structures, classical potentials are a must. Here, in this work, we derive a new classical potential for Mg-bearing spinels in which the B cations are Al, Ga, and/or In. The potential does well in describing the DFT energetics of various spinel structures as a function of chemistry and inversion. In particular, it reproduces the thermodynamically favorable MgAlGaO 4 structure while correctly predicting that neither MgAlInO 4 nor MgGaInO 4 are stable. Further, it reproduces physical trends in elastic properties as compared against experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Effects of interstitial oxygen on ω transformations and twin formation in bcc NbTaTiHf multi-principal element alloy from first-principles

Transformation- and twinning-induced plasticity (TRIP and TWIP) have been reported to contribute to the low-temperature deformation of some body-centered cubic (bcc) multi-principal element alloys (MPEAs) containing large fractions of group IV transition metals. The influence of interstitial solutes on the mechanisms underlying these forms of plasticity, however, remains unclear. Using first-principles calculations, we study the effects of interstitial O atoms on the relative stability of bcc and ω phases and on unstable and twin boundary stacking fault energy profiles in a representative bcc MPEA with high group-IV elemental fraction: NbTaTiHf. We find that O additions generally promote the relaxation of ω configurations back to their parent bcc structure, therefore inhibiting ω transformation. Calculations of the Rice parameter for bulk bcc and phases, as well as bcc-ω interfaces, further show that ω formation is a potent embrittlement factor, an effect that is enhanced by O additions, suggesting that the formation of bcc-ω interfaces is energetically preferred over the formation of the bulk ω phase. By contrast, the Rice parameter for twin boundaries indicates that these interfaces do not embrittle the material, even with O atoms at twin boundaries, providing a more favorable pathway for plastic deformation compared to ω transformation.

Density functional theory↗
DOE OSTI2020

Efficient discontinuous Galerkin implementations and preconditioners for implicit unsteady compressible flow simulations

This work presents and compares efficient implementations of high-order discontinuous Galerkin methods: a modal matrix-free discontinuous Galerkin (DG) method, a hybridizable discontinuous Galerkin (HDG) method, and a primal formulation of HDG, applied to the implicit solution of unsteady compressible flows. The matrix-free implementation allows for a reduction of the memory footprint of the solver when dealing with implicit time-accurate discretizations. HDG reduces the number of globally-coupled degrees of freedom relative to DG, at high order, by statically condensing element-interior degrees of freedom from the system in favor of face unknowns. The primal formulation further reduces the element-interior degrees of freedom by eliminating the gradient as a separate unknown. This paper introduces a p-multigrid preconditioner implementation for these discretizations and presents results for various flow problems. Benefits of the p-multigrid strategy relative to simpler, less expensive, preconditioners are observed for stiff systems, such as those arising from low-Mach number flows at high-order approximation. The p-multigrid preconditioner also shows excellent scalability for parallel computations. Finally, additional savings in both speed and memory occur with a matrix-free/reduced version of the preconditioner.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Sensitivity of puff dynamics and airborne droplet nuclei distribution to variations in violent expiration events

Here, we use large eddy simulations to investigate the puff and droplet dynamics from violent expiratory events such as coughs and sneezes in the first few seconds following an ejection. For each of the eleven simulations considered, over 60,000 droplets are ejected and individually tracked using the point-particle Euler–Lagrange approach. We test the sensitivity of the puff and droplet dynamics to various parameters including the ejection volume, momentum, and orientation. We also explore the effect of the mouth shape on the aforementioned dynamics by considering elliptical and circular inlet cross-sections. The results from the simulations compare favorably with a recent theoretical framework put forth by Balachandar et al. in terms of the puff size and propagation velocity. More importantly however, the theory is able to accurately predict the number and size spectra of the potentially virus-laden droplet nuclei that remain airborne within the puff. We observe that the ejection angle and mouth shape do not significantly affect the puff and droplet dynamics. Additionally, we quantify the carrying capacity of the detached puff portions in terms of the number and size spectra of droplets/droplet nuclei suspended within.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Multirate partitioned Runge–Kutta methods for coupled Navier–Stokes equations

Earth system models are complex integrated models of atmosphere, ocean, sea ice, and land surface. Coupling the components can be a significant challenge due to the difference in physics, temporal, and spatial scales. Further, this study explores multirate partitioned Runge-Kutta methods for the fluid-fluid interaction problem and demonstrates its parallel performance by using the PETSc library. We consider compressible Navier-Stokes equations with gravity coupled through a rigid-lid interface. Our large-scale numerical experiments reveal that multirate partitioned Runge-Kutta coupling schemes (1) can conserve total mass; (2) have second-order accuracy in time; and (3) provide favorable strong- and weak-scaling performance on modern computing architectures. We also show that the speedup factors of multirate partitioned Runge-Kutta methods match theoretical expectations over their base (single-rate) method.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Surviving nanoscale interfacial stability in extreme thermal expansion contrast Zn(CN) 2 -epoxy resin matrix composites

Here, we report the remarkable stability at the nanoscale matrix-filler interface in a series of overall low coefficient of thermal expansion (CTE) zinc cyanide (ZC)/epoxy resin composites. These interfaces demonstrate the stability of highly contrasting materials. These composites are designed to optimize the properties of both materials: epoxy resins are valued for their high strength, superb chemical resistance, low dielectric properties, and adhesive abilities. ZC, a model flexible framework, exhibits negative thermal expansion (NTE) behavior and compensates for the epoxy resin’s detrimentally large thermal expansion values. Filled resins minimize the native 65 ppm/°C CTE values to 19 ppm/° C over a wide temperature range while maintaining linearity, yielding composites compatible with typical metal and ceramic substrates’—commonly bonded to the resin— 0–20 ppm/° C CTE values. Additionally, with significant intrinsic differences in the materials, we use wide-angle X-ray diffraction and atomic force microscopy to elucidate thermal movement. The ZC/epoxy interface is surprisingly resilient over 1,000 cycles between -55° C and 75° C, extremes typical to aerospace applications. Further, we show that the resin viscosity remains less than 21 Pa*s at up to 30 vol% ZC, an important consideration for workability. We also discuss settling distribution, glass transition temperature, flexural strength, and longevity. All results point favorably towards practical implementation in low thermal expansion application needs: filling the gap in component design.

42 ENGINEERING↗
DOE OSTI2024

Supersonic hot jet ablative testing and analysis of boron nitride nanotube hybrid composites

Boron nitride nanotubes (BNNTs) are high-strength, high-modulus nanotubes with high thermal and oxidative stabilities. Two hybrid composites were prepared with satin weave carbon fiber (CF) and resole-type phenolic resin: one with surface layers of BNNTs and one with alternating interlayers of BNNTs. The samples were subjected to hot jet tests that simulate realistic high-pressure-temperature conditions to understand the behavior of BNNTs under high-pressure erosion. Adding BNNTs to CF/phenolic laminates enhanced the ablation resistance by reinforcing the char material and mitigated localized thermal damage. Hybrid laminates exhibited up to 14% lower weight loss, 55% increase in flexural modulus, higher thermal diffusivity, and improved char yield and microstructure compared to CF/phenolic samples. The surface layer hybrid had many surviving nanotubes reinforcing the char and crystalline oxide structures that could mitigate further oxygen diffusion. Further, various characterization methods were used to deduce possible mechanisms and their products, indicating that BNNTs could serve as growth templates for direct crystalline boron oxide formation. Overall, hybrid BNNT/CF/phenolic laminates displayed better ablation resistance and favorable microstructure evolution under high-pressure conditions.

42 ENGINEERING↗
DOE OSTI2025

Initial assessment of alternative carbon fiber geometries for design of cost-effective compressive performance: Size effect studies

Carbon fiber provides opportunity to reduce weight in structural composites, including wind turbine blades, due to the material's superior specific stiffness and specific strength compared to alternatives. Despite these advantages, cost and compressive performance are considered weaknesses for carbon fiber products available today. Studies to produce low-cost carbon fiber alternatives, including the use of textile-derived precursor systems, have shown progress and merit through the DOE/ORNL low-cost carbon fiber initiatives. Here, this work focuses on enabling increases in compressive strength through design of the carbon fiber geometry, applicable to both textile and conventional precursor systems, while also providing opportunities to reduce carbon fiber processing costs. Fiber-resin interface and fiber alignment are among the most frequently cited factors controlling composite compressive performance. However, it is believed that there is opportunity in traditionally unexplored routes to increasing compressive strength through alteration of the carbon fiber geometry by increasing the fiber area moment of inertia and/or the fiber perimeter and interfacial area. This paper presents initial results from manufacturing carbon fiber materials to assess the impacts of carbon fiber size on tested composite compressive performance with projected neutral or even beneficial impact on fiber and composite manufacturing economics. Carbon fiber systems with increasing size illustrate a favorable correlation for compressive performance greater than predicted from a micromechanical failure model. The manufacturing and mechanical test results support the hypothesis of this work that alterations to fiber geometry can be used to produce improvements of the compressive strength of carbon fiber reinforced polymers and provide incentive for related work in designing alternative shapes to further enhance compressive performance.

42 ENGINEERING↗
DOE OSTI2026

Carbon fiber design improvements based on economic models

Circular fiber geometries are predominant in commercial carbon fiber material systems, but the use of this fiber shape has numerous limitations. Circular geometries have continuous symmetry, which is helpful for various processing considerations, but also have the largest possible maximum diffusion thickness for a given fiber area. This characteristic means that circular carbon fibers always have the highest material processing cost and lowest production throughput compared to any other fiber shape with the same area and tow count. To quantify material cost and other benefits for non-circular carbon fiber geometries, process models for polyacrylonitrile based carbon fiber production are developed in relationship to the carbon fiber shape, size, and tow count. For a given fiber shape, precursor production costs are shown to favor maximizing fiber size while conversion costs are minimized by the smallest fiber size. These competing cost trends result in a numerically optimal fiber size for a given shape and tow count, while both cost components are decreased by increasing tow count. Cost–performance tradeoffs for three lobe fiber geometries are studied by supplementing the cost trends with a numerical failure model to predict compressive strength for discrete shape variants. The shape selection is shown to be more sensitive to variations in cost than compressive strength while suboptimal shape designs can improve manufacturing robustness and achievable fiber volume fractions. Finally, an optimal three lobe carbon fiber is identified that balances the set of considerations while reducing costs and embodied energy and increasing production throughput compared to a commercial carbon fiber.

Carbon fiber↗
DOE OSTI2022

Optimization of matrix viscosity improves polypropylene fiber dispersion and properties of engineered cementitious composites

Engineered cementitious composites (ECC) is a durable cementitious material with high tensile ductility and strain-hardening characteristics. Although considered as a cost-effective fiber for ECC, polypropylene (PP) fiber is reportedly difficult to disperse in mortar matrix due to its high aspect ratio and hydrophobicity. In this study, the matrix viscosity of a low-carbon ECC based on limestone calcined clay cement was tailored as a variable to improve PP fiber dispersion under a pre-determined mixing protocol. The effect of matrix viscosity on the composite fresh and hardened properties was investigated experimentally. Results suggested an optimal range of matrix viscosity (10.3–11.5 Pa . s) favors the composite tensile strength and strain capacity at 28 days. At the optimal state with a 0.1 % viscosity modifying admixture (VMA)-to-binder mass ratio, PP-ECC achieved 7.0 % tensile strain capacity and 3.5 MPa ultimate tensile strength. When matrix viscosity falls outside the desired range, both ultimate tensile strength and strain capacity were diminished. By tailoring the VMA dosage, the matrix viscosity can be adjusted for desired fiber dispersion, workability, and mechanical properties. Finally, the findings of this study provide a technical reference for the practical design and application of PP-ECC.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Genetically pliable green algae for bioproduction of modified fatty acids, nutritional therapeutic oils, and biopharmaceuticals

Homologous recombination (HR) is an essential tool for complex metabolic engineering in yeast, but transgene integration into plant and green algal nuclear genomes predominantly occurs by non-homologous end-joining. Species of the closely related, oleaginous trebouxiophytes Auxenochlorella and Prototheca, are unusual among the green algae in that HR is the favored mechanism for DNA integration into the nuclear genome. This property enables locus-specific targeting of gene cassettes encoding multiple enzymes for manipulating existing biochemical pathways or introducing new functions. Genetic malleability, and regulatory approval for human consumption, coupled with robust fermentation performance at industrial scale, establishes Auxenochlorella and Prototheca as prime candidates for algal production of biochemicals and biomaterials. The examples presented here highlight strain improvement and engineering for synthesis of hydroxylated fatty acids for biomaterials, structured triglycerides resembling human milk fat for infant nutrition, very-long-chain mono- and polyunsaturated fatty acids with nutraceutical or therapeutic potential, and cannabinoids for pharmacological applications.

Moseley, Jeffrey L. [University of California, Ber↗
DOE OSTI2025

Chromium segregation-induced oxide evolution in Ni-10Cr alloys during high-temperature oxidation

The oxidation behavior of a Ni-10(wt%)Cr alloy under high-temperature O 2 conditions is investigated using transmission electron microscopy and first-principles calculations. Results reveal that chromium segregation plays a central role in driving the evolution of complex oxide phase structures during oxidation. At low Cr concentrations, Cr preferentially segregates to NiO grain boundaries or internal pores, substituting for Ni atoms and forming Ni(Cr)O solid solutions. As Cr content increases, enhanced diffusion promotes Cr penetration into the NiO lattice, leading to the formation of multiphase oxide structures. First-principles modeling corroborates these findings: at low Cr concentrations, Cr atoms favor surface and grain-boundary segregation, while higher concentrations lead to Cr aggregation within the NiO bulk. Furthermore, the integrated experimental-theoretical approach provides atomistic insights into Cr-mediated mass transport mechanisms during alloy oxidation and offers valuable guidance for controlling oxide growth kinetics and phase stability in Ni-Cr alloys, with implications for improving oxidation resistance in high-temperature structural applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A conservative phase-space moving-grid strategy for a 1D-2V Vlasov–Fokker–Planck Solver

In this work, we develop a conservative configuration- and velocity-space (i.e., phase-space) moving-grid strategy for the Vlasov–Fokker–Planck (VFP) equation in a planar geometry. The velocity-space grid is normalized and shifted in terms of the thermal speed and the bulk-fluid velocity, respectively. The configuration-space grid is moved according to a mesh-motion-partial-differential equation (MMPDE), which equidistributes a monitor function that is inversely proportional to the gradient-length scales of the macroscopic plasma quantities. The resulting inertial terms in the transformed VFP equations are discretized to ensure the discrete conservation of mass, momentum, and energy. To satisfy the discrete conservation theorems in the presence of phase-space mesh motion, we employ the method of discrete nonlinear constraints – explored in previous studies – but the underlying symmetries are determined in a much more efficient manner than before. The conservative grid-adaptivity strategy provides an efficient scheme that resolves important physical structures in the phase-space while controlling the computational complexity at all times. We demonstrate the favorable features of the proposed algorithm through a set of test cases of increasing complexity. The problems test independent components of the algorithms, as well as the integrated capability on settings relevant to inertial confinement fusion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Causal explicit algorithm for heat conduction in a plasma

Hyperbolic heat conduction extends standard Spitzer-Harm heat conduction by including a term proportional to the time derivative of the heat flux. The new term arises from a kinetic derivation of the heat flux that includes higher order corrections. Here we present a causal explicit numerical algorithm for solving the nonlinear hyperbolic heat conduction equation in an unmagnetized plasma. The maximum stable timestep for the causal explicit algorithm scales linearly with the cell size, owing to the hyperbolic nature of the problem. This is in contrast to the quadratic scaling of the maximum stable timestep with the cell size for the parabolic forward time centered space algorithm. The favorable scaling of the timestep with the cell size enables a practical explicit implementation of heat conduction in high-performance massively parallel plasma codes. In particular, we have implemented the causal explicit algorithm in the laser plasma interaction code pF3D. We verify the CE algorithm and analyze its convergence rate by simulating a harmonic mode, which has an analytic solution within the context of the HHC model. We also compare simulations using the CE algorithm to those using the forward time centered space algorithm on a pair of test problems: evolution in time of a Gaussian temperature perturbation in a uniform plasma and heat transport in the presence of inverse bremsstrahlung heating by a Gaussian laser speckle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

GPU-enabled extreme-scale turbulence simulations: Fourier pseudo-spectral algorithms at the exascale using OpenMP offloading

Fourier pseudo-spectral methods for nonlinear partial differential equations are of wide interest in many areas of advanced computational science, including direct numerical simulation of three-dimensional (3-D) turbulence governed by the Navier-Stokes equations in fluid dynamics. This paper presents a new capability for simulating turbulence at a new record resolution up to 35 trillion grid points, on the world's first exascale computer, Frontier, comprising AMD MI250x GPUs with HPE's Slingshot interconnect and operated by the US Department of Energy's Oak Ridge Leadership Computing Facility (OLCF). Key programming strategies designed to take maximum advantage of the machine architecture involve performing almost all computations on the GPU which has the same memory capacity as the CPU, performing all-to-all communication among sets of parallel processes directly on the GPU, and targeting GPUs efficiently using OpenMP offloading for intensive number-crunching including 1-D Fast Fourier Transforms (FFT) performed using AMD ROCm library calls. With 99% of computing power on Frontier being on the GPU, leaving the CPU idle leads to a net performance gain via avoiding the overhead of data movement between host and device except when needed for some I/O purposes. Memory footprint including the size of communication buffers for MPI_ALLTOALL is managed carefully to maximize the largest problem size possible for a given node count. Detailed performance data including separate contributions from different categories of operations to the elapsed wall time per step are reported for five grid resolutions, from 2048 3 on a single node to 32768 3 on 4096 or 8192 nodes out of 9408 on the system. Both 1D and 2D domain decompositions which divide a 3D periodic domain into slabs and pencils respectively are implemented. The present code suite (labeled by the acronym GESTS, GPUs for Extreme Scale Turbulence Simulations) achieves a figure of merit (in grid points per second) exceeding goals set in the Center for Accelerated Application Readiness (CAAR) program for Frontier. The performance attained is highly favorable in both weak scaling and strong scaling, with notable departures only for 2048 3 where communication is entirely intra-node, and for 32768 3 , where a challenge due to small message sizes does arise. Communication performance is addressed further using a lightweight test code that performs all-to-all communication in a manner matching the full turbulence simulation code. Performance at large problem sizes is affected by both small message size due to high node counts as well as dragonfly network topology features on the machine, but is consistent with official expectations of sustained performance on Frontier. Overall, although not perfect, the scalability achieved at the extreme problem size of 32768 3 (and up to 8192 nodes — which corresponds to hardware rated at just under 1 exaflop/sec of theoretical peak computational performance) is arguably better than the scalability observed using prior state-of-the-art algorithms on Frontier's predecessor machine (Summit) at OLCF. New science results for the study of intermittency in turbulence enabled by this code and its extensions are to be reported separately in the near future.

3D fast Fourier transform↗