Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FORMULATION”

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

Modeling the role of local crystallographic correlations in microstructures of Ti-6Al-4V using a correlated structure visco-plastic self-consistent polycrystal plasticity formulation

This paper presents a multi-level crystal plasticity-based simulation framework for modeling mechanical response and microstructure evolution of Ti-6Al-4V with α-lath/lamellar microstructures. Here, the model is a correlated structure visco-plastic self-consistent (CS-VPSC) formulation linking three scales: a single crystals micro-scale, a lath/lamellar colony meso-scale, and a lath/lamellar aggregate macro-scale. A selected hardening law for the evolution of critical resolved shear stress per slip system used in CS-VPSC is phenomenological. However, it adjusts the resistances of basal and prismatic slip systems based on the geometry of slip transfer between adjacent lamellae. Consistent with experimental evidences, the resolved shear stress on the pyramidal slip planes is dependent not only on the stress in the direction of slip but also on the two orthogonal shear stress components and the three normal stress components (non-Schmid effects). Electron backscatter diffraction (EBSD) data in conjunction with a procedure relying on α→β phase transformation is used to construct paired variants of α-laths/lamellae satisfying their local crystallographic correlations. The procedure fits volume fractions of individual laths/lamellae with the experimental EBSD data and selects a distribution of habit planes between adjacent variants with respect to the loading direction. The simulation framework is applied to interpret the deformation behavior in tension and compression along two sample directions of Ti-6Al-4V fabricated via laser powder bed fusion. Moreover, the model is used to simulate texture evolution during rolling of the material to large strains. It is demonstrated that the model is capable of predicting plastic anisotropy/asymmetry and the concomitant texture evolution. While the model reveals a significant effect of habit plane inclination with respect to the loading direction on yield stress, the comparison of the data and model predictions shows that a random distribution of habit planes fits the flow response. It is further inferred that the tension-compression asymmetry arises from the non-Schmid effects.

36 MATERIALS SCIENCE↗
DOE OSTI2014

A dynamical formulation of one-dimensional scattering theory and its applications in optics

We develop a dynamical formulation of one-dimensional scattering theory where the reflection and transmission amplitudes for a general, possibly complex and energy-dependent, scattering potential are given as solutions of a set of dynamical equations. By decoupling and partially integrating these equations, we reduce the scattering problem to a second order linear differential equation with universal initial conditions that is equivalent to an initial-value time-independent Schrödinger equation. We give explicit formulas for the reflection and transmission amplitudes in terms of the solution of either of these equations and employ them to outline an inverse-scattering method for constructing finite-range potentials with desirable scattering properties at any prescribed wavelength. In particular, we construct optical potentials displaying threshold lasing, antilasing, and unidirectional invisibility. -- Highlights: • Proposes a dynamical theory of scattering in one dimension. • Derives and solves dynamical equations for scattering data. • Gives a new inverse scattering prescription. • Constructs optical potentials with desired scattering properties.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Fierz bilinear formulation of the Maxwell–Dirac equations and symmetry reductions

We study the Maxwell–Dirac equations in a manifestly gauge invariant presentation using only the spinor bilinear scalar and pseudoscalar densities, and the vector and pseudovector currents, together with their quadratic Fierz relations. The internally produced vector potential is expressed via algebraic manipulation of the Dirac equation, as a rational function of the Fierz bilinears and first derivatives (valid on the support of the scalar density), which allows a gauge invariant vector potential to be defined. This leads to a Fierz bilinear formulation of the Maxwell tensor and of the Maxwell–Dirac equations, without any reference to gauge dependent quantities. We show how demanding invariance of tensor fields under the action of a fixed (but arbitrary) Lie subgroup of the Poincaré group leads to symmetry reduced equations. The procedure is illustrated, and the reduced equations worked out explicitly for standard spherical and cylindrical cases, which are coupled third order nonlinear PDEs. Spherical symmetry necessitates the existence of magnetic monopoles, which do not affect the coupled Maxwell–Dirac system due to magnetic terms cancelling. In this paper we do not take up numerical computations. As a demonstration of the power of our approach, we also work out the symmetry reduced equations for two distinct classes of dimension 4 one-parameter families of Poincaré subgroups, one splitting and one non-splitting. The splitting class yields no solutions, whereas for the non-splitting class we find a family of formal exact solutions in closed form. - Highlights: • Maxwell–Dirac equations derived in manifestly gauge invariant tensor form. • Invariant scalar and four vector fields for four Poincaré subgroups derived, including two unusual cases. • Symmetry reduction imposed on Maxwell–Dirac equations under example subgroups. • Magnetic monopole arises for spherically symmetric case, consistent with charge quantization condition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Beyond the Shannon–Khinchin formulation: The composability axiom and the universal-group entropy

