Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iterative limiters”

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 541 records · Page 30

Learning local equivariant representations for large-scale atomistic dynamics

Abstract A simultaneously accurate and computationally efficient parametrization of the potential energy surface of molecules and materials is a long-standing goal in the natural sciences. While atom-centered message passing neural networks (MPNNs) have shown remarkable accuracy, their information propagation has limited the accessible length-scales. Local methods, conversely, scale to large simulations but have suffered from inferior accuracy. This work introduces Allegro, a strictly local equivariant deep neural network interatomic potential architecture that simultaneously exhibits excellent accuracy and scalability. Allegro represents a many-body potential using iterated tensor products of learned equivariant representations without atom-centered message passing. Allegro obtains improvements over state-of-the-art methods on QM9 and revMD17. A single tensor product layer outperforms existing deep MPNNs and transformers on QM9. Furthermore, Allegro displays remarkable generalization to out-of-distribution data. Molecular simulations using Allegro recover structural and kinetic properties of an amorphous electrolyte in excellent agreement with ab-initio simulations. Finally, we demonstrate parallelization with a simulation of 100 million atoms.

74 ATOMIC AND MOLECULAR PHYSICS↗

Applying Constrained Bayesian Optimization to the Design of Critical Experiments

Often when planning a criticality experiment, many design configurations are iteratively investigated with a Monte Carlo transport code. The goal is that the experiment will be optimal with respect to some variable, like the fraction of fissions occurring at a certain energy range, while simultaneously being critical. Unfortunately, the Monte Carlo transport simulations are expensive, which can ultimately limit the number of configurations that can be explored. In this work, we present how Gaussian processes (GPs) can be used as a reduced-order model in a constrained Bayesian optimization (CBO) algorithm to design a criticality experiment. The GPs replace the Monte Carlo transport simulations that explore the design space. The CBO algorithm efficiently identifies new points in the design space to run the Monte Carlo transport code while respecting the criticality constraint. It does so in a manner that both improves the accuracy of the GP and finds the approximate global optimum. We demonstrate the performance of CBO with the design of a Thermal Epithermal eXperiment (TEX) for the criticality safety validation of nuclear waste models of the Hanford Tank Farm.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Forward marching procedure for separated boundary-layer flows

A forward-marching procedure for separated boundary-layer flows which permits the rapid and accurate solution of flows of limited extent is presented. The streamwise convection of vorticity in the reversed flow region is neglected, and this approximation is incorporated into a previously developed (Carter, 1974) inverse boundary-layer procedure. The equations are solved by the Crank-Nicolson finite-difference scheme in which column iteration is carried out at each streamwise station. Instabilities encountered in the column iterations are removed by introducing timelike terms in the finite-difference equations. This provides both unconditional diagonal dominance and a column iterative scheme, found to be stable using the von Neumann stability analysis.

Carter, J. E.↗

The UTRC wind energy conversion system performance analysis for horizontal axis wind turbines (WECSPER)

The theory for the UTRC Energy Conversion System Performance Analysis (WECSPER) for the prediction of horizontal axis wind turbine performance is presented. Major features of the analysis are the ability to: (1) treat the wind turbine blades as lifting lines with a prescribed wake model; (2) solve for the wake-induced inflow and blade circulation using real nonlinear airfoil data; and (3) iterate internally to obtain a compatible wake transport velocity and blade loading solution. This analysis also provides an approximate treatment of wake distortions due to tower shadow or wind shear profiles. Finally, selected results of internal UTRC application of the analysis to existing wind turbines and correlation with limited test data are described.

Egolf, T. A.↗

Development of an autonomous video rendezvous and docking system, phase 3

Field-of-view limitations proved troublesome. Higher resolution was required. Side thrusters were too weak. The strategy logic was improved and the Kalman filter was augmented to estimate target attitude and tumble rate. Two separate filters were used. The new filter estimates target attitude and angular momentum. The Newton-Raphson iteration improves image interpretation.

Tietz, J. C.↗

SYS-5620-Project - Achieving Robust Laser Performance utilizing Historical Shot Experiments - A Systems Engineering Proposal

