Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “boundary effects”

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 127 records · Page 7

Effect of finite ion mass on relativistic self-induced transparency of plasma layers with a sharp boundary

We consider the effect of ion mobility on the threshold of relativistic self-induced transparency under irradiation of thin plasma layers by circularly polarised laser radiation. An analytical model of the motion of ions during the removal of electrons from the layer surface by the ponderomotive force of laser radiation is constructed. The model is used to analyse the motion of probe electrons in the resulting electromagnetic field. It is shown that the higher the ion mobility and the longer the laser pulse, the more stable the plasma layer to longitudinal perturbations of a finite value and the higher the threshold of self-induced transparency. These conclusions are verified by one-dimensional numerical simulation of a complete system of kinetic equations for plasma and Maxwell’s equations. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A multiscale anisotropic polymer network model coupled with phase field fracture

The study of polymers has continued to gain substantial attention due to their expanding range of applications, spanning essential engineering fields to emerging domains like stretchable electronics, soft robotics, and implantable sensors. These materials exhibit remarkable properties, primarily stemming from their intricate polymer chain network, which, in turn, increases the complexity of precisely modeling their behavior. Especially for modeling elastomers and their fracture behavior, accurately accounting for the deformations of the polymer chains is vital for predicting the rupture in highly stretched chains. Despite the importance, many robust multiscale continuum frameworks for modeling elastomer fracture tend to simplify network deformations by assuming uniform behavior among chains in all directions. Recognizing this limitation, our study proposes a multiscale fracture model that accounts for the anisotropic nature of elastomer network responses. At the microscale, damage in the chains is assumed to be driven by both the chain's entropy and the internal energy due to molecular bond distortions. In order to bridge the stretching in the chains to the macroscale deformation, we employ the maximal advance path constraint network model, inherently accommodating anisotropic network responses. As a result, chains oriented differently can be predicted to exhibit varying stretch and, consequently, different damage levels. To drive macroscale fracture based on damages in these chains, we utilize the micromorphic regularization theory, which involves the introduction of dual local-global damage variables at the macroscale. The macroscale local damage variable is obtained through the homogenization of the chain damage values, resulting in the prediction of an isotropic material response. The macroscale global damage variable is subjected to nonlocal effects and boundary conditions in a thermodynamically consistent phase field continuum formulation. Moreover, the total dissipation in the system is considered to be mainly due to the breaking of the molecular bonds at the microscale. To validate our model, we employ the double-edge notched tensile test as a benchmark, comparing simulation predictions with existing experimental data. Additionally, to enhance our understanding of the fracturing process, we conduct uniaxial tensile experiments on a square film made up of polydimethylsiloxane (PDMS) rubber embedded with a hole and notches and then compare our simulation predictions with the experimental observations. Furthermore, we visualize the evolution of stretch and damage values in chains oriented along different directions to assess the predictive capacity of the model. In conclusion, the results are also compared with another existing model to evaluate the utility of our model in accurately simulating the fracture behavior of rubber-like materials.

42 ENGINEERING↗

On observers in holographic maps

A straightforward gravitational path integral calculation implies that closed universes are trivial, described by a one dimensional Hilbert space. Two recent papers by Harlow-Usatyuk-Zhao and Abdalla-Antonini-Iliesiu-Levine have sought to ameliorate this issue by defining special rules to incorporate observers into the path integral. However, the proposed rules are different, leading to differing results for the Hilbert space dimension. Moreover, the former work offers a holographic map realized using a non-isometric code construction to complement their path integral result and clarify its physics. In this work, we propose a non-isometric code that implements the second construction, allowing thorough comparison. Our prescription may be thought of as simply removing the portion of the map that acts on the observer, while preserving the rest, creating an effective holographic boundary at the observer-environment interface. This proposal can be directly applied to general holographic maps for both open and closed universes of any dimension.

AdS-CFT correspondence↗

Effective data sampling strategies and boundary condition constraints of physics-informed neural networks for identifying material properties in solid mechanics

Material identification is critical for understanding the relationship between mechanical properties and the associated mechanical functions. However, material identification is a challenging task, especially when the characteristic of the material is highly nonlinear in nature, as is common in biological tissue. In this work, we identify unknown material properties in continuum solid mechanics via physics-informed neural networks (PINNs). To improve the accuracy and efficiency of PINNs, we develop efficient strategies to nonuniformly sample observational data. We also investigate different approaches to enforce Dirichlet-type boundary conditions (BCs) as soft or hard constraints. Finally, we apply the proposed methods to a diverse set of time-dependent and time-independent solid mechanic examples that span linear elastic and hyperelastic material space. The estimated material parameters achieve relative errors of less than 1%. As such, this work is relevant to diverse applications, including optimizing structural integrity and developing novel materials.

