Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite differences”

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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Matching Crystal Structures Atom-to-Atom

Finding an optimal match between two different crystal structures underpins many important materials science problems, including describing solid-solid phase transitions and developing models for interface and grain boundary structures. Here, we formulate the matching of crystals as an optimization problem where the goal is to find the alignment and the atom-to-atom map that minimize a given cost function such as the Euclidean distance between the atoms. We construct an algorithm that directly solves this problem for large finite portions of the crystals and retrieves the periodicity of the match subsequently. We demonstrate its capacity to describe transformation pathways between known polymorphs and to reproduce experimentally realized structures of semi-coherent interfaces. Additionally, from our findings, we define a rigorous metric for measuring distances between crystal structures that can be used to properly quantify their geometric (Euclidean) closeness.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Verification of the kinetic electron role in the microinstabilities in a negative triangularity model equilibrium

Effect of kinetic electrons on negative triangularity plasmas has been investigated and compared against the corresponding positive triangularity plasmas, using the global gyrokinetic code X-point Gyrokinetic Code with scale-separated delta-f option without Coulomb collisions. Our model magnetic equilibria have strong positive and negative triangularities and weak magnetic shear. However, unusually large ρ i > a and low density plasmas are chosen to maximize the nonlocal effect to investigate the finite ρ i effect and to be clearly away from kinetic ballooning modes. Similar conclusions to previous flux tube and global simulations have been obtained in this highly nonlocal model plasma: it is essential to include kinetic electrons in the micro-instability study of negative triangularity plasmas. Most physics findings agree with existing reports, with some disagreement. We offer a new “effective trapping fraction” concept that can add to the explanation of the growth rate difference between NT and PT plasmas, pointing to the significant variation in trapped particle fractions that have turning points in the mode growth regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Finite ion temperature effects on electrostatic instabilities in partially magnetized plasmas

In this paper, the effects of warm ions on electrostatic, kinetic instabilities observed in partially magnetized plasmas, such as the electron cyclotron drift instability (ECDI) and modified two-stream instability (MTSI), are investigated. The kinetic dispersion relation that accounts for warm, non-magnetized ions and warm, magnetized electrons is solved using a complex root-finding technique [A. C. Denig and K. Hara, Phys. Plasmas 30, 032108 (2023)]. It is observed that finite ion temperature introduces ion Landau damping effects, decreasing the growth rate of most unstable modes. The high-frequency, short-wavelength ECDI modes are more strongly affected by ion Landau damping than the low-frequency, long-wavelength MTSI modes. The solutions of the kinetic dispersion relation for MTSI are compared with the fluid dispersion relation of MTSI, which was originally proposed by McBride et al. [Phys. Fluids 15 2367–2383 (1972)]. The differences in the dispersion characteristics between the kinetic and fluid MTSI limits are assessed, illustrating the effects of electron and ion Landau damping on kinetic instabilities in partially magnetized plasmas.

Denig, A. C. (ORCID:0000000253339572)↗
DOE OSTI2023

High Resolution Polar Kerr Effect Studies of CsV 3 Sb 5 : Tests for Time-Reversal Symmetry Breaking below the Charge-Order Transition

We report high resolution polar Kerr effect measurements on CsV 3 Sb 5 single crystals in search for signatures of spontaneous time reversal symmetry breaking below the charge order transition at T * = 94 K. Utilizing two different versions of zero-area loop Sagnac interferometers operating at 1550 nm wavelength, each with the fundamental attribute that without a time reversal symmetry breaking sample at its path, the interferometer is perfectly reciprocal, we find no observable Kerr effect to within the noise floor limit of the apparatus at 30 nanoradians. Simultaneous coherent reflection ratio measurements confirm the sharpness of the charge order transition in the same optical volume as the Kerr measurements. At finite magnetic field we observe a sharp onset of a diamagnetic shift in the Kerr signal at T *, which persists down to the lowest temperature without change in trend. Since 1550 nm is an energy that was shown to capture all features of the optical properties of the material that interact with the charge order transition, we are led to conclude that it is highly unlikely that time reversal symmetry is broken in the charge ordered state in CsV 3 Sb 5 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

