Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “alternating minimization”

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 379 records · Page 21

Investigation of long-term light stability of negative charge injected into oxide-nitride-oxide passivation stack of crystalline silicon solar cells

A negatively charged oxide-nitride-oxide stack for field-effect passivation of crystalline silicon solar cells is discussed. The negative charge was injected into the stack by a plasma charge injection technology. Charge stability was studied by exposing samples to AM1.5 simulation visible light and full-spectrum light at temperatures ranging from 55 to 78 °C for up to 300 h. Charge injection and loss were quantified based on shifts in the flatband voltage of capacitance–voltage curves measured with a mercury probe. The most probable mechanism of charge loss was found to be diffusion of negative charged hydrogen atoms through nitride and bottom oxide. The optimum recipe for each layer of the stack was investigated to minimize the loss of injected charge. The flatband voltage decay of the optimized stack was found to fit a power-law trend, suggesting the dispersive transport of hydrogen atoms with a dispersion parameter of ~0.06–0.07. The optimized stack is projected to maintain a negative charge density of about 3.6 × 10 12 cm –2 or more after 25 years of field operation in an environment such as Arizona, which would be sufficient for field-effect passivation under one-sun illumination. Furthermore, the high stability of the negative injected charge makes the plasma charging technology a safer and lower cost alternative to Al 2 O 3 -passivation technology commonly used to passivate p-type surfaces.

14 SOLAR ENERGY↗

Quarkyonic matter equation of state in beta-equilibrium

Quark matter may appear due to a hadronic-quark transition in the core of a hybrid star. Quarkyonic matter is an approach in which both quarks and nucleons appear as quasiparticles in a crossover transition, and provides an explicit realization of early ideas concerning quark matter (e.g., the MIT bag model). This description has recently been employed by McLerran and Reddy to model chargeless (pure neutron) matter with an approach that has the virtue that the speed of sound rises quickly at a neutron-quark transition so as to satisfy observational constraints on the neutron star maximum mass (≳ 2 $M_⊙$) and the radius of a 1.4 $M_⊙$ star ($R_{1.4}$ ≲ 13.5 km). Traditional models involving first-order transitions result in softer pressure-energy density relations that have difficulty satisfying these constraints except with very narrow choices of parameters. Here, we propose a variation of quarkyonic matter involving protons and leptons whose energy can be explicitly minimized to achieve both chemical and beta equilibrium, which cannot be done in the chargeless formulation. Quarkyonic stellar models are able to satisfy observed mass and radius constraints with a wide range of model parameters, avoiding the obligatory fine-tuning of conventional hybrid star models, including requiring the transition density to be very close to the nuclear saturation density. Our formulation fits experimental and theoretical properties of the nuclear symmetry energy and pure neutron matter, and contains as few as three free parameters. This makes it an ideal tool for the study of high-density matter that is an efficient alternative to piecewise polytrope or spectral decomposition methods.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Improving Gas Reactor Design with Complex Non-Standard Reaction Mechanisms in a Reactive Flow Model

Fluorinated gases are critical to a number of high-technology industries including semiconductor manufacturing and related processes such as flat panel display, LED, and photovoltaic manufacturing. These fluorinated gases have very high global warming potentials (GWPs) and, unfortunately are chemically very stable, in some cases lasting for thousands of years in the atmosphere. For environmental reasons, it is imperative that the release of these gases be minimized, but existing methods to destroy them requires a significant amount of energy. The goal of this project was to develop a more energy efficient method of destroying these high GWP gases using computational fluid dynamics, but the computational requirements are significant. In this project, we bring together resources available at the national labs, including high-performance computing hardware and modern computational methods, in order to explore an industrially important problem with key environmental and economic impact. In this project, we examined the state-of-the-art models for the physical processes involved, and made subtle changes to arrive at alternative formulations considerably better suited for computation and study. We examined several representative scenarios and developed a set of analyses techniques and tools that could be used to explore the systems in much greater detail. This exploratory project has demonstrated to ALZETA that HPC can indeed provide a unique view into their complex engineering problem - one that is simply inaccessible by any other means. The tools allow us to explore the details of the chemical process, and allow one to construct "what-if" scenarios to address shortcomings of the physical device in its current operating modes. However, as with many sources of completely new information, many new questions are illuminated with this capability. It will take considerable time, effort and experience before we learn to probe the solutions, pose operating scenarios and confirm modifications to the device that lead to real engineering advances.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