The notion of entropy is ubiquitous both in natural and social sciences. In the last two decades, a considerable effort has been devoted to the study of new entropic forms, which generalize the standard Boltzmann–Gibbs (BG) entropy and could be applicable in thermodynamics, quantum mechanics and information theory. In Khinchin (1957), by extending previous ideas of Shannon (1948) and Shannon and Weaver (1949), Khinchin proposed a characterization of the BG entropy, based on four requirements, nowadays known as the Shannon–Khinchin (SK) axioms. The purpose of this paper is twofold. First, we show that there exists an intrinsic group-theoretical structure behind the notion of entropy. It comes from the requirement of composability of an entropy with respect to the union of two statistically independent systems, that we propose in an axiomatic formulation. Second, we show that there exists a simple universal family of trace-form entropies. This class contains many well known examples of entropies and infinitely many new ones, a priori multi-parametric. Due to its specific relation with Lazard’s universal formal group of algebraic topology, the new general entropy introduced in this work will be called the universal-group entropy. A new example of multi-parametric entropy is explicitly constructed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Construction of distinct discrete time scattering quantum walk formulations on the honeycomb lattice

Here we address the scattering construction of discrete time quantum walks on the honeycomb lattice. We write the system general (unitary) one step time evolution operator in terms of: (i) complete arbitrary scattering matrices Γ{sup (j,k)}, defined on the sites (j,k) of the lattice and; and (ii) topological directional functions Φ, which represent the distinct ways we can keep track of the propagation directions along the three bonds attached to each (j,k). By imposing the Φ’s to comply with all the honeycomb translational and point group symmetries (in the case of Γ{sup (j,k)}=Γ), we obtain in total ten independent model versions, presenting different dynamical features. To study some of their traits, we consider the idea of ‘characteristic paths’ CPs (closely related to classical random walks), determining the CPs for each one of the ten formulations. We then discuss many numerical examples of time evolution dynamics for our scattering quantum walks. For the calculations we use few Γ’s, as the Discrete Fourier Transform (DFT) and Grover matrices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

First order formulation of the Yang–Mills theory in a background field

The background gauge renormalization of the first order formulation of the Yang–Mills theory is studied by using the BRST identities. Together with the background symmetry, these identities allow for an iterative proof of renormalizability to all orders in perturbation theory. However, due to the fact that certain improper diagrams which violate the BRST symmetry should be removed, the renormalizability must be deduced indirectly. The recursive method involves rescalings and mixings of the fields, which lead to a renormalized effective action for the background field theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Helicity, spin, and infra-zilch of light: A Lorentz covariant formulation

Highlights: • Spin and orbital parts of angular momentum. • The helicity tensor contains the helicity, spin and infra-zilch currents. • The helicity tensor is a Noether current for a duality-symmetric action. • Symmetry generator for the helicity tensor. • The helicity tensor is conserved in Lorenz gauge. In this paper, a novel conserved Lorentz covariant tensor, termed the helicity tensor, is introduced in Maxwell theory. The conservation of the helicity tensor expresses the conservation laws contained in the helicity array, introduced by Cameron et al. (2012), including helicity, spin, as well as the spin-flux or infra-zilch. The Lorentz covariance of the helicity tensor is in contrast to previous formulations of the helicity hierarchy of conservation laws, which required the non-Lorentz covariant transverse gauge. The helicity tensor is shown to arise as a Noether current for a variational symmetry of a duality-symmetric Lagrangian for Maxwell theory. This symmetry transformation generalizes the duality symmetry and includes the symmetry underlying the conservation of the spin part of the angular momentum.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Effect of Aging on the NOx Storage and Regeneration Characteristics of Fully Formulated Lean NOx Trap Catalysts

