Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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

Interfacial Effects on Transport Coefficient Measurements in Li-ion Battery Electrolytes

Development of Li + -containing electrolytes with improved transport properties requires reliable, reproducible, and ideally low volume techniques to rigorously understand ion-transport with varying composition. Precisely measuring the complete set of transport coefficients in liquid electrolytes under battery-relevant operating conditions is difficult and the reliability of these methods are sparsely described in electrolyte transport literature. In this work, we apply a potentiostatic polarization-based transport characterization approach typically used for polymer electrolytes to liquid electrolyte systems in an attempt to fully measure all transport coefficients (conductivity, total salt diffusion coefficient, thermodynamic factor and transference number) for the model system of LiPF 6 in an ethylene carbonate—ethyl methyl carbonate (EC:EMC) mixture. Using systematic timescale and statistical analyses, we find that transport coefficients measured using potentiostatic polarization of Li-Li symmetric cells exhibit strong correlation to Li electrode interfacial resistance, indicating that such methods are probing both bulk and interfacial phenomena. This reveals a major roadblock in characterizing electrolyte systems where the interfacial resistance is significantly larger than ohmic electrolyte resistance. As a result, we find that methods that rely on potentiostatic Li metal stripping/plating do not readily result in reliable liquid electrolyte transport coefficients, unlike similar methods for solid polymer electrolytes, where interfacial resistances are typically smaller than electrolyte resistances at the elevated temperatures typically of interest for such electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2018

Measurement and modeling of the multiwavelength optical properties of uncoated flame-generated soot

Optical properties of flame-generated black carbon (BC) containing soot particles were quantified at multiple wavelengths for particles produced using two different flames: a methane diffusion flame and an ethylene premixed flame. Measurements were made for (i) nascent soot particles,(ii) thermally denuded nascent particles, and (iii) particles that were coated and then thermally denuded, leading to the collapse of the initiallylacy, fractal-like morphology. The measured mass absorption coefficients(MACs) depended on soot maturity and generation but were similar between flames for similar conditions. For mature soot, here corresponding to particles with volume-equivalent diameters > ~ 160 nm, the MAC and absorption Ångström exponent (AAE) values were independent of particle collapse while the single-scatter albedo increased. The MAC values for these larger particles were also size-independent. The mean MAC value at 532 nm for larger particles was 9.1±1.1 m 2 g -1 , about 17 % higher than that recommended by Bond and Bergstrom (2006), and the AAE was close to unity. Effective, theory-specific complex refractive index(RI) values are derived from the observations with two widely used methods:Lorenz–Mie theory and the Rayleigh–Debye–Gans (RDG) approximation. Mietheory systematically underpredicts the observed absorption cross sections at all wavelengths for larger particles (with x >0.9 ) independent of the complex RI used, while RDG provides good agreement. (The dimensionless size parameter x = π d p / λ , where d p is particle diameter and λ is wavelength.) Importantly, this implies that the use of Mie theory within air quality and climate models, as is common, likely leads to under predictions in the absorption by BC, with the extent of underprediction depending on the assumed BC size distribution and complex RI used. We suggest that it is more appropriate to assume a constant,size-independent (but wavelength-specific) MAC to represent absorption by uncoated BC particles within models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Artificial ground reflector size and position effects on energy yield and economics of single‐axis‐tracked bifacial photovoltaics

Abstract Artificial ground reflectors improve bifacial energy yield by increasing both front and rear‐incident irradiance. Studies have demonstrated an increase in energy yield due to the addition of artificial reflectors; however, they have not addressed the effect of varying reflector dimensions and placement on system performance and the impact of these parameters on the reflectors' financial viability. We studied the effect of high albedo (70% reflective) artificial reflectors on single‐axis‐tracked bifacial photovoltaic systems through ray‐trace modeling and field measurements. In the field, we tested a range of reflector configurations by varying reflector size and placement and demonstrated that reflectors increased daily energy yield up to 6.2% relative to natural albedo for PERC modules. To confirm the accuracy of our model, we compared modeled and measured power and found a root mean square error (RMSE) of 5.4% on an hourly basis. We modeled a typical meteorological year in Golden, Colorado, to demonstrate the effects of artificial reflectors under a wide range of operating conditions. Seventy percent reflective material can increase total incident irradiance by 1.9%–8.6% and total energy yield by 0.9%–4.5% annually after clipping is considered with a DC–AC ratio of 1.2. Clipping has a significant effect on reflector impact and must be included when assessing reflector viability because it reduces reflector energy gain. We calculated a maximum viable cost for these improvements of up to $2.50–4.60/m 2 , including both material and installation, in Golden. We expanded our analysis to cover a latitude range of 32–48°N and demonstrated that higher‐latitude installations with lower energy yield and higher diffuse irradiance content can support higher reflector costs. In both modeling and field tests, and for all locations, the ideal placement of the reflectors was found to be directly underneath the module due to the optimized rear irradiance increase.