In situ Study on the Evolution of Atomic and Electronic Structure of LaTiO 3 /SrTiO 3 System

As with other transition-metal oxide interfaces, LaTiO 3 /SrTiO 3 interfaces exhibit interesting electronic properties, such as a two-dimensional electron gas. While the available Mott states and multiple pathways to metallicity in LaTiO 3 can lead to a variety of applications in oxide electronics, the origin of the different properties remains uncertain. Herein, utilizing real-time, in situ x-ray diffraction at the synchrotron, robust but thin LaTiO 3 /SrTiO 3 heterostructures were fabricated using oxide molecular beam epitaxy. A thickness of six unit cells was determined to be critical, at which point the LaTiO 3 /SrTiO 3 heterostructures begin to sustain sharp, well-defined interfaces. The layers of LaTiO 3 and SrTiO 3 were sufficiently thin to benefit from the finite escape length of electrons at resonant soft x-ray photoemission, and the heterostructures were transferred in vacuo after growth to another chamber for synchrotron-based x-ray photoemission spectroscopy studies. With the x-ray energy tuned to the Ti-2p resonance, the LaTiO 3 /SrTiO 3 interfaces could be probed with enough penetration depth and selectivity. It is shown that all of the heterointerfaces exhibit the two-dimensional electron gas.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The rigorous upscaling of advection-dominated transport in heterogeneous porous media via the Method of Finite Averages

Systems involving advection-dominated transport through heterogeneous porous and fractured media are ubiquitous in subsurface engineering applications. However, upscaling such systems continues to challenge rigorous modeling efforts, particularly when advection is stronger than diffusion at fine spatial scales (i.e., when the Péclet number is greater than one at length scales that characterize a system’s unit-cells, representative elementary volumes, or averaging regions). Here, in this work, we propose and validate a strategy for extending the Method of Finite Averages (MoFA), a rigorous upscaling methodology for heterogeneous porous media, to upscale transport systems experiencing stronger advection than diffusion at fine scales (i.e., fine-scale Péclet numbers greater than one). We detail the strategy, the physical conditions under which it can be applied while retaining a priori modeling error guarantees, and implement the strategy to obtain a MoFA model for advective-diffusive transport that accommodates advective physics at fine spatial scales. We then perform two numerical experiments considering systems with system-scale Péclet numbers of 300 and 1000 — which correspond to fine-scale Péclet numbers of 30 and 100, respectively — to verify that the error guarantees are met under the strategy. After, we conduct a numerical study to demonstrate the strategy’s advantages over the original MoFA methodology. The results suggest that rigorously-upscaled transport models for heterogeneous porous media experiencing advective physics at finer spatial scales can be derived through MoFA and resolved orders of magnitude faster than their pore-scale counterparts. The results also suggest that the presented strategy is limited to modeling shallow concentration gradients when there are large differences between the time scales related to advection and a system’s temporally-varying boundary conditions. This limitation hinders the strategy’s practicality in modeling more advective systems, and as such, opportunity exists for developing additional strategies that accommodate rapidly-varying boundary conditions — and consequentially, steeper concentration gradients — while modeling advective systems with MoFA.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Measurement-driven, model-based estimation of residual stress and its effects on fatigue crack growth. Part 1: Validation of an eigenstrain model

The objective of this paper is to validate a measurement-driven, model-based approach to estimate residual stress (RS) in samples machined from quenched aluminum stock. Model input is derived from measurement of RS in the parent stock. Validation is performed for prismatic T-sections removed from bars at different locations. We find RS predicted agrees with RS measured, by contour and neutron diffraction methods, with root-mean-square model-measurement difference of 22 MPa. Follow-on work (in Part 2) applies the RS estimation to samples representative of aircraft structures and examines the effects of RS on fatigue crack growth in the RS-bearing samples.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A sensitivity analysis of twinning crystal plasticity finite element model using single crystal and poly crystal Zircaloy

