Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “difference scaling”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

N 5 -Scaling Excited-State-Specific Perturbation Theory

In this paper, we show that by working in a basis similar to that of the natural transition orbitals and using a modified zeroth-order Hamiltonian, the cost of a recently introduced perturbative correction to excited-state mean field theory can be reduced from seventh to fifth order in the system size. The (occupied) 2 (virtual) 3 asymptotic scaling matches that of ground-state second-order Møller-Plesset theory but with a significantly higher prefactor because the bottleneck is iterative: it appears in the Krylov-subspace-based solution of the linear equation that yields the first-order wave function. Here, we discuss the details of the modified zeroth-order Hamiltonian we use to reduce the cost and the automatic code generation process we used to derive and verify the cost scaling of the different terms. Overall, we find that our modifications have little impact on the method's accuracy, which remains competitive with singles and doubles equation-of-motion coupled cluster.

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

The Role of Phase Separation on Rayleigh-Plateau Type Instabilities in Alloys

Classical molecular dynamics (MD) simulations are used to investigate the role of phase separation (PS) on the Rayleigh-Plateau (RP) instability. Ni–Ag bulk structures are created at temperatures (2000 K and 1400 K) that generate different PS length scales, λ PS , relative to the RP instability length scale, λRP. Rectanguloids are then cut from the bulk structures and patterned with a perturbation of certain amplitude and wavelength, λ RP . It is found that when λ PS << λ RP (2000 K), the patterned rectanguloids break up into nanoparticles in a manner consistent with classical RP theory, whereas when λ PS << λ RP (1400 K), soluto-capillarity affects the RP instability significantly. Specifically, since Ag has a lower surface energy than Ni, Ag migrates to cover neighboring Ni regions, therefore modifying the RP instability. Thus, we demonstrate that the phase separation length scale of an immiscible alloy can be exploited to direct the assembly of functional bimetallic alloys.

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

Anisotropic Hot Spot Formation at a Grain Boundary in Shock-Compressed TATB High Explosive Crystal

Secondary high explosives (HEs) exhibit rich microstructure that promotes the formation of hot spots responsible for detonation initiation, but the role of microstructural interfaces remains poorly quantified. To this end, we develop extensions for the generalized crystal-cutting method (GCCM) to prepare molecular dynamics (MD) simulation cells containing grain boundaries (GBs) and other crystal–crystal interfaces with prescribed tilt and twist orientations. Using the GCCM, we perform MD simulations of shock interactions with a GB between the (001) and (100) crystal facets in the secondary HE TATB (1,3,5-triamino-2,4,6-trinitrobenzene). Our MD simulations reveal a strong directional dependence to the formation of a hot spot at the GB interface. In particular, transmission of the shock from the (001) grain to the (100) grain yields a hot spot in the (100) grain at the GB interface, whereas no hot spot is produced when an equivalent shock transits the GB in the opposite direction. We trace the origin of this GB anisotropy to three dominant factors: (1) the intrinsic differences in shock-deformation mechanisms and wave structures for the bulk (100) and (001) grains, which leads to distinct geometries and mechanical impedances upon shock arrival to the GB depending on which grains are donor or acceptor for the transmitted shock; (2) the different time intervals separating the initial shock rise and the formation of steady wave structures in the respective donor–acceptor configurations; and (3) the differences in time scales required to re-establish local thermal equilibrium. Interfacial hot spots form when these factors combine to impede development of a steady two-wave structure and instead induce a localized, pseudosingly shocked region that undergoes a higher rate of work production (resulting in a higher temperature) compared to when the steady two-wave structure develops further from the interface. The extensions to the GCCM approach presented here are anticipated to facilitate a wide range of MD studies that focus on understanding the role of crystal–crystal interfaces in molecular materials.

organic↗
DOE OSTI2022

Discrete Orientations of Interfacial Waters Direct Crystallization of Mica-Binding Proteins