In order to elucidate the effect of washcoat composition on lean NO{sub x} trap (LNT) aging characteristics, fully formulated monolithic LNT catalysts containing varying amounts of Pt, Rh and BaO were subjected to accelerated aging on a bench reactor. Subsequent catalyst evaluation revealed that in all cases aging resulted in deterioration of the NO{sub x} conversion as a consequence of impaired NO{sub x} storage and NO{sub x} reduction functions, while increased selectivity to NH{sub 3} was observed in the temperature range 250--450 C. Elemental analysis, H{sub 2} chemisorption and TEM data revealed two main changes which account for the degradation in LNT performance. First, residual sulfur in the catalysts, associated with the Ba phase, decreased catalyst NO{sub x} storage capacity. Second, sintering of the precious metals in the washcoat occurred, resulting in decreased contact between the Pt and Ba, and hence in less efficient NO{sub x} spillover from Pt to Ba during NO{sub x} adsorption, as well as decreased rates of reductant spillover from Pt to Ba and reverse NO{sub x} spillover during catalyst regeneration. For the aged catalysts, halving the Pt loading from 100 to 50 g/ft{sup 3} was found to result in a significant decrease in overall NO{sub x} conversion, while for catalysts with the same 100 g/ft{sup 3} Pt loading, increasing the relative amount of Pt on the NO{sub x} storage components (BaO and La-stabilized CeO{sub 2}), as opposed to an Al{sub 2}O{sub 3} support material (where it was co-located with Rh), was found to be beneficial. The effect of Rh loading on aged catalyst performance was found to be marginal within the range studied (10--20 g/ft{sup 3}), as was the effect of BaO loading in the range 30--45 g/L.

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

Formulation of low clinker blended cements and concrete with enhanced fresh and hardened properties

Blended cements appear as the most promising way to massively produce cements with reduced environmental footprints. However, considering the generally lower early mechanical performance of highly blended cements, chemical admixtures are increasingly relied upon to increase early strength. However, in many cases this can negatively impact other properties as rheology and durability. This paper proposes an approach to formulate low clinker cements and concrete using moderate alkali activation or gypsum to enhance early and long-term compressive strength, without inducing excessive alkali-silica reaction (ASR) expansion. It examines blended cements containing burnt oil shale along with limestone and fly ash or slag as supplementary cementitious materials. It optimizes the combination of polycarboxylate ether (PCE) superplasticizers with the strength enhancers. Along with the use of Ca(NO{sub 3}){sub 2} to drive PCE adsorption, this provides good rheology control, increased early and long-term compressive strength, as well as acceptable ASR expansion in cements with only 50% clinker.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A data-driven linear formulation of the optimal demand response scheduling problem for an industrial air separation unit

Demand response (DR) has become a key element in balancing the power grid as the contribution of time-varying renewable power generation increases. Chemical plants are appealing candidates for DR programs as they offer large, concentrated and flexible loads. DR participation calls for frequent production rate changes over time scales that overlap with the dominant dynamics of the plant. Production scheduling should therefore consider the process dynamics explicitly. Here we present a data-driven approach for modelling the scheduling-relevant dynamics based on historical closed-loop operating data using autoregressive with extra inputs (ARX) models. We introduce a new, linear scheduling problem formulation based on the ARX representation, and demonstrate its implementation on an industrial air separation unit.

42 ENGINEERING↗
DOE OSTI2019

A phase-field formulation for dynamic cohesive fracture

We extend a phase-field/gradient damage formulation for cohesive fracture to the dynamic case. The model is characterized by a regularized fracture energy that is linear in the damage field, as well as non-polynomial degradation functions. Two categories of degradation functions are examined, and a process to derive a given degradation function based on a local stress-strain response in the cohesive zone is presented. The resulting model is characterized by a linear elastic regime prior to the onset of damage, and controlled strain-softening thereafter. The governing equations are derived according to macro- and microforce balance theories, naturally accounting for the irreversible nature of the fracture process by introducing suitable constraints for the kinetics of the underlying microstructural changes. The model is complemented by an efficient staggered solution scheme based on an augmented Lagrangian method. Numerical examples demonstrate that the proposed model is a robust and effective method for simulating cohesive crack propagation, with particular emphasis on dynamic fracture.

42 ENGINEERING↗
DOE OSTI2016

Compositional effects on the ignition of FACE gasolines [Compositional effects on the ignition of FACE gasoline fuels: experiments, surrogate fuel formulation, and chemical kinetic modeling]