The National Ignition Facility (NIF) at Lawrence Livermore National Laboratory precisely guides, amplifies, reflects, and focuses 192 powerful laser beams into a target about the size of a pencil eraser in a few billionths of a second, delivering more than 2 million joules of ultraviolet energy and 500 trillion watts of peak power. A crucial goal of the system is to trigger precise implosions of fuel capsules. This is achieved by delivering all 192 beams at user-specified times and locations on the target, minimizing any deviation from the requested performance. Power requirements can vary substantially on each experiment and the facility supports numerous amplifier pumping configurations and their attendant nonlinear effects. To achieve the tight performance required across such a broad array of configurations, constant comparison of measured and requested power delivery are tracked and long-term trends analyzed as a guide to understanding future performance. A very common question asked is given the current state of the laser today – how well would a similar experiment from the past perform today? Additionally, if one were to specify an alternate amplifier configuration using today’s model would and damage limits be exceeded and would there be an increase in performance. The physics model used to make these predictions and equipment protection checks is called the Virtual Beamline (VBL). VBL is used with an incoming desired pulse shape and energy to be delivered on target, and then does an iterative solve to predict the needed injected pulse in the front-end of the system to achieve this result. To effectively guide and predict future NIF experiment performance, laser scientists explore current amplifier configurations and compare them with historical data utilizing a tool called the Reverify Toolbox.

42 ENGINEERING↗

Feasible Actuator Range Modifier (FARM), a Tool Aiding the Solution of Unit Dispatch Problems for Advanced Energy Systems

Integrated energy systems (IESs) seek to minimize power generating costs in future power grids through the coupling of different energy technologies. To accommodate fluctuations in load demand due to the penetration of renewable energy sources, flexible operation capabilities must be fully exploited, and even power plants that are traditionally considered as base-load units need to be operated according to unconventional paradigms. Thermomechanical loads induced by frequent power adjustments can accelerate the wear and tear. If a unit is flexibly operated without respecting limits on materials, the risk of failures of expensive components will eventually increase, nullifying the additional profits ensured by flexible operation. In addition to the bounds on power variations (explicit constraints),the solution of the unit dispatch problem needs to meet the limits on the variation of key process variables, including temperature, pressure and flow rate (implicit constraints).The FARM (Feasible Actuator Range Modifier) module was developed to enable existing optimization algorithms to identify solutions to the unit dispatch problem that are both economically favorable and technologically sustainable. Thanks to the iterative dispatcher–validator scheme, FARM permits addressing all the imposed constraints without excessively increasing the computational costs. In this work, the algorithms constituting the module are described, and the performance was assessed by solving the unit dispatch problem for an IES composed of three units, i.e., balance of plant, gas turbine, and high-temperature steam electrolysis. Finally, the FARM module provides dedicated tools for visualizing the response of the constrained variables of interest during operational transients and a tool aiding the operator at making decisions. These techniques might represent the first step towards the deployment of an ecological interface design (EID) for IES units.

47 OTHER INSTRUMENTATION↗

Multigrid Acceleration of Time-Accurate Navier-Stokes Calculations

A numerical scheme to solve the unsteady Navier-Stokes equations is described. The scheme is fully implicit in time and is unconditionally stable (at least for first- and second-order discretizations of the physical time derivatives). With unconditional stability, the choice of the time step is based on the physical phenomena to be resolved rather than limited by numerical stability. This is especially important for high Reynolds number viscous flows, where the spatial variation of grid cell size can be as much as six orders of magnitude. A multigrid-multiblock, steady-state, three-dimensional Navier-Stokes solver, TLNS3D, was modified to iteratively invert the equations at each physical time step. The implementation of this procedure in TLNS3D is discussed. The implications of applying several popular turbulence models to unsteady flow are also considered. Numerical results are presented to show the application of the scheme to various two-dimensional turbulent flows. The results of a three-dimensional laminar flow calculation are also given.

Melson, N. Duane↗

Spacecraft Orbits for Earth/Mars-Lander Radio Relay