PWM for Simultaneous Minimization of Switching Losses and CM Voltage in Current Source Inverters

This article proposes a new pulsewidth modulation (PWM) scheme that reduces common-mode (CM) voltage and switching losses simultaneously in CSI. The proposed PWM [i.e., reduced voltages PWM (RVM)] reduces CM voltages and switching losses by selecting the optimal zero current vector and PWM sequence based on the CSI output voltage conditions. This article also proves the proposed RVM achieves both objectives (i.e., reduced CM voltage and switching losses) simultaneously for all CSI operating conditions without compromising the other objective. The proposed RVM inherently achieves zero-voltage switching (ZVS) conditions during the zero current vector transitions, further reducing the switching losses. Experimental results demonstrate significant reductions in CM voltage, reaching up to approximately 95%, accompanied by an increase of up to about 5.5% in CSI efficiency across a wide operating range compared to alternative PWM methods, confirming the predicted advantages of the proposed RVM.

42 ENGINEERING↗

Life cycle cost, energy, and carbon emissions of molds for precast concrete: Exploring the impacts of material choices and additive manufacturing

Molds for precast concrete are commonly used to create simple- to complex-shaped concrete products away from construction sites. These molds are often handmade from wood; however, additively manufacturing (AM, or 3D printing) fiber-reinforced polymer composites is an advantageous alternative, producing significantly more durable, highly complex molds faster, but likely at a higher cost. Here, this study explores the impact of material and production variables on the cost, energy, and carbon emissions of employing composite AM molds over the full lifecycle. The case study employed techno-economic and life cycle assessments to show that using wood flour–poly(lactic acid) or recycled carbon fiber–acrylonitrile butadiene styrene for AM molds can be less expensive than conventional wood molds, especially when considering use phase costs. While wood molds have the least environmental impacts due to wood's higher biogenic carbon sequestration and minimal processing, optimizing AM designs could reduce energy demand, carbon emissions, and cost.

36 MATERIALS SCIENCE↗

Shadow energy functionals and potentials in Born–Oppenheimer molecular dynamics

In Born–Oppenheimer molecular dynamics (BOMD) simulations based on the density functional theory (DFT), the potential energy and the interatomic forces are calculated from an electronic ground state density that is determined by an iterative self-consistent field optimization procedure, which, in practice, never is fully converged. The calculated energies and forces are, therefore, only approximate, which may lead to an unphysical energy drift and instabilities. Here, we discuss an alternative shadow BOMD approach that is based on backward error analysis. Instead of calculating approximate solutions for an underlying exact regular Born–Oppenheimer potential, we do the opposite. Instead, we calculate the exact electron density, energies, and forces, but for an underlying approximate shadow Born–Oppenheimer potential energy surface. In this way, the calculated forces are conservative with respect to the approximate shadow potential and generate accurate molecular trajectories with long-term energy stabilities. We show how such shadow Born–Oppenheimer potentials can be constructed at different levels of accuracy as a function of the integration time step, δt, from the constrained minimization of a sequence of systematically improvable, but approximate, shadow energy density functionals. For each energy functional, there is a corresponding ground state Born–Oppenheimer potential. These pairs of shadow energy functionals and potentials are higher-level generalizations of the original “zeroth-level” shadow energy functionals and potentials used in extended Lagrangian BOMD. The proposed shadow energy functionals and potentials are useful only within this extended dynamical framework, where also the electronic degrees of freedom are propagated as dynamical field variables together with the atomic positions and velocities. The theory is quite general and can be applied to MD simulations using approximate DFT, Hartree–Fock, or semi-empirical methods, as well as to coarse-grained flexible charge models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Energy efficiency and renewable energy utilization in desalination systems

