Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dependence”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Mineralogy-dependent sulfide oxidation via polysulfide and thiosulfate pathways during weathering of mixed-sulfide bearing mine waste rock

Tracking the S oxidation pathway in sulfide-bearing mine waste-rock piles is complicated by variations in water content, O 2 and Fe 3+ concentrations, microbial diversity, mineralogy, the occurrence of a variety of S species associated with incomplete oxidation, and nonlinear coupling between physicochemical processes. Synchrotron-based S K-edge X-ray absorption near edge structure (XANES) spectroscopy facilitates the identification and quantification of S species (i.e., S 2– , $S_{2}O_{2}^{2–}$, S 8 , $S_{2}O_{3}^{2–}$, and $S_{2}O_{4}^{2–}$) in a weathered mine waste-rock pile containing pyrite, sphalerite, galena, pyrrhotite, and chalcopyrite. Mineralogy-dependent polysulfide and thiosulfate pathways both affect the S oxidation within the waste-rock pile. Currently, the polysulfide pathway, with S 8 as the intermediate S species, generates high concentrations of dissolved metals (from the rapid oxidation of monosulfides including sphalerite, galena, pyrrhotite, and chalcopyrite). As these monosulfides are depleted, the thiosulfate pathway of pyrite oxidation, including intermediate oxidation products (i.e., S 8 and $S_{2}O_{3}^{2–}$), will become the dominant oxidation pathway, resulting in the potential for generation of additional acidic drainage. This information highlights the importance of identifying the sulfide-mineral oxidation reaction pathways when attempting to estimate acid generation potential, and when developing remediation strategies for the storage and management of sulfide-bearing waste rock.

58 GEOSCIENCES↗
DOE OSTI2023

The physics of desiccation cracks 1: Ductile fracturing and dependence on relative humidity

Ductile deformation is ubiquitously found in the shrinkage of geomaterials. The existence of ductility requires elastoplastic mechanics when analyzing the structure deformation under external stress. Here, in this work, we explore the physics of ductile fracturing based on the results from desiccation experiments and triaxial tests. By using the digital image correlation (DIC) method to generate the strain maps of samples undergoing desiccation cracking under different relative humidities, we obtain results showing the previously postulated Cnoidal Wave ductile failure patterns propagating under atmospheric condition-controlled crack velocities. We then correlate these observations with rate-dependent plasticity models calibrated through triaxial tests undergoing several unloading–reloading cycles and velocity stepping. This work demonstrates the necessity to consider ductility in soil cracking, indicating that the formation of crack patterns in soil desiccation is a slow and predictable process.

58 GEOSCIENCES↗
DOE OSTI2023

Dynamical responses of constrained pipe conveying fluids and its dependence on the modeling of the contact force

Accurately modeling the impact force used in the analysis of loosely constrained cantilevered pipes conveying fluid is imperative. If little information is known of the motion-limiting constraints used in experiments, the analysis of the system may yield inaccurate predictions. Here in this work, multiple forcing representations of the impact force are defined and analyzed for a cantilevered pipe that conveys fluid. Depending on the representation of the impact force, the dynamics of the pipe can vary greatly when only the stiffness of the constraints is known from experiments. Three gap sizes of the constraints are analyzed, and the representation of the impact force used to analyze the system is found to significantly affect the response of the pipe at each gap size. An investigation on the effects of the vibro-impact force representation is performed through using basin of attraction analysis and nonlinear characterization of the system’s response.

42 ENGINEERING↗
DOE OSTI2022

Understanding syringeability and injectability of high molecular weight PEO solution through time-dependent force-distance profiles

Opioid medications play a vital role in treating moderate to severe pain. Unfortunately, many drug misusers and abusers attempt to alter the formulations or properties of these drugs by manipulation, (e.g., crushing, chewing, smoking, snorting, injecting). The intravenous (IV) route is most dangerous to abusers, as the drugs directly enter the circulatory system and produce intense euphoria. To obtain a full understanding of the impact of syringe factors (e.g., needle gauge size, needle length, syringe barrel size), on the ease of injection, we undertook a comprehensive assessment of syringeability and injectability of manipulated abuse-deterrent formulations (ADFs). A texture analyzer-based testing method was developed for the measurement of the resistance force of pulling, holding, and pushing phases of injections. Results showed that the finer needle gauge sizes required higher injection force to withdraw drug solutions. In addition, the syringed liquid volume was highly dependent on needle gauge size, holding time, and sample viscosity. In most cases, a lower needle gauge number and a longer holding time increased the syringed volume. Needle length was highly correlated to injection force (R 2 = 0.99). Using longer needles to inject drug solution requires greater force. Furthermore, large barrel size was correlated to pushing force (R 2 = 0.99); thus, increasing the difficulty of pushing the plunger of a large syringe with one hand. Lastly, relationships between injection force, sample viscosity, and testing conditions were elucidated using a mathematical model, which could be used in the future to assess and predict injection force of solution samples.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Temperature-dependent deformation behavior of a CuZr-based bulk metallic glass composite

