Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Step”

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 811 records · Page 45

Multirate linearly-implicit GARK schemes

Many complex applications require the solution of initial-value problems where some components change fast, while others vary slowly. Multirate schemes apply different step sizes to resolve different components of the system, according to their dynamics, in order to achieve increased computational efficiency. The stiff components of the system, fast or slow, are best discretized with implicit base methods in order to ensure numerical stability. To this end, linearly implicit methods are particularly attractive as they solve only linear systems of equations at each step. This paper develops the Multirate GARK-ROS/ROW (MR-GARK-ROS/ROW) framework for linearly-implicit multirate time integration. The order conditions theory considers both exact and approximative Jacobians. The effectiveness of implicit multirate methods depends on the coupling between the slow and fast computations; an array of efficient coupling strategies and the resulting numerical schemes are analyzed. Multirate infinitesimal step linearly-implicit methods, that allow arbitrarily small micro-steps and offer extreme computational flexibility, are constructed. The new unifying framework includes existing multirate Rosenbrock(-W) methods as particular cases, and opens the possibility to develop new classes of highly effective linearly implicit multirate integrators.

97 MATHEMATICS AND COMPUTING↗

Influence of interrupted cooling on the development of bimodal γ' precipitate distributions in ATI 718Plus

The influence of aging treatments on the evolution of the γ' precipitates in ATI 718Plus was studied using scanning electron microscopy, energy dispersive spectroscopy, transmission electron microscopy, and atomic probe tomography. A set of aging treatments was performed on solution-treated samples to obtain various γ' size distributions. The aging temperatures ranged from 720 to 900°C with holding times of either two hours (for the 900°C aging temperature) or 10 h (for the 720°C and 750°C aging temperatures) followed by either air cooling, furnace cooling, or water quenching. The cooling method and aging temperature influenced the distribution of the solute elements. This work discussed the evolution of the γ' precipitation during single-step aging and two-step aging involving continuous cooling or interrupted cooling. Finally, a bimodal size distribution of the γ' precipitates was obtained through interrupted cooling two-step aging treatments, while single-step aging resulted in unimodal γ' precipitate distributions.

36 MATERIALS SCIENCE↗

Microstructure and Temperature Dependent Indentation Response of Additively Manufactured Precipitation-Strengthened Al 0.3 Ti 0.2 Co 0.7 CrFeNi 1.7 High Entropy Alloy

The temperature-dependent [from room temperature (RT) to 500°C] nanoindentation behavior of a precipitation-strengthened Al 0.3 Ti 0.2 Co 0.7 CrFeNi 1.7 high-entropy alloy (HEA) processed via two different additive manufacturing (AM) techniques was investigated in the as-deposited and annealed conditions. The hierarchically heterogeneous microstructures were achieved via simple one-step annealing treatments, exploiting the residual stresses in the AM-processed HEA to partially recrystallize the microstructure, performed remarkably better than the nearly homogeneous microstructures in the as-deposited state. Further, the one-step annealed conditions revealed < 6.6% reduction in hardness values at 500°C compared to RT, while the as-deposited conditions showed a > 18% reduction in the hardness. The one-step annealed conditions also exhibited significantly higher hardness than the as-deposited conditions owing to their L1 2 -strengthened FCC microstructures. Furthermore, serrated yielding or the Portevin-Le Chatlier effect indicative of microstructural instability was observed during nanoindentation deformation (at 500°C) for the as-deposited conditions but not for the one-step annealed conditions. This, therefore, signifies the robustness of the hierarchically heterogeneous microstructures at elevated temperatures presenting a strong avenue for tuning the HEAs for future nuclear reactor applications.

36 MATERIALS SCIENCE↗

Mechanistic Insights into the Catalytic Condensation of Methyl Ketones on MgO Surfaces

