Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “relaxation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Evaluation of Minimum NOx Emission From Ammonia Combustion

Abstract Ammonia (NH3) is being explored as a hydrogen carrier with no carbon emissions. However, if burned directly as NH3, rather than being completely decomposed back to N2/H2, the fuel-bound nitrogen comes with a potentially significant NOx emissions penalty. Indeed, several existing studies are showing ammonia combustion NOx emissions that exceed current natural gas fueled, DLN technologies by one to two orders of magnitude. Therefore, it is important to establish the theoretical minimum NOx emissions for an ammonia combustor, to determine how much NOx levels can be reduced via further technology development. In other words, the purpose of this work is not to analyze the performance of a specific combustor but, rather, the fundamental limits of what is achievable. This study quantifies this minimum NOx level for a two-stage combustor system for a given combustor exit temperature and residence time, with a constraint on unburned fuel levels. As expected, the optimum configuration is a rich front end combustor to burn and crack ammonia with significant H2 production, followed by an NO relaxation reactor, followed by a lean stage that consumes the remaining H2. The optimum residence time and stoichiometry of each zone are determined in the fast mixing limit, which essentially balances between NOx production in the primary and secondary zones. These results show minimum NOx levels are in 200–400 ppm range at 1 bar, but drop to levels of ∼25 ppm at 20 bar. These NOx emissions are dominated by NOx production in the primary stage which relaxes to equilibrium levels quite slowly. As processes controlling NOx relaxation to equilibrium in the primary stage dominate overall NO emission levels, combustor NOx sensitivities are essentially opposite that of natural gas fired, DLN systems. Specifically, NOx values drop with increased combustor residence time, increased pressure, and increased combustor exit temperature. These results also suggest that the most important strategy for NOx minimization is to provide sufficient relaxation time after the primary zone for NOx to approach equilibrium—this can be done via kinetic means to accelerate this relaxation rate, such as enhancing pressure or temperature, or increasing residence times. Indeed, this work shows that low pressure combustors specifically optimized for ammonia will have residence times that are one to two orders of magnitude larger than current natural gas systems. By doing so, NOx levels below 10 ppm may be achievable. Finally, we discuss the sensitivity of these values to uncertainties in ammonia kinetics.

Engineering↗
DOE OSTI2023

Microscopic activated dynamics theory of the shear rheology and stress overshoot in ultradense glass-forming fluids and colloidal suspensions

We formulate a particle and force level, activated dynamics-based statistical mechanical theory for the continuous startup nonlinear shear rheology of ultradense glass-forming hard sphere fluids and colloidal suspensions in the context of the elastically collective nonlinear Langevin equation approach and a generalized Maxwell model constitutive equation. Activated structural relaxation is described as a coupled local-nonlocal event involving caging and longer range collective elasticity which controls the characteristic stress relaxation time. Theoretical predictions for the deformation-induced enhancement of mobility, the onset of relaxation acceleration at remarkably low values of stress, strain, or shear rate, apparent power law thinning of the steady-state structural relaxation time and viscosity, a nonvanishing activation barrier in the shear thinning regime, an apparent Herschel–Buckley form of the shear rate dependence of the steady-state shear stress, exponential growth of different measures of a yield or flow stress with packing fraction, and reduced fragility and dynamic heterogeneity under deformation were previously shown to be in good agreement with experiments. The central new question we address here is the defining feature of the transient response—the stress overshoot. In contrast to the steady-state flow regime, understanding the transient response requires an explicit treatment of the coupled nonequilibrium evolution of structure, elastic modulus, and stress relaxation time. We formulate a new quantitative model for this aspect in a physically motivated and computationally tractable manner. Theoretical predictions for the stress overshoot are shown to be in good agreement with experimental observations in the metastable ultradense regime of hard sphere colloidal suspensions as a function of shear rate and packing fraction, and accounting for deformation-assisted activated motion appears to be crucial for both the transient and steady-state responses.

Mechanics↗
DOE OSTI2021

A Bilevel Approach for Identifying the Worst Contingencies for Nonconvex Alternating Current Power Systems