In this study, the effect of cryogenic testing temperatures on the deformation behavior of a CuZr-based bulk metallic glass composite (BMGC) has been investigated using in situ high energy X-ray diffraction. The yield strength and work-hardenability of the studied BMGC show an increasing trend from 298 K to 153 K. A higher yield strength at a lower temperature can be attributed to the larger difference in the volume shrinkage between the crystalline phase and the amorphous matrix. The increasing martensitic transformation degree of the crystalline phase from 298 K to 153 K results in the enhanced work-hardenability. These results lead to a better understanding of the temperature dependence of phase-transformation induced plasticity in CuZr-based BMGCs. (C) 2020 Elsevier B.V. All rights reserved.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Insight Into the High-Potential Branch of the Alternate Nad+-Dependent NADPH:Ferredoxin Oxidoreductase II (NfnII) from Pyrococcus Furiosus. Evidence for the Gating Step in Electron Bifurcation

We have investigated the rapid-reaction kinetics of the NAD+-dependent NADPH:ferredoxin oxidoreductase II (NfnII) from Pyrococcus furiosus, permitting a comparison with recent work done with the paralog NfnI from the same organism. The half-potentials of the electron-bifurcating L-FAD are highly crossed in both NfnI, meaning the potential for the quinone/semiquinone couple is significantly lower than that for the semiquinone/hydroquinone couple so that the semiquinone oxidation state is thermodynamically unstable. The same appears to be the case with NfnII on the basis of its similar behavior in transient absorption spectroscopy experiments and the absence of any evidence for FAD*- accumulation in the course of reductive titrations (which would be manifested as a transient increase in absorbance at ~380 nm)1. Reductive titrations with both the one-electron donor sodium dithionite and the obligate two-electron donor NADPH demonstrate that, in contrast to NfnI, little FADH* accumulates in the high-potential pathway of NfnII in the course of reduction, as reflected in the absence of a transient increase in absorbance in the 500-600 nm region. This indicates that the half-potentials of the S-FAD are crossed in NfnII by a minimum of 120 mV. Rapid-reaction experiments mixing oxidized NfnII with NADPH also show no evidence of S-FADH* accumulation. Furthermore, the enzyme is only partially reduced at the end of the reaction with NADPH, indicating that there is little electron transfer into the high-potential pathway of NfnII. When the reaction is carried out in the presence of the one-electron carrier ferredoxin, only approximately one equivalent of ferredoxin becomes reduced, in contrast to the maximum of three equivalents seen with NfnI. This is consistent with only a single electron bifurcation event taking place under these conditions with NfnII, with only a single high-potential electron passing into the [2Fe-2S] cluster of the high-potential pathway of NfnII but not progressing further to the S-FAD. This accounts for the observation that NfnII has greatly reduced bifurcating activity compared to NfnI. The lack of electron transfer into the S-FAD of NfnII due to its crossed half-potentials prevents successive bifurcation activity and indicates that electron transfer within the high-potential branch gates bifurcation activity in NfnI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Distance-dependent dielectric constant at the calcite/electrolyte interface: Implication for surface complexation modeling