Understanding the basis of templated molecular assembly on a solid surface requires a fundamental comprehension of both short- and long-range aqueous response to the surface under a variety of solution conditions. Herein we provide a detailed picture of how the molecular-scale response to different mica surfaces yields distinct solvent orientations that produce quasi-static directional potentials onto which macromolecules can adsorb. We connect this directionality to observed (a)symmetric epitaxial alignment of designed proteins onto these surfaces, corroborate our findings with 3D atomic force microscopy experiments, and identify slight differences in surface structure as the origin of this effect. Furthermore, our work provides a detailed picture of the intrinsic electrolyte response in the vicinity of mineral interfaces, with clear predictions for experiment, and highlights the role of solvent on the predictive assembly of hierarchical materials on mineral surfaces.

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

Distinct Surface and Bulk Thermal Behaviors of LiNi 0.6 Mn 0.2 Co 0.2 O 2 Cathode Materials as a Function of State of Charge

Understanding how structural and chemical transformations take place in particles under thermal conditions can inform designing thermally robust electrode materials. Such a study necessitates the use of diagnostic techniques that are capable of probing the transformations at multiple length scales and at different states of charge (SOC). In this study, the thermal behavior ofLiNi 0.6 Mn 0.2 Co 0.2 O 2 (NMC-622) was examined as a function of SOC, using an array of bulk and surface-sensitive techniques. In general, thermal stability decreases as lithium content is lowered and conversion in the bulk to progressively reduced metal oxides (spinels, rock salt) occurs as the temperature is raised. Hard X-ray absorption spectroscopy (XAS) and X-ray Raman spectroscopy (XRS) experiments, which probe the bulk, reveal that Ni and Co are eventually reduced when partially delithiated samples (regardless of the SOC) are heated, although Mn is not. Surface-sensitive synchrotron techniques, such as soft XAS and transmission X-ray microscopy (TXM), however, reveal that for 50% delithiated samples, apparent oxidation of nickel occurs at particle surfaces under some circumstances. This is partially compensated by reduction of cobalt but may also be a consequence of redistribution of lithium ions upon heating. TXM results indicate the movement of reduced nickel ions into particle interiors or oxidized nickel ions to the surface or both. These experiments illustrate the complexity of the thermal behavior of NMC cathode materials. The study also informs the importance of investigating the surface and bulk difference as a function of SOC when studying the thermal behaviors of battery materials.

25 ENERGY STORAGE↗
DOE OSTI2020

Permeability Hysteresis From Microchannels Opening During Dissolution/Reprecipitation Cycle

Abstract Permeability is a critical parameter for geological resources characterization. Its evolution with respect to porosity is particularly interesting and many research initiatives focus on deriving such relationships, to understand some hydraulic impacts of microstructure alteration. Permeability evolution from chemical reactions, for instance, can become complex as flow channels may open during rock dissolution. In this contribution, we show that this phenomenon can lead to irregular porosity‐permeability curves and permeability hysteresis after reprecipitation. Current approaches describing permeability as a simple function of porosity can therefore not capture this behavior, and we advocate instead the use of dynamic modeling for such scenarios. We demonstrate our approach by modeling a dissolution/precipitation cycle for a unit cell pore channel and quantify the process at larger scale on three different rock samples, whose microstructures are reconstructed from segmented micro‐Computerised Tomography scans.

Lesueur, Martin↗
DOE OSTI2024

Increased Asian Sulfate Aerosol Emissions Remarkably Enhance Sahel Summer Precipitation

Observational evidence shows that Sahel summer precipitation has experienced a considerable increase since the 1980s, coinciding with significant diverging trends of increased sulfate emissions in Asia and decreased emissions in Europe (dipole pattern of aerosols between Asia and Europe). The decrease in European sulfate aerosols has substantial effects on the Sahel summer precipitation increase, but the corresponding effect of increased Asian sulfate is unknown. Multi-model simulations in the Precipitation Driver and Response Model Intercomparison Project (PDRMIP) show, compared to decreased European aerosols, that increased Asian aerosols similarly enhance the Sahel summer precipitation but with different large-scale atmospheric circulation changes. Further analysis of the Sixth Coupled Model Intercomparison Project (CMIP6) simulations under historical attribution and various emission scenarios reinforces the results about the climate impacts of anthropogenic aerosols and suggests that in future scenarios with strong international cooperation and rapid climate mitigations (SSP2-45), the Sahel drought will be intensified likely due to the decline in Asian aerosol emissions. Our results suggest that Asian anthropogenic aerosols are likely a non-negligible driver of the recent recovery in Sahel precipitation amounts

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Q-dependent collective relaxation dynamics of glass-forming liquid Ca0.4K0.6(NO3)1.4 investigated by wide-angle neutron spin-echo