We address the bilevel optimization problem of identifying the most critical attacks to an alternating current (AC) power flow network. The upper-level binary maximization problem consists of choosing an attack that is treated as a parameter in the lower-level defender minimization problem. Instances of the lower-level global minimization problem by themselves are NP-hard due to the nonconvex AC power flow constraints, and bilevel solution approaches commonly apply a convex relaxation or approximation to allow for tractable bilevel reformulations at the cost of underestimating some power system vulnerabilities. Our main contribution is to provide an alternative branch-and-bound algorithm whose upper bounding mechanism (in a maximization context) is based on a reformulation that avoids relaxation of the AC power flow constraints in the lower-level defender problem. Lower bounding is provided with semidefinite programming (SDP) relaxed solutions to the lower-level problem. We establish finite termination with guarantees of either a globally optimal solution to the original bilevel problem, or a globally optimal solution to the SDP-relaxed bilevel problem which is included in a vetted list of upper-level attack solutions, at least one of which is a globally optimal solution to the bilevel problem. We demonstrate through computational experiments applied to IEEE case instances both the relevance of our contribution, and the effectiveness of our contributed algorithm for identifying power system vulnerabilities without resorting to convex relaxations of the lower-level problem. We conclude with a discussion of future extensions and improvements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Composing preconditioners for multiphysics PDE systems with applications to Generalized MHD

New patch smoothers or relaxation techniques are developed for solving linear matrix equations coming from systems of discretized partial differential equations (PDEs). One key linear solver challenge for many PDE systems arises when the resulting discretization matrix has a near null space that has a large dimension, which can occur in generalized magnetohydrodynamic (GMHD) systems. Patch-based relaxation is highly effective for problems when the null space can be spanned by a basis of locally supported vectors. The patch-based relaxation methods that we develop can be used either within an algebraic multigrid (AMG) hierarchy or as stand-alone preconditioners. These patch-based relaxation techniques are a form of well-known overlapping Schwarz methods where the computational domain is covered with a series of overlapping sub-domains (or patches). Patch relaxation then corresponds to solving a set of independent linear systems associated with each patch. In the context of GMHD, we also reformulate the underlying discrete representation used to generate a suitable set of matrix equations. In general, deriving a discretization that accurately approximates the curl operator and the Hall term while also producing linear systems with physically meaningful near null space properties can be challenging. Unfortunately, many natural discretization choices lead to a near null space that includes non-physical oscillatory modes and where it is not possible to span the near null space with a minimal set of locally supported basis vectors. Further discretization research is needed to understand the resulting trade-offs between accuracy, stability, and ease in solving the associated linear systems.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

Adaptive Inverse Control for Rotorcraft Vibration Reduction

This thesis extends the Least Mean Square (LMS) algorithm to solve the mult!ple-input, multiple-output problem of alleviating N/Rev (revolutions per minute by number of blades) helicopter fuselage vibration by means of adaptive inverse control. A frequency domain locally linear model is used to represent the transfer matrix relating the higher harmonic pitch control inputs to the harmonic vibration outputs to be controlled. By using the inverse matrix as the controller gain matrix, an adaptive inverse regulator is formed to alleviate the N/Rev vibration. The stability and rate of convergence properties of the extended LMS algorithm are discussed. It is shown that the stability ranges for the elements of the stability gain matrix are directly related to the eigenvalues of the vibration signal information matrix for the learning phase, but not for the control phase. The overall conclusion is that the LMS adaptive inverse control method can form a robust vibration control system, but will require some tuning of the input sensor gains, the stability gain matrix, and the amount of control relaxation to be used. The learning curve of the controller during the learning phase is shown to be quantitatively close to that predicted by averaging the learning curves of the normal modes. For higher order transfer matrices, a rough estimate of the inverse is needed to start the algorithm efficiently. The simulation results indicate that the factor which most influences LMS adaptive inverse control is the product of the control relaxation and the the stability gain matrix. A small stability gain matrix makes the controller less sensitive to relaxation selection, and permits faster and more stable vibration reduction, than by choosing the stability gain matrix large and the control relaxation term small. It is shown that the best selections of the stability gain matrix elements and the amount of control relaxation is basically a compromise between slow, stable convergence and fast convergence with increased possibility of unstable identification. In the simulation studies, the LMS adaptive inverse control algorithm is shown to be capable of adapting the inverse (controller) matrix to track changes in the flight conditions. The algorithm converges quickly for moderate disturbances, while taking longer for larger disturbances. Perfect knowledge of the inverse matrix is not required for good control of the N/Rev vibration. However it is shown that measurement noise will prevent the LMS adaptive inverse control technique from controlling the vibration, unless the signal averaging method presented is incorporated into the algorithm.