14 SOLAR ENERGY↗
DOE OSTI2023

Effects of surface vacancy-adatom pair formation on PFC tungsten surface morphological response

Abstract We report a simulation study of the effect of He-irradiation-induced surface vacancy-adatom pair formation on the surface morphological evolution of plasma-facing component (PFC) tungsten and examine a number of factors that impact such evolution. Our analysis is based on self-consistent dynamical simulations according to an atomistically-informed, continuum-scale surface evolution model that has been developed following a hierarchical multiscale modeling strategy and can access the spatiotemporal scales of relevance to fuzz formation. The model accounts for the flux of surface adatoms generated as a result of the surface vacancy-adatom pair formation effect upon He implantation, which contributes to the anisotropic growth of surface nanostructural features due to the different rates of adatom diffusion along and across step edges of islands on the tungsten surface. We have carried out atomic-scale computations of optimal diffusion pathways along and across island step edges on the W(110) surface and calculated Ehrlich–Schwoebel (ES) barriers in adatom diffusion along and across such step edges. This aspect of surface adatom diffusion contributes to anisotropic surface atomic fluxes, terrace and step diffusive currents, and has been incorporated into our PFC surface evolution model, which predicts the formation of preferentially aligned nanoridge stripe patterns on the PFC surface. We establish that these anisotropic diffusive currents accelerate nanotendril growth on the PFC surface and the onset of surface nanostructure pattern formation. We also explore systematically the dependence of the PFC surface morphological response on the surface temperature and He ion incident flux, characterize in detail the resulting surface topographies and growth kinetics, and compare the predicted surface morphologies with experimental observations. Our simulation predictions for the emerging surface nanostructure patterns under certain plasma exposure conditions are consistent with experimental findings in the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Kinetic Model for Predicting Trace Gas Uptake and Reaction

A model is developed to describe trace gas uptake and reaction with applications to aerosols and microdroplets. Gas uptake by the liquid is formulated as a coupled equilibria that links gas, surface, and bulk regions of the droplet or solution. Previously, this framework was used in explicit stochastic reaction-diffusion simulations to predict the reactive uptake kinetics of ozone with droplets containing aqueous aconitic acid, maleic acid, and sodium nitrite. With the use of prior data and simulation results, a new equation for the uptake coefficient is derived, which accounts for both surface and bulk reactions. Lambert W functions are used to obtain closed form solutions to the integrated rate laws for the multiphase kinetics; similar to previous expressions that describe Michaelis-Menten enzyme kinetics. Together these equations couple interface and bulk processes over a wide range of conditions and do not require many of the limiting assumptions needed to apply resistor model formulations to explain trace gas uptake and reaction.

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

Entropy modeling on serrated flows in carburized steels

Samples comprised of carburized steel alloys, 4615, 4720, and 10B22, were tension tested at room temperature and at a strain rate of 2 × 10 -4 s -1 . The serrated-flow behavior was modeled and analyzed, using the refined composite multiscale entropy (RCMSE) method. High energy X-ray diffraction revealed the body-centered-cubic (BCC) structure of the three alloys in the as-received condition. Moreover, Vickers-hardness experiments were performed, which showed an increase in the hardness of the alloy that corresponded to increasing impurity concentrations. Here the diffusion coefficient of carbon in the three alloys were estimated, respectively, to be 2.5 × 10 -11 m 2 /s, 5.1 × 10 -11 m 2 /s, and 2.8 × 10 -11 m 2 /s for alloys, 10B22, 4615, and 4720, respectively, which agree with the values found in the literature for α-iron. Moreover, the sample entropy, on average, increased with respect to the heating time, and, consequently, the amount of carbon that had diffused into the alloy. Finally, the above results show that the complexity of the serration behavior in the alloys studied increases with respect to the amount of carbon that has diffused into the material.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Novel vertical silicon photodiodes based on salicided polysilicon trenched contacts

