Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “argument”

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

Impact of the adhesion-force lever-arm “a” on the Rock ‘n’ Roll resuspension model and how to compute it from contact mechanics

The paper proposes a methodology to fill a gap in the context of mechanistic particle resuspension models related to the adhesion and lift force lever arms ( a ): to date, only one experimental observation exists from a centrifuge experiment with a specific particle/surface pair which found r P / a = 100 and assumed that the contact points were symmetric with respect to the particle center. However, the Reeks and Hall’s (2001) model is highly sensitive to this parameter, which therefore cannot be used as-is for different particle diameters or materials. This paper proposes a methodology to compute the adhesion and lift force lever arms ( a A and a L , respectively) based on Atomic Force Microcopy (AFM) surface roughness measurements and contact mechanics arguments. In practice, a spherical, perfectly smooth particle is brought in contact (computationally) with the measured surface roughness and is then rotated by the point of contact to determine the two or more asperities over which the particle would sit on the surface. A method to compute the asperity radius of curvature is devised to compute the asperity deformation and the correct partitioning of the total deformation between the particle and the asperity is presented. For the glass surface under consideration, the asperity deformation was negligible, i.e., it was below the van der Waals radius, which is the limit of applicability of contact mechanics. The Reeks and Hall’s (2001) model is then extended to account for the joint probability of a A and a L . Significant differences in particle resuspended after a day in common atmospheric conditions are found between the approach that uses the lever-arm distributions or their mean values. Also, important differences are found when comparing the lever-arm distribution and r P / a = 100. Hence, one should be cautious of using r P / a = 100 for surface/particle pairs differing from those in the original study. Finally, limits of applicability of the theories used are presented together with considerations regarding how to extend the method to non-spherical particles and different materials.

42 ENGINEERING↗
DOE OSTI2025

A Markov chain Monte Carlo (MCMC) Bayesian inference approach to analyze apparent activation barriers and reaction orders from microreactor data

Statistical analysis of steady-state catalytic kinetic data is often limited by data sparsity due to the slow pace at which the data is collected. Data sparsity and limitations in statistical analysis make it difficult to differentiate between mechanistic models and catalytic sites. A Bayesian inference tool is reported for catalysis researchers to estimate error in the determination of reaction orders from steady state microreactor data. The benefits of a Bayesian inference approach are discussed, as an alternative to the more common frequentist approach. The approach incorporates prior knowledge of the system and the data collected to form an error estimate on reaction orders. We investigated the effects of three distinct data treatments—individual fitting of trials, pooled analysis, and constrained regression methods—on the precision and uncertainty of reaction order determinations. To assess the robustness of our findings, we conducted sensitivity analyses to evaluate the influence of Bayesian parameters on uncertainty estimation. Additionally, we utilized synthetic data to illustrate how data quality impacts the precision of uncertainty assessments. We show Bayesian analysis can obtain a more precise estimation of error with a sparse data set than a frequentist analysis. Finally, this work provides strong evidence that the adoption of Bayesian analysis of kinetic data may help researchers make more precise arguments as to the strength of their evidence for a particular mechanistic hypothesis, or in comparing across different catalysts.

42 ENGINEERING↗
DOE OSTI2010

Probabilistic models for stochastic elliptic partial differential equations

Mathematical requirements that the random coefficients of stochastic elliptical partial differential equations must satisfy such that they have unique solutions have been studied extensively. Yet, additional constraints that these coefficients must satisfy to provide realistic representations for physical quantities, referred to as physical requirements, have not been examined systematically. It is shown that current models for random coefficients constructed solely by mathematical considerations can violate physical constraints and, consequently, be of limited practical use. We develop alternative models for the random coefficients of stochastic differential equations that satisfy both mathematical and physical constraints. Theoretical arguments are presented to show potential limitations of current models and establish properties of the models developed in this study. Numerical examples are used to illustrate the construction of the proposed models, assess the performance of these models, and demonstrate the sensitivity of the solutions of stochastic differential equations to probabilistic characteristics of their random coefficients.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

A method for solving stochastic equations by reduced order models and local approximations

A method is proposed for solving equations with random entries, referred to as stochastic equations (SEs). The method is based on two recent developments. The first approximates the response surface giving the solution of a stochastic equation as a function of its random parameters by a finite set of hyperplanes tangent to it at expansion points selected by geometrical arguments. The second approximates the vector of random parameters in the definition of a stochastic equation by a simple random vector, referred to as stochastic reduced order model (SROM), and uses it to construct a SROM for the solution of this equation. The proposed method is a direct extension of these two methods. It uses SROMs to select expansion points, rather than selecting these points by geometrical considerations, and represents the solution by linear and/or higher order local approximations. The implementation and the performance of the method are illustrated by numerical examples involving random eigenvalue problems and stochastic algebraic/differential equations. The method is conceptually simple, non-intrusive, efficient relative to classical Monte Carlo simulation, accurate, and guaranteed to converge to the exact solution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Finite Volume schemes on unstructured grids for non-local models: Application to the simulation of heat transport in plasmas