Hypothesis: The electrical double layer formed at the mineral/electrolyte interface is often modeled using mean-field approaches based on a continuum description of the solvent whose dielectric constant is assumed to decrease monotonically with decreasing distance to the surface. In contrast, molecular simulations show that the solvent polarizability oscillates near the surface similar to the water density profile - as shown previously, for example, by Bonthuis et al. (D.J. Bonthuis, S. Gekle, R.R. Netz, Dielectric Profile of Interfacial Water and its Effect on Double-Layer Capacitance, Phys Rev Lett 107(16) (2011) 166102). Here, we showed that molecular and mesoscale pictures agree by spatially averaging the dielectric constant obtained from molecular dynamics simulations over the distances relevant to the mean-field representation. In addition, the values of capacitances used to describe the electrical double layer in Surface Complexation Models (SCMs) of the mineral/electrolyte interface can be estimated using molecularly informed spatially averaged dielectric constants and positions of hydration layers. Experiments: First, we used molecular dynamics simulations to model the calcite 101¯4/electrolyte interface. Next, by using atomistic trajectories, we calculated the distance-dependent static dielectric constant and water density in the direction normal to the. Finally, we applied spatial compartmentalization consistent with the model of parallel-plate capacitors connected in series to estimate SCM capacitances. Findings: Computationally expensive simulations are required to determine the dielectric constant profile of interfacial water near the mineral surface. On the other hand, water density profiles are readily assessable from much shorter simulation trajectories. Our simulations confirmed that dielectric and water density oscillations at the interface are correlated. Here, we parametrized linear regression models to estimate the dielectric constant directly from the local water density. This is a significant computational shortcut compared to slowly converging calculations relying on total dipole moment fluctuations. The amplitude of the interfacial dielectric constant oscillation can exceed the dielectric constant of the bulk water, suggesting an ice-like frozen state, but only if there are no electrolyte ions. The interfacial accumulation of electrolyte ions causes a decrease in the dielectric constant due to the reduction of water density and re-orientation of water dipoles in ion hydration shells. Finally, we show how to use the computed dielectric properties to estimate SCM's capacitances.

58 GEOSCIENCES↗
DOE OSTI2026

Quantitative SANS and multi-model analysis of spacer-dependent micellization of urea-based gemini surfactants

The micellization behavior of urea-based cationic gemini surfactants was investigated using small-angle neutron scattering (SANS) with multi-model form factor analysis. A homologous series of surfactants with urea group included in the hydrophobic tail and polymethylene spacers consisting of two to ten methylene units was analyzed using three form factor models: a core–shell ellipsoid and two variants of homogeneous ellipsoids. The results from all models show a consistent trend of the micelle structures, confirming that the spacer length critically influences micellar geometry, aggregation number, and hydration. The surfactant with four CH 2 groups in the spacer formed the largest micelles with the highest aggregation number, while longer spacers led to progressively smaller, more compact aggregates. The shell hydration—quantified as the volume fraction of heavy water within the hydrophilic region—decreased systematically with increasing spacer length due to enhanced hydrophobicity of the headgroup-spacer region. Intermicellar interactions, modeled as screened Coulomb interaction using the rescaled mean spherical approximation (RMSA), revealed the strongest electrostatic repulsion for the case of four methylene groups in the spacer, corresponding to the highest micellar charge and largest interparticle spacing. The observed spacer-dependent trends were robust across all modeling approaches, demonstrating that the spacer length serves as a key structural determinant of self-assembly in this type of urea-based gemini systems. These findings provide insight into the design of gemini surfactants with tailored aggregation behavior for applications in drug delivery, nanostructure templating, and solubilization technologies.

Core–shell ellipsoid model↗
DOE OSTI2020

Time-step dependent force interpolation scheme for suppressing numerical Cherenkov instability in relativistic particle-in-cell simulations

The WT scheme, a piecewise polynomial force interpolation scheme with time-step dependency, is proposed in this paper for relativistic particle-in-cell (PIC) simulations. The WT scheme removes the lowest order numerical Cherenkov instability (NCI) growth rate for arbitrary time steps allowed by the Courant condition. While NCI from higher order resonances is still present, the numerical tests show that for smaller time steps, the numerical instability grows much slower than using the optimal time step found in previous studies. The WT scheme is efficient for improving the quality and flexibility of relativistic PIC simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Parallel transport dynamics for mixed quantum states with applications to time-dependent density functional theory

Direct simulation of the von Neumann dynamics for a general (pure or mixed) quantum state can often be expensive. One prominent example is the real-time time-dependent density functional theory (rt-TDDFT), a widely used framework for the first principle description of many-electron dynamics in chemical and materials systems. Practical rt-TDDFT calculations often avoid the direct simulation of the von Neumann equation, and solve instead a set of Schrödinger equations, of which the dynamics is equivalent to that of the von Neumann equation. However, the time step size employed by the Schrödinger dynamics is often much smaller. Here, in order to improve the time step size and the overall efficiency of the simulation, we generalize a recent work of the parallel transport (PT) dynamics for simulating pure states [An, Lin, Multiscale Model. Simul. 18, 612, 2020] to general quantum states. The PT dynamics provides the optimal gauge choice, and can employ a time step size comparable to that of the von Neumann dynamics. Going beyond the linear and near adiabatic regime in previous studies, we find that the error of the PT dynamics can be bounded by certain commutators between Hamiltonians, density matrices, and their derived quantities. Such a commutator structure is not present in the Schrödinger dynamics. We demonstrate that the parallel transport-implicit midpoint (PT-IM) method is a suitable method for simulating the PT dynamics, especially when the spectral radius of the Hamiltonian is large. The commutator structure of the error bound, and numerical results for model rt-TDDFT calculations in both linear and nonlinear regimes, confirm the advantage of the PT dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Conservative remapping of material-dependent fields between possibly misaligned material regions