The classical concept of silicon photodiodes comprises of a planar design characterized by heavily doped emitters. Such geometry has low collection efficiency of the photons absorbed close to the surface. An alternative, promising, approach is to use a vertical design. Nevertheless, realization of such design is technologically challenged, hence hardly explored. Herein, a novel type of silicon photodiodes, based on salicided polysilicon trenched contacts, is presented. These contacts can be prepared up to 10 μm in depth, without showing any leakage current associated with the increase in the contact area. Consequently, the trenched photodiodes revealed better performance than no-trench photodiodes. A simple two dimensional model was developed, allowing to estimate the conditions under which a vertical design has the potential to have better performance than that of a planar design. At large, the deeper the trench is, the better is the vertical design relative to the planar (up to 10 μm for silicon). The vertical design is more advantageous for materials characterized by short diffusion lengths of the carriers. Salicided polysilicon trenched contacts open new opportunities for the design of solar cells and image sensors. For example, these contacts may passivate high contact area buried contacts, by virtue of the conformity of polysilicon interlayer, thus lowering the via resistance induced recombination enhancement effect.

14 SOLAR ENERGY↗
DOE OSTI2019

Formulation analysis and computation of an optimization-based local-to-nonlocal coupling method

We present an optimization-based coupling method for local and nonlocal continuum models. Our approach couches the coupling of the models into a control problem where the states are the solutions of the nonlocal and local equations, the objective is to minimize their mismatch on the overlap of the local and nonlocal problem domains, and the virtual controls are the nonlocal volume constraint and the local boundary condition. We present the method in the context of Local-to-Nonlocal diffusion coupling. Numerical examples illustrate the theoretical properties of the approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Mass transfer coefficients, kL, and air-CO2 ingassing rates in 3.4 m2 and 1-acre raceway ponds.

The US DOE 2016 Billion Ton Update (Langholtz 2016) projected that CO2 in flue-gas from power plants or similar sources would limit U.S. algal biofuel potential to under 5 billion gallons gasoline equivalent (gge) per year, due to lack of sufficient land, water and other requirements near such flue gas sources. 2nd generation carbon capture technologies are proposed to overcome this constraint (Davis 2017), and could expand this resource potential nearly 10-fold, assuming CO2 costs near the future, 2025 DOE NETL, goal of $40/metric ton CO2. An alternative pathway is the direct uptake of air CO2 into the algal ponds through an increase in the air-CO2 transfer rate resulting from the reaction of CO2 with hydroxide ions. However, such a ‘chemical enhancement’ in mass transfer depends on a high pH in the culture, and, by extension, the ability of microalgae to maintain high rates of carbon assimilation under such conditions. Abiotic experiments characterized the air-CO2 mass-transfer rate in 1-acre and 3.4 m2 raceway ponds as a function of pH, with rates approaching 10 g C/m2-day at pH 12 in the brackish water tested, equivalent to a biomass productivity potential of near 20 g AFDW/m2-day (@ 0.47 g C/g AFDW). The mass transfer rate was found to be nearly independent of the mass-transfer coefficient, indicating that at turbulence levels achievable in typical raceway ponds, the system is reaction rate, rather than diffusion limited. In biotic trials with an unknown green microalga, biomass productivity in ponds receiving CO2 from only direct surface air-CO2 exchange averaged 5 +/- 0.5 g AFDW/m2-day in August - September central California conditions, about half-that of experimental controls fertilized with supplemental CO2. pH in the air-only treatment reached a maximum of 10.5, supporting a model predicted air-CO2 ingassing rate between 2-3 g C/m2-day, consistent with the observed biomass carbon assimilation rate. Results suggest that unique alkaliphilic strains are required if reliant on direct in-pond air-CO2 transfer for algal inorganic carbon supply, and that reaching an economically viable biomass productivity will require strains that thrive at pH 11 and above.