Ketone coupling via aldol condensation is one of the promising routes to produce cyclic and value-added precursors for renewable hydrocarbon biofuels. A first-principles-based microkinetic modeling is performed to evaluate the surface-mediated reaction mechanisms and the role of water molecules in the observed activities for 2-pentanone and 3-pentanone aldol condensation on dehydroxylated MgO(111) surface and hydroxylated terminated surface[OH-MgO(111)]. Here, we have identified the enhancement of the surface OH group to MgO(111) surface catalytic activity by destabilizing the binding strength of reaction intermediates and reducing the energy barriers of rate-determining steps(proton transfer and dehydration steps). The 2-pentanone has one elementary step less in the complete reaction mechanism of aldol condensation and preferable energy barrier for proton transfer and dehydration steps, revealing 2-pentanone as terminal ketone is more reactive than 3-pentanone as central ketone. The water molecules dominated the OH-MgO( 111) surface after further addition of water, leading to the reduction of turnover frequency of the aldol condensation dimer product as the loss of aldol condensation reaction intermediates in competitive adsorption with water molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of a Griffin model of the advanced test reactor

In the pursuit of a higher fidelity deterministic simulation capability of the Advanced Test Reactor, it is important to have a fast yet accurate deterministic neutronics model. Here, to achieve this, we employed an advanced two-step method. The first step involves generating homogenized cross sections using OpenMC, a cutting-edge Monte Carlo neutron transport code. OpenMC offers excellent modular capabilities, allowing for easy component integration and flexibility in incorporating new designs into the model. The second step involves deterministic transport calculations, which are performed using Griffin, a reactor physics application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE). To ensure the accurate spatial resolution and assignment of material cross sections, a Cubit-generated mesh for the Advanced Test Reactor is utilized as an intermediate step between the OpenMC and Griffin models; Griffin utilizes the mesh for its finite element solution, while OpenMC material identifications are written to the mesh file to be used in Griffin material assignments. Additionally, a Python-based script converts the cross sections generated by OpenMC into the ISOXML format required by Griffin. Initial comparisons using the Griffin diffusion solver indicated good agreement between the neutron multiplication factors obtained from the standalone OpenMC model and the Griffin model, with differences of less than 10 pcm in the 2D geometry configuration; it was later determined that this agreement was likely due to compensating effect and was more likely on the order of –700 pcm relative to the OpenMC solution. However, in three-dimensional calculations, an unacceptably large error (almost 8,000 pcm) was found in the Griffin solution with the diffusion solver. Subsequent calculations using Griffin’s discrete ordinates solver demonstrated substantially improved agreement, within 116 pcm of the OpenMC solution used to generate the cross sections for Griffin. Building on this capability, future work will seek to perform more detailed validation calculations. The ultimate goal is to evaluate both transient and multiphysics simulations of the reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Adaptive primal–dual control for distributed energy resource management

With the increased adoption of distributed energy resources (DERs) in distribution networks, their coordinated control with a DER management system (DERMS) that provides grid services (e.g., voltage regulation, virtual power plant) is becoming more necessary. One particular type of DERMS using primal–dual control has recently been found to be very effective at providing multiple grid services among an aggregation of DERs; however, the main parameter, the primal–dual step size, must be manually tuned for the DERMS to be effective, which can take a considerable amount of engineering time and labor. To this end, we design a simple method that self-tunes the step size(s) and adapts it to changing system conditions. Additionally, it gives the DER management operator the ability to prioritize among possibly competing grid services. Here we evaluate the automatic tuning method on a simulation model of a real-world feeder in Colorado with data obtained from an electric utility. Through a variety of scenarios, we demonstrate that the DERMS with automatically and adaptively tuned step sizes provides higher-quality grid services than a DERMS with a manually tuned step size.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Experimental demonstration of a 5 kW th granular-flow reactor for solar thermochemical energy storage with aluminum-doped calcium manganite particles