97 MATHEMATICS AND COMPUTING↗

Half-closed discontinuous Galerkin discretisations

Here we introduce the concept of half-closed nodes for nodal discontinuous Galerkin (DG) discretisations. Unlike more commonly used closed nodes in DG, where on every element nodes are placed on all of its boundaries, half-closed nodes only require nodes to be placed on a subset of the element's boundaries. The effect of using different nodes on DG operator sparsity is studied and we find in particular for there to be no difference in the sparsity pattern of the Laplace operator whether closed or half-closed nodes are used. On quadrilateral/hexahedral elements we use the Gauss-Radau points as the half-closed nodes of choice, which we demonstrate is able to speed up DG operator assembly in addition to leverage previously known superconvergence results. We also discuss in this work some linear solver techniques commonly used for Finite Element or discontinuous Galerkin methods such as static condensation and block-based methods, and how they can be applied to half-closed DG discretisations.

97 MATHEMATICS AND COMPUTING↗

Atomistic determination of Peierls barriers of dislocation glide in nickel

The Peierls barrier measures the lattice resistance to dislocation glide in crystalline solids. We use the nudged elastic band (NEB) method to calculate the Peierls barriers for screw and edge dislocation glide in a face-centered cubic (FCC) metal of Ni. The minimum energy paths (MEPs) across single or sequential Peierls barriers are determined under shear loading. The NEB results show the decreasing Peierls barrier with increasing shear stress, giving the Peierls stress at which the Peierls barrier vanishes. The effects of boundary condition and system size on Peierls barriers are studied by comparing strain- and stress-controlled NEB results. Furthermore, the free-end NEB methods are applied to determine MEPs with improved computational efficiency. The NEB results are also used to evaluate the energetic driving force of dislocation glide, which is consistent with that determined from the Peach-Koehler force. The accuracy of the present NEB results based on an empirical interatomic potential is assessed by comparison with a machine-learning potential. This work demonstrates the robust and efficient quantification of Peierls barriers to dislocation glide in an FCC metal, and it lays a solid foundation for the atomistic determination of Peierls barriers in compositionally complex alloys with the FCC structure in future studies.

42 ENGINEERING↗

New scaling law for turbulent boundary layers with high surface mass transfer

A fundamental understanding of the phenomena occurring in the turbulent boundary layer in the presence of surface blowing is limited, and considerable disagreements persist even in describing primary quantities, such as the boundary layer profile. The theories based on the linear boundary layer equations show that the thickness of the sublayer increases in the presence of surface blowing; therefore, the viscous sublayer and law of the wall modify. In this study, direct numerical simulations (DNS) of turbulent boundary layers with uniform surface mass transfers are carried out in order to scale the velocity profile. Emphasis is placed on moderate to high mass transfer rates, which are relevant to the most common hybrid rockets configuration. DNS data are used to establish a functional law of the wall and a law of wake by means of the relation between the wall shear stress and surface mass transfer. Analysis of the mean kinetic energy budget shows that the magnitude of turbulent kinetic energy increases by surface mass transfer, and the production rate extends significantly in the inner layer as the injection rate increases. DNS data of various surface blowing are used to complete the closure of turbulence kinetic energy equation and develop an eddy viscosity model. The predicted turbulent kinetic energy and eddy viscosity agree with DNS data for moderate to high blowing rates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Experimental and modeling studies of the plasma chemistry in a humid Ar radiofrequency atmospheric pressure plasma jet

Abstract While humid atmospheric pressure plasmas are extensively modeled, reaction set validation for these conditions remains limited. We present a detailed comparison of a modelling and experimental study with a goal to elucidate the plasma chemistry in a humid Ar radiofrequency-driven atmospheric pressure plasma jet. A large group of species including radicals (H, OH, O, HO 2 ) and long-lived species (H 2 , O 2 and H 2 O 2 ) in the jet effluent was experimentally quantified by molecular beam mass spectroscopy (MBMS). MBMS measurements of H 2 O 2 , OH and H were validated by direct comparison with a liquid phase colorimetric measurement, laser-induced fluorescence (LIF) and two-photon absorption LIF respectively. While an excellent agreement was found for OH and H 2 O 2 by both techniques, a significant difference was found for H and shown to be due to boundary layer effects at the MBMS sampling substrate. The measured O, OH, HO 2 and H 2 are in good agreement with the plug model while H and O 2 were underestimated and H 2 O 2 was overestimated by the model. The accuracy of both the used reaction set and the diagnostics, as well as the observed discrepancies between the modeling and experimental results, are critically assessed. The results presented in this work enable us to identify further data needs for describing H 2 O vapor chemistry in low-temperature plasmas.