Desalination technologies and industry have advanced significantly in the past two decades to meet the growing freshwater demands stimulated by the compounding issues of both water quality and quantity in many regions of the world. As desalination processes are energy demanding, there have been many efforts dedicated to improve energy-efficiency of the process units, enhance energy conservation and recovery, and increase renewable energy integration in desalination plants. This review article highlights recent key advances and discusses possible venues for further development in desalination energy portfolio to reduce specific energy consumption and, via integration with solar energy, to minimize the environmental footprint associated with freshwater production. First, an overview of current desalination technologies and their energy requirements are presented followed by a discussion on opportunities for improving energy efficiency and energy recovery in both membrane and thermal desalination technologies. Then, various combinations of renewable energy driven desalination plants are discussed with some recent highlights in solar energy driven membrane, thermal and hybrid desalination processes. Technological readiness levels for novel desalination processes, their perceived impact and expected near-future developments in renewable energy integrated desalination technologies are presented. Lastly, the potential for solar driven desalination as a cost-competitive freshwater supply alternative is discussed

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Distributed Optimization in Distribution Systems: Use Cases, Limitations, and Research Needs

We report electric distribution grid operations typically rely on both centralized optimization and local non-optimal control techniques. As an alternative, distribution system operational practices can consider distributed optimization techniques that leverage communications among various neighboring agents to achieve optimal operation. With the rapidly increasing integration of distributed energy resources (DERs), distributed optimization algorithms are growing in importance due to their potential advantages in scalability, flexibility, privacy, and robustness relative to centralized optimization. Implementation of distributed optimization offers multiple challenges and also opportunities. This paper provides a comprehensive review of the recent advancements in distributed optimization for electric distribution systems and classifications using key attributes. Problem formulations and distributed optimization algorithms are provided for example use cases, including volt/var control, market clearing process, loss minimization, and conservation voltage reduction. Finally, this paper also presents future research needs for the applicability of distributed optimization algorithms in the distribution system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Inverse problem in the large momentum effective theory framework

One proposal to compute parton distributions from first principles is the large momentum effective theory (LaMET), which requires the Fourier transform of matrix elements computed nonperturbatively. Lattice quantum chromodynamics (QCD) provides calculations of these matrix elements over a finite range of Fourier harmonics that are often noisy or unreliable in the largest computed harmonics. It has been suggested that enforcing an exponential decay of the missing harmonics helps alleviate this issue. Using nonperturbative data, we show that the uncertainty introduced by this inverse problem in a realistic setup remains significant without very restrictive assumptions, and that the importance of the exact asymptotic behavior is minimal for values of 𝑥 where the framework is currently applicable. We show that the crux of the inverse problem lies in harmonics of the order of 𝜆 = 𝑧⁢𝑃 𝑧 ∼ 5–15, where the signal in the lattice data is often barely existent in current studies, and the asymptotic behavior is not firmly established. We stress the need for more sophisticated techniques to account for this inverse problem, whether in the LaMET or related frameworks like the short-distance factorization. We also address a misconception that, with available lattice methods, the LaMET framework allows a “direct” computation of the 𝑥-dependence, whereas the alternative short-distance factorization only gives access to moments or fits of the 𝑥-dependence.