Jacklin, Stephen A.↗
NASA NTRS1991

Viscoelastic properties of addition-cured polyimides used in high temperature polymer matrix composites

Viscoelastic properties of the addition cured polyimide, PMR-15, were studied using dynamic mechanical and stress relaxation tests. For temperatures below the glass transition temperature, T sub g, the dynamic mechanical properties measured using a temperature scan rate of 10 C/min were strongly affected by the presence of absorbed moisture in the resin. Dynamic mechanical properties measured as a function of time during an isothermal hold provided an indication of chemical changes occurring in the resin. For temperatures above (T sub g + 20 C), the storage modulus increased continuously as a function of time indicating that additional crosslinking is occurring in the resin. Because of these changes in chemical structures, the stress relaxation modulus could not be measured over any useful time interval for temperatures above T sub g. For temperatures below T sub g, dynamic mechanical properties appeared to be unaffected by chemical changes for times exceeding 1 hr. Since the duration of the stress relaxation tests was less than 1 hr, the stress relaxation modulus could be measured. As long as the moisture content of the resin was less than 2 pct, stress relaxation curves measured at different temperatures could be superimposed using horizontal shifts along the log(time) axis with only small shifts along the vertical axis.

Roberts, Gary D↗
NASA NTRS1993

Mapping robust parallel multigrid algorithms to scalable memory architectures

The convergence rate of standard multigrid algorithms degenerates on problems with stretched grids or anisotropic operators. The usual cure for this is the use of line or plane relaxation. However, multigrid algorithms based on line and plane relaxation have limited and awkward parallelism and are quite difficult to map effectively to highly parallel architectures. Newer multigrid algorithms that overcome anisotropy through the use of multiple coarse grids rather than line relaxation are better suited to massively parallel architectures because they require only simple point-relaxation smoothers. The parallel implementation of a V-cycle multiple semi-coarsened grid (MSG) algorithm or distributed-memory architectures such as the Intel iPSC/860 and Paragon computers is addressed. The MSG algorithms provide two levels of parallelism: parallelism within the relaxation or interpolation on each grid and across the grids on each multigrid level. Both levels of parallelism must be exploited to map these algorithms effectively to parallel architectures. A mapping of an MSG algorithm to distributed-memory architectures that demonstrate how both levels of parallelism can be exploited is described. The results is a robust and effective multigrid algorithm for distributed-memory machines.

Overman, Andrea↗
NASA NTRS1993

Mapping robust parallel multigrid algorithms to scalable memory architectures

The convergence rate of standard multigrid algorithms degenerates on problems with stretched grids or anisotropic operators. The usual cure for this is the use of line or plane relaxation. However, multigrid algorithms based on line and plane relaxation have limited and awkward parallelism and are quite difficult to map effectively to highly parallel architectures. Newer multigrid algorithms that overcome anisotropy through the use of multiple coarse grids rather than relaxation are better suited to massively parallel architectures because they require only simple point-relaxation smoothers. In this paper, we look at the parallel implementation of a V-cycle multiple semicoarsened grid (MSG) algorithm on distributed-memory architectures such as the Intel iPSC/860 and Paragon computers. The MSG algorithms provide two levels of parallelism: parallelism within the relaxation or interpolation on each grid and across the grids on each multigrid level. Both levels of parallelism must be exploited to map these algorithms effectively to parallel architectures. This paper describes a mapping of an MSG algorithm to distributed-memory architectures that demonstrates how both levels of parallelism can be exploited. The result is a robust and effective multigrid algorithm for distributed-memory machines.