As regulatory measures for improved fuel economy and decreased emissions are pushing gasoline engine combustion technologies towards extreme conditions (i.e., boosted and intercooled intake with exhaust gas recirculation), fuel ignition characteristics become increasingly important for enabling stable operation. Here, this study explores the effects of chemical composition on the fundamental ignition behavior of gasoline fuels. Two well-characterized, high-octane, non-oxygenated FACE (Fuels for Advanced Combustion Engines) gasolines, FACE F and FACE G, having similar antiknock indices but different octane sensitivities and chemical compositions are studied. Ignition experiments were conducted in shock tubes and a rapid compression machine (RCM) at nominal pressures of 20 and 40 atm, equivalence ratios of 0.5 and 1.0, and temperatures ranging from 650 to 1270 K. Results at temperatures above 900 K indicate that ignition delay time is similar for these fuels. However, RCM measurements below 900 K demonstrate a stronger negative temperature coefficient behavior for FACE F gasoline having lower octane sensitivity. In addition, RCM pressure profiles under two-stage ignition conditions illustrate that the magnitude of low-temperature heat release (LTHR) increases with decreasing fuel octane sensitivity. However, intermediate-temperature heat release is shown to increase as fuel octane sensitivity increases. Various surrogate fuel mixtures were formulated to conduct chemical kinetic modeling, and complex multicomponent surrogate mixtures were shown to reproduce experimentally observed trends better than simpler two- and three-component mixtures composed of n-heptane, iso-octane, and toluene. Measurements in a Cooperative Fuels Research (CFR) engine demonstrated that the multicomponent surrogates accurately captured the antiknock quality of the FACE gasolines. Simulations were performed using multicomponent surrogates for FACE F and G to reveal the underlying chemical kinetics linking fuel composition with ignition characteristics. Finally, a key discovery of this work is the kinetic coupling between aromatics and naphthenes, which affects the radical pool population and thereby controls ignition.

02 PETROLEUM↗
DOE OSTI2020

Front propagation formulation for large eddy simulation of turbulent premixed flames

The front propagation formulation (FPF) for large eddy simulation (LES) of turbulent premixed flames is presented and discussed. The FPF method adopts a regularized Dirac delta function to minimize the spurious propagation of under-resolved filtered flame fronts in LES, while reproducing the reaction characteristics of premixed flame fronts. Here, the extension of the FPF method to LES and the relation to existing sub-filter combustion models are discussed. An algebraic model for the regularized Dirac delta function is presented and validated using direct numerical simulation (DNS) data for statistically planar premixed flames propagating in homogeneous isotropic turbulence. A posteriori validation for a filtered 1-D planar laminar flame is performed to demonstrate the characteristics of FPF in reproducing the propagation speed and the thickness of the under-resolved flame front. LES of a turbulent Bunsen flame is also performed to demonstrate the performance of the proposed method.

42 ENGINEERING↗
DOE OSTI2019

An MINLP formulation for the optimization of multicomponent distillation configurations

Designing configurations for multicomponent distillation, a ubiquitous process in chemical and petrochemical industries, is often challenging. This is because, as the number of components increases, the number of admissible distillation configurations grows rapidly and these configurations vary substantially in their energy needs. Consequently, if a method could identify a few energy-efficient choices from this large set of alternatives, it would be extremely attractive to process designers. This paper develops here such a method by solving a Mixed Integer Nonlinear Program (MINLP) that is formulated to pick, among the regular-column configurations of Shah and Agrawal (2010b), those configurations that have a low vapor-duty requirement. To compute the minimum vapor-duty requirement for each column within the configuration, we use techniques that rely on the Underwood’s method. The combined difficulty arising from the nonlinearity of Underwood equations and the combinatorial explosion of the choice-set of alternatives poses unmistakable challenges for the branch-and-bound algorithm, the current method of choice to globally solve MINLPs. To address this difficulty, we exploit the structure of Underwood equations and derive valid cuts that expedite the convergence of branch-and-bound by enabling global solvers, such as BARON, infer tighter bounds on Underwood roots. This provides a quick way to identify a few lucrative alternative configurations for separation of a given non-azeotropic mixture. We illustrate the practicality of our approach on a case-study concerning heavy-crude distillation and on various other examples from the literature.

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

A matrix-free hyperviscosity formulation for high-order ALE hydrodynamics

