Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analytical model”

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 757 records · Page 42

Highly selective Si 3 N 4 /SiO 2 etching using an NF 3 /N 2 /O 2 /H 2 remote plasma. II. Surface reaction mechanism

Developing processes for highly selective etching of silicon nitride (Si 3 N 4 ) with respect to silicon dioxide (SiO 2 ) is a major priority for semiconductor fabrication processing. In this paper and in Paper I [Volynets et al., J. Vac. Sci. Technol. A 38, 023007 (2020)], mechanisms are discussed for highly selective Si3N4 etching in a remote plasma based on experimental and theoretical investigations. The Si 3 N 4 /SiO 2 etch selectivity of up to 380 was experimentally produced using a remote plasma sustained in NF 3 /N 2 /O 2 /H 2 mixtures. A selectivity strongly depends on the flow rate of H 2 , an effect attributed to the formation of HF molecules in vibrationally excited states that accelerate etching reactions. Based on experimental measurements and zero-dimensional plasma simulations, an analytical etching model was developed for etch rates as a function of process parameters. Reaction rates and sticking coefficients were provided by quantum chemistry models and also fitted to the experimental results. Etch rates from the analytical model show good agreement with the experimental results and demonstrate why certain etchants accelerate or inhibit the etch process. In particular, the modeling shows the important role of HF molecules in the first vibrationally excited state [HF( v = 1)] in achieving high Si 3 N 4 /SiO 2 selectivity.

42 ENGINEERING↗

Spatial Resolution of a microCT Scanner

We have measured the spatial resolution of the 450 kV test-bed microCT scanner. We have also modeled the resolution analytically. The model and measurement show correspondence within 5%. Our measurement method was to use a thin tungsten wire to act as an impulse function. Therefore, the image is the impulse response, also called the point spread function (PSF) of the system. The Fourier transform (FT) of the impulse response is the modulation transfer function (MTF) of the system. We report the MTF here. By moving the wire at different positions in the scan field of view (FOV), we can measure MTF at different positions in the FOV. Since our imaging system is radially symmetric, we report MTF in two directions, radial and tangential as shown in Figure 1. The radial direction is along the line from isocenter to the wire, and the tangential is perpendicular to it.

42 ENGINEERING↗

Comparing the MELCOR Aerosol Deposition Model with Exact Analytical Solutions

Exact analytical solutions are presented for the evolution of the aerosol particle mass density function in a control volume for particle deposition due to gravitational settling, thermophoresis, and diffusion. The solutions are for arbitrary initial mass density functions and are applied for an initial lognormal density function. Integration of these solutions provides the suspended mass in the control volume as a function of time. These solutions serve as an exact benchmark to assess the accuracy of numerical methods. For the numerical algorithm used in MELCOR, excellent agreement is obtained for gravitational settling, diffusive deposition, and thermophoretic deposition for the suspended aerosol mass. In all cases, the default number of discrete particle size bins of 10 is shown to converge, with hardly any advantage to using 20 size bins.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Black and gray box learning of amplitude equations: Application to phase field systems

Here, we present a data-driven approach to learning surrogate models for amplitude equations and illustrate its application to interfacial dynamics of phase field systems. In particular, we demonstrate learning effective partial differential equations describing the evolution of phase field interfaces from full phase field data. We illustrate this on a model phase field system, where analytical approximate equations for the dynamics of the phase field interface (a higher-order eikonal equation and its approximation, the Kardar-Parisi-Zhang equation) are known. For this system, we discuss data-driven approaches for the identification of equations that accurately describe the front interface dynamics. When the analytical approximate models mentioned above become inaccurate, as we move beyond the region of validity of the underlying assumptions, the data-driven equations outperform them. In these regimes, going beyond black box identification, we explore different approaches to learning data-driven corrections to the analytically approximate models, leading to effective gray box partial differential equations.

42 ENGINEERING↗

A Three-Dimensional, Analytical Wind Turbine Wake Model: Flow Acceleration, Empirical Correlations, and Continuity

A new, three-dimensional, analytical, steady-state wake model is presented that includes local flow acceleration near the rotor, improving the wake description compared to existing models. Wake structures such as the momentum deficit and regions of accelerated flow are concisely described with compound and normal Gaussian functions. Large-eddy simulations (LES) are used as training data to develop the model using two in-line turbines under various inflow conditions parameterized by hub-height wind speed and turbulence intensity. Mass conservation is considered by fixing two components of the wake velocity model and optimizing the third to best satisfy continuity; after which, the model performs comparably if not better than existing work with regards to both relative error and mass consistency. The final model demonstrates a high degree of flexibility, making use of empirical correlations to scale across different inflow conditions. The inclusion of these effects is capable of revealing unused opportunities for enhanced power generation by aligning wake trajectories with these regions of accelerated flow.