Overman, Andrea↗
NASA NTRS1994

New constitutive equation for the volume viscosity in fluids

The traditional volume viscosity, Stokes' hypothesis, and acoustical relaxation are reviewed. The lossy Navier-Stokes Equation is applied to periodic (acoustic) flow, and it is shown that the traditional volume viscosity leads to a result which contradicts that describing acoustical relaxation. It is demonstrated that the addition of a second volume viscosity term to the constitutive equation, to account for pressure relaxation, resolves the conflict, and leads to a direct correspondence between the volume viscosity parameters and the acoustical relaxation parameters. The representation of volume viscosity is formulated for the case of multiple relaxations, as occur in air. Finally, an application of the new constitutive equation to a simple convective compressible flow, namely a linearly accelerating flow, demonstrates the impact of volume viscosity upon the flow and the physical conditions for which it is important.

Zuckerwar, Allan J.↗
NASA NTRS1994

CO2 cooling in terrestrial planet thermospheres

We examine the recent progress in the debate on the CO2-O relaxation rate, its temperature dependence, and its corresponding impact on the thermospheric heat budgets of Venus, Earth, and Mars. This comparative approach provides the broadest range of conditions under which a common CO2-O relaxation rate should provide consistent results. New global mean calculations are presented for the heat budgets of these three planets using large CO2-O relaxation rates that have been inferred recently from Earth CO2 radiance measurements and laboratory studies. Results indicate that available Venus and Mars data constrain the CO2-O relaxation rate to be 2-4 x 10(exp -12)/cu cm/s at 300 K. For Venus, this strong cooling serves as an effective thermostat that gives rise to a small variation of thermospheric temperatures over the solar cycle, just as observed. Conversely, CO2 cooling does not appear to be dominant in the dayside heat budget of the Mars thermosphere over most of the solar cycle. For the Earth, this strong cooling implies that the lower thermosphere does not typically require significant eddy diffusion or heat conduction. However, global-scale dynamics or an additional heating mechanism may be needed to restore calculated temperatures to observed values when relaxation rates exceeding 2 x 10(exp -12)/cu cm/s are employed.

Bougher, S. W.↗
NASA NTRS1999

Inspection of Space Station Cold Plate Using Visual and Automated Holographic Techniques

Real-time holography has been used to confirm the presence of non-uniformity in the construction of an International Space Station cold plate. Ultrasonic C-scans have previously shown suspected areas of cooling fin disbonds. But both neural-net processed and visual holography did not evidence any progressive permanent changes resulting from 3000 pressurization and relaxation cycles of a Dash 8 cold plate. Neural-net and visual inspections were performed of characteristic patterns generated from electronic time-average holograms of the vibrating cold plate. Normal modes of vibration were excited at very low amplitudes for this purpose, The neural nets were trained to flag very small changes in the mode shapes as encoded in the characteristic patterns. Both the whole cold plate and a zoomed region were inspected. The inspections were conducted before, after, and during pressurization and relaxation cycles of the cold plate. A water-filled cold plate was pressurized to 120 psig (827 kPa) and relaxed for each cycle. Each cycle required 5 seconds. Both the artificial neural networks and the inspectors were unable to detect changes in the mode shapes of the relaxed cold plate. The cold plate was also inspected visually using real-time holography and double-exposure holography. Regions of non-uniformity correlating with the C-scans were apparent, but the interference patterns did not change after 3000 pressurization and relaxation cycles. These tests constituted the first practical application of a neural-net inspection technique developed originally with support from the Director's Discretionary Fund at the Glenn Research Center at Lewis Field.

Decker, Arthur J.↗
NASA NTRS2002

A Thermodynamic Theory Of Solid Viscoelasticity. Part 1: Linear Viscoelasticity.

