Search NASA⌕ Search

SEARCH · Search NASA

Results for “convention equation”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Coupling a Lagrangian–Eulerian Spark-Ignition (LESI) model with LES combustion models for engine simulations

In the United States transportation sector, Light-Duty Vehicles (LDVs) are the largest energy consumers and CO 2 emitters. Electrification of LDVs is posed as a potential solution, but SI engines can still contribute to decarbonization. Car manufacturers have turned to unconventional engine operation to increase the efficiency of Spark-Ignition (SI) engines and reduce the carbon emissions of their fleets. Dilute, lean, and stratified-charge engine operation has the potential for engine efficiency improvements at the expense of increased cyclic variability and combustion instability. At such demanding engine conditions, the spark ignition event is key for flame initiation and propagation and for enhanced combustion stability. Reliable and accurate spark ignition models can help design ignition systems that reduce cyclic variability. Multiple computational spark-ignition models exist that perform well under conventional conditions, but the underlying physics needs to be expanded, for unconventional engine operation. In this paper, a hybrid Lagrangian–Eulerian Spark-Ignition (LESI) model is coupled with different turbulent flame propagation models for engine simulations. LESI relies on Lagrangian arc tracking and Eulerian energy deposition. The LESI model is coupled with the Well-Stirred Reactor (WSR), Thickened Flame Model (TFM), and g-equation model and used to simulate several cycles of a Direct-Injection Spark-Ignition (DISI) engine using a commercial Computational Fluid Dynamics (CFD) engine solver. The results showcase the successful coupling of LESI with the combustion models. Global engine metrics, such as pressure and Apparent Heat Release Rate (AHRR), for each simulation setup are compared to experimental engine results, for validation. In addition, results highlight the successful prediction of spark channel movement by comparing simulation images to experimental optical engine images. Finally, the successful coupling of LESI to combustion models, making it a usable model in the engine modeling community, is emphasized and future development details are discussed.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Multi-scale model based design of membrane reactor/separator processes for intensified hydrogen production through the water gas shift reaction

This paper aims to first quantify the impact of various diffusion models (Maxwell-Stefan, Wilke, Dusty-Gas) on the predictions of a multi-scale membrane reactor/separator mathematical model, and to then demonstrate this model's use for the design and process intensification of membrane reactor/separator systems for hydrogen production. This multi-scale model captures velocity, temperature and species' concentration profiles along the catalyst pellet's radial direction, and along the reactor's axial direction, by solving the momentum, energy, and species transport equations, accounting for convection, conduction, reaction, and diffusion mechanisms. In the first part of work, the effect of pellet-scale design parameters (mean pore diameter, volumetric porosity, tortuosity factor, etc.) and various species' flux models on the model predictions is studied. In the second part, the study focuses on the comparison, in terms of their process intensification characteristics, of various hydrogen production processes. These include a conventional high-temperature shift reactor (HTSR)/low-temperature shift reactor (LTSR) sequence, a novel HTSR/membrane separator (MS)/LTSR/MS sequence, and a process that involves low-temperature shift membrane reactors-LTSMR in a series.

08 HYDROGEN↗
DOE OSTI2026

Geometrical optics without singularities: using the ray time as the coordinate space

Geometrical optics (GO) is widely used for reduced modelling of waves in plasmas, but it fails near reflection points, where it predicts a spurious singularity of the wave amplitude. We show how to avoid this singularity by adopting a different representation of the wave equation. Instead of the physical coordinate 𝑥 and the wavevector 𝑘, we use the ray time 𝜏 as the new canonical coordinate and the ray energy ℎ as the associated canonical momentum. To derive the envelope equation in the 𝜏-representation, we construct the Weyl symbol calculus on the (𝜏,ℎ) space and show that the corresponding Weyl symbols are related to their (𝑥,𝑘) counterparts by the Airy transform. This allows us to express the coefficients in the envelope equation through the known properties of the original dispersion operator. When necessary, solutions of this equation can be mapped to the 𝑥-space using a generalised metaplectic transform. However, the field per se might not even be needed in practice. Instead, knowing the corresponding Wigner function usually suffices for linear and quasilinear calculations. As a Weyl symbol itself, the Wigner function can be mapped analytically, using the aforementioned Airy transform. We show that the standard Airy patterns that form in regions where conventional GO fails are successfully reproduced within metaplectic GO (MGO) simply by remapping the field from the 𝜏-space to the 𝑥-space. An extension to mode-converting waves is also presented. This formulation, which we call generalised MGO, can be particularly useful, for example, for reduced modelling of the O–X conversion in inhomogeneous plasma near the critical density, an effect that is important for fusion applications and also occurs in the ionosphere. Overall, MGO can replace GO for any practical purposes, because it better handles cutoffs and is similar otherwise.