The popularity of crystal plasticity finite element method (CPFEM) models is increasing due to their ability to predict the mechanical response of crystalline materials such as metals and metal alloys more accurately than traditional continuum mechanics models. This is since the crystal plasticity models consider the effect of atomic structure, microstructural morphology, and properties of individual grains. These CPFEM models use a large number of material parameters in order to capture the mesoscale physics which comes with the downside of the tedious calibration process. In this paper, a CPFEM code was developed to include the twinning induced grain reorientation and subsequent crystallographic slip for HPC material. The developed code is incorporated in a large-scale, parallelized nonlinear solver WARP3D. Further, a sensitivity analysis with respect to 22 material parameters was then conducted using single crystal and polycrystal representative volume element (RVE) of Zircaloy material. Loading was applied along five different crystallographic orientations for single crystal RVE and along three directions namely, rolling (RD), transverse (TD), and normal (ND) direction for polycrystal RVE. Results obtained from the sensitivity analysis were used for the calibration of material parameters for Zircaloy. Finally, developed code along with calibrated material parameters was used to investigate the effect of the hydride phase formation in Zircaloy which is a typical case observed for nuclear applications. It was found that the volume fraction of the hydride phase has a significant impact on the mechanical properties of Zircaloy.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Positive geometry, local triangulations, and the dual of the Amplituhedron

We initiate the systematic study of local positive spaces which arise in the context of the Amplituhedron construction for scattering amplitudes in planar maximally supersymmetric Yang-Mills theory. We show that all local positive spaces relevant for one-loop MHV amplitudes are characterized by certain sign-flip conditions and are associated with surprisingly simple logarithmic forms. In the maximal sign-flip case they are finite one-loop octagons. Particular combinations of sign-flip spaces can be glued into new local positive geometries. These correspond to local pentagon integrands that appear in the local expansion of the MHV one-loop amplitude. We show that, geometrically, these pentagons do not triangulate the original Amplituhedron space but rather its twin “Amplituhedron-Prime”. This new geometry has the same boundary structure as the Amplituhedron (and therefore the same logarithmic form) but differs in the bulk as a geometric space. On certain two-dimensional boundaries, where the Amplituhedron geometry reduces to a polygon, we check that both spaces map to the same dual polygon. Interestingly, we find that the pentagons internally triangulate that dual space. This gives a direct evidence that the chiral pentagons are natural building blocks for a yet-to-be discovered dual Amplituhedron.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Rupture Jumping and Seismic Complexity in Models of Earthquake Cycles for Fault Stepovers with Off-Fault Plasticity

ABSTRACT Fault stepovers are prime examples of geometric complexity in natural fault zones that may affect seismic hazard by determining whether an earthquake rupture continues propagating or abruptly stops. However, the long-term pattern of seismicity near-fault stepovers and underlying mechanisms of rupture jumping in the context of earthquake cycles are rarely studied. Leveraging a hybrid numerical scheme combining the finite element and the spectral boundary integral methods, FEBE, we carry out fully dynamic simulations of sequences of earthquakes and aseismic slip for both compressive and tensile stepovers with off-fault plasticity. We consider a rate-and-state friction law for the fault friction and pressure-sensitive Drucker–Prager plasticity for the off-fault bulk response. We observe that the accumulation of plastic deformation, an indication of off-fault damage, is significantly different in the two cases, with more plastic deformation projected in the overlapping region for the tensile stepover. The seismic pattern for a tensile stepover is more complex than for a compressive stepover, and incorporating plasticity also increases complexity, relative to the elastic case. A tensile stepover with off-fault plasticity shows rupture segmentation, temporal clustering, and frequent rupture jumping from one fault to another. These results shed light on possible mechanisms of rupture jumping in fault stepovers as well as the long-term evolution of the fault zone.