The present series of three consecutive papers develops a general theory for linear and finite solid viscoelasticity. Because the most important object for nonlinear studies are rubber-like materials, the general approach is specified in a form convenient for solving problems important for many industries that involve rubber-like materials. General linear and nonlinear theories for non-isothermal deformations of viscoelastic solids are developed based on the quasi-linear approach of non-equilibrium thermodynamics. In this, the first paper of the series, we analyze non-isothermal linear viscoelasticity, which is applicable in a range of small strains not only to all synthetic polymers and bio-polymers but also to some non-polymeric materials. Although the linear case seems to be well developed, there still are some reasons to implement a thermodynamic derivation of constitutive equations for solid-like, non-isothermal, linear viscoelasticity. The most important is the thermodynamic modeling of thermo-rheological complexity , i.e. different temperature dependences of relaxation parameters in various parts of relaxation spectrum. A special structure of interaction matrices is established for different physical mechanisms contributed to the normal relaxation modes. This structure seems to be in accord with observations, and creates a simple mathematical framework for both continuum and molecular theories of the thermo-rheological complex relaxation phenomena. Finally, a unified approach is briefly discussed that, in principle, allows combining both the long time (discrete) and short time (continuous) descriptions of relaxation behaviors for polymers in the rubbery and glassy regions.

Freed, Alan D.↗
NASA NTRS1993

Evaluation of muscle injury using magnetic resonance imaging

The objective of this study was to investigate spin echo T2 relaxation time changes in thigh muscles after intense eccentric exercise in healthy men. Spin echo and calculated T2 relaxation time images of the thighs were obtained on several occasions after exercise of one limb; the contralateral limb served as control. Muscle damage was verified by elevated levels of serum creatine kinase (CK). Thirty percent of the time no exercise effect was discernible on the magnetic resonance (MR) images. In all positive MR images (70%) the semitendinosus muscle was positive, while the biceps femoris, short head, and gracilis muscles were also positive in 50% and 25% of the total cases, respectively. The peak T2 relaxation time and serum CK were correlated (r = 0.94, p<0.01); temporal changes in muscle T2 relaxation time and serum CK were similar, although T2 relaxation time remained positive after serum CK returned to background levels. We conclude that magnetic resonance imaging can serve as a useful tool in the evaluation of eccentric exercise muscle damage by providing a quantitative indicator of damage and its resolution as well as the specific areas and muscles.

NASA Discipline Number 26-10↗
NASA NTRS2011

Cognitive Behavioral Training and Education for Spaceflight Operations

Cognitive behavioral-training (CBT) is an evidence-based practice commonly used to help treat insomnia, and is part of NASA's countermeasure regimen for Fatigue Management. CBT addresses the life style and habits of individuals that are maladaptive to managing stress and fatigue. This includes addressing learned behavioral responses that may cause stress and lead to an increased sense of fatigue. While the initial cause of onset of fatigue in the individual may be no longer present, the perception and engrained anticipation of fatigue persist and cause an exaggerated state of tension. CBT combined with relaxation training allows the individual to unlearn the maladaptive beliefs and behaviors and replace them with routines and techniques that allow cognitive restructuring and resultant relief from stress. CBT allows for elimination in individuals of unwanted ruminating thoughts and anticipatory anxiety by, for example, training the individuals to practice stressful situations in a relaxed state. As a result of CBT, relaxation can be accomplished in many ways, such as progressive muscle relaxation, meditation and guided imagery. CBT is not therapy, but rather the synthesis of behavioral countermeasures. CBT utilizes progressive relaxation as a means of reinforcing educational and cognitive countermeasures. These countermeasures include: masking, elimination of distracting thoughts, anxiety control, split attention, cognitive restructuring and other advanced psychological techniques.

Moonmaw, Ronald↗
DOE OSTI2024

Instrumentation for the In-Core Real-Time Mechanical Testing of Structural Materials (INCREASE) Project