Dutrieux, Hervé [Aix-Marseille Université, Marseil↗

Technical Report on Waveform Fit Metrics for Global Models

The new WAVEFORMS Initiative in the Ground-based Nuclear Detonation Detection (GNDD) program includes an increased emphasis on the development of Earth models and methods to predict entire seismic and acoustic waveforms more accurately. In general, this increased emphasis is predicated on the need to better characterize seismic events and provide improved model-based discrimination between event types including earthquakes and explosions. More specifically, while current moment tensor inversion methods tend to work well for larger events (M>~4) using tuned 1-D Earth models, the development of state-of-the-art 3-D models and methods is required for the prediction of shorter period waves over large areas for discrimination of smaller events. There is no standard metric for model-based waveform prediction accuracy used in the waveform modeling/inversion community. However, there are several popular waveform misfit definitions; and minimizing the corresponding objective functions is the goal of waveform inversion. Some example misfit definitions employed for adjoint waveform tomography include measures of simple travel time differences (e.g. Tape et al., 2010), cross-correlation travel time differences (e.g. Luo and Schuster, 1991), multi-taper frequency dependent methods (e.g. Lei et al., 2020), time-frequency phase misfit functions (e.g. Fichtner 2010; Rodgers et al., 2022), normalized cross-correlation methods (e.g. Tao et al., 2018), and others. In some cases, these misfit definitions also involve complicated weighting schemes and summations over multiple frequency bands making it difficult to duplicate the misfit measurement with alternative models and datasets. Although each of the misfit definitions mentioned above are useful for developing waveform models, the actual misfit values are not usually meaningful outside of a given project, model, and/or dataset. Therefore, it is difficult to understand and communicate model performance for predicting waveforms and comparing to other models and/or new model iterations with a different dataset. Therefore, there is a need for a generalized method for evaluating overall model performance that is independent from the specific misfit chosen to develop the waveform models that is also intuitive and meaningful. In this report, we describe a new metric we refer to as ‘Percent of Correlated Signal’. The following sections describe and demonstrate the metric with a case study event and a more rigorous test using a random selection of globally distributed events. While the focus here is on global tomography models, the metric is meant to applicable to regional ‘wiggle-for-wiggle’ waveform models/studies as well.

58 GEOSCIENCES↗

Opportunities for Hybrid Wind and Solar PV Plants in India

The share of variable renewable energy (VRE) on India's grid has surpassed 100 GW, and the government has ambitious plans reach 450 GW by 2030. One strategy to increase wind and solar PV deployment is through the co-location of wind and solar PV plants to form a single hybrid power plant. Hybrid plants have the potential to reduce transmission infrastructure costs and variability in the output power profile compared to a standalone plant with a single technology, and this resource analysis aims to take a first step towards quantifying the potential savings from hybridizing wind and solar PV plants in India and the size of this opportunity. We utilize a brute-force optimization to minimize the levelized cost of energy (LCOE) for standalone wind, standalone solar PV, and hybrid wind/solar PV plants across all of India. By comparing these LCOEs, we determine that locations where hybrid plants exhibit potential cost savings and grid benefits exhibit both; a high interconnection cost and; a wind capacity factor between roughly 34% and 38%. However, because our work does not capture the value of the electricity generated by looking at energy prices, nor does it quantify the potential of hybrids to provide other value streams such as firm capacity and reserves. Further, because the work does not compare solar PV and wind hybrids to alternative generation technologies or storage systems, it cannot be considered a holistic cost-benefit analysis.

14 SOLAR ENERGY↗

Constrained variational optimization of counting-time allocation in sequential scattering measurements: Application to Bonse–Hart USANS

Sequential scattering measurements are often performed under a fixed experimental-time budget, even though the expected count rate varies strongly across the measured coordinate. When the dwell time at each measurement position can be controlled independently, this variation creates a general resource-allocation problem: how should the available time be distributed to minimize the uncertainty of the reconstructed profile? We formulate this problem as a constrained variational optimization for measurements governed by Poisson counting statistics. When each measurement is treated independently, minimizing the averaged squared relative uncertainty yields an inverse-square-root intensity allocation. The formulation is then generalized to include correlations between neighboring measurements and an instrumental resolution operator, leading to an allocation criterion that equalizes the marginal reduction in posterior uncertainty per unit measurement time. Bonse–Hart ultra-small-angle neutron scattering (USANS), in which reciprocal space is sampled sequentially through analyzer-angle stepping, provides an experimentally grounded application. Computational benchmarking shows that the optimized allocation outperforms uniform-time and constant-relative-error strategies, while application to an experimentally measured graphite USANS profile from the Spallation Neutron Source, using Poisson resampling under alternative schedules, demonstrates how counting time should be redistributed toward weak-intensity regions under an identical total duration. The resulting framework applies to sequential scattering and related scanning measurements whenever local dwell times are adjustable and directly determine the measurement uncertainties, and when the relevant correlation and instrumental-response models are available.

Tung, Chi-Huan [ORNL] (ORCID:0000000221972074)↗

Efficient 3-D velocity model building using joint inline and crossline plane-wave wave-equation migration velocity analyses

SUMMARY Wave-equation migration velocity analysis (WEMVA) is an image-domain inversion method for velocity model building. Automatic plane-wave WEMVA (PWEMVA) calculates the moveouts of plane-wave common-image gathers (CIGs) by searching a best-fitting parabola with semblance analysis and backprojects residual CIG moveouts into wavefield wave paths with a reflection tomographic kernel. However, 3-D PWEMVA is very computationally expensive because 3-D reflection tomographic inversion requires at least five 3-D reverse-time migrations per iteration and stores two types of source wavefields at model boundaries. We develop a joint inline and crossline PWEMVA method for efficient 3-D velocity model building. We alternatively implement the inline and crossline PWEMVAs with a constraint for each other, in which we iteratively construct the 3-D velocity model update through 1-D spline interpolation of 2-D gradients. The inline and crossline joint inversion is practical since PWEMVA only inverts for low-wavenumber velocity perturbations along wave paths, and the method can take less than 1 per cent of the computational cost of full 3-D PWEMVA. To construct unaliased plane waves for our joint inline and crossline PWEMVA, we develop a 3-D data interpolation method in the frequency–wavenumber (FK) domain to recover regularly and randomly missing traces. The method minimizes the misfit on sufficiently localized data subsets with iterative optimal step lengths and a gradient preconditioner that iteratively selects dominant dips along different azimuths. In numerical experiments, we use a 3-D synthetic seismic data set and a land 3-D field seismic data set acquired at the Farnsworth CO2-EOR (enhanced oil recovery) field to demonstrate the efficacy of our velocity model building and data interpolation methods.

Liu, Xuejian↗

A Business Case Evaluation of Gas Switching Reforming (GSR) Technology: A Promising Technology for Natural Gas Reforming with Integrated CO2 Capture

Hydrogen is essential in the transition to sustainable energy, and developing low-carbon production methods is a key research focus. Traditional steam methane reforming (SMR) dominates the hydrogen industry but contributes substantially to CO2 emissions. In response, Gas Switching Reforming (GSR) has emerged as a novel process that integrates carbon capture and utilizes process heat more efficiently. Unlike other reforming methods, GSR consolidates oxidation and reduction reactions within a single reactor, which minimizes external energy inputs and simplifies scaling. Like conventional steam methane reforming (SMR), GSR can be integrated with water-gas shift and pressure swing adsorption units for pure hydrogen production. This work presents a comprehensive business case analysis of GSR technology based on experimental results in Technology Readiness Level 3, Life Cycle Assessment (LCA) and Techno-Economic (TEA) evaluation incorporating ASPEN Plus process modeling considering different configurations and energy scenarios. The TEA incorporates data from kinetic experiments from various catalysts to evaluate the GSR process under various conditions. The goal of this work is to evaluate GSR’s potential to serve as a low-carbon alternative to SMR, focusing on global warming potential and additional impact categories to evaluate a wide spectrum of environmental impacts. Comparative assessments were conducted with SMR, chemical loop reforming (CLR), and proton exchange membrane (PEM) electrolysis to explore trade-offs across environmental metrics. The environmental impact assessment of this work encompasses the entire hydrogen production lifecycle from raw material extraction to plant decommissioning, using a cradle-to-gate boundary. Preliminary findings highlight that GSR, when integrated with low-carbon energy sources, could significantly reduce environmental impacts, making it a promising candidate for low-carbon hydrogen infrastructure. The insights from this business case evaluation aim to guide industry in scale-up and commercialization of this promising clean energy technology.

03 NATURAL GAS↗

Graphite Intercalation Compounds Derived by Green Chemistry as Oxygen Reduction Reaction Catalysts

Precious group metal (PGM) catalysts such as Pt supported on carbon supports are expensive catalysts utilized for the oxygen reduction reaction (ORR) due to their unmatched catalytic activity and durability. As an alternative, PGM-free ORR electrocatalysts that offer respectable catalytic activity are being pursued. Most of the notable PGM-free catalysts are obtained either from a bottom-up approach synthesis utilizing nitrogen-rich polymers as building blocks, or from a top down approach, where nitrogen and metal moieties are incorporated to carbonaceous matrixes. The systematic understanding of the origin of catalytic activity for either case is speculative and currently employed synthesis techniques typically generate large amounts of hazardous waste such as acids, oxidizing agents, and solvents. Herein, for the first time, we investigate the catalytic activity of graphite-based materials obtained via intercalation strategies that minimally perturb the graphitic backbone. Here, our outlined approaches demonstrate initial efforts to not only elucidate the role of each element but also significantly reduce the use of hazardous chemicals, which remains a pressing challenge. Graphite intercalation compounds (GIC) were obtained using fewer steps and solvent-free processes. X-ray diffraction and Raman results confirm the successful intercalation of FeCl 3 between graphite layers. Electrochemical data shows that the ORR performance of FeCl 3 -intercalated GIC displays slight improvement where the onset potential reaches 0.77 V vs RHE in alkaline environments. However, expansion of the graphite and solvent-free incorporation of iron and nitrogen moieties resulted in a significant increase in ORR activity with onset potential to 0.89 V vs RHE, a maximum half-wave of 0.72 V vs RHE, and a limiting current of about 2.5 mA cm –2 . We anticipate that the use of near solvent-free processes that result in a high yield of catalysts along with the fundamental insight into the origin of electrochemical activity will tremendously impact the methodologies for developing next-generation ORR catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An improved guess for the variational calculation of charge-transfer excitations in large systems

Ab initio quantum-chemical methods that perform well for computing the electronic ground state are not straightforwardly transferable to electronically excited states, particularly in large molecular systems. Wave function theory offers high accuracy, but is often prohibitively expensive. Methods based on time-dependent density functional theory (TD-DFT) are crucially sensitive to the chosen exchange-correlation functional (XCF) parameterization, and system-specific tuning protocols were therefore proposed to address the method's robustness. Methods based on the variational relaxation of the excited-state electron density showcased promising results for the calculation of charge-transfer excitations, but the complex shape of the electronic hypersurface makes convergence to a specific excited state much more difficult than for the ground state when standard variational techniques are applied. We address the latter aspect by providing suitable initial guesses, which we obtain by two separate constrained algorithms. Combined with the squared-gradient minimization algorithm for all-electrons relaxation in a freeze-and-release scheme (FRZ-SGM), we demonstrate that orbital-optimized density functional theory (OO-DFT) calculations can reliably converge to the charge-transfer states of interest even for large molecular systems. We test the FRZ-SGM method on a phenothiazine-anthraquinone CT excitation in a supramolecular Pd(II) coordination cage complex as a function of the cage conformation. This compound has been studied experimentally prior to our work. We compare this freeze-and-release scheme to two XCF reparameterizations, which were recently proposed as low-cost TD-DFT-based alternatives to variational methods. Two dye-semiconductor complexes, which were previously investigated in the context of photovoltaic applications, serve as a second example to investigate the convergence and stability of the FRZ-SGM approach. Our results demonstrate that FRZ-SGM provides reliable convergence for charge-transfer excited states and avoids variational collapse to lower-lying electronic states, whereas time-dependent DFT calculations with an adequate tuning procedure for the range-separation parameter provide a computationally efficient initial estimate of the corresponding energies, with a computational cost comparable to that of configuration-interaction singles (CIS) calculations.

Bogo, Nicola↗

Hybrid Quantum-Classical Neural Networks

Deep learning is one of the most successful and far-reaching strategies used in machine learning today. However, the scale and utility of neural networks is still greatly limited by the current hardware used to train them. These concerns have become increasingly pressing as conventional computers are soon expected to approach the physical limitations that will slow their performance improvements in the near future. For these reasons, scientists have begun to explore alternative computing platforms, like quantum computers, for training neural networks. In recent years, variational quantum circuits have emerged as one of the most successful approaches to quantum deep learning on noisy intermediate scale quantum devices. We propose a hybrid quantum-classical neural network architecture where each neuron is a variational quantum circuit. We empirically analyze the performance of this hybrid neural network on a series of binary classification data sets using a simulated IBM universal quantum computer and a state-of-the-art IBM universal quantum computer. On the simulated hardware, we observe that the hybrid neural network achieves around 10% higher classification accuracy and 20% better minimization of the cost function than an individual variational quantum circuit. On the quantum hardware, we observe that each model only performs well when the qubit and gate count is sufficiently small.

Arthur, Davis↗