Geochemistry & Geophysics↗
DOE OSTI2022

EEL: A MOOSE-based Application for 3D Electro-Chemo-Thermo-Mechanical Modeling of Solid-State Batteries

This report describes a tool for three dimensional, high fidelity, coupled electro-chemo-thermo-mechanical modeling of solid-state batteries. A complete version of the tool is available as open source software at https://github.com/hugary1995/eel.git. The theoretical framework of the tool revolves around an inf-sup statement on a total potential, comprising the Helmholtz free energy, the electrical kinetic potential, the chemical potential, the Fourier potential, the chemical reaction potential, and the external power expenditure. The tool uses the finite element framework of Multiphysics Object-Oriented Simulation Environment (MOOSE) and a variational formulation to solve the boundary value problem for solid-state battery incorporating the full set of multiphysics couplings. The variational formulation also enables a modular software architecture for the tool so additional physics can be easily included by specifying the new contribution to the total potential. The report discusses several method of manufactured solutions that were used to verify the implementation of different physics and the Butler-Volmer reaction kinetics at the electrolyte-electrode interface. Finally, the report discusses results from complete charge/discharge simulations of a quasi-1D and a 3D solid-state battery.

25 ENERGY STORAGE↗
DOE OSTI2020

Technology Demonstration of a High-Pressure Swirl Oxy-Coal Combustor

This technical report presents the exploration of the design and prototyping of a High-Pressure Swirl Oxy-coal Combustor. Pressurized oxy-coal combustion systems have the potential to improve efficiency along with an increased carbon capture rate. Reduction of flue gas at higher pressure, smaller system size, and capital cost reductions render high-pressure oxy-coal systems particularly attractive as next-generation energy-producing systems. High-pressure oxy-coal combustion systems are a recent concept, and thus operability issues of combustor designs for such systems are not fully understood. Significant challenges exist to maintain oxy-coal combustion stability at elevated pressure and a high CO 2 diluent environment. Although a body of knowledge exists for high-pressure oxygen combustion in rocket engines (or similar applications), it is yet to be strategized how these fundamental concepts can be translated to low-temperature CO 2 diluent combustion regimes. The realization of the pressurized oxy-coal based systems requires combustor components to be designed and demonstrated for an operating pressure over 10 bar. However, pressurized oxy-coal combustor design information at this pressure range and scale relevant to validate those proposed systems is currently limited. Experimental data from MWth scale oxy-coal combustors are needed to identify the optimal trade-off between net efficiency and systems size. The proposed effort is aimed at demonstrating a 1 MWth down-fired swirl Oxy-Coal combustor and investigate the interrelation between combustor operating conditions (pressure; flame stability; flue gas recirculation ratio) and conversion efficiencies to minimize oxygen requirements. One of the key challenges is to configure burner design (i.e., swirl number and injector) and operating conditions for high-pressure oxy-coal combustion systems. These experiments differ from current systems partly due to the high theoretical flame temperature and related burner operability issues associated with oxy-combustion. An ASPEN PLUS® model study for 550 MWe TIPS and ENEL pressurized oxy-coal systems with CO 2 recirculation was performed to evaluate system design, subsystems sizing, and operating condition determination. The system analysis effort included TRL and technology gap determination of subsystems and critical components. This information was scaled to develop design requirements (design pressure and flue gas recirculation: RR Flue Gas = $\frac{m_{flue}}{m_{total}}$) for the 1 MWth combustor. The effects of a wide range of carbon dioxide recirculation ratios on the thermal efficiency of ENEL and TIPS cycles are studied. The pressure of 10 bar and 80 bar are used for ENEL and TIPS cycles, respectively. The thermal efficiency of ENEL is significantly higher than the efficiency of TIPS at a pressure of less than 10 bar. The insights from system analysis were then used to design a 1 MWth swirl oxy-coal combustor. Flame temperature analysis and material strength analysis was performed to determine the combustor thickness. The structural integrity of the combustor was validated by finite element analysis using Abacus® and Hypermesh®. Feasibility of igniters and secondary burners are investigated in successful high-pressure oxy-methane combustion. The secondary burners are designed in such a way that it can operate between 100 to 500 kW firing input. Three generations of the pintle injector were designed based on swirl numbers (S=0, 0.9, and 1.2). Key pintle injector parameters such as pintle size, pintle orifice size, spray pattern were investigated by cold flow tests. Information from these tests was used to modify injector design for smooth and successful operation. A 5 mm pintle orifice size was decided upon as the optimum size for oxy-coal operation for the combustor. Shadow sizing experiments were performed to identify the atomization rate of each injector. Different coal water slurry mixtures (30 – 50% coal by wt% in the mixture) at various total momentum ratios (TMR) were investigated for this purpose. These experiments provided decisive information to choose the best design of the injector. The injector with 1.2 swirl provided higher atomization in all cases than other designs. The mean equivalent droplet size of the jet was similar at different TMR and mixture ratios using this injector, thus making it suitable for use in most cases. Therefore, the 1.2 swirl-pintle injector was chosen for the shakedown test. The combustor and other sub-systems, including feed systems and control and data acquisition, have been manufactured, assembled, and integrated. The total system integration and installation began on July 1, 2020. The shake-down tests and initial operational capability demonstration are expected to be completed by September 30, 2020.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