Physics↗

Effect of geometric and magnetic boundary conditions on magnetic islands in 3D force-free ideal MHD equilibria

Three-dimensional boundary features and magnetic perturbations affect the magnetic topology of the force-free ideal MHD equilibria formed within a given volume. The use of the PSI-Tet 3D finite-element equilibrium code reveals the existence of magnetic islands inside such equilibria within boundaries with three-dimensional features in configurations that directly depend on the natures of those features. Here, the nature and limits of the influence of such features are explored, as is the influence of flux boundary conditions to qualitatively identify the structure of magnetic islands in Taylor state equilibria with three-dimensional boundary conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Integrability technique for fluid flow induced deformation of a boundary hair

The deformation of a dense carpet of hair due to Stokes flow in a channel can be described by a nonlinear integrodifferential equation for the shape of a single hair, which possesses several solutions for a given choice of parameters. Although it was posed in a previous study and it bears a resemblance to the pendulum problem from mechanics, this equation has not been analytically solved until now. Despite the presence on an integral with a nonlinear functional dependence on the dependent variable, the system is integrable. We compare the analytically obtained solution to a finite-difference numerical approach, identify the physically realizable solution branch, and briefly study the solution structure through a conserved energylike quantity. Time-dependent fluid-structure interactions are a rich and complex subject to investigate, and we argue that the solution discussed herein can be used as a basis for understanding these systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Numerical Study of Liquid Piston Compression Using Large-Eddy Simulation and Volume-of-Fluid Approach

Efforts to increase the efficiency of residential and commercial air conditioners and heat pumps have demonstrated that the compressor accounts for most of the system’s electrical energy usage. Therefore, the efficiency of this component should be improved to reduce its energy usage. The US Department of Energy’s Oak Ridge National Laboratory developed a near-isothermal liquid piston compressor (LPC) that uses propylene glycol (PG) to compress CO2. This report presents numerical studies of the LPC in which the compression chamber fills with injected PG from the bottom inlet. Numerical simulations were performed using the large-eddy simulation (LES) with the wall-adapting local eddy-viscosity (WALE) subgrid-scale model coupled with the multiphase volume-of-fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfer within the compression chamber. In this effort, effects of boundary conditions applied to the LES-VOF calculations (e.g., no wall, an adiabatic wall, and a wall with a heat flux subscribed) to the overall pressure and temperature of CO2 gas as well as the transient evolution of flow and heat transfer evolution within the compression chamber are investigated and discussed. It was found that the LES calculation with no wall have shown no dynamical flow patterns and the volume-averaged temperature of CO2 increased from 305 to 392.7 K, while LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO2 temperatures. Results of LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO2 gas layers and circulation flow patterns within the chamber that contribute to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor.

Nguyen, Thien D.↗

A Comparison of the Optoelectronic Properties of High-Efficiency Polycrystalline and Epitaxial Cu(In,Ga)Se2 Photovoltaic Films

Over the last several decades, champion photovoltaic (PV) devices using CuInGaSe2 (CIGS) as the absorber material have been achieved using polycrystalline films exclusively. This has led to the assumption that polycrystalline CIGS generally outperform single-crystal CIGS in PV devices. However, recently, very high-quality epitaxial CIGS has been grown on GaAs substrates producing PV device efficiencies of 20.0%. These results have revived the debate on what effects grain boundaries have on PV device efficiencies. In this contribution, we compare the optoelectronic properties of polycrystalline CIGS films with those of high-efficiency epitaxial CIGS films. This comparison reveals that grain boundaries are associated with properties that negatively impact PV device efficiency. Additionally, we find that the grain interiors in polycrystalline films exhibit properties that are similar to the high-performance epitaxial films. Our results suggest that it may be possible to achieve higher device efficiencies with epitaxial CIGS than with polycrystalline films.

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

Rheological Properties of Enzymatically Hydrolyzed Corn Stover Pretreated via Deacetylation and Mechanical Refining