In the so-called Spitzer-Haerm regime, equations of plasma physics reduce to a nonlinear parabolic equation for the electronic temperature. Coming back to the derivation of this limiting equation through hydrodynamic regime arguments, one is led to construct a hierarchy of models where the heat fluxes are defined through a non-local relation which can be reinterpreted as well by introducing coupled diffusion equations. We address the question of designing numerical methods to simulate these equations. The basic requirement for the scheme is to be asymptotically consistent with the Spitzer-Haerm regime. Furthermore, the constraints of physically realistic simulations make the use of unstructured meshes unavoidable. We develop a Finite Volume scheme, based on Vertex-Based discretization, which reaches these objectives. We discuss on numerical grounds the efficiency of the method, and the ability of the generalized models in capturing relevant phenomena missed by the asymptotic problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

A novel adaptive time stepping variant of the Boris–Buneman integrator for the simulation of particle accelerators with space charge

We show that adaptive time stepping in particle accelerator simulation is an enhancement for certain problems. The new algorithm has been implemented in the OPAL (Object Oriented Parallel Accelerator Library) framework. The idea is to adjust the frequency of costly self-field calculations, which are needed to model Coulomb interaction (space charge) effects. In analogy to a Kepler orbit simulation that requires a higher time step resolution at the close encounter, we propose to choose the time step based on the magnitude of the space charge forces. Inspired by geometric integration techniques, our algorithm chooses the time step proportional to a function of the current phase space state instead of calculating a local error estimate like a conventional adaptive procedure. Building on recent work, a more profound argument is given on how exactly the time step should be chosen. An intermediate algorithm, initially built to allow a clearer analysis by introducing separate time steps for external field and self-field integration, turned out to be useful by its own, for a large class of problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Absorbing layers for the Dirac equation

This work is devoted to the construction of perfectly matched layers (PML) for the Dirac equation, that not only arises in relativistic quantum mechanics but also in the dynamics of electrons in graphene or in topological insulators. While the resulting equations are stable at the continuous level, some care is necessary in order to obtain a stable scheme at the discrete level. This is related to the so-called fermion doubling problem. For this matter, we consider the numerical scheme introduced by Hammer et al. [19], and combine it with the discretized PML equations. We state some arguments for the stability of the resulting scheme, and perform simulations in two dimensions. The perfectly matched layers are shown to exhibit, in various configurations, superior absorption than the absorbing potential method and the so-called transport-like boundary conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Stable evaluation of Green's functions in cylindrically stratified regions with uniaxial anisotropic layers

We present a robust algorithm for the computation of electromagnetic fields radiated by point sources (Hertzian dipoles) in cylindrically stratified media where each layer may exhibit material properties (permittivity, permeability, and conductivity) with uniaxial anisotropy. Analytical expressions are obtained based on the spectral representation of the tensor Green's function based on cylindrical Bessel and Hankel eigenfunctions, and extended for layered uniaxial media. Due to the poor scaling of these eigenfunctions for extreme arguments and/or orders, direct numerical evaluation of such expressions can produce numerical instability, i.e., underflow, overflow, and/or round-off errors under finite precision arithmetic. To circumvent these problems, we develop a numerically stable formulation through suitable rescaling of various expressions involved in the computational chain, to yield a robust algorithm for all parameter ranges. Numerical results are presented to illustrate the robustness of the formulation including cases of practical interest.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Scalable algorithms for three-field mixed finite element coupled poromechanics

We introduce a class of block preconditioners for accelerating the iterative solution of coupled poromechanics equations based on a three-field formulation. The use of a displacement/velocity/pressure mixed finite-element method combined with a first order backward difference formula for the approximation of time derivatives produces a sequence of linear systems with a 3 x 3 unsymmetric and indefinite block matrix. The preconditioners are obtained by approximating the two-level Schur complement with the aid of physically-based arguments that can be also generalized in a purely algebraic approach. A theoretical and experimental analysis is presented that provides evidence of the robustness, efficiency and scalability of the proposed algorithm. In conclusion, the performance is also assessed for a real-world challenging consolidation experiment of a shallow formation.

42 ENGINEERING↗
DOE OSTI2017

Thermodynamical effects and high resolution methods for compressible fluid flows