A two-step cycle was discussed for solar thermochemical energy storage based on aluminum-doped calcium manganite reduction/oxidation reactions for direct integration into Air Brayton cycles. The two steps encompassed 1) the storage of concentrated solar direct irradiation via the thermal reduction of aluminum-doped calcium manganite and 2) the delivery of heat to an Air-Brayton cycle via re-oxidation of oxygen-deficient aluminum-doped calcium manganite. The re-oxidized aluminum-doped calcium manganite was fed back to the first step to complete the cycle. A 5 kWth solar thermochemical reactor operating under vacuum was fabricated and tested to examine the first cycle reduction step. Reactor operating conditions and high-flux solar simulator control were tuned for continuous reactor operation with particle temperatures > 1073 K. Continuous operation was achieved using intermittent, dense granular flows. A maximum absorption efficiency of 64.7% was shown, accounting for both sensible and chemical heat storage.

25 ENERGY STORAGE↗

Origins of wear-induced tungsten corrosion defects in semiconductor manufacturing during tungsten chemical mechanical polishing

During chemical mechanical polishing (CMP) process optimization, a unique corrosion defect was identified in the tungsten (W) contact plugs, which were investigated to develop practical solutions to improve CMP mass productivity. The corrosion defect occurred in the absence of external incident energy and was observed to occur only at the W plugs connected to the P + /n-well close to the wafer edge. To examine the root cause of the corrosion, a series of experiments was carried out at each W CMP step, using in-line defect inspection. Additional experiments were conducted using on-line inductively coupled plasma mass spectrometry (ICP-MS) to understand pH- and potential-dependent electrochemical dissolution behaviors of W. Experimentally, W corrosion was confirmed to be exacerbated if the surface oxide of the W plugs is removed by the deionized water (DIW) rinsing step under wafer pressurization, then exposed to DIW and alkaline cleaning solution. Although the DIW rinsing step is commonly performed to clean polishing pads right after the main polishing, a small amount of applied pressure (0.5 psi) without slurry supply can lead to corrosion defects. Finally, this wear-induced corrosion of the W plugs can be prevented by removing the applied pressure during the post DIW rinsing step of W CMP.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A low-rank solver for the stochastic unsteady Navier–Stokes problem

Here we study a low-rank iterative solver for the unsteady Navier–Stokes equations for incompressible flows with a stochastic viscosity. The equations are discretized using the stochastic Galerkin method, and we consider an all-at-once formulation where the algebraic systems at all the time steps are collected and solved simultaneously. The problem is linearized with Picard’s method. To efficiently solve the linear systems at each step, we use low-rank tensor representations within the Krylov subspace method, which leads to significant reductions in storage requirements and computational costs. Combined with effective mean-based preconditioners and the idea of inexact solve, we show that only a small number of linear iterations are needed at each Picard step. The proposed algorithm is tested with a model of flow in a two-dimensional symmetric step domain with different settings to demonstrate the computational efficiency.

97 MATHEMATICS AND COMPUTING↗

Data-driven variational multiscale reduced order models

We propose a new data-driven reduced order model (ROM) framework that centers around the hierarchical structure of the variational multiscale (VMS) methodology and utilizes data to increase the ROM accuracy at a modest computational cost. The VMS methodology is a natural fit for the hierarchical structure of the ROM basis: In the first step, we use the ROM projection to separate the scales into three categories: (i) resolved large scales, (ii) resolved small scales, and (iii) unresolved scales. In the second step, we explicitly identify the VMS–ROM closure terms, i.e., the terms representing the interactions among the three types of scales. In the third step, we use available data to model the VMS–ROM closure terms. Thus, instead of phenomenological models used in VMS for standard numerical discretizations (e.g., eddy viscosity models), we utilize available data to construct new structural VMS–ROM closure models. Specifically, we build ROM operators (vectors, matrices, and tensors) that are closest to the true ROM closure terms evaluated with the available data. We test the new data-driven VMS–ROM in the numerical simulation of four test cases: (i) the 1D Burgers equation with viscosity coefficient v = 10 -3 ; (ii) a 2D flow past a circular cylinder at Reynolds numbers, Re = 100, Re = 500, and Re= 1000; (iii) the quasi-geostrophic equations at Reynolds number Re = 450 and Rossby number Ro = 0.0036; and (iv) a 2D flow over a backward facing step at Reynolds number Re = 1000. The numerical results show that the data-driven VMS–ROM is significantly more accurate than standard ROMs.

42 ENGINEERING↗