Abstract The relaxation behavior of glass formers exhibits spatial heterogeneity and dramatically changes upon cooling towards the glass transition. However, the underlying mechanisms of the dynamics at different microscopic length scales are not fully understood. Employing the recently developed wide-angle neutron spin-echo spectroscopy technique, we measured the Q -dependent coherent intermediate scattering function of a prototypical ionic glass former Ca 0.4 K 0.6 (NO 3 ) 1.4 , in the highly viscous liquid state. In contrast to the structure modulated dynamics for Q < 2.4 Å −1 , i.e., at and below the structure factor main peak, for Q > 2.4 Å −1 , beyond the first minimum above the structure factor main peak, the stretching exponent exhibits no temperature dependence and concomitantly the relaxation time shows smaller deviations from Arrhenius behavior. This finding indicates a change in the dominant relaxation mechanisms around a characteristic length of 2π/(2.4 Å −1 ) ≈ 2.6 Å, below which the relaxation process exhibits a temperature independent distribution and more Arrhenius-like behavior.

42 ENGINEERING↗
DOE OSTI2025

Effect of multiple calcination cycles on CO 2 capture efficiency during carbonation of MgO in a mineral looping process

Magnesium oxide (MgO) is considered as a potential sorbent for direct air capture of carbon dioxide in a looping process. Previous research on mineral looping for carbon capture from flue gas using MgO has shown deactivation of the sorbent with repeated cycles but repeated cycling for carbonation at ambient conditions has not been yet investigated. Here, we tested three cycles of carbonation for MgO nano-scale powders with different starting surface area. We find that carbonation efficiency is higher at higher surface area. No decrease in carbonation efficiency was observed with cycling, instead we find that carbonation efficiency and surface area are stable or even increasing with cycling. Based on our experimental data, we hypothesize that the carbonation in presence of relative humidity leads to first hydration of MgO, leading to formation of brucite (Mg(OH) 2 ). This formation of Mg(OH) 2 is a volume-increasing reaction, which leads to fracturing of MgO particles and results in an increase in surface area. We observed formation of amorphous and crystalline hydrated carbonates. The crystalline phase observed was nesquehonite (MgCO 3 ⋅ 3H 2 O). Our results show that MgO does not show any deactivation with repeated cycling for carbon capture at ambient conditions and in the presence of humidity. These findings therefore indicate that MgO is a suitable candidate as sorbent material for direct air capture of carbon dioxide.

Ambient weathering↗
DOE OSTI2020

Chapter 15: Summary and Outlook

In the various chapters of this book, numerous characterization techniques are presented that can be applied to thin-film solar cells to determine (micro)structural, compositional, electrical, and optoelectronic properties. What has yet to be more fully elucidated is to what extent these characterization techniques can be combined, in a correlative way, to enhance the information gathered on materials and devices. Indeed, it is valuable to consider combining techniques to verify the relevance of the measured materials properties or to obtain them on different length scales-to compare surface with bulk properties, or to correlate structure and composition of materials with electrical and optoelectronic properties.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2021

Test methodology for validation of multi-frequency models of renewable energy generators using small-signal perturbations

Providing trustworthy and accurate multi-frequency (or harmonic) models for renewable energy generators (REG) is an ongoing challenge for harmonic studies. There have been effective attempts to propose and design a test device to validate the harmonic models, mainly based on shunt current perturbations. However, using additional devices for perturbations is costly for converter-based test sites. This paper provides the test specifications to extend the application of the grid emulators for voltage perturbations and appropriate harmonic model validation. Besides, the effects of the sequence couplings, initial emissions, and power set-points on the test results have been overlooked in the literature. Considering these effects, this paper proposes a generic test methodology to obtain more accurate models in the sequence domain. The experimental verification of the proposed methodology is demonstrated using a 7 MVA grid emulator for testing of a 2 MVA photo-voltaic converter and a 2 MVA Type 3 wind turbine. This way, the test challenges, specifications, and recommendations are presented using the MW-scale experiments on different REGs. Furthermore, the effects of sequence couplings and initial emissions on the calculation results are investigated and compared. The proposed methodology is applicable for harmonic model validation as well as empirical modelling.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