plasma waves↗
DOE OSTI2023

Synthetically non-Hermitian nonlinear wave-like behavior in a topological mechanical metamaterial

Topological mechanical metamaterials have enabled new ways to control stress and deformation propagation. Exemplified by Maxwell lattices, they have been studied extensively using a linearized formalism. Herein, we study a two-dimensional topological Maxwell lattice by exploring its large deformation quasi-static response using geometric numerical simulations and experiments. We observe spatial nonlinear wave-like phenomena such as harmonic generation, localized domain switching, amplification-enhanced frequency conversion, and solitary waves. We further map our linearized, homogenized system to a non-Hermitian, nonreciprocal, one-dimensional wave equation, revealing an equivalence between the deformation fields of two-dimensional topological Maxwell lattices and nonlinear dynamical phenomena in one-dimensional active systems. Our study opens a regime for topological mechanical metamaterials and expands their application potential in areas including adaptive and smart materials and mechanical logic, wherein concepts from nonlinear dynamics may be used to create intricate, tailored spatial deformation and stress fields greatly transcending conventional elasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Multidimensional Rietveld refinement of high-pressure neutron diffraction data of PbNCN

High-pressure neutron powder diffraction data from PbNCN were collected on the high-pressure diffraction beamline SNAP located at the Spallation Neutron Source (SNS) of Oak Ridge National Laboratory (Tennessee, USA). The diffraction data were analyzed using the novel method of multidimensional (two dimensions for now, potentially more in the future) Rietveld refinement and, for comparison, employing the conventional Rietveld method. To achieve two-dimensional analysis, a detailed description of the SNAP instrument characteristics was created, serving as an instrument parameter file, and then yielding both cell and spatial parameters as refined under pressure for the first time for solid-state cyanamides/carbodiimides. The bulk modulus B 0 = 25.1 (15) GPa and its derivative B ′ 0 = 11.1 (8) were extracted for PbNCN following the Vinet equation of state. Surprisingly, an internal transition was observed beyond 2.0 (2) GPa, resulting from switching the bond multiplicities (and bending direction) of the NCN 2− complex anion. The results were corroborated using electronic structure calculation from first principles, highlighting both local structural and chemical bonding details.

Meinerzhagen, Yannick↗
DOE OSTI2022

Smoothed boundary method for simulating incompressible flow in complex geometries

Simulating flow through porous media with explicit considerations of complex microstructures is very challenging using conventional sharp-interface methods because of the difficulties in generating meshes conformal to complex geometries. In this work, a diffuse interface embedded boundary method known as the Smoothed Boundary Method (SBM) is utilized to facilitate simulations of fluid dynamics involving complex geometries. In diffuse-interface methods, the geometry is described by a domain parameter. The SBM allows the straightforward reformulation of the time-dependent Navier-Stokes equations in terms of this domain parameter, using only algebraic identities. Thus, enforcing the appropriate boundary conditions at the irregular embedded boundary is greatly simplified. Adaptive mesh refinement is used to increase the accuracy of the diffuse interface method by allowing thinner interfacial thickness to be used in the domain parameter. Here, the SBM-formulated Navier-Stokes equations are solved with the Finite Difference Method on refined mesh systems. Sharp-interface Finite Element Method simulations using the commercial software COMSOL on body-conforming meshes are also provided for comparison. Favorable agreement between the two methods is observed. Since it is no longer necessary for the mesh to conform to the complex geometry, the grid system for the SBM simulations can be generated rapidly and without additional manual interventions, making the entire simulation process more expedient.