09 BIOMASS FUELS↗
DOE OSTI2016

Analysis of non-equilibrium phenomena in inductively coupled plasma generators

This article addresses the modeling of non-equilibrium phenomena in inductively coupled plasma discharges. In the proposed computational model, the electromagnetic induction equation is solved together with the set of Navier-Stokes equations in order to compute the electromagnetic and flow fields, accounting for their mutual interaction. Semi-classical statistical thermodynamics is used to determine the plasma thermodynamic properties, while transport properties are obtained from kinetic principles, with the method of Chapman and Enskog. Particle ambipolar diffusive fluxes are found by solving the Stefan-Maxwell equations with a simple iterative method. Two physico-mathematical formulations are used to model the chemical reaction processes: (1) A Local Thermodynamics Equilibrium (LTE) formulation and (2) a thermo-chemical non-equilibrium (TCNEQ) formulation. In the TCNEQ model, thermal non-equilibrium between the translational energy mode of the gas and the vibrational energy mode of individual molecules is accounted for. The electronic states of the chemical species are assumed in equilibrium with the vibrational temperature, whereas the rotational energy mode is assumed to be equilibrated with translation. Three different physical models are used to account for the coupling of chemistry and energy transfer processes. Numerical simulations obtained with the LTE and TCNEQ formulations are used to characterize the extent of non-equilibrium of the flow inside the Plasmatron facility at the von Karman Institute. Each model was tested using different kinetic mechanisms to assess the sensitivity of the results to variations in the reaction parameters. A comparison of temperatures and composition profiles at the outlet of the torch demonstrates that the flow is in non-equilibrium for operating conditions characterized by pressures below 30 000 Pa, frequency 0.37 MHz, input power 80 kW, and mass flow 8 g/s.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Trapped electron plasma formation and equilibrium with a low-power radio-frequency drive

Penning-Malmberg traps confining electron plasmas usually rely on external sources like thermo- and photocathodes. It has been already demonstrated that electron plasmas of comparable densities can be produced by applying a radio-frequency (RF) power to any inner electrode of the trap. Such excitation may result in significant electron heating and ionization of the residual gas with the formation of a plasma column when the RF frequency is of the order or larger than the typical axial bounce frequencies of few-eV electrons, even at RF amplitude of few volts. While discharges are common in plasma generation at higher pressures and RF power, this mechanism is not yet well explored in our working conditions, namely ultra-high vacuum and very low RF power. This plasma production mechanism is very sensitive to the experimental conditions. Interesting phenomena can be observed: transition from a diffuse to a narrow-section, denser plasma column; presence of low-order diocotron modes in transient and steady-state plasmas; modulation of the m=1 diocotron mode and suppression of its instability despite the presence of positive ions and resistive loads. These observations are reported here, and possible explanations are discussed. In addition, a possible electron heating mechanism is investigated with a single-particle, one-dimensional model described by an area-preserving map where an electron bounces within a square potential well and the RF excitation is modelled by a time-oscillating square barrier. The low-energy part of the Poincaré plot includes both quasi-periodic and chaotic regions, where heating up to ionization energies is achievable. Results of a systematic analysis of the map extracting its chaotic properties and scaling laws as a function of the control parameters are reported.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Analytical population dynamics underlying harmful algal blooms triggered by prey avoidance

A modified version of the NPZD ecosystem model is used to analytically examine the effects of predation avoidance, a possible mechanism for triggering harmful algal blooms (HAB). To resolve HAB development caused by predation avoidance, an additional phytoplankton functional group is considered, one that has slower nutrient uptake and better predation avoidance characteristics than the non-harmful phytoplankton group used in traditional NPZD models. Because the two phytoplankton groups (one non-harmful and one HAB) compete for only one resource within the same system, steady state (equilibrium) conditions cannot occur without the presence of zooplankton; only the non-harmful phytoplankton group, which defeats the HAB group in the resource competition, can survive in the equilibrium. The presence of sufficient zooplankton effectively acts to replenish the nutrient pool by consuming the non-harmful phytoplankton. When this occurs, two equilibrium states are found: one with both phytoplankton groups coexisting, and one that only includes the HAB group. Further, the condition required for equilibrium is that the total nitrogen within the system should be larger than a threshold determined by model coefficients. The threshold and feasibility of the equilibrium are sensitive to the relative HAB predation avoidance coefficient. If the coefficient is larger than the ratio of net growth rates between the HAB and non-harmful phytoplankton group, the threshold becomes infinite, and an equilibrium is not feasible. The time scale for the system to reach an equilibrium state that includes a HAB group is determined asymptotically. The dependence of a threshold condition as a controlling factor may explain the regime shift of dominant species causing HABs. The ecosystem model is fully implemented into the Regional Ocean Modeling System and applied to an idealized coastal embayment (with depths and geometry taken from San Francisco Bay) to show numerically the dominance of prey avoidance dynamics in a natural shallow water environment that includes advection and diffusion. The analytical results improve strategies for HAB modeling and provide guidance for setting model coefficients necessary to resolve a HAB event.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