In situ x-ray and electron scattering studies of oxide molecular beam epitaxial growth

We perform the first in situ synchrotron X-ray ray diffraction (SXRD) / reflection high energy electron diffraction (RHEED) studies on the growth of complex oxide thin films by molecular beam epitaxy. The unique deposition chamber, located at the Advanced Photon Source, allows the preparation of complex oxide samples with monolayer precision and facilitates the formation of direct correlations between in situ X-ray studies and the more prevalent RHEED investigations. Importantly, because SXRD and RHEED probe different atomic-scale processes during thin film synthesis, their concomitant use enables the extraction of details concerning growth behavior than one cannot determine from either probe alone. We describe the results of such in situ studies on the epitaxial growth of perovskite LaNiO3on (La 0.18 Sr 0.82 )(Al 0.59 Ta 0.41 )O 3 (001). We find that during the earliest stages of growth, the RHEED and X-ray signals do not agree with each other, demonstrating that while regular RHEED oscillations may imply high quality growth, the film-substrate interface can undergo significant changes during deposition due to the occurrence of interdiffusion at the growth temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Efficient numerical algorithm for multi-level ionization of high-atomic-number gases

An efficient numerical algorithm for laser driven multi-level ionization of high-atomic-number gases is proposed and implemented in an electromagnetic particle-in-cell code SPACE. The algorithm is based on analytical solutions to the system of differential equations describing ionization evolution. Using analytical solutions resolves the multiscale issue of ionization due to different characteristic time scales of ionization processes and the main code time step. Algorithm efficiency and memory requirements are significantly improved by using a locally reduced system of differential equations. The algorithm also assigns proper orbital quantum numbers and their projections to ionization states. The algorithm is verified and validated using experimental data.

Cheng, A. (ORCID:000000021945282X)↗
DOE OSTI2024

Estimating the carbon footprint of computational fluid dynamics

In this study, computational resources have grown exponentially in the past few decades. These machines make possible research and design in fields as diverse as medicine, astronomy, and engineering. Despite ever-increasing computational capabilities, direct simulation of complex systems has remained challenging owing to the degrees of freedom involved. At the cusp of exascale computing, high-resolution simulation of practical problems with minimal model assumptions may soon experience a renaissance. However, growing reliance on modern computers comes at the cost of a growing carbon footprint. To illustrate this, we examine historic computations in fluid dynamics where larger computers have afforded the opportunity to simulate flows at increasingly relevant Reynolds numbers. Under a variety of flow configurations, the carbon footprint of such simulations is found to scale roughly with the fourth power of Reynolds number. This is primarily explained by the computation cost in core-hours, which is also described by similar scaling, though regional differences in renewable energy use also play a role. Using the established correlation, we examine a large database of simulations to develop estimates for the carbon footprint of computational fluid dynamics in a given year. Collectively, the analysis provides an additional benchmark for new computations where, in addition to balancing considerations of model fidelity, carbon footprint should also be considered.

42 ENGINEERING↗
DOE OSTI2022

Single Gaussian process method for arbitrary tokamak regimes with a statistical analysis

Abstract Gaussian process regression is a Bayesian method for inferring profiles based on input data. The technique is increasing in popularity in the fusion community due to its many advantages over traditional fitting techniques including intrinsic uncertainty quantification and robustness to over-fitting. This work investigates the use of a new method, the change-point method, for handling the varying length scales found in different tokamak regimes. The use of the Student’s t-distribution for the Bayesian likelihood probability is also investigated and shown to be advantageous in providing good fits in profiles with many outliers. To compare different methods, synthetic data generated from analytic profiles is used to create a database enabling a quantitative statistical comparison of which methods perform the best. Using a full Bayesian approach with the change-point method, Matérn kernel for the prior probability, and Student’s t-distribution for the likelihood is shown to give the best results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Simulations of the fault current limiting operation of a long-length REBCO CORC ® superconducting cable cooled by helium gas