36 MATERIALS SCIENCE↗
DOE OSTI2023

On the design of stable, consistent, and conservative high-order methods for multi-material hydrodynamics

Obtaining stable and high-order numerical solutions for multi-material hydrodynamics is an open challenge. Although slope limiters are widely used to maintain monotonicity near discontinuities, typical limiting procedures violate closure laws at the discrete level when applied to multi-material hydrodynamics equations. Due to this, the high-order expansions of quantities related by the closure laws are no longer consistent. The commonly observed symptom of this consistency-violation is that the numerical method fails to maintain constant pressure and velocity across material interfaces. This leads to sub-optimal convergence rates for smooth multi-material problems as well. Specialized limiting procedures that satisfy consistency while maintaining conservation need to be developed for such equations. A novel procedure that re-instates consistency into slope-limited high-order discretizations applied to the multi-material hydrodynamics equations is presented here. Using simple examples, it is demonstrated that the presented method satisfies closure laws at the discrete level, while maintaining conservative properties of the high-order method. Furthermore, this procedure involves a projection step which relies on the compact basis of the underlying spatial discretization, i.e. for discontinuous schemes (viz. DG and FV) the projection is local, and does not involve global matrix solves. Comparisons with conventional approaches emphasizes the necessity of the consistent closure-law preserving limiting approach, in order to maintain design order of accuracy for smooth multi-material problems.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Particle production by time-varying dark energy and the end of cosmic expansion

We consider various possible consequences of time-varying dark energy due to a quintessence scalar field whose energy density is partially converted to particles as the field evolves down its potential. This particle production acts as a source of thermal friction on the field that can make it difficult to distinguish whether dark energy is due to a radiating field rolling down a steep potential, a purely self-interacting field moving down a flatter potential, or a cosmological constant. By reducing the acceleration of the scalar field, thermal friction increases the amount of accelerated expansion and can cause a sizable bump in the quintessence equation of state. We take special interest in the case where a steep potential rapidly changes from positive to negative as the field evolves, resulting in the end of cosmic expansion and the beginning of contraction. Even in this case, we find that thermal friction lengthens the period of accelerated expansion and consequently delays the end of cosmic expansion, making it challenging to detect the impending transition to contraction using conventional cosmological tests. However, particle production can also provide alternative avenues for detection by generating a background of thermal dark radiation, partly comprised of neutrinos or other particles, whose energy density exceeds the remnant photon energy density.

cosmological neutrinos↗
DOE OSTI2025

Particle Production by Time-Varying Dark Energy and the End of Cosmic Expansion

We consider various possible consequences of time-varying dark energy due to a quintessence scalar field whose energy density is partially converted to particles as the field evolves down its potential. This particle production acts as a source of thermal friction on the field that can make it difficult to distinguish whether dark energy is due to a radiating field rolling down a steep potential, a purely self-interacting field moving down a flatter potential, or a cosmological constant. By reducing the acceleration of the scalar field, thermal friction increases the amount of accelerated expansion and can cause a sizable bump in the quintessence equation of state. We take special interest in the case where a steep potential rapidly changes from positive to negative as the field evolves, resulting in the end of cosmic expansion and the beginning of contraction. Even in this case, we find that thermal friction lengthens the period of accelerated expansion and consequently delays the end of cosmic expansion, making it challenging to detect the impending transition to contraction using conventional cosmological tests. However, particle production can also provide alternative avenues for detection by generating a background of thermal dark radiation, partly comprised of neutrinos or other particles, whose energy density exceeds the remnant photon energy density.