modeling↗

Cutting Corners: Curvilinear-Surface-Based Gravity Models for Asteroids and Comets

We report that asteroids and comets are often highly irregular in shape and can contain large density inhomogenities. For such bodies, the simple point-mass gravity model does not capture the dynamics of the environment and this has driven a need for higher fidelity alternatives. A number of approaches have been developed with one of the most popular being the analytic polyhedral model. The analytic polyhedral model approximates the irregular shaped body as a polyhedron, typically consisting of triangular facets, with some assumed internal density distribution. Several analytic formulations have been derived for the gravitational fields of a constant-density polyhedron. A thorough history of the model is provided in Ref. [6] as well as a comparison of the implementations of Refs. [2–4]. The three versions provide near identical accuracy with the implementation of Werner requiring the fewest transcendental function evaluations suggesting its superior efficiency. Implementations with linear density contrasts and arbitrary polynomial density contrasts have also been developed.

79 ASTRONOMY AND ASTROPHYSICS↗

Twisted bilayer graphene as topological heavy fermion: II. Analytical approximations of the model parameters

The recently-introduced topological heavy fermion (THF) model [1] of twisted bilayer graphene (TBG) aims to reconcile the quantum-dot-like electronic structure of the latter observed by scanning tunneling microscopy, with its electron delocalization seen in transport measurements. The THF model achieves this by coupling localized (heavy) fermions with anomalous conduction electrons. Originally, the parameters of the THF model were obtained numerically from the Bistritzer–Macdonald (BM) model of TBG [1]. In this work, we derive analytical expressions for the THF model parameters as a function of the twist angle, the ratio between the tunneling amplitudes at the AA and AB regions (w 0 /w 1 ), and the screening length of the interaction potential. By numerically computing the THF model parameters across an extensive experimentally-relevant parameter space, we show that the resulting approximations are remarkably good, i.e., within the 30% relative error for almost the entire parameter space. At the single-particle level, the THF model accurately captures the energy spectrum of the BM model over a large phase space of angles and tunneling amplitude ratios. When interactions are included, we also show that the THF description of TBG is good around the magic angle for realistic values of the tunneling amplitude ratios (0.6 ≤ w 0 /w 1 ≤ 1.0), for which the hybridization between the localized and conduction fermions γ is smaller than the onsite repulsion of the heavy fermions U 1 (i.e., |γ| < U 1 ).

Approximation methods↗

Fast Surface Oxygen Release Kinetics Accelerate Nanoparticle Exsolution in Perovskite Oxides

Exsolution is a recent advancement for fabricating oxide-supported metal nanoparticle catalysts via phase precipitation out of a host oxide. A fundamental understanding and control of the exsolution kinetics are needed to engineer exsolved nanoparticles to obtain higher catalytic activity toward clean energy and fuel conversion. Since oxygen release via oxygen vacancy formation in the host oxide is behind oxide reduction and metal exsolution, we hypothesize that the kinetics of metal exsolution should depend on the kinetics of oxygen release, in addition to the kinetics of metal cation diffusion. In this report we probe the surface exsolution kinetics both experimentally and theoretically using thin-film perovskite SrTi 0.65 Fe 0.35 O 3 (STF) as a model system. We quantitatively demonstrated that in this system the surface oxygen release governs the metal nanoparticle exsolution kinetics. As a result, by increasing the oxygen release rate in STF, either by reducing the sample thickness or by increasing the surface reactivity, one can effectively accelerate the Fe 0 exsolution kinetics. Fast oxygen release kinetics in STF not only shortened the prereduction time prior to the exsolution onset, but also increased the total quantity of exsolved Fe 0 over time, which agrees well with the predictions from our analytical kinetic modeling. The consistency between the results obtained from in situ experiments and analytical modeling provides a predictive capability for tailoring exsolution, and highlights the importance of engineering host oxide surface oxygen release kinetics in designing exsolved nanocatalysts.

36 MATERIALS SCIENCE↗

A machine learning approach for clinker quality prediction and nonlinear model predictive control design for a rotary cement kiln