A report discusses a network of spacecraft, in orbit around Mars, used to relay radio communications between Earth stations and mobile exploratory robots (rovers) as well as stationary scientific instruments that have been landed on the Mars surface. The relay spacecraft include two already in orbit plus several others planned to arrive at Mars in the years 2004 through 2008. A major portion of the report is devoted to the orbit of the G. Macroni Orbiter, which is in the midst of an iterative design process and is intended to be the first Mars orbiter designed primarily for radio relay. Candidate orbits are analyzed with a view toward choosing one that maximizes the amount of time available for communication with surface units, taking account of visibility as a function of position, the limit on communication distance at low power, and the fact that surface units can transmit more easily when they are in sunlight. Two promising new orbits for Mars relay satellites are identified: a 1/2-sol apoapsis-at-constant-time-of-day equatorial orbit and a 1/4-sol apoapsis-at-constant-time-of-day, critical-inclination orbit.

Noreen, Gary↗

CFD Assessment of Forward Booster Separation Motor Ignition Overpressure on ET XT 718 Ice/Frost Ramp

Computational fluid dynamics assessment of the forward booster separation motor ignition over-pressure was performed on the space shuttle external tank X(sub T) 718 ice/frost ramp using the flow solver OVERFLOW. The main objective of this study was the investigation of the over-pressure during solid rocket booster separation and its affect on the local pressure and air-load environments. Delta pressure and plume impingement were investigated as a possible contributing factor to the cause of the debris loss on shuttle missions STS-125 and STS-127. A simplified computational model of the Space Shuttle Launch Vehicle was developed consisting of just the external tank and the solid rocket boosters with separation motor nozzles and plumes. The simplified model was validated by comparison to full fidelity computational model of the Space Shuttle without the separation motors. Quasi steady-state plume solutions were used to calibrate the thrust of the separation motors. Time-accurate simulations of the firing of the booster-separation motors were performed. Parametric studies of the time-step size and the number of sub-iterations were used to find the best converged solution. The computed solutions were compared to previous OVERFLOW steady-state runs of the separation motors with reaction control system jets and to ground test data. The results indicated that delta pressure from the overpressure was small and within design limits, and thus was unlikely to have contributed to the foam losses.

Tejnil, Edward↗

Iterative redox activation promotes interfacial synergy in an Ag/Cu x O catalyst for oxygen reduction

In this study, an iterative electrochemical activation process is demonstrated to maximize the interface between Ag and Cu x O. Morphological reconstruction and phase transformation of Cu x O occurs through the iterative electrochemical redox reaction, as indicated by energy-dispersive spectroscopy and X-ray photoelectron spectroscopy. The activated Ag/Cu x O/C catalysts exhibits enhanced oxygen reduction reaction performance, with an onset potential of 0.86 V vs. the reversible hydrogen electrode, a Tafel slope of 46mV dec -1 , and a high stability as compared to Ag/C. A mechanistic study using density functional theory shows that the weakened binding energy of the OH intermediate, originating from the charge transfer from Ag to Cu x O, improves the ORR activity at the interface sites of the Ag/Cu x O electrocatalysts. Ag/Cu x O shows a higher limiting potential of 0.14V for interface sites than isolated Ag nanoparticles. The interfacial charge transfer between Ag and CuO is verified experimentally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A new tubular ironless linear induction motor using non-interrupted winding method

The ironless linear induction motor device can be used for any application which requires a large motion span or rotation restrictedly by the ironless feature due to the strong magnetic field environment. One of the possible applications is ITER diagnostics shutters. The tokamak reactor systems have a severe radiation issue during the operation. The traditional Kapton-insulated wires have significantly lower radiation resistance than an expected large radiation dose at ITER Motion Stark Effect (MSE) shutter location. Furthermore, the traditional single-sided linear induction motor winds the wire to make the N-turn pitch coils independently and then connect them in series to make a one-phase winding. Therefore, it is a big challenge to make so many radiation-protection connectors in such a limited space as inside the vacuum vessel of fusion reactors. The approved mineral-insulated kinds of cables need a large bend radius for the end turns, so the traditional single-sided linear motor with two-layer three-phase windings need to be upgraded to the tubular linear motors. The new proposed tubular ironless linear induction motor can not only overcome the large bend radius issue with the mineral-insulated cables but also allow the motor coils wound continuously around the tubular-shaped stator, forward and backward N times to make a N-turn winding for each phase, in order to resolve the radiation protection issue with the winding connection. The 3D transient Maxwell electromagnetic model with a tubular three-phase linear induction motor has been analyzed. Furthermore, this paper will present the detailed design and winding method of ironless linear induction motor actuators, and results of the 3D transient numerical simulation of the actuator performance.