In this work, we propose an interpolation or remapping algorithm of material-dependent fields on polyhedral meshes where any source or target cell contains only one material. It is conservative and it preserves sharp material boundaries on the target mesh, even if the source and target regions delineating the same material are slightly misaligned. If those material regions are aligned, then the algorithm is also linearity-preserving and bounds-preserving. For a given material, it consists of a conservative field reconstruction on a target mesh part from a source mesh part associated with that material, followed by a repair step in case of misaligned boundaries. No assumption is made regarding the topology of the input meshes

36 MATERIALS SCIENCE↗
DOE OSTI2024

A fully implicit, asymptotic-preserving, semi-Lagrangian algorithm for the time dependent anisotropic heat transport equation

In this paper, we extend the operator-split asymptotic-preserving, semi-Lagrangian algorithm for time dependent anisotropic heat transport equation proposed in Chacón et al. (2014) [18] to use a fully implicit time integration with backward differentiation formulas. The proposed implicit method can deal with arbitrary heat-transport anisotropy ratios $\mathcal{X}$∥ /$ \mathcal{X}$⟂ $\ggg$ 1 (with $\mathcal{X}$∥, $ \mathcal{X}$⟂ the parallel and perpendicular heat diffusivities, respectively) in complicated magnetic field topologies in an accurate and efficient manner. Further, the implicit algorithm is second-order accurate temporally and demonstrates an accurate treatment at boundary layers (e.g., island separatrices), which was not ensured by the operator-split implementation. The condition number of the resulting algebraic system is independent of the anisotropy ratio, and is inverted with preconditioned GMRES. We propose a simple preconditioner that renders the finite-dimensional linear operator compact, resulting in mesh-independent convergence rates for topologically simple magnetic fields, and convergence rates scaling as ~ (NΔt) 1/4 (with N the total mesh size and Δt the timestep) in topologically complex magnetic-field configurations. We demonstrate the accuracy and performance of the approach with test problems of varying complexity, including an analytically tractable boundary-layer problem in a straight magnetic field, and a topologically complex magnetic field featuring magnetic islands with extreme anisotropy ratios $\mathcal{X}$∥ /$ \mathcal{X}$⟂ = 10 10 ) .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A review of low-rank methods for time-dependent kinetic simulations

Time-dependent kinetic models are ubiquitous in computational science and engineering. The underlying integro-differential equations in these models are high-dimensional, comprised of a six–dimensional phase space, making simulations of such phenomena extremely expensive. In this article we demonstrate that in many situations, the solution to kinetics problems lives on a low dimensional manifold that can be described by a low-rank matrix or tensor approximation. We then review the recent development of so-called low-rank methods that evolve the solution on this manifold. The two classes of methods we review are the dynamical low-rank (DLR) method, which derives differential equations for the low-rank factors, and a Step-and-Truncate (SAT) approach, which projects the solution onto the low-rank representation after each time step. Thorough discussions of time integrators, tensor decompositions, and method properties such as structure preservation and computational efficiency are included. We further show examples of low-rank methods as applied to particle transport and plasma dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Crystallization-induced emission enhancement of alkyl chain-dependent pyrene-based luminogens: Visual detection of nitro-explosives

In view of the widespread applications and urgent need for efficient blue emitting materials, the synthesis of such materials has seen renewed interest. In the present work, three alkyl chain-dependent blue emitting pyrene-based luminogens were synthesized and characterized. Here, the experimental results show that these luminogens exhibited high photoluminescence efficiency both in solution (>86%) and in the solid state (>35%). Interestingly, 1,3,6,8-tetrakis (pentylphenyl)pyrene (TPPy), as a photoelectric functional material, was endowed with higher quantum yield (72%) and competitiveness due to a remarkable crystallization-induced emission enhancement (CIEE) property. More importantly, TPPy exhibits high sensitivity and selectivity for ortho-nitroaniline (o-NA) with low limit of detection (9.99 × 10 -8 M), which is mainly due to the distinctive odd-even effects of the terminal alkyl chain. The results demonstrate that this type of luminogen can potentially be utilized for practical applications in the field of contaminant or explosive detection as a sensitive and portable fluorimeter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Antiferromagnet thickness dependence and rotatable spins in exchange biased CoO/Fe films