Conductor-on-round-core (CORC®) cables composed of rare-earth-barium-copper-oxide high-temperature superconducting (HTS) tapes are of great interest for power transmission applications due to their many advantages such as high power density, light weight, and low loss. Closed circulation loops of cryogenic helium gas can be used to cool HTS cables down to low temperatures to significantly improve their current-carrying capacity. Coupled circuit-electromagnetic-thermal finite element simulations implemented in the COMSOL Multiphysics package were developed, validated, and then used for simulating the fault current limiting (FCL) performance and the cooling processes of an 8-layer CORC® cable cooled with a flow of cryogenic helium gas. In the simulations, the temperature dependence of the electrical and thermal properties of all component materials is implemented for improved accuracy. To overcome computational challenges caused by the considerable difference in geometrical scales (i.e. few-µm-thick HTS layers versus 10 m-long HTS cable), the model is divided into two separate simulations. The first simulation is performed on the transverse cross-section of the cable to calculate the electric field, heating power and temperature rise in each component of a CORC® cable during FCL operation. The heating power calculated in the first simulation is transferred to the second model to simulate the cooling of a 10 m-long cable after the fault is cleared. The effect of the helium gas flow rate on the cooling process is also investigated to develop strategic approaches for optimizing cooling systems for HTS cables with FCL capability. The simulations indicated that a 40 ms fault with a voltage drop of 20 V m −1 along the cable can result in a temperature increase from 60 K to about 165 K inside the cable, and it takes about 500 s to cool the cable back to nearly 60 K with a flow of cold helium gas at a rate of 5 g s −1 .

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Numerical modeling of pedestal stability and broadband turbulence of wide-pedestal QH-mode plasmas on DIII-D

The wide-pedestal quiescent high confinement mode discovered on DIII-D in recent years is a stationary and quiescent H-mode (QH-mode) with the pedestal width exceeding EPED prediction by at least 25%. Its characteristics, such as low rotation, high energy confinement and edge localized mode-free operation, make it an attractive operation mode for future reactors. Linear and nonlinear simulations using BOUT++ reduced two fluid MHD models and awere carried out to investigate the bursty broadband turbulence often observed in the edge of wide-pedestal QH-mode plasmas. Two kinds of MHD-scale instabilities in different spatial locations within the pedestal were found in the simulations: one mild peeling–ballooning (PB) mode γ PB < 0.04ω A ) located near the minimum in E r</:sub> well propagating in ion diamagnetic drift direction; and one drift-Alfvén wave locates at smaller radius compared to E r</:sub> well propagating in the electron diamagnetic drift direction and unstable only when the parallel electron dynamics is included in the simulation. The coupling between drift wave and shear Alfvén wave provides a possible cause of the experimentally observed local profile flattening in the upper-pedestal. The rotation direction, mode location, as well as the wavenumber of these two modes from BOUT++ simulations agree reasonably well with the experimental measurements, while the lack of quantitative agreement is likely due to the lack of trapped electron physics in current fluid model. This work presents improved physics understanding of the pedestal stability and turbulence dynamics for wide-pedestal QH-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Watch and learn—a generalized approach for transferrable learning in deep neural networks via physical principles

Transfer learning refers to the use of knowledge gained while solving a machine learning task and applying it to the solution of a closely related problem. Such an approach has enabled scientific breakthroughs in computer vision and natural language processing where the weights learned in state-of-the-art models can be used to initialize models for other tasks which dramatically improve their performance and save computational time. Here we demonstrate an unsupervised learning approach augmented with basic physical principles that achieves fully transferrable learning for problems in statistical physics across different physical regimes. By coupling a sequence model based on a recurrent neural network to an extensive deep neural network, we are able to learn the equilibrium probability distributions and inter-particle interaction models of classical statistical mechanical systems. Our approach, distribution-consistent learning, DCL, is a general strategy that works for a variety of canonical statistical mechanical models (Ising and Potts) as well as disordered interaction potentials. Using data collected from a single set of observation conditions, DCL successfully extrapolates across all temperatures, thermodynamic phases, and can be applied to different length-scales. This constitutes a fully transferrable physics-based learning in a generalizable approach.

97 MATHEMATICS AND COMPUTING↗