Ironless↗

Impact of collisionality on turbulence in the edge of tokamak plasma using 3D global simulations

Collisionality is one of the key parameters in determining turbulent transport in the plasma edge, regulating phenomena such as 'shoulder formation', separation of scale lengths in the scrape-off layer (SOL), turbulence damping and zonal flow dynamics. Understanding its role is therefore of primary importance for future reactors like ITER. Getting reliable predictions and a better characterization of plasma flow properties when varying collisionality remains, however, a critical challenge for the simulations. This paper focuses on the impact of varying collisionality in a non-isothermal three-dimensional fluid model of the plasma edge. A high field side limited configuration encompassing open and closed magnetic field lines with parameters typical of a medium-sized tokamak is considered. The present model can consistently account for the variations of collisionality and its impact on both the parallel resistivity η ∥ and the ion and electron parallel thermal conductivities χ ∥e,i . Details on mean flow and turbulence properties are given. Changing collisionality leads to significant changes in the flow properties both on the mean and fluctuating quantities. In particular, lowering collisionality decreases the size of coherent structures, the fluctuation levels of turbulence, and steepens the density and temperature equilibrium profiles around the separatrix leading to a global reduction of the turbulent transport. The SOL width is observed to increase with collisionality, eventually resulting in the disappearance of the scale lengths separation between near and far SOL, consistently with previous experimental observations. At low collisionality, where the presence of narrow feature is well-established, a contribution of heat conduction increases up to compete with heat convection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

De Novo Design of a Highly Stable Ovoid TIM Barrel: Unlocking Pocket Shape towards Functional Design

The ability to finely control the structure of protein folds is an important prerequisite to functional protein design. The TIM barrel fold is an important target for these efforts as it is highly enriched for diverse functions in nature. Although a TIM barrel protein has been designed de novo, the ability to finely alter the curvature of the central beta barrel and the overall architecture of the fold remains elusive, limiting its utility for functional design. Here, we report the de novo design of a TIM barrel with ovoid (twofold) symmetry, drawing inspiration from natural beta and TIM barrels with ovoid curvature. We use an autoregressive backbone sampling strategy to implement our hypothesis for elongated barrel curvature, followed by an iterative enrichment sequence design protocol to obtain sequences which yield a high proportion of successfully folding designs. Designed sequences are highly stable and fold to the designed barrel curvature as determined by a 2.1 Å resolution crystal structure. The designs show robustness to drastic mutations, retaining high melting temperatures even when multiple charged residues are buried in the hydrophobic core or when the hydrophobic core is ablated to alanine. As a scaffold with a greater capacity for hosting diverse hydrogen bonding networks and installation of binding pockets or active sites, the ovoid TIM barrel represents a major step towards the de novo design of functional TIM barrels.

Chu, Alexander E.↗

Dislocation dynamics during the growth of silicon ribbon

The thermal viscoplastic stresses and the dislocation densities in silicon ribbon are computed for an axially changing thermal profile by using an iterative finite difference method. A material constitutive equation (Haasen-Sumino model) which involves an internal variable (mobile dislocation density) is used. The results are interpreted as showing that there is a maximum width of silicon ribbon that can be grown when viscoplasticity and dislocations are considered. This maximum width limitation does not exist if the material behavior is elastic.

Dillon, O. W., Jr.↗

Batched Sparse Linear Algebra (Final Report for Subcontract B648960)