The numerical approximation of compressible hydrodynamics is at the core of high-energy density (HED) multiphysics simulations as shocks are the driving force in experiments like inertial confinement fusion (ICF). In this work, we describe our extension of the hyperviscosity technique, originally developed for shock treatment in finite difference simulations, for use in arbitrarily high-order finite element methods for Lagrangian hydrodynamics. Hyperviscosity enables shock capturing while preserving the high-order properties of the underlying discretization away from the shock region. Specifically, we compute a high-order term based on a product of the mesh length scale to a high power scaled by a hyper-Laplacian operator applied to a scalar field. We then form the total artificial viscosity by taking a non-linear blend of this term and a traditional artificial viscosity term. We also present a matrix-free formulation for computing the finite element based hyper-Laplacian operator. Such matrix-free methods have superior performance characteristics compared to traditional full matrix assembly approaches and offer advantages for GPU based HPC hardware. We demonstrate the numerical convergence of our method and its application to complex, multi-material ALE simulations on high-order (curved) meshes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Reducing impacts of artificial ponding in modeling salt marshes using a conductivity-formulated subgrid model

The landforms of salt marsh systems are rather complex, with meandering channel networks cutting through low-lying, extensive marsh platforms. Some small-scale topographic details are critical to flooding and draining processes but cannot be reflected in the Digital Elevation Model (DEM) due to limited measurement resolution and vegetation bias. Numerical models need a grid with high enough resolution to capture the small-scale connectivity to correctly model the small-scale flow processes, making the model computationally unaffordable. Here, in this study, we derived a subgrid model, which couples two flow components, one is the regular subgrid model for simulating the subgrid-resolved flows, and the other is a conductivity-formulated porous flow model for modeling small-scale gut flows unresolved at the subgrid level. The model validation was performed using an idealized case with a narrow slot that was unresolved in the coupled flow model. Model results were compared with a slot-resolved regular subgrid model where only surface flows are solved and suggested that, with a proper model parameterization, the coupled flow model is capable of reproducing flooding and draining processes similar to the slot-resolved subgrid model of surface flow. The model was applied to simulating flooding and draining processes in a meso-tidal salt marsh. The model results demonstrate that the coupled subgrid model can be used to reduce the impacts of artificial ponding that usually occur in numerical simulations of salt marshes with numerous narrow gullies and creeks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A general non-Fourier Stefan problem formulation that accounts for memory effects

The Stefan problem is the classical model of a melting phase change. In heterogeneous systems, such phase changes can exhibit non-Fourier (anomalous) behaviors, where the advance of the melt interface does not follow the expected time scaling. These situations can be modeled by replacing the derivatives, in the governing partial differential equations, with fractional order derivatives. In particular, replacing the time derivatives leads to non-Fourier models that account for memory effects in the system. In this work, by using appropriate time convolution integrals, a general thermodynamic balance statement for melting phase problems, explicitly accounting for memory effects, is developed. From this balance, a general model formulation applicable to problems involving melting over a temperature range (i.e., a mushy region) is derived. A key component in this model is the representation of memory effects through the use of fractional derivative based constitutive models of the enthalpy and heat flux. Further, on shrinking the mushy region to a single isotherm, a general sharp interface melting model is obtained. Here, in contrast to the classic Stefan problem, the fractional derivatives induce a natural regularization, such that the constitutive models for enthalpy and heat flux are continuous at the melt interface; a result confirmed through numerical simulation. To further support the theoretical findings, a physical example of a non-Fourier Stefan problem is presented. Overall the development and results in this paper underscore the importance of explicitly relating the development of fractional calculus models to the appropriate thermodynamic balance statements.

42 ENGINEERING↗
DOE OSTI2020

Modeling of the thermo-mechanical response and texture evolution of WE43 Mg alloy in the dynamic recrystallization regime using a viscoplastic self-consistent formulation

This paper presents a microstructure sensitive model for predicting mechanical response and texture evolution of metals in the dynamic recrystallization regime. A recently proposed viscoplastic self-consistent (VPSC) formulation for the prediction of recrystallization driven by strain energy and intragranular misorientation is extended to hexagonal close-packed (hcp) metals. The model is applied to the dynamic recrystallization of magnesium alloy WE43 at different temperatures and strain rates. Model predictions in terms of stress-strain response and texture evolution are compared to the experimental measurements and acceptable agreement is achieved. According to the model predictions, superplastic behavior of nuclei was found to be the dominant softening mechanism at high temperatures and low strain rates. High concentration of precipitates at the grain boundaries and presence of alloying elements are the likely causes of low boundary mobility, resulting in nucleation dominated dynamic recrystallization. Relatively strong basal compression textures indicate dominant activity of basal slip, which can be achieved only through large difference in slip resistance between soft basal and hard prismatic and pyramidal modes.

36 MATERIALS SCIENCE↗