Thermomechanical conversion in metals: dislocation plasticity model evaluation of the Taylor-Quinney coefficient

Using a partitioned-energy thermodynamic framework which assigns energy to that of atomic configurational stored energy of cold work and kinetic-vibrational, in this study we derive an important constraint on the Taylor-Quinney coefficient, which quantifies the fraction of plastic work that is converted into heat during plastic deformation. Associated with the two energy contributions are two separate temperatures – the ordinary temperature for the thermal energy and the effective temperature for the configurational energy. We show that the Taylor-Quinney coefficient is a function of the thermodynamically defined effective temperature that measures the atomic configurational disorder in the material. Finite-element analysis of recently published experiments on the aluminum alloy 6016-T4 [1], using the thermodynamic dislocation theory (TDT), shows good agreement between theory and experiment for both stress-strain behavior and temporal evolution of the temperature. The simulations include both conductive and convective thermal energy loss during the experiments, and significant thermal gradients exist within the simulation results. Computed values of the differential Taylor-Quinney coefficient are also presented and suggest a value which differs between materials and increases with increasing strain.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Two-body weak currents in heavy nuclei

In light and medium-mass nuclei, two-body weak currents from chiral effective field theory account for a significant portion of the phenomenological quenching of Gamow-Teller transition matrix elements. Here we examine the systematic effects of two-body axial currents on Gamow-Teller strength and β-decay rates in heavy nuclei within energy-density functional theory. Using a Skyrme functional and the charge-changing finite amplitude method, we add the contributions of two-body currents to the usual one-body linear response in the Gamow-Teller channel, both exactly and though a densitymatrix expansion. The two-body currents, as expected, usually quench both summed Gamow-Teller strength and decay rates, but by an amount that decreases as the neutron excess grows. They can enhance individual low-lying transitions, however, leading to decay rates that are quite different from those that an energy-independent quenching would produce, particularly in neutron-rich nuclei. We show that both these effects are related to changes in the total nucleon density as the number of neutrons increases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Study of recirculating liquid fuel in a 1D critical stationary system