This report finalizes design specifications for developing batched kernels for small tensor operations for unassembled matrix-free iterative solvers, batched solvers for partially assembled operators, and batched solvers with support for various sparse formats. The outcome of the project milestones is a set of interfaces to Batched Sparse LA solvers running on hardware accelerators for use in ECP Libraries and Applications. It is part of the development of sparse batched kernels, solvers/preconditioners as well as creating interoperability in xSDK libraries with sparse and dense batched functions to benefit ECP applications. The participants included representatives from ECP libraries (not limited to the xSDK project), applications, and vendors (AMD, Intel, and NVIDIA). Batched sparse linear algebra solvers form the new frontier for algorithmic development and performance engineering. Many applications (ECP and non-ECP alike) require simultaneous solutions of small linear systems of equations that are structurally sparse. To move towards high hardware utilization, it is important to provide these applications with appropriate interfaces to efficient batched sparse solvers running on modern hardware accelerators. We present interface designs in use by HPC software libraries supporting batched sparse linear algebra and the development of sparse batched kernel codes for solvers and preconditioners. We also address the potential interoperability opportunities to keep the software portable between the major hardware accelerators from AMD, Intel, and NVIDIA. The presented interface specifications includes batched band, sparse iterative, and sparse direct solvers. This report summarizes progress in Kokkos Kernels and the xSDK libraries MAGMA, Ginkgo, hypre, SUNDIALS, and SuperLU_dist.

97 MATHEMATICS AND COMPUTING↗

Bifurcation-like transition of divertor conditions induced by X-point radiation in KSTAR L-mode plasmas *

Abstract Density ramps with ion grad B drift directed into lower single null KSTAR L-mode plasmas are associated with a simultaneous and abrupt reduction of the divertor particle flux on both low- and high-field-side targets when the mid-plane line averaged electron density reaches a given level. Target embedded Langmuir probe signals show a clear ‘cliff edge’ behavior similar to that observed in the divertor target electron temperature in DIII-D H-mode plasmas (Eldon et al 2017 Nucl. Fusion 57 066039; McLean et al 2015 J. Nucl. Mater. 463 533–6). The collapse of the particle flux is observed along the whole divertor target area (from private flux region to the far scrape-off layer (SOL)). The critical upstream density of this target flux cliff is invariant under fuel gas throughput modulation. The transition along the cliff occurs in tens of milliseconds. With the cliff, carbon impurities and deuterium neutrals transported through the X-point to the core produce a strong radiation spot near the X-point, seen on bolometric signals, and increase the upstream density. The experimental observations are consistent with time-dependent SOLPS-ITER simulations, which also demonstrate an abrupt transition of the target flux and upstream density with the increase in X-point radiation. The timescale of the cliff predicted by SOLPS-ITER is consistent with the experiment, although, it is influenced by gas throughput or time-dependent numerical methods. In the L-mode phase space of separatrix electron density and temperature, branches are divided based on target temperature, because the latter is strongly coupled to the radiation front and ionization front due to the monotonic characteristic of the parallel electron temperature distribution. Since the H-mode condition operates at a much higher upstream density and electron temperature in phase space, dissipation from sputtered carbon alone leads to the density limit before reaching the X-point radiation condition. This is therefore consistent with the fact that cliffs have never been observed in H-mode KSTAR experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A universal inequality on the unitary 2D CFT partition function

We prove the conjecture proposed by Hartman, Keller and Stoica (HKS) [1]: the grand-canonical free energy of a unitary 2D CFT with a sparse spectrum below the scaling dimension $\frac{c}{12}$ + ϵ and below the twist $\frac{c}{12}$ is universal in the large c limit for all β L β R ≠ 4π 2 . The technique of the proof allows us to derive a one-parameter (with parameter α ∈ (0, 1]) family of universal inequalities on the unitary 2D CFT partition function with general central charge c ⩾ 0, using analytical modular bootstrap. We derive an iterative equation for the domain of validity of the inequality on the (β L , β R ) plane. The infinite iteration of this equation gives the boundary of maximal-validity domain, which depends on the parameter α in the inequality.

AdS-CFT Correspondence↗