One of the fundamental differences of compressible fluid flows from incompressible fluid flows is the involvement of thermodynamics. This difference should be manifested in the design of numerical schemes. Unfortunately, the role of entropy, expressing irreversibility, is often neglected even though the entropy inequality, as a conceptual derivative, is verified for some first order schemes. In this paper, we refine the GRP solver to illustrate how the thermodynamical variation is integrated into the design of high resolution methods for compressible fluid flows and demonstrate numerically the importance of thermodynamic effects in the resolution of strong waves. As a by-product, we show that the GRP solver works for generic equations of state, and is independent of technical arguments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Adaptive spectral solution method for the Landau and Lenard-Balescu equations

In this paper, we present an adaptive spectral method for solving the Landau/Fokker-Planck equation for electron-ion systems. The heart of the algorithm is an expansion in Laguerre polynomials, which has several advantages, including automatic conservation of both energy and particles without the need for any special discretization or time-stepping schemes. One drawback of such an expansion is the O(N 3 ) memory requirement, where N is the number of polynomials used. This can impose an inconvenient limit in cases of practical interest, such as when two particle species have widely separated temperatures. The algorithm we describe here addresses this problem by periodically re-projecting the solution onto a judicious choice of new basis functions that are still Laguerre polynomials but have arguments adapted to the current physical conditions. This results in a reduction in the number of polynomials needed, at the expense of increased solution time. Because the equations are solved with little difficulty, this added time is not of much concern compared to the savings in memory. To demonstrate the algorithm, we solve several relaxation problems that could not be computed with the spectral method without re-projection. Another major advantage of this method is that it can be used for collision operators more complicated than that of the Landau equation, and we demonstrate this here by using it to solve the non-degenerate quantum Lenard-Balescu (QLB) equation for a hydrogen plasma. We conclude with some comparisons of temperature relaxation problems solved with the latter equation and the Landau equation with a Coulomb logarithm inspired by the properties of the QLB operator. We find that with this choice of Coulomb logarithm, there is little difference between using the two equations for these particular systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Permutation-adapted complete and independent basis for atomic cluster expansion descriptors

Atomic cluster expansion (ACE) methods provide a systematic way to describe particle local environments of arbitrary body order. For practical applications it is often required that the basis of cluster functions be symmetrized with respect to rotations and permutations. Existing methodologies yield sets of symmetrized functions that are over-complete. These methodologies thus require an additional numerical procedure, such as singular value decomposition (SVD), to eliminate redundant functions. In this work, it is shown that analytical linear relationships for subsets of cluster functions may be derived using recursion and permutation properties of generalized Wigner symbols. From these relationships, subsets (blocks) of cluster functions can be selected such that, within each block, functions are guaranteed to be linearly independent. It is conjectured that this block-wise independent set of permutation-adapted rotation and permutation invariant (PA-RPI) functions forms a complete, independent basis for ACE. Along with the first analytical proofs of block-wise linear dependence of ACE cluster functions and other theoretical arguments, numerical results are offered to demonstrate this. The utility of the method is demonstrated in the development of an ACE interatomic potential for tantalum. Using the new basis functions in combination with Bayesian compressive sensing sparse regression, some high degree descriptors are observed to persist and help achieve high-accuracy models.

Angular momentum↗
DOE OSTI2010

Capillary induced self-assembly of thin foils into 3Dstructures

Self-assembly of complex structures is common in nature. Self-assembly principles provide a promising way to fabricate three-dimensional, micro- or millimeter scale devices. In the present paper, we present a generalized analytical study of the self-folding of thin plates into deterministic 3D shapes through fluid–solid interactions. Based on the beam theory, a mechanics model is developed, incorporating the two competing components—a capillary force promoting folding and the bending rigidity of the foil that resists folding into a 3D structure. Through an equivalence argument of thin foils of different geometry, an effective folding parameter, which uniquely characterizes the driving force for folding, has been identified. A criterion for spontaneous folding of any shaped 2D patterned foil based on the effective folding parameter is thus established. Furthermore, the model predictions show excellent agreement with experimental measurements made on a variety of materials, indicating that the assumptions used in the analysis arevalid.

42 ENGINEERING↗
DOE OSTI2026

A parametric study of slow dynamic nonlinear elasticity with comparisons to models

Several phenomenological models that aspire to quantitative description of anomalous nonlinear mesoscopic elasticity are reviewed and compared with laboratory measurements. This class of nonlinearity, best known perhaps for slow dynamics and aging, is seen widely in imperfectly consolidated granular solids but is not well understood. Typical slow dynamic tests show that a modest conditioning oscillatory "pump" strain depresses material stiffness, which then recovers like the logarithm of time after conditioning ceases. Several phenomenological models based on physical arguments have been proposed that predict the material stiffness response to arbitrary pump strain histories during conditioning and recovery. Approximate closed form and numerical solutions to the models are presented that predict the quantitative influence of three key pump parameters: the pump's strain amplitude, the pump's strain rate, and the pump’s duration. Laboratory measurements on Berea sandstone, concrete and a confined single aluminum bead find that slow dynamic responses are linear in pump strain and independent of pump frequency. Measurements also show that, after pump-off, stiffness recovers over times far longer than the pump duration. These observations and others are compared to model predictions. One of the considered models, based on a picture of fast brittle damage and slow healing, successfully matches all these behaviors.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Revisiting the relationship between demand growth and electricity prices