Abstract Cement manufacturing is energy‐intensive (5Gj/t) and comprises a significant portion of the energy footprint of concrete systems. Incorporating modern monitoring, simulation and control systems will allow lower energy use, lower environmental impact, and lower costs of this widely used construction material. One of the goals of the CESMII roadmap project on the Smart Manufacturing of Cement included developing an analytical process model for clinker quality that includes the chemistry of the kiln feed and accounts for critical process variables. This predictive model will be used in nonlinear model predictive control system designed to significantly reduce process energy use while maintaining or improving product quality. In the cement manufacturing plant used in this study, the kiln feed (meal) is tested every 12 h and used to estimate the mineral composition of the cement kiln output (clinker) using the stoichiometry‐based Bogue's model and the expertise of the plant operators. During kiln operation, kiln output (clinker) is sampled and tested every 2 h to measure its chemical and mineral composition. The predicted and measured values of the clinker composition are used by the plant operators to adjust the kiln input stream and the production process characteristics to maintain stable operation and uniform product quality. However, the time delay between prediction and testing, along with inaccuracies inherent in the Bogue's model have made any process changes designed to minimize energy use problematic, especially in‐light of potential clinker quality issues that process changes often pose. A new analytical model that integrates quality information and process operation information has been developed from data collected from 2 years of production from an operating cement facility. To make the model fuel‐type‐independent, consumed heat energy was computed in the model instead of fuel type and amount. A Feedforward Network was trained and tailored from collected data. Many data‐based simulations were conducted to quantitatively evaluate the proposed model and the 5‐fold cross‐validation procedure was used to test the models. The resulting predictive model was shown to have a low root mean square error (MSE) with respect to the estimated clinker mineral composition compared to that using the industry standard “Bogue’ model”. The end goal of this work was to develop a single machine learning tool that allows the use of quality control data and process control variables to improve energy efficiency of the process in a continuous fashion. The proposed nonlinear model predictive control system (NMPC) can generate predicted kiln production characteristics based on manipulated variables in manner that accurately follows the target product quality values. Simulation results also show that the proposed model produced accurate predictions of kiln outputs that fell within the required constraints, while manipulating control variables within typical operational ranges.

Ali, Asem M.↗

An analytic asymmetric-piston model for the impact of mode-1 shell asymmetry on ICF implosions

For many years, low mode asymmetry in inertially confined fusion (ICF) implosions has been recognized as a potential performance limiting factor, but analysis has been limited to using simulations and searching for data correlations. In this paper, an analytically solvable model based upon the simple picture of an asymmetric piston is presented. Asymmetry of the shell driving the implosion, as opposed to asymmetry in the hot-spot, is key to the model. The model provides a unifying framework for the action of mode-1 shell asymmetry and the resulting connections between various diagnostic signatures. A key variable in the model is the shell asymmetry fraction, f, which is related to the areal density variation of the shell surrounding the hot-spot. It is shown that f is simply related to the observed hot-spot mode-1 velocity and to the concept of residual energy in an implosion. The model presented in this paper yields explicit expressions for the hot-spot diameter, stagnation pressure, hot-spot energy, inertial confinement-time, Lawson parameter, hot-spot temperature, and fusion yield under the action of mode-1 asymmetry. Agreement is found between the theory scalings when compared to ICF implosion data from the National Ignition Facility and to large ensembles of detailed simulations, making the theory a useful tool for interpreting data. The theory provides a basis for setting tolerable limits on asymmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Considerations for Distributed Energy Resource Integration in Puerto Rico: DOE Multi-Lab Grid Modeling Support for Puerto Rico; Analytical Support for Interconnection and IEEE Std 1547-2018 National Renewable Energy Laboratory (Task 3.0)

To ensure the sustainable long-term recovery of Puerto Rico's electric power grid from hurricanes María and Irma and to build capacity to manage future potential natural disasters in the most secure and resilient way, the U.S. Department of Energy (DOE) convened experts from multiple national laboratories to develop a comprehensive set of data, models, analytic tools, and studies, considering inputs from a wide variety of stakeholder groups, to support technically sound recommendations for Puerto Rico's energy investment decisions. A resilient electric grid is vital to Puerto Rico's security, economy, and way of life, and it will provide the foundation for essential services that people and businesses on the island rely on every day. This report shows progress for a grid modeling task under the DOE-sponsored project that is a collaboration among the National Renewable Energy Laboratory (NREL) and other national laboratories. In Phase 1 of the multilab effort to support Puerto Rico's recovery, NREL provided the utility company Puerto Rico Electric Power Authority (PREPA) recommendations for a new framework of interconnection standards to accelerate the integration of utility-scale, transmission-connected renewable electrical generation and energy storage that ensure cross-technology compatibility and enable high deployment levels without compromising grid reliably, safety, or security. This Phase 2 report focuses on the interconnection of distributed energy resources (DERs) to the electric distribution system in Puerto Rico. This report is intended to familiarize the reader with Puerto Rico's distribution infrastructure and operational practices and procedures that are relevant to DER interconnection. The report also provides considerations for streamlining the interconnection process given the expected increase in deployments resulting from Puerto Rico's renewable portfolio standard (RPS) goal of 100% renewables by 2050. Accordingly, the report identifies considerations and concerns associated with the increase in intermittent generation, strategies for DER interconnection best practices, and the potential use of the latest technological solutions identified in the latest revision of the Institute of Electrical and Electronics Engineers (IEEE) 1547-2018 interconnection standard.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Flattening the Radial Temperature Profile across the Transformational Challenge Reactor Core