Cosmology and Nongalactic Astrophysics (astro-ph.C↗
DOE OSTI2024

The advancing wave front on a sloping channel covered by a rod canopy following an instantaneous dam break

The drag coefficient Cd for a rigid and uniformly distributed rod canopy covering a sloping channel following the instantaneous collapse of a dam was examined using flume experiments. The measurements included space x and time t high resolution images of the water surface h(x, t) for multiple channel bed slopes So and water depths behind the dam Ho along with drag estimates provided by sequential load cells. Using these data, an analysis of the Saint-Venant equation (SVE) for the front speed was conducted using the diffusive wave approximation. An inferred Cd=0.4 from the h(x, t) data near the advancing front region, also confirmed by load cell measurements, is much reduced relative to its independently measured steady-uniform flow case. This finding suggests that drag reduction mechanisms associated with transients and flow disturbances are more likely to play a dominant role when compared to conventional sheltering or blocking effects on Cd examined in uniform flow. The increased air volume entrained into the advancing wave front region as determined from an inflow–outflow volume balance partly explains the Cd reduction from unity.

Mechanics↗
DOE OSTI2022

Phase transformation and growth mechanism of RF sputtered ferroelectric lead scandium tantalate (PbSc 0.5 Ta 0.5 O 3 ) films

Lead scandium tantalate (PbSc 0.5 Ta 0.5 O 3 , PST), an order/disorder ferroelectric, is a potential candidate for electrocaloric cooling and pyroelectric infrared (IR) detector. In this work, we report the phase transformation kinetics from two series of samples containing pure amorphous and mixture of amorphous and pyrochlore to desired perovskite phase using postdeposition rapid thermal processing (RTP) as well as growth mechanism of RF sputtered PST thin films using excess lead target on platinized silicon (Pt/Ti/SiO 2 /Si) substrates. We find that small changes in the temperature ramp have a large effect on the degree of perovskite conversion (ferroelectric phase), orientation (crystallographic texture), and long-range order parameter (< S 111 >). Through isothermal annealing, we obtained optimal perovskite phase at ≥700°C temperature. The phase transformation is characterized by spontaneous formation of center-type in-plane radial rosette-like structures revealed by scanning electron microscopy. The PST perovskite crystallites were found to coexist with pyrochlore in RTP annealed films. The volume fractions for perovskite and pyrochlore phase were obtained from the analysis of “rosettes” and respective X-ray diffraction intensities which helped to determine various parameters associated with phase kinetics (n, k, and activation energy, Ea) and accompanying growth. The effective activation energies of perovskite transition and growth were found to be 332 ± 11 kJ/mol (345 ± 11 kJ/mol) and 114 ± 10 kJ/mol (122 ± 10 kJ/mol), respectively, for pure amorphous only (and mixed amorphous and pyrochlore) phase following nucleation-growth controlled Avrami's equation. A linear growth rate (n~1) for the perovskite phase indicates predominant interface-controlled process and diffusion-limited phenomena thus inhibiting rosette size owing to reactant depletion and soft impingement at the grain boundary. However, the growth behavior is isotropic in two-dimension parallel to the plane of the substrates for both sample series. Furthermore, lead loss was severe for in-situ growth and RTP combined with conventional furnace annealing than those of RTP only films, which were closer to stoichiometric albeit with excess lead and marginal oxygen vacancies (V o ).

36 MATERIALS SCIENCE↗
DOE OSTI2022

Application of the Continuum Damage Mechanics Wilshire-Cano-Stewart (WCS) Model

In this study, the applications of the continuum damage mechanics-based Wilshire-Cano-Stewart (WCS) model are explored to predict rupture time, minimum-creep-strain-rate (MCSR), damage, damage evolution, and creep deformation. Increase knowledge in manufacturing methods has pushed the limit of material science and the development of new materials. Conventional testing is required to qualify materials against creep which according to the ASME B&PV III code, 10,000+ hours of experiments are necessary for each heat before materials are put into service. This process is costly and not feasible for new materials. As an alternative, models have been employed to predict creep behaviors and reduce the amount of time necessary for material qualification. Many models have been developed to predict distinct creep behaviors and the question of which model is best remains. Amongst current models, the WCS model has emerge with the ability to predict multiple behaviors using an explicit analytical approach with the ability to predict long-term creep. In this study the novel continuum damage mechanics WCS is employed in multiple applications. The goals of the study are (a) to discuss and determine the framework of the WCS model and validated it mathematically and using parametric simulations, (b) applied the model to accelerated creep data to show the capabilities of the model with non-conventional data, and (c) applied a novel numerical method, the datum temperature method (DTM) to show the model extrapolations and interpolations capabilities with limited and reduce data sets. To accomplish these goals, data is gathered for alloy P91 and Inconel 718 to develop and post-audit validate the model. The benefits of using the WCS model is that it provides an explicit stress and temperature dependency ideal for extrapolations, the coupled equations are suitable for finite element analysis (FEA) implementation, and it follows an explicit calibration approach. The model also proves that it can be applied to accelerated testing data and using the DTM.

Cano, Jaime A↗
DOE OSTI2023

Integral boundary conditions in phase field models

Modeling the chemical, electric and thermal transport as well as phase transitions and the accompanying mesoscale microstructure evolution within a material in an electronic device setting involves the solution of partial differential equations often with integral boundary conditions. Employing the familiar Poisson equation describing the electric potential evolution in a material exhibiting insulator to metal transitions, we exploit a special property of such an integral boundary condition, and we properly formulate the variational problem and establish its well-posedness. Next, we compare our method with the commonly-used Lagrange multiplier method that can also handle such boundary conditions. Numerical experiments demonstrate that our new method achieves optimal convergence rate in contrast to the conventional Lagrange multiplier method. Furthermore, the linear system derived from our method is symmetric positive definite, and can be efficiently solved by Conjugate Gradient method with algebraic multigrid preconditioning.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Physics-informed neural network with transfer learning (TL-PINN) based on domain similarity measure for prediction of nuclear reactor transients

Nuclear reactor safety and efficiency can be enhanced through the development of accurate and fast methods for prediction of reactor transient (RT) states. Physics informed neural networks (PINNs) leverage deep learning methods to provide an alternative approach to RT modeling. Applications of PINNs in monitoring of RTs for operator support requires near real-time model performance. However, as with all machine learning models, development of a PINN involves time-consuming model training. Here, we show that a transfer learning (TL-PINN) approach achieves significant performance gain, as measured by reduction of the number of iterations for model training. Using point kinetic equations (PKEs) model with six neutron precursor groups, constructed with experimental parameters of the Purdue University Reactor One (PUR-1) research reactor, we generated different RTs with experimentally relevant range of variables. The RTs were characterized using Hausdorff and Fréchet distance. We have demonstrated that pre-training TL-PINN on one RT results in up to two orders of magnitude acceleration in prediction of a different RT. The mean error for conventional PINN and TL-PINN models prediction of neutron densities is smaller than 1%. We have developed a correlation between TL-PINN performance acceleration and similarity measure of RTs, which can be used as a guide for application of TL-PINNs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Mechanisms of Zn removal from water by amorphous geopolymer: Molecular-level insights from X-ray absorption spectroscopy, isotope fractionation, and surface complexation modeling

Porous geopolymers have attracted widespread attention as promising heavy metal adsorbents that can be synthesized from aluminosilicate solid wastes. However, the precise microstructural evidence for adsorbed heavy metals on geopolymers remains unclear due to the insensitivity of conventional characterization techniques on minerals with amorphous structure and surface disorder. Batch adsorption and column experiments coupled with X-ray absorption spectroscopy (XAS), Zn stable isotope, and surface complexation model (SCM) were employed to reveal the Zn removal mechanisms with coal fly ash porous geopolymer (CFAPG) at a molecular scale. The macroscopic kinetic and isothermal adsorption of Zn on CFAPG were well described by the pseudo-second-order model and Bi_Langmuir equation, respectively, indicating the presence of abundant heterogeneous active sites on the CFAPG surface. Further, two types of active sites on the CFAPG surface were identified by XAS coupled with Zn isotopes in batch experiments at pH <= 6.0: one is pH-dependent and associated with tetrahedral zinc coordi-nation, and the other is pH-insensitive and associated with octahedral zinc coordination; these sites were confirmed by the bidentate SCM as the variable charge site (surface complexation, >S-OH) and the permanent negative charge site (cation exchange, >X - ), respectively. Furthermore, the important contribution of surface co -precipitation besides surface complexation and cation exchange to the Zn adsorption on CFAPG was identified by XAS coupled with SCM in a flow-through column experiment at pH >6.0. These investigations provide a systemic understanding of the Zn adsorption mechanisms on CFAPG and an SCM reference for the application and prediction of geopolymers in heavy metal-contaminated water remediations.

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

Adiabatic Evolution of Cloud Droplet Spectral Width: A New Look at an Old Problem

Abstract Spectral width of the cloud droplet spectrum is important for radiative properties and drizzle/rain development in warm ice‐free clouds. We use an adiabatic rising parcel model to study activation and diffusional growth of cloud droplets, focusing on the spectral width evolution, and contrasting clean and polluted environments. A comprehensive droplet growth equation is used that includes kinetic, solute, and surface tension effects. We show that those effects have an appreciable impact on the spectral width evolution above the cloud base. Without those effects, the droplet area standard deviation should not change once activation is completed. In contrast, simulation results show that the area standard deviation does increase with height, especially for weak updrafts and polluted environments. Implications of those results for cloud modeling, especially applying conventional bin microphysics, are discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Demystifying the Stern layer at a metal–electrolyte interface: Local dielectric constant, specific ion adsorption, and partial charge transfer

Electric double layer (EDL) represents one of the most basic concepts in electrochemistry and is pertinent to diverse engineering applications ranging from electrocatalysis to energy storage. Whereas phenomenological and coarse-grained models have been long established to describe ionic distributions in the diffuse layer, a faithful prediction of the physicochemical properties of the electrode–electrolyte interface from a molecular perspective remains a daunting challenge. In this work, we investigate the charging behavior of an Ag (111) electrode in NaF aqueous solutions leveraging experimental results and theoretical calculations based on the classical density functional theory for ion distributions in the diffuse layer and on the joint density functional theory (JDFT) for the electronic structure. When the Ag electrode is applied with a negative voltage, the surface charge density can be reasonably described by assuming a neutral Stern layer with the dielectric constant dependent on the local electric field as predicted by the Kirkwood equation. However, the specific adsorption of F− ions must be considered when the electrode is positively charged and the fluoride adsorption can be attributed to both physical and chemical interactions. Qualitatively, F− binding and partial charge transfer are supported by JDFT calculations, which predict an increased binding energy as the voltage increases. Our findings shed insight on the molecular characteristics of the Stern layer and the charge behavior of adsorbed species not specified by conventional EDL models.

Wang, Xuepeng (ORCID:0000000285704546)↗
DOE OSTI2026

Microscopic constraints for the equation of state and structure of neutron stars: A Bayesian model mixing framework

Bayesian model mixing (BMM) is a statistical technique that can combine constraints from different regions of an input space in a principled way. Here we extend our BMM framework for the equation of state (EOS) of strongly interacting matter from symmetric nuclear matter to asymmetric matter, specifically focusing on zero-temperature, charge-neutral, 𝛽-equilibrated matter. We use Gaussian processes (GPs) to infer constraints on the neutron-star matter EOS at intermediate densities from two different microscopic theories: chiral effective-field theory (𝜒⁢EFT) at baryon densities around nuclear saturation, 𝑛 𝐵 ∼ 𝑛 0 , and perturbative QCD at asymptotically high baryon densities, 𝑛 𝐵 ⩾ 20⁢𝑛 0 . The uncertainties of the 𝜒⁢EFT and pQCD EOSs are obtained using the BUQEYE truncation error model. We demonstrate the flexibility of our framework through the use of two categories of GP kernels: conventional stationary kernels and a nonstationary changepoint kernel. We use the latter to explore potential constraints on the dense matter EOS by including exogenous data representing theory predictions and heavy-ion collision measurements at densities ⩾ 2⁢𝑛 0 . We also use our EOSs to obtain neutron-star mass-radius relations and their uncertainties. Finally, our framework, whose implementation will be available through a GitHub repository, provides a prior distribution for the EOS that can be used in large-scale neutron-star inference frameworks.

Bayesian methods↗