Lignocellulosic biomass is a feedstock for fuels and chemicals that does not compete with food resources and has less contaminants than refuse-derived biomass feedstocks. To convert lignocellulosics to biofuels or value-added products, multiple processing steps are typically necessary. One method of producing biofuels from lignocellulosic biomass utilizes a deacetylation and mechanical refining pretreatment and an enzymatic hydrolysis reaction to produce fermentable sugars from cellulose and hemicellulose. The rheological properties of biomass, such as yield stress and plastic viscosity, change during enzymatic hydrolysis and alter the energy requirements for pumping and mixing, an important consideration for the design of processing equipment. The dynamic changes in rheological properties that occur in a corn stover feedstock undergoing enzymatic hydrolysis are characterized in this work, and the influence on pressure losses in piping systems is estimated. Two rheometer geometries were fabricated with stereolithography 3D printing to reduce wall slip and sample ejection. The slurries have complex rheological behaviors that include shear-thinning behavior. Shear stress ramps were performed on samples at 20 and 50 degrees C using the custom geometries, and the Herschel-Bulkley model was fit to the data. The dynamic nature of the rheological properties is correlated with changes in the average fiber length at various extents of reaction, and the influence of solids concentration on the observed rheology and piping pressure losses is discussed.

09 BIOMASS FUELS↗

Elucidating the mechanical and thermal response of nanotwinned Ni alloys

Transformative advances in nanoscale materials synthesis, characterization and modeling are enabling the synthesis of materials with unprecedent properties and greater understanding of the nanoscale mechanisms that underpin these properties. Nanotwinned Cu alloys have received considerable attention due to their impressive balance of strength and ductility, but they have limited microstructural stability. Recent instantiation of nanotwins in sputter deposited Ni-Mo-W, and several commercial Ni-based superalloys, point to the potential development of a new class of high temperature materials with a very beneficial suite of mechanical and physical properties. The experimental study described in this report was undertaken to elucidate the nanoscale origins of the thermal and mechanical behavior of nanotwinned Ni alloys. Micropillar compression experiments have shown nanotwinned Ni 85 Mo 15-x W x alloys to possess unusually high strengths above 3.5GPa and enhanced microstructural stability. The strength of these nanotwinned alloys is determined by the abrupt formation of shear bands. This study focused on identifying the nanoscale trigger, or triggers, for shear banding in nanotwinned materials using in situ experiments and atomic-scale postmortem analysis. Contrasting and comparing the mechanical response and postmortem nanostructures of “Mo-rich” and “W-rich” nanotwinned specimens elucidated the importance of the lateral motion of easy glide defects that results in erasure of twins and local coarsening of the nanotwinned microstructure. This twin coarsening was then associated with very localized plasticity, strain softening, and shear band formation. The difference between alloys was further studied by considering the role of various material factors: stacking fault energy, coherent twin boundary spacing and flatness, grain size, and the interactions of solute atoms with twin and grain boundaries. The effect of alloy composition and sputtering power and temperature were also investigated and used to control twin spacing and geometry. Combined with atomic-scale simulations, these experiments insights should allow us to identify strategies for achieving concomitant ultrahigh strength and ductility in nanotwinned Ni alloys. In parallel but synergistic studies, our work was buoyed by collaborative state-of-the-art nanoscale orientation and strain mapping at the DOE National Center for Electron Microscopy (NCEM). For example, novel 4D-STEM techniques were used to acquire nanoscale strain maps. The extremely fine twin spacing limited our measurements of atomic-scale thermal expansion within twins, but recent results from NCEM suggest that it will soon be possible to conduct such measurements and to unravel the origins of the novel thermal expansion characteristics of nanotwinned Ni alloys.

36 MATERIALS SCIENCE↗

Liquid Piston with Spray Cooling Near-Isothermal Compressor