The Transformational Challenge Reactor (TCR) program is demonstrating an agile development approach to advanced nuclear reactor design, which has traditionally utilized a linear design process. In leveraging artificial intelligence, additive manufacturing, advanced materials, and cutting-edge modeling and simulation, the TCR program aims to minimize the high cost and lengthy deployment timelines now standard in the nuclear industry. Within a relatively short period of time, a robust and mature advanced gas-cooled reactor was iteratively designed under the TCR program, using these cutting-edge technologies. The TCR is a 3 MWt gas-cooled microreactor fueled with uranium nitride (UN) tristructural isotropic (TRISO) fuel particles. Though manufactured via traditional means, these UN TRISO particles are loaded into additively manufactured silicon carbide (SiC) cans [4]. Once loaded with TRISO particles, the SiC cans are densified using a chemical vapor infiltration process. The additively manufactured SiC enables significantly more freedom in the design of the fuel form than could ever be achieved using traditionally manufactured SiC. The helium coolant, pressurized to 5 MPa, enters the core at 300°C and nominally exits it at 500°C. Typically, the most thermally limiting components in any reactor design are the fuel assemblies in the core center. To provide a wide thermal margin in these central fuel assemblies, the flow may be biased toward the center of the core to more effectively cool these fuel assemblies with more power deposition and flatten the core’s radial temperature distribution. An analytical fluid model of the TCR core was developed to explore methods for biasing the flow away from the cooler outer fuel assemblies and towards the hotter inner ones. Higher-fidelity models developed in STAR-CCM+ 2020.3.1, a computational fluid dynamics code, were then utilized to verify the analytical model’s findings.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Breakdown of order fractionalization in the CPT model

Here, we present an analysis of the half-filled CPT model, an analytically tractable Kondo lattice model with Yao-Lee spin-spin interactions on a 3D hyperoctagon lattice, proposed by Coleman, Panigrahi, and Tsvelik. Previous studies have established that the CPT model exhibits odd-frequency triplet superconductivity and order fractionalization. Through asymptotic analyses in the small- and large- Kondo coupling limits, we identify a quantum critical point at , marking a transition from a superconductor to a Kondo insulator. By estimating the vison gap energy to account for thermal gauge fluctuations, we determine the energy scales governing the thermal breakdown of order fractionalization. Moreover, at large the Kondo insulator undergoes orbital decoupling, leading to the formation of a decoupled Kitaev orbital liquid. These findings and analogies with the -gauged model lead us to propose a tentative phase diagram for the CPT model at half-filling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Gravitation and Mesh Adaption

The Gravitation and Mesh Adaptation (GaMA) toolbox consists of a set of Matlab classes for modeling the environment around asteroids and comets. A variety of gravitational models and supporting algorithms from are consolidated alongside a custom meshing utility tailored for the application. This allows the user to work within a single streamlined environment to import and manipulate surface definitions, probe the dynamical environment, integrate trajectories, and post-process results. For other applications, modified surface meshes can be exported. Gravity Models: 1) Werner's analytic polyhedron model 2) Mascon model 3) Gottlieb's spherical harmonic model 4) Approximate polyhedron models 5) Curvilinear surface models 6) Custom composite models Meshing Features: 1) Array-based half-edge data structure 2) Supports curvilinear surface definitions up to degree 4 3) Ray-tracing 4) Coarsening 5) Feature-based refinement 6) Projection 7) Smoothing 8) Mesh quality validity tests 9) Mesh repair Additional Features: 1) Solar radiation pressure model 2) Distant 3rd body model 3) Collision detection 4) Trajectory integration and post-processing 5) Visualization of surface fields

Pearl, Jason↗

Comparative investigations of multi-fidelity modeling on performance of electrostatically-actuated cracked micro-beams

Silicon is a commonly used material for the fabrication of beams for use in micro-electrical-mechanical systems (MEMS). Although silicon is a brittle material, it has been shown to accumulate fatigue damage at the micro-scale. Understanding the effect this has on the overall device performance is critical to the design of reliable devices. Analytical methods for modeling damage provide expedient results but are limited by broad modeling assumptions. Numerical models account for more detailed physical phenomena but can be computationally intensive. In this work, two different crack scenarios are modeled using both analytical techniques and 3D computational simulations. First, the effects of a single surface crack on the static deflection and natural frequency of an electrostatically actuated micro-beam are formulated and compared. Then, a new method for approximating damage associated with realistic distributed crack networks is formulated for use in an analytical model and numerical simulations. A method for utilizing experimentally derived crack statistics to inform the analytical and numerical distributed crack models is developed. Good agreement between the analytical and numerical models is obtained for both crack scenarios. Altogether, these models can be used to effectively simulate a variety of damage and fatigue behaviors in silicon-based MEMS devices.

42 ENGINEERING↗