Here, in this Commentary, we aim to clarify the relationship between electricity demand and prices in the short- and longer-term. We outline three dimensions that shape this relationship: system capacity utilization, system expansion costs, and cost allocation in the rate design process. We explain why demand growth has historically been largely associated with falling electricity prices and explore arguments that future C&I demand growth may increase prices.

O'Shaughnessy, Eric [Lawrence Berkeley National La↗
DOE OSTI2010

Synthesis, structures and magnetic properties of pseudo-hollandite chromium sulfides

The pseudo-hollandite chromium sulfides, ACr{sub 5}S{sub 8} (A=K, Rb, Cs) and A'{sub 0.5}Cr{sub 5}S{sub 8} (A'=Sr, Ba), have been synthesized and investigated in structural and magnetic properties. All the compounds crystallize in the isostructure with a monoclinic C2/m. Its crystal structure has triangular lattices and double chains made of Cr{sup 3+} (d{sup 3}; S=3/2) triangles. The magnetic susceptibilities of all compounds behave as Curie-Weiss types in high temperature region. From magnetic susceptibility and specific heat measurements, all the compounds have antiferromagnetic ground states. The Neel temperatures are rather low compared to Weiss temperatures, reflecting magnetic frustration in the triangular lattices and double chains. The magnetic transitions in KCr{sub 5}S{sub 8} and Ba{sub 0.5}Cr{sub 5}S{sub 8} are a two-step transition around 50 and 60 K, respectively, while RbCr{sub 5}S{sub 8} shows a sharp magnetic transition at 42 K, accompanied by magnetoelastic behavior. CsCr{sub 5}S{sub 8} shows a structural transition around 100 K, followed by a magnetic transition at 10 K. In Sr{sub 0.5}Cr{sub 5}S{sub 8}, ferromagnetic interaction develops below 100 K and then a three-dimensional antiferromagnetic order takes place at 30 K. These magnetic properties are discussed from A-cation dependences of local structures and a model of magnetic structure for RbCr{sub 5}S{sub 8} is proposed on the basis of the arguments of magnetic interactions and neutron diffraction data. It is different from the known magnetic structure of TlCr{sub 5}Se{sub 8}. - Graphical abstract: Temperature dependences of magnetic susceptibilities for pseudo-hollandite sulfides, ACr{sub 5}S{sub 8} (A=K, Rb, Cs) and A{sub 0.5}'Cr{sub 5}S{sub 8} (A'=Sr, Ba). The observed magnetic frustration effect originates in the triangular lattice.

36 MATERIALS SCIENCE↗
DOE OSTI2011

Local ordering and magnetism in Ga{sub 0.9}Fe{sub 3.1}N

Prior investigations of the ternary nitride series Ga{sub 1-x}Fe{sub 3+x}N (0{<=}x{<=}1) have indicated a transition from ferromagnetic {gamma}'-Fe{sub 4}N to antiferromagnetic 'GaFe{sub 3}N'. The ternary nitride 'GaFe{sub 3}N' has been magnetically and spectroscopically reinvestigated in order to explore the weakening of the ferromagnetic interactions through the gradual incorporation of gallium into {gamma}'-Fe{sub 4}N. A hysteretic loop at RT reveals the presence of a minority phase of only 0.1-0.2 at%, in accord with the sound two-step synthesis. The composition of the gallium-richest phase 'GaFe{sub 3}N' was clarified by Prompt Gamma-ray Activation Analysis and leads to the berthollide formula Ga{sub 0.91(1)}Fe{sub 3.09(10)}N{sub 1.05(7)}. Magnetic measurements indicate a transition around 8 K, further supported by Moessbauer spectral data. The weakening of the ferromagnetic coupling through an increasing gallium concentration is explained by a simple Stoner argument. In Ga{sub 0.9}Fe{sub 3.1}N, the presence of iron on the gallium site affects the magnetism by the formation of 13-atom iron clusters. - Graphical Abstract: The crystal structure of GaFe{sub 3}N with green nitrogen atoms in the very center, red iron atoms at the face centers, and gray gallium atoms at the corner positions. Highlights: > Almost phase-pure synthesis of Ga{sub 0.9}Fe{sub 3.1}N. > Prompt gamma-ray activation analysis yields precise composition. > Magnetic ordering of the facial Fe atoms at the lowest temperatures. > Moessbauer spectroscopy suggests percolation or RKKY-type interaction. > Fe{sub 13} clusters due to berthollide character.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