The emergence of exchange bias and coercivity enhancement has been investigated in epitaxial CoO/Fe films with varied antiferromagnet (AF) thicknesses, even smaller than the critical value where the frozen CoO spins are detectable. Vector magnetometry and first-order reversal curve (FORC) measurements reveal different CoO thickness dependence of the exchange bias and coercivity enhancement, including the evolution of magnetization reversal from a high coercivity, low bias phase due to rotatable CoO moments to a high bias, low coercivity phase due to frozen CoO moments. Furthermore, the AF domain state is found to be metastable, which can be reoriented by external and exchange fields prior to the appearance of frozen spins, pointing to a generic origin of the training effect. Monte Carlo simulations show that the AF anisotropy energy barrier and the rotatable spins induced by magnetic field and exchange interaction at the interface are responsible for the observed effects.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Configurational stability of a crack propagating in a material with mode-dependent fracture energy – Part II: Drift of fracture facets in mixed-mode I+II+III

In earlier papers (Leblond et al., 2011; 2019), we presented linear stability analyses of the coplanar propagation of a crack loaded in mixed-mode I+III, based on a “double” propagation criterion combining (Griffith, 1920)’s energetic condition and (Goldstein and Salganik, 1974)’s principle of local symmetry. The difference between the two papers was that in the more recent one, the local value of the critical energy-release-rate was no longer considered as a constant, but heuristically allowed to depend upon the ratio of the local mode III to mode I stress intensity factors. This led to a much improved, qualitatively acceptable agreement of theory and experiments, for the “threshold” value of the ratio of the unperturbed mode III to mode I stress intensity factors, above which coplanar propagation becomes unstable. In this paper, the analysis is extended to the case where a small additional mode II loading component is present in the initially planar configuration of the crack, generating a small, general kink of this crack from the moment it is applied. The main new effect resulting from presence of such a loading component is that the instability modes present above the threshold must drift along the crack front during its propagation. This prediction may be useful for future theoretical interpretations of a number of experiments where such a drifting motion was indeed observed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Grain-boundary fracture mechanisms in Li 7 La 3 Zr 2 O 12 (LLZO) solid electrolytes: When phase transformation acts as a temperature-dependent toughening mechanism

Garnet-type, solid electrolytes, such as Li 7 La 3 Zr 2 O 12 (LLZO), are a promising alternative to liquid electrolytes for lithium-metal batteries. However, such solid-electrolyte materials frequently exhibit undesirable lithium (Li) metal plating and fracture along grain boundaries. In this work, we employ atomistic simulations to investigate the mechanisms and key fracture properties associated with intergranular fracture along one such boundary. Our results show that, in the case of a Σ5 (310) grain boundary, this boundary exhibits brittle fracture behavior, i.e. the absence of dislocation activity ahead of the propagating crack tip, accompanied with a decrease in work of separation, peak stress, and maximum stress intensity factor as the temperature increases from 300 K to 1500 K. As the crack propagates, we predict two temperature-dependent Li clustering regimes. For temperatures at or below 900 K, Li tends to cluster in the bulk region away from the crack plane driven by a void-coalescence mechanism concomitant a simultaneous cubic-to-tetragonal phase transition. The tetragonalization of LLZO in this temperature regime acts as an emerging toughening mechanism. At higher temperatures, this phase transition mechanism is suppressed leading to a more uniform distribution of Li throughout the grain-boundary system and lower fracture properties as compared to lower temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Application of pH-dependent environmental protection agency test method 1313 to low-activity nuclear waste glass

Here, in this study, we present results of a room temperature, pH-dependent leach test established by United States Environmental Protection Agency (EPA) Leaching Environmental Assessment Framework (LEAF) (EPA Method 1313 and further abbreviated EPA 1313) to study the alteration behavior of low-activity nuclear waste (LAW) glasses. Implementation of EPA 1313, which is closer to LAW disposal temperatures, could allow for more flexibility in the formulation of Hanford LAW glasses, specifically in terms of increased waste loading. A set of 48 LAW glasses were selected to test with EPA 1313. Briefly, the test method consists of adding glass to a number of solutions at different initial pH values and at a fixed mass-to-solution-volume ratio. The glass is then allowed to corrode for two days and the final solution pH and concentration of different glass elements in solution is measured. The test was modified from the published procedure to constrain particle size to provide more reproducible results, reduce the studied pH range to a more applicable range, and to reduce solution volume to reduce costs associated with material preparation. Compositional modeling was used to predict the measured elemental releases at various pH. The model predictions improved in tests where there was a relatively small difference between initial and final solution pH indicating that pH controls would result in better predictive models.

Accelerated aging↗