DECA: Discrete Event inspired Cellular Automata for grain structure prediction in additive manufacturing

Microstructure largely dictates macroscopic material properties and is strongly affected by processing. Therefore, the simulation of microstructure evolution in response to thermal fields during processing is of significant interest within the computational materials science community. Additive manufacturing (AM) has emerged as a technique for producing complex geometries and unique microstructures. Yet, complex and rapid thermal cycles in AM pose computational challenges for existing microstructure models. This work proposes a discrete event inspired cellular automata (CA) approach, titled DECA, to accelerate simulation of grain structure evolution in AM. In contrast to conventional time-stepped CA models, this model directly solves the times capture events would take place allowing for stepping in events rather than time (a technique also found in the field of discrete-event simulation). In comparison to purely serial discrete-event models, DECA allows for temporary violation of the causality constraint, but detects and corrects these violations, leading to an emergent phenomenon dubbed causality rippling, in which previously calculated capture events are overwritten. The amount of repeated calculations, defined by the capture ratio, is taken as a measure of computational inefficiency, and the model parameters that affect this ratio are evaluated. The new DECA approach was found to be more computationally efficient than conventional time-stepped CA models while guaranteeing an accurate solution, which can only be achieved in the conventional models for vanishingly small time steps. Finally, opportunities for parallelization and scaling of the new approach are discussed.

36 MATERIALS SCIENCE↗

Power quality disturbances diagnosis: A 2D densely connected convolutional network framework

The fast and accurate diagnosis of power quality disturbances (PQD) aids in avoiding shutdowns and unnecessary procedures, concerning electric energy distribution systems. As such, a number of techniques have been tested and applied in order to reach this objective. Majority of the techniques applied are two-step based. On the first step, power quality disturbances features are extracted. Second step, considering features extracted, disturbance classification is implemented. Recently, relevant literature has presented data-driven signal processing-based approaches, as deep convolutional neural networks (DCNN), which can implement both processing steps while providing automated recognition of patterns and outliers in data. However, not considered by state-of-art, power quality disturbances are evolving in nature, while all possible regularities might not be represented in the dataset. In this work a 2 Dimension Densely Connected Convolutional Network (2D-DenseNet) framework is presented. Further, a case study with synthetic disturbance events are analyzed. Easy-to-implement formulation, built on the 2D-DenseNet, without hard-to-design parameters, highlight potential aspects for real-life implementation.

42 ENGINEERING↗

Corrosion-resistant polyurethane coatings from structure-homogenized biorefinery lignin through fractionation and oxypropylation

Using lignin for polyurethane (PU) synthesis is challenging because of its heterogeneous structure, i.e., broad molecule weight (MW) distribution and diverse hydroxyl (OH) functionality. This complexity influences the crosslinking reaction between lignin and isocyanate and affects the quality of PU products. Here, we propose a two-step strategy for lignin structure homogenization by lowering its MW and uniforming OH functionality (step 2). Step 1 involves a bioethanol fractionation process, yielding fractionated lignin with one order lower MW than the raw lignin while exhibiting better processability for the subsequent process. Step 2 is performed by an oxyalkylation process using fractionated lignin that converts phenolic and carboxylic OHs to aliphatic OHs. As a result, we obtained structure-homogenized lignin with good compatibility and crosslinking with hexamethylene diisocyanate for PU coating synthesis. Moreover, PU coatings prepared from structure-homogenized lignin display a defect-free microscopic structure and show better corrosion resistance than coatings fabricated from its raw lignin alternative. Overall, the structure homogenization process is a beneficial way of using lignin as biopolyol in PU synthesis toward the application-oriented design of corrosion-resistant lignin-based PU coatings.

36 MATERIALS SCIENCE↗

DFT and microkinetic comparison of Pt, Pd and Rh-catalyzed ammonia oxidation