Pseudomorphic growth of InAs on misoriented GaAs for extending quantum cascade laser wavelength

The authors have studied the impact of epilayer strain on the deposition of InAs/GaAs on (100) and (111)B with 2° offset toward<2-1-1> surfaces. Consequences of a 7% lattice mismatch between these orientations in the form of three-dimensional growth are less apparent for (111)B with 2° offset toward<2-1-1> surfaces compared to (100). By exploring a range of molecular beam epitaxy process parameters for InAs/GaAs growth and utilizing scanning electron microscopy, atomic force microscopy, and Raman spectroscopy to evaluate the quality of these strained layers, the authors develop empirical models that describe the influence of the process conditions in regards to surface roughness with >92% accuracy. The smoothest InAs/GaAs samples demonstrated average surface roughness of 0.08 nm for 10 μm{sup 2} areas, albeit at very low deposition rates. The authors have found the most important process conditions to be substrate temperature and deposition rate, leading us to believe that controlling diffusion length may be the key to reducing defects in severely strained structures. InGaAs/AlGaAs quantum cascade laser structures were also produced on (111)B with 2° offset toward<2-1-1> to take advantage of the piezoelectric effect, and the modified laser transitions due to these effects were observed.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Study on hydrogen isotopes permeation in fluidized state of liquid lithium-lead

Lithium-lead (Li-Pb) is one of the most promising candidate materials for the liquid blanket of fusion reactors. Hydrogen transfer under a fluidized condition of Li-Pb is investigated experimentally to design a Li-Pb blanket system. Li-Pb eutectic alloy flows inside a Ni tube in the experimental system, where H{sub 2} permeates into and out of the forced Li-Pb flow. The overall H{sub 2} permeation rate is analyzed using a mass balance model. Hydrogen atoms diffuse in Ni and Li-Pb. The steady-state H{sub 2} permeation rate obtained by this experiment is smaller than the result of the calculation model. A resistance factor is introduced to the present analysis in order to evaluate the influence of other H{sub 2} transfer mechanisms, such as diffusion in Li-Pb and dissolution reaction between Ni and Li-Pb. The contribution of the resistance to the overall H{sub 2} permeation rate becomes large when the flow rate of Li-Pb is low. This is because the boundary layer thickness between Ni and Li-Pb affects the overall H{sub 2} permeation rate. When the flow velocity of Li-Pb increases, the thickness of the boundary layer becomes thin, and the driving force of H{sub 2} permeation through the Ni wall becomes bigger. (authors)

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

Analytical homogenization techniques applied to the Fickian diffusion: Effective diffusivity coefficient

For multiple applications in nuclear energy, the ability to accurately represent material behavior with a simplified model is important to facilitate practical engineering-scale simulations. In this work, we focus on the homogenized thermal response of a medium containing spherical inclusions, similar to a fuel form (compact or pebble) containing TRISO particles. An extensive survey on effective thermal conductivity modeling was performed in our previous study, considering a random distribution of mono-sized spherical inclusions in a continuous matrix. Using the analogy between heat conduction and the simplified Fickian diffusion (or fission product species conservation), we can use the same analytical homogenization methods to obtain ETC as for the effective diffusivity coefficient (EDC). We performed several numerical experiments at varying conditions to assess the validity of our hypothesis for EDC calculations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Connecting the failure of K-theory inside and above vegetation canopies and ejection-sweep cycles by a large eddy simulation