Idaho National Laboratory (INL), in collaboration with the Electric Power Research Institute (EPRI), the Nuclear Regulatory Commission (NRC), the French Atomic and Alternative Energies Commission (CEA), the Joint Research Center (JRC), the Nuclear Research and Consultancy Group (NRG), and the Research Center Rez (CVR), started a Joint Experimental Program (JEEP) project that operates within the Nuclear Energy Agency’s Framework for Irradiation Experiments (FIDES II) program in order to develop capabilities for the in-core real-time mechanical testing of structural materials. This effort will focus on designing a shared capsule capable of housing a variety of in-core mechanical testing instrumentation allowing enhanced experiments for the material science community. The outcome of the project would be high-priority, stress relaxation data for stainless-steel-based materials provided by EPRI and CEA. Stress relaxation is a major phenomenon that contributes to material degradation in nuclear reactor components. Currently, nuclear material stress relaxation is assessed both before and after irradiation, using complex and costly post-irradiation examination (PIE) activities. In-situ data would support the development of precision modeling and simulation of this degradation phenomena and would provide validation and benchmarking for existing models using the PIE data. As part of the U.S. Department of Energy Advanced Sensor and Instrumentation (ASI) program, INL has fabricated and tested out-of-core mechanical test instrumentation. This instrumentation was designed for easy adaptation to the irradiation capsule proposed under this JEEP, and can be deployed to measure real-time stress relaxation under pressurized-water reactor (PWR) conditions. This initial effort will partially serve to replace the testing capabilities lost because of shutting down the Halden Boiling Water Reactor (HBWR). The project would provide these capabilities to the international community via a shared capsule design that is easily adaptable to additional material test reactors. The design features will incorporate expansion to PWR and non-light-water reactor (LWR) environments that will be developed in future work. The capsule and instrumentation will be demonstrated in the Massachusetts Institute of Technology Reactor (MITR) for phase I and the Petten High Flux Reactor (HFR) for phase II irradiations to deliver real-time stress relaxation data on high priority stainless steel structural materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Unveiling the Structural Origins of Dynamic Diversity in Pd-Based Metallic Glasses

The β -relaxation is one of the major dynamic behaviors in metallic glasses (MGs) and exhibits diverse features. Despite decades of efforts, the understanding of its structural origin and contribution to the overall dynamics of MG systems is still unclear. Here two palladium-based Pd—Cu—P and Pd—Ni—P MGs are reported with distinct different β -relaxation behaviors and reveal the structural origins for the difference using the advanced X-ray photon correlation spectroscopy and absorption fine structure techniques together with the first-principles calculations. The pronounced β -relaxation and fast atomic dynamics in the Pd—Cu—P MG mainly come from the strong mobility of Cu atoms and their locally favored structures. In contrast, the motion of Ni atoms is constrained by P atoms in the Pd—Ni—P MG, leading to the weakened β -relaxation peak and sluggish dynamics. Finally, the correlation of atomic dynamics with microscopic structures provides a way to understand the structural origins of different dynamic behaviors as well as the nature of aging in disordered materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A simple and effective five-equation two-phase numerical model for liquid-vapor phase transition in cavitating flows