Ammonia oxidation is the heart of the Ostwald process and is important in emissions control. Catalytic behaviors are a function of conditions and are observed to vary across the platinum group metals (PGMs) Pt, Pd, Rh. Here, we combine density functional theory computations and microkinetic modeling to rationalize these dependencies. We compute reactions over model (2 1 1) and (1 1 1) surfaces of PGMs. Binding energies are similar on Pd and Pt and generally greater on Rh, while activation energies vary across all metals. Rates on (2 1 1) surfaces are greater than (1 1 1) surfaces. Additionally, the stepped Pt is most active and stepped Rh most selective to N 2 at ammonia slip conditions, while at Ostwald process conditions, stepped Pd is most active and stepped Pt most selective to NO. Degree of rate and selectivity control analysis provides insights into the reactions limiting performance of PGMs. Both activation barriers and surface coverages influence rates and selectivities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Catalytic site ensembles: A context to reexamine the Langmuir-Hinshelwood kinetic description

The Langmuir-Hinshelwood formalism describes catalytic reactions of Langmuirian surface species under the assumption that all adsorbates are randomly-distributed – enabling adjacency of surface-bound intermediates to be determined solely by coverages of single-site occupants. We demonstrate herein that this approximation is inappropriate even for simple catalytic reactions (e.g. A + A → A 2 ) and manifestly neglects islanding of slowly-consumed species and partitioning of highly-reactive species inherently engendered by ≥ two-site elementary steps (e.g. A*–A* → A 2(g) + *–*). Rigorous description of kinetically-consequential islanding/partitioning phenomena requires explicit description of the coverage and chemical dynamics of all multi-site ensembles. Higher-order, ensemble-specific rate terms identify the particular microscopic events relevant to each ensemble, and, in doing so, reveal that each elementary step (e.g. A (g) adsorption) describes not one event (e.g. A (g) + * → A*) , but a sum over all ensemble-specific paths (e.g. A (g) + *–* → A*–* and A (g) + A*–* → A*–A*). De-convoluting each elementary step into its constituent multi-site paths proffers kinetic detail otherwise inaccessible – enabling (i) identification of rate- and selectivity-determining site ensembles, (ii) calculation of rates and degrees of rate control of ensemble-specific elementary steps, (iii) incorporation of adsorbate surface diffusion, (iv) incorporation of lateral adsorbate interactions, and (v) quantitative description of catalysis of multi-site-occupying intermediates (e.g. *C n H m * species in hydrocarbon (de-)hydrogenation and C-C bond coupling/cleavage reactions) which we demonstrate here is inaccessible to the Langmuir-Hinshelwood formalism even if adsorbate surface diffusion is infinitely-fast.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comparative life cycle analysis on ethylene production from electrocatalytic reduction of carbon dioxide

Ethylene is one of the largest greenhouse gas emitters and the most diversly used commodity chemicals globally. Electrocatalytic reduction of CO 2 to ethylene received great attention from the research society to decarbonize the ethylene production. Here, in this study, a life-cycle analysis is conducted using the Greenhouse Gases, Regulated Emissions, and Energy use in Technologies (GREET) model on the three electrocatalytic CO 2 -reduction pathways (or “e-ethylene” pathways): i) cascade conversion via carbon monoxide intermediate; ii) single-step conversion in membrane electrode assembly (MEA); and iii) single-step conversion in alkaline flow cell. The results showed that the electricity consumption is the lowest for the cascade pathway (164 MJ/kg), thus resulting in the lowest cradle-to-gate carbon intensity [18 kgCO 2 e/kg with United States (US) average grid)] among the three pathways followed by the single-step MEA (32 kgCO 2 e/kg) and then by the single-step alkaline (56 kgCO 2 e/kg). However, all three e-ethylene pathways were significantly more carbon-intensive than their fossil-based counterpart (1.1 kgCO 2 e/kg) due to their excessive energy consumption with the current state of technology. With renewable electricity, all three pathways yielded negative carbon intensity: from -3.1 kgCO 2 e/kg to -1.6 kgCO 2 e/kg depending on the source of CO 2 . The threshold carbon intensity of electricity (TCIE), defined as the upper bound of the carbon intensity of electricity to achieve lower carbon intensity for e-ethylene compared to fossil-based ethylene, is calculated for both current and future state of e-ethylene technologies. The cascade pathway had the highest TCIE out of the three e-ethylene pathways for both current (92 gCO 2 e/kWh) and future (124 gCO 2 e/kWh) state of technologies. However, the carbon intensity of average US grid (i.e., 467 and 303 gCO 2 e/kWh for current and future projections) were higher than the TCIEs of the corresponding timeline. Thus, reducing electricity requirement for e-ethylene pathways and bringing low-carbon generation mix in the United States (US) grid faster than the current projection are both essential to decarbonize ethylene and its downstream chemicals/polymers.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