Several studies have been conducted to investigate the physics of liquid fuel reactors, showing also applications with molten salts. A liquid nuclear fuel implies changes in the neutron balance equation to take into account the precursors' displacement and the emission of delayed neutrons in a different position than at the original fission. This requires to upgrade the computer codes normally used to calculate nuclear reactors using only solid fuel. In this work, we revisit a simple problem with liquid fuel, which is proposed for the verification of the numerical solutions obtained by advanced computer codes. This problem studies criticality with constant coefficients, thus neglecting thermal feedback. We elaborate on the analytical solution of the problem, deriving also a generalized eigenvalue problem by finite-volume integration over the cells of a discretized mesh to study the evolution of the dominance ratio with fuel velocity. Finally, we investigate the influence of the fuel velocity on the reactivity of the system. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Finite N indices and the giant graviton expansion

The superconformal index of N = 4 super-Yang Mills theory with U(N) gauge group can be written as a matrix integral over the gauge group. Recently, Murthy demonstrated that this integral can be reexpressed as a sum of terms corresponding to a giant graviton expansion of the index, and provided an explicit formula for the case of a single giant graviton. Here we give similar explicit formulae for an arbitrary number, m ≥ 1, of giant gravitons. We provide 1/2 and 1/16 BPS index examples up to the order where three giant gravitons contribute and demonstrate that the expansion of the matrix integral differs from the giant graviton expansion computed in the supergravity dual. This shows that the giant graviton expansion is not necessarily unique once two or more giant gravitons start appearing.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Numerical Analysis of Fuel Effects on Advanced Compression Ignition Using a Cooperative Fuel Research Engine Computational Fluid Dynamics Model

Growing environmental concerns and demand for a better fuel economy are driving forces that motivate the research for more advanced engines. Multi-mode combustion strategies have gained attention for their potential to provide high thermal efficiency and low emissions for light-duty applications. These strategies target optimizing the engine performance by correlating different combustion modes to load operating conditions. The extension from boosted spark ignition (SI) mode at high loads to advanced compression ignition (ACI) mode at low loads can be achieved by increasing the compression ratio and utilizing intake air heating. Further, in order to enable an accurate control of intake charge condition for ACI mode and rapid mode-switches, it is essential to gain fundamental insights into the autoignition process. Within the scope of ACI, homogeneous charge compression ignition (HCCI) mode is of significant interest. It is known for its potential benefits, operation at low fuel consumption, low NOx, and particulate matter (PM) emissions. In the present work, a virtual Cooperative Fuel Research (CFR) engine model is used to analyze fuel effects on ACI combustion. In particular, the effect of fuel octane sensitivity (S) (at constant Research Octane Number (RON)) on autoignition propensity is assessed under beyond-RON (BRON) and beyond-MON (BMON) ACI conditions. The three-dimensional CFR engine computational fluid dynamics (CFD) model employs a finite-rate chemistry approach with a multi-zone binning strategy to capture autoignition. Two binary blends with Research Octane Number (RON) of 90 are chosen for this study: primary reference fuel (PRF) with S = 0 and toluene–heptane (TH) blend with S = 10.8, representing paraffinic and aromatic gasoline surrogates. Reduced mechanisms for these blends are generated from a detailed gasoline surrogate kinetic mechanism. Simulation results with the reduced mechanisms are validated against experimental data from an in-house CFR engine, with respect to in-cylinder pressure, heat release rate, and combustion phasing. Thereafter, the sensitivity of combustion behavior to ACI operating condition (BRON versus BMON), air-fuel ratio (λ = 2 and 3), and engine speed (600 and 900 rpm) is analyzed for both fuels. It is shown that the sensitivity of a fuel’s autoignition characteristics to λ and engine speed significantly differs at BRON and BMON conditions. Moreover, this sensitivity is found to vary among fuels, despite the same RON. It is also observed that the presence of low-temperature heat release (LTHR) under BRON condition leads to more sequential autoignition and longer combustion duration than BMON condition. Finally, the study indicates that the octane index (OI) fails to capture the trend in the variation of autoignition propensity with S under the BMON condition.

33 ADVANCED PROPULSION SYSTEMS↗