Numerical difficulties, notably the non-monotonic behavior of the Wood speed of sound and the volume fraction positivity, associated with the reduced five-equation two-phase flow model of Kapila et al. (2001) [A.K. Kapila, R. Menikoff, J.B. Bdzil, S.F. Son, D.S. Stewart, 2001. Two-phase modeling of deflagration-to-detonation transition in granular materials: reduced equations, Physics of Fluids 13(10), 3002–3024] have been resolved in the past through the introduction of a frozen speed of sound and an algebraic approach for mechanical relaxation afforded by a pressure non-equilibrium six-equation model proposed by [R. Saurel, F. Petitpas, R.A. Berry, 2009. Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixture, J. Comput. Phys. 228, 1678–1712]. By contrast, it is explored and demonstrated in this work that these difficulties can in fact still be resolved within the numerical scheme for solving the reduced five-equation model by numerically replacing the Wood speed of sound for the estimates of wave speeds in the approximate Riemann solver HLLC with the monotonic mixture speed of sound for a transport five-equation model. For shock interface (artificial mixture separating pure or nearly pure fluids) interaction problems, with the apparent advantage of monotonic behavior of the speed of sound in the interface, the effect of the numerical replacement is also confined to the interface. Differences other than the behavior of the speed of sound within the interface in the solutions due to the replacement diminish with increasing resolution when reasonable solution can be obtained with Wood speed of sound. For cavitating/expansion problems in physical fluid mixture, it is pointed out and explained why acoustics in the numerical solutions still propagate at the Wood speed of sound (therefore consistent with the reduced five-equation model) even though in some cases a much higher speed of sound like the numerical replacement above for solving the reduced five-equation model or the frozen speed of sound for solving a six-equation model is used for the estimates of wave speeds in the HLLC scheme. A variant of the five-equation two-phase flow model by Saurel et al. (2008) [R. Saurel, F. Petitpas, R. Abgrall, 2008. Modelling phase transition in metastable liquids: application to cavitating and flashing flows, J. Fluid Mech. 607, 313–350] is then constructed for liquid-vapor phase transition in cavitating flows. The relaxation toward thermo-chemical equilibrium during phase transition is achieved by solving a simple system of algebraic equations for the equilibrium state variables for better efficiency, following Pelanti and Shyue (2014) [M. Pelanti, K.-M. Shyue, 2014. A mixture-energy-consistent six-equation two-phase numerical model for fluid with interfaces, cavitation and evaporation waves. J. Comput. Phys. 259, 331–357]. Therefore, the current model retains both the simplicity afforded by the five-equation model and the efficiency of the algebraic relaxation solver. An alternative algebraic approach for handling the non-conservative term (the so-called K∇ · u term) in the reduced five-equation model for mechanical equilibrium of a liquid-vapor mixture is also explored by enforcing the thermal equilibrium at the same time. Finally, numerical results of sample tests in both one and two dimensions in the literature as well as that in three dimensions demonstrate the effectiveness and ability of the proposed model to simulate cavitating flows. An interesting mechanism of shock generation by acoustics in water due to phase transition is then found by the numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Bicrystallography-informed Frenkel–Kontorova model for interlayer dislocations in strained 2D heterostructures

In recent years, van der Waals (vdW) heterostructures and homostructures, which consist of stacks of two-dimensional (2D) materials, have risen to prominence due to their association with exotic quantum phenomena originating from correlated electronic states harbored by them. Atomistic scale relaxation effects play an extremely important role in the electronic scale quantum physics of these systems, providing means of manipulation of these materials and allowing them to be tailored for emergent technologies. We investigate such structural relaxation effects in this work using atomistic and mesoscale models, within the context of twisted bilayer graphene — a well-known heterostructure system that features moiré patterns arising from the lattices of the two graphene layers. For small twist angles, atomic relaxation effects in this system are associated with the natural emergence of interface dislocations or strain solitons, which result from the cyclic nature of the generalized stacking fault energy (GSFE), that measures the interface energy based on the relative movement of the two layers. Here, in this work, we first demonstrate using atomistic simulations that atomic reconstruction in bilayer graphene under a large twist also results from interface dislocations, although the Burgers vectors of such dislocations are considerably smaller than those observed in small-twist systems. To reveal the translational invariance of the heterointerface responsible for the formation of such dislocations, we derive the translational symmetry of the GSFE of a 2D heterostructure using the notions of coincident site lattices (CSLs) and displacement shift complete lattices (DSCLs). The workhorse for this exercise is a recently developed Smith normal form bicrystallography framework. Next, we construct a bicrystallography-informed and frame-invariant Frenkel–Kontorova model, which can predict the formation of strain solitons in arbitrary 2D heterostructures, and apply it to study a heterostrained, large-twist bilayer graphene system. Our mesoscale model is found to produce results consistent with atomistic simulations. We anticipate that the model will be invaluable in predicting structural relaxation and for providing insights into various heterostructure systems, especially in cases where the fundamental unit cell is large and therefore, atomistic simulations are computationally expensive.

2D heterostructures↗