The goal of this project was to prototype and characterize the performance of a liquid-piston spray-cooled gas compressor. The working principle of the compressor enables optimized high-efficiency operation over a very wide range of operating conditions, unlike conventional compressors that are optimized for a narrow range of operating conditions. The compressor technology is suitable for many applications, such as gas pipeline transport, gas storage, and commercial and residential heat pumps. Both physical testing and computational fluid dynamics (CFD) modeling of the processes using the Oak Ridge National Laboratory high-performance computing center were completed. The experimental and CFD studies focused on a near-isothermal liquid-piston compressor (LPC) that uses propylene glycol to compress CO 2 . The first prototype demonstrated isothermal operation during several sequentially executed cycles of CO 2 compression and raised the temperature of the compressed CO 2 by only 2 K, compared with approximately 6 K when the gas was compressed non-isothermally. Isothermal operation was demonstrated at CO 2 flow rates of up to 2 L/min. The second prototype was designed with two compression chambers to allow continuous flow of high-pressure CO 2 . However, the design of the valve train to direct flow between the compression chambers was not sufficient to allow demonstration of CO 2 compression. Numerical simulations of the LPC in which the compression chamber was filled with propylene glycol injected from the bottom inlet were performed using large eddy simulation (LES) with the wall-adapting local eddy-viscosity subgrid-scale model coupled with the multiphase volume of fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfers within the compression chamber. In this effort, the effects of boundary conditions applied to the LES-VOF calculations (i.e., no wall, an adiabatic wall, and a wall with a heat flux subscribed) on the overall pressure and temperature of the CO 2 gas as well as the transient evolution of flow and heat transfer within the compression chamber were investigated and are discussed in this report. The LES calculation with no wall showed no dynamical flow patterns, and the volume-averaged temperature of CO 2 increased from 305 to 392.7 K, whereas LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO 2 temperatures. The results of the LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO 2 gas layers and circulation flow patterns within the chamber that contributed to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor. Though an industrial partnership for commercializing the compressor was not secured, the technology was attractive for an industrial partner to use in two research proposals in response to US Department of Energy funding opportunity announcements.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Applications of Nickelate perovskites for neuromorphic computing from electronic structure and Machine Learning

While the limit of Moore's law is presently being reached with current microelectronic technologies, we need to develop new paradigms that overcome this limitation. In that respect, neuromorphic computing is a concept that emulates the neural behavior and response of the human brain, and it has been recognized as a promising alternative approach. In this research project, we will perform multi-fidelity scale bridging to explore the potential use of materials with metal to insulator transition for neuromorphic applications. In particular, rare earth nickelates are promising for such purposes, as the transition in these materials is quite sensitive to a broad set of different external stimuli. Our multi-fidelity approach will bridge the high-fidelity electronic structure calculations with classical potentials. We will bridge dynamical mean field theory with a classical atomistic representation via a deep learning force field. The neural network is trained with energies, charges, and forces obtained by accurate electronic structure theories based on Dynamical Mean Field Theory. The configurational space is generated from known crystal phases, ab initio molecular dynamics with exchange-correlation functionals corrected with the Hubbard model, disordered phases with different concentrations of oxygen vacancies, and nonsymmetrical positions and induced strain by grain interfaces or contact with a substrate. Strategies to train the model with a reduced number of training examples are obtained from active learning methods, and new structures for improving the learning process are generated by using machine learning autoencoders. This classical potential will be validated through a diversity of electronic structure methods and represents an important step to combine the flexibility and accuracy of first-principles with the speed of classical potentials. The generated multi-fidelity surrogate model will be used to understand the role of strain, oxygen vacancies, proton doping, the variation of the crystal phase, substrate effects, vibrational effects as the octahedral rotation, grain boundaries and defect effects on the response of a Metal to Insulator Transition (MIT) in correlated materials. Long time and large-scale simulations will help understand the role of different stimuli to control the hysteresis of the MIT, as it has been experimentally suggested. Selected configurations will be analyzed with higher-level theories to provide an accurate electronic description and to study how the orbitals and charges are rearranged under different conditions.

36 MATERIALS SCIENCE↗

Production of Creep-Resistant FeCrAl-ODS

Four ODS FeCrAl alloys with base alloy composition Fe-16Cr-4Al-2Mo-0.2Y (wt.%) were mechanically alloyed with minor additions of CeO 2 and Y 2 O 3 additions and were successfully extruded into rectangular-shaped bars. Two oxide compositions were explored, which were designated NC with oxide addition of 0.5CeO 2 plus 0.2Yin solution of the base alloy powder and PC with additions of 0.5CeO 2 and 0.7Y 2 O 3 to elevate the O level. The microstructure of both alloys revealed by SEM and EBSD to consist of very small grains with sizes of ~450-520 nm and weak texture components. Sub-size tensile specimens fabricated from both alloys were tensile tested from room temperature to 800ºC. The results for both alloys from room temperature to 400ºC showed very high strengths, but low ductilities, especially uniform strain. From 500ºC to 800ºC, the strengths decreased rapidly, and the ductilities increased significantly. The reason for the rapid decreases in tensile strengths above 500ºC for both alloys is not clear, but will require advanced characterization techniques including atom probe tomography, to determine the size and number density of the nano-size oxide particles, and TEM/EFTEM, for investigating the dispersion of nano-size oxide particles on grain boundaries to assess the Zener grain boundary pinning effect, which will hinder grain coarsening at high temperatures that may enhance the Hall Petch strengthening mechanism.

36 MATERIALS SCIENCE↗