Parameterizations of biosphere-atmosphere interaction processes in climate models and other hydrological applications require characterization of turbulent transport of momentum and scalars between vegetation canopies and the atmosphere, which is often modeled using a turbulent analogy to molecular diffusion processes. However, simple flux-gradient approaches (K-theory) fail for canopy turbulence. One cause is turbulent transport by large coherent eddies at the canopy scale, which can be linked to sweep-ejection events, and bear signatures of non-local organized eddy motions. K-theory, that parameterizes the turbulent flux or stress proportional to the local concentration or velocity gradient, fails to account for these non-local organized motions. The connection to sweep-ejection cycles and the local turbulent flux can be traced back to the turbulence triple moment $\overline{C'W'W'}$. In this work, we use large-eddy simulation to investigate the diagnostic connection between the failure of K-theory and sweep-ejection motions. Analyzed schemes are quadrant analysis (QA) and a complete and incomplete cumulant expansion (CEM and ICEM) method. The latter approaches introduce a turbulence timescale in the modeling. Furthermore, we find that the momentum flux needs a different formulation for the turbulence timescale than the sensible heat flux. In conclusion, accounting for buoyancy in stratified conditions is also deemed to be important in addition to accounting for non-local events to predict the correct momentum or scalar fluxes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Slip velocity boundary conditions for the lattice Boltzmann modeling of microchannel flows

Abstract Slip flows in ducts are important in numerous engineering applications, most notably in microchannel flows. Compared to the standard no‐slip Dirichlet condition, the case of slip formulates as a Robin‐type condition for the fluid tangential velocity. Such an increase in mathematical complexity is accompanied by a more challenging numerical transcription. The present work concerns with this topic, addressing the modeling of the slip velocity boundary condition in the lattice Boltzmann method (LBM) applied to steady slow viscous flows inside ducts of nontrivial shapes. As novelty, we extend the newly revised local second‐order boundary (LSOB) Dirichlet fluid flow method [ Philos. Trans. R. Soc. A 378, 20190404 (2020)] to implement the slip velocity condition within the two‐relaxation‐time (TRT) framework. The LSOB follows an in‐node philosophy where its operation principle seeks to explicitly reconstruct the unknown boundary populations in the form of a third‐order accurate Chapman–Enskog expansion, where the wall slip condition is built‐in as a normal Taylor‐type condition. The key point of this approach is that the required first‐ and second‐order momentum derivatives, rather than computed through nonlocal finite difference approximations, are locally determined through a simple local linear algebra procedure, whose formulation is particularly aided by the TRT symmetry argument. To express the obtained derivatives, two approaches are considered, called and , which operate with node and wall variables, respectively. These two formulations are developed to prescribe the physical slip condition over plane and curved walls, including the corners. Their consistency and accuracy characteristics are examined against alternative linkwise strategies to impose the wall slip velocity, such as the kinetic‐based diffusive bounce‐back scheme, the central linear interpolation slip scheme, and the multireflection slip scheme. The several slip schemes are tested over different 3D microchannel configurations, with walls not conforming with the LBM uniform mesh. Numerical tests confirm the advanced accuracy characteristics of the proposed LSOB slip boundary scheme, revealing the added challenge of the wall slip modeling, and that parabolic accuracy is a necessary requirement to reach second‐order accuracy within this problem class.

Silva, Goncalo↗
DOE OSTI2020

First-principles equations of state and structures of liquid metals in multi-megabar conditions

Liquid metals at extreme pressures and temperatures are widely interested in the high-pressure community. Based on density functional theory molecular dynamics, we conduct first-principles investigations on the equation of state ( EOS) and structures of four metals (Cu, Fe, Pb, and Sn) at 1.5–5 megabar conditions and 5 × 10 3 –4 × 10 4 K. Our first-principles EOS data enable evaluating the performance of four EOS models in predicting Hugoniot densities and temperatures of the four systems. We find the melting temperature of Cu is 1000–2000 K higher and shows a similar Clapeyron slope, in comparison to those of Fe. Our structure, coordination number, and diffusivity analysis indicates all the four liquid metals form similar simple close-packed structures. Our results set theoretical benchmarks for EOS development and structures of metals in their liquid states and under dynamic compression.

36 MATERIALS SCIENCE↗