An analysis of the spatio-temporal resolution of the immersed boundary method with direct forcing

The immersed boundary method (IBM) with direct forcing is very popular in the simulation of rigid particulate flows. In the IBM, an interaction force is introduced at the interface between fluid and particle in order to approximate the no-slip boundary condition. The interaction force is calculated through dividing the velocity difference (or error) between fluid and particle at the interface by the time step. Here, a dynamic equation for the velocity difference is derived. Additionally, analyses on the dynamic equation provide a few new findings: (i) The interaction force is the solution of a least-squares error problem, with the direct implication that the Lagrangian marker distribution has no effect on the large scale flow structure once the distribution of Lagrangian markers become saturated along the interface (i.e., each marker remains properly correlated with all its neighbors); (ii) The Lagrangian volume-weight is a relaxation factor to control how fast the velocity error decays to the ideal value of zero; (iii) The optimal choice of the Lagrangian volume-weight is the largest value permissible by a stability condition. A comprehensive convergence analysis with regard to the spatial and temporal resolution is presented for the velocity error and also for the shear-stress and surface pressure. In three simple canonical problems, it is analytically and numerically shown that the IBM results converge to the theoretical solutions obtained with precise imposition of no-slip and no-penetration boundary conditions. It is observed that it is not necessary to match the Lagrangian marker volume-weight to that of the local Eulerian cell volume and in fact this matching leads to lower than optimal computational efficiency. However, it is found that extremely high Eulerian grid resolution and small time step have to be used to obtain high precision simulation results. Especially, the time step should be inversely proportional to the particle Reynolds number for low Reynolds number flows. For high frequency oscillation problems, the grid size needs to be reduced by a factor of the square root of the frequency, and the time step to be reduced by a factor of the frequency. The theoretical findings here can be used to alleviate the technical difficulties in simulating non-spherical particles by not requiring the Lagrangian marker distribution to match the Eulerian grids and also in the implementation of IBM on non-uniform Eulerian grids. The present work also provides simple practical guidance on the choice of temporal and spatial resolution so as to control the simulation error a priori.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Reprint of "Selective conversion of CO into ethanol on Cu(511) surface reconstructed from Cu(pc): Operando studies by electrochemical scanning tunneling microscopy, mass spectrometry, quartz crystal nanobalance, and infrared spectroscopy"

A polycrystalline copper, surface-terminated by a well-defined (511)-oriented facet, was electrochemically generated by a series of step-wise surface reconstruction and iterations of mild oxidative-reductive processes in 0.1 M KOH. The electrochemical reduction of CO on the resultant stepped surface was investigated by four surface-sensitive operando methodologies: electrochemical scanning tunneling microscopy (STM), electrochemical quartz crystal nanobalance (EQCN), differential electrochemical mass spectrometry (DEMS), and polarization-modulation infrared spectroscopy (PMIRS). The stepped surface catalyzed the facile conversion of CO into ethanol, the exclusive alcohol product at a low overpotential of -1.06 V (SHE) or - 0.3 V (RHE). Additionally, the chemisorption of CO was found to be a necessary prelude to ethanol production; i.e. the surface coverages, rather than solution concentrations, of CO and its surface-bound intermediates primarily dictate the reaction rates (current densities). Contrary to the expected predominance of undercoordinated step-site reactivity over the coordination chemistry of vicinal surfaces, vibrational spectroscopic evidence reveals the involvement of terrace-bound CO adsorbates during the multi-atomic transformations associated with the production of ethanol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