Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Radiative transfer simulations”

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 343 records · Page 19

Machine-learning-based inversion of nuclear responses

A microscopic description of the interaction of atomic nuclei with external electroweak probes is required for elucidating aspects of short-range nuclear dynamics and for the correct interpretation of neutrino oscillation experiments. Nuclear quantum Monte Carlo methods infer the nuclear electroweak response functions from their Laplace transforms. Inverting the Laplace transform is a notoriously ill-posed problem; and Bayesian techniques, such as maximum entropy, are typically used to reconstruct the original response functions in the quasielastic region. In this work, we present a physics-informed artificial neural network architecture suitable for approximating the inverse of the Laplace transform. Utilizing simulated, albeit realistic, electromagnetic response functions, we show that this physics-informed artificial neural network outperforms maximum entropy in both the low-energy transfer and the quasielastic regions, thereby allowing for robust calculations of electron scattering and neutrino scattering on nuclei and inclusive muon capture rates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

How Can Floor Covering Influence Buildings’ Demand Flexibility?

Although the thermal mass of floors in buildings has been demonstrated to help shift cooling load, there is still a lack of information about how floor covering can influence the floor’s load shifting capability and buildings’ demand flexibility. To fill this gap, we estimated demand flexibility based on the daily peak cooling load reduction for different floor configurations and regions, using EnergyPlus simulations. As a demand response strategy, we used precooling and global temperature adjustment. The result demonstrated an adverse impact of floor covering on the building’s demand flexibility. Specifically, under the same demand response strategy, the daily peak cooling load reductions were up to 20–34% for a concrete floor whereas they were only 17–29% for a carpet-covered concrete floor. This is because floor covering hinders convective coupling between the concrete floor surface and the zone air and reduces radiative heat transfer between the concrete floor surface and the surrounding environment. In hot climates such as Phoenix, floor covering almost negated the concrete floor’s load shifting capability and yielded low demand flexibility as a wood floor, representing low thermal mass. Sensitivity analyses showed that floor covering’s effects can be more profound with a larger carpet-covered area, a greater temperature adjustment depth, or a higher radiant heat gain. With this effect ignored for a given building, its demand flexibility would be overestimated, which could prevent grid operators from obtaining sufficient demand flexibility to maintain a grid. Our findings also imply that for more efficient grid-interactive buildings, a traditional standard for floor design could be modified with increasing renewable penetration.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Demonstration of low-mode shape control in indirect-drive double shell implosions at the NIF

Double shell inertial confinement fusion is a concept for achieving robust thermonuclear burn that uses dense metal shells to compress deuterium-tritium (DT) fuel to fusion conditions. Double shell implosions are typically indirectly driven and involve a target that consists of a low-Z ablator, a foam layer, and a high-Z pusher surrounding the DT fuel. The goal of the campaign is to achieve a volumetric burn as radiation losses from the DT fuel are trapped by the opaque high-Z shell. The overall performance of double shell implosions relies on the efficient collisional transfer of kinetic energy between layers. The efficiency of this transfer (and therefore the overall performance of a given implosion) is degraded by the presence of low-mode asymmetries. P2 asymmetries are often observed in spatially resolved 2D radiographs of nominal double shell implosions. This work discusses three such experiments: one with an oblate P2 asymmetry, one with a prolate P2 asymmetry, and one with an approximate spherical symmetry. After performing a shape analysis of the oblate and prolate implosions to quantify asymmetries, these experimental results were compared with the results of hydrodynamic simulations for the two experiments. Differences between the experiment and simulation were then used to design an approximately spherical implosion by altering the incident laser cone fraction. Radiographs from the experiment that implemented the modified cone fraction show evidence of an implosion that is approximately spherical until bang time. This design is intended to serve as a point design for future studies that will seek to optimize various aspects of the double shell target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of turbulence-radiation coupling in pressurized oxy-combustion

Concerns over climate change have led to numerous efforts in developing low-carbon energy technologies. Pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emission in power generation. Due to the strong impact of pressure on thermal radiation, heat transfer in POC differs significantly from the situation of conventional atmospheric pressure combustion. Thus thermal radiation in POC needs to be investigated to aid new combustor development. The present computational work is a step in this direction, initiating a systematic analysis of thermal radiation and heat transfer in a pilot-scale POC combustor, which has been developed at Washington University in St. Louis (WUSTL). In a POC process, pulverized coal is burned under elevated pressure and O2-CO2 environment. While most previous computational works on flame radiation focused on atmospheric pressure condition, this work considered a pressurized flame. Specifically, a 15-bar POC combustor of power 50 kW is modeled employing the Ansys FLUENT commercial platform, using both Reynolds-averaged Navier-Stokes (RANS) modeling and large-eddy simulation (LES). A recently published and validated global radiation model is used to predict the radiative property of the flue gas. The discrete ordinates (DO) radiation model is chosen to solve the radiative transport equation. Incident radiation on the walls of the combustor is identified and investigated. It is revealed that for this pilot-scale, pressurized combustor, thermal radiation exhibits similar patterns in the RANS and LES models. However, the magnitudes of radiation are different in both models, presumably because interactions between radiation and turbulence are embedded in the LES model but not in the RANS model. This difference thereby underlines the major goal of this research: investigating the impact of such turbulence-to-radiation coupling on thermal radiation in a pilot-scale POC setting. Most previous studies on turbulence-to-radiation interactions focused on bench-scale, atmospheric pressure flames. This work is the first effort to extend this study to a pilot-scale, pressurized flame.

Li, Lei↗

A Parameterization of Sub-grid Topographical Effects on Solar Radiation in the E3SM Land Model (Version 1.0): Implementation and Evaluation Over the Tibetan Plateau

Topography exerts significant influences on the incoming solar radiation at the land surface. A few stand-alone regional and global atmospheric models have included parameterizations for sub-grid topographic effects on solar radiation. However, nearly all Earth System Models (ESMs) that participated in the Coupled Model Intercomparison Project (CMIP6) use a plane-parallel (PP) radiative transfer scheme that assumes the terrain is flat. In this study, we incorporated a well-validated sub-grid topographic (TOP) parameterization in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) version 1.0 to quantify the effects of sub-grid topography on solar radiation flux, including the shadow effects and multi-scattering between adjacent terrain. The Moderate Resolution Imaging Spectroradiometer (MODIS) data was used to evaluate the performance of ELM. We studied the role of sub-grid topography by performing ELM simulations with the PP and TOP schemes over the Tibetan Plateau (TP). Additional ELM simulations were performed at multiple spatial resolutions to investigate the role of spatial scale on sub-grid topographic effects on solar radiation. When compared to MODIS data, incorporating the sub-grid topographic effects overall reduces the biases of ELM in simulating surface energy balance, snow cover and surface temperature especially in the high-elevation and snow-cover regions over the TP. Topography has non-negligible effects on surface energy budget, snow cover, and surface temperature over the TP. The absolute differences in surface energy fluxes for net solar radiation, latent heat flux, and sensible heat flux between TOP and PP exceed 20 W/m 2 , 10 W/m 2 , and 5 W/m 2 , respectively. The differences in land surface albedo, snow cover fraction, and surface temperature between TOP and PP exceed 0.1, 20%, and 1K, respectively. The magnitude of the sub-grid topographic effects is dependent on seasons and elevations, and is also sensitive to the spatial scales. Although the sub-grid topographic effects on solar radiation are more significant with more spatial details at finer spatial scales, they cannot be simply neglected at coarse spatial scales. The inclusion of sub-grid topographic effects on solar radiation parameterization in ELM will contribute to advancing our understanding of the role of the surface topography on terrestrial processes over complex terrain.

58 GEOSCIENCES↗

Improved Particle Heat Transfer by way of Bimodal Particle Distributions for High Temperature Solar Thermal Energy

High temperature solar thermal facilities are looking to increase operating temperatures through novel heat transfer media, one such being solid particles. These particles operating at high temperatures will require transferring their thermal energy into another working fluid like supercritical carbon dioxide which can be used in advanced power cycles. Achieving high heat transfer between the particles and supercritical carbon dioxide is essential to high efficiency and low-cost operation. Therefore, optimizing the thermal conductivity of these particles is one potential way to ensure high performance. Traditionally, unimodal particle distributions have been employed in high temperature particle solar power plants. However, ambient temperature testing of bimodal particle distributions has revealed a superior thermal conductivity when compared to its unimodal counterpart at the same temperature. This data was obtained by certified, off-the-shelf instruments that can effectively simulate the conditions a particle would be exposed to in a high temperature solar thermal system. Data obtained in this way suggests that the increased thermal conductively imputed by a bimodal particle distribution is significant at working temperatures in solar facilities. Furthermore, the thermal conductivity of these bimodal particle distributions peaks when the best combination of large and small particles is applied. At high temperatures, binary particle distributions are compared to monodispersed distributions of larger particles where heat transfer is more prolific due to the increased surface radiation. Various thermal conductivity, porosity and heat exchanger models are explored in conjunction with data acquired up to 700 C.

Stout, Dallin (ORCID:0009000294586091)↗

WRF-ELM v1.0: a regional climate model to study land–atmosphere interactions over heterogeneous land use regions

Abstract. The Energy Exascale Earth System Model (E3SM) Land Model (ELM) is a state-of-the-art land surface model that simulates the intricate interactions between the terrestrial land surface and other components of the Earth system. Originating from the Community Land Model (CLM) version 4.5, ELM has been under active development, with added new features and functionality, including plant hydraulics, radiation–topography interaction, subsurface multiphase flow, and more explicit land use and management practices. This study integrates ELM v2.1 with the Weather Research and Forecasting (WRF; WRF-ELM) model through a modified Lightweight Infrastructure for Land Atmosphere Coupling (LILAC) framework, enabling affordable high-resolution regional modeling by leveraging ELM's innovative features alongside WRF's diverse atmospheric parameterization options. This framework includes a top-level driver for variable communication between WRF and ELM and Earth System Modeling Framework (ESMF) caps for the WRF atmospheric component and ELM workflow control, encompassing initialization, execution, and finalization. Importantly, this LILAC–ESMF framework demonstrates a more modular approach compared to previous coupling efforts between WRF and land surface models. It maintains the integrity of ELM's source code structure and facilitates the transfer of future developments in ELM to WRF-ELM. To test the ability of the coupled model to capture land–atmosphere interactions over regions with a variety of land uses and land covers, we conducted high-resolution (4 km) WRF-ELM ensemble simulations over the Great Lakes region (GLR) in the summer of 2018 and systematically compared the results against observations, reanalysis data, and WRF-CTSM (WRF coupled with the Community Terrestrial Systems Model). In general, the coupled WRF-ELM model has reasonably captured the spatial distribution of surface state variables and fluxes across the GLR, particularly over the natural vegetation areas. The evaluation results provide a baseline reference for further improvements in ELM in the regional application of high-resolution weather and climate predictions. Our work serves as an example to the model development community for expanding an advanced land surface model's capability to represent fully-coupled land–atmosphere interactions at fine spatial scales. The development and release of WRF-ELM marks a significant advancement for the ELM user community, providing opportunities for fine-scale regional representation, parameter calibration in coupled mode, and examination of new schemes with atmospheric feedback.

54 ENVIRONMENTAL SCIENCES↗

Physicochemical and Molecular Insights into the Boundary Layer and Free Troposphere Aerosol Interactions over the Southern Great Plains

Ambient aerosols’ vertical profiles are critical for evaluating the role of aerosols in atmospheric chemistry and radiative transfer, but limited data on these profiles hinders our ability to fully assess their impact on the Earth's radiative balance. Here, in this study, we investigated the size-, time-, and altitude resolved composition of individual particles and bulk molecular composition of particle samples collected by an uncrewed aerial system–ArcticShark over the Southern Great Plains. Single particle microanalysis shows that, the free tropospheric (FT) samples are dominated (56-66%) by carbonaceous sulfate particles, while boundary layer (BL) samples are dominated (57-74%) by carbonaceous particles. Back trajectory simulations suggest that FT particles are likely influenced by long-range transport and have undergone aqueous-phase processing. Conversely, in-situ size distribution data shows evidence of particle growth in the upper BL and just below the FT. This observation may indicate vertical transport of particles from an elevated aerosol layer in the FT, possibly linked to a new particle formation event. This observation is further supported by high resolution molecular composition data, which reveals particle volatility increasing with increasing size, which aligns with the growth event. This study aids in fundamental understanding of the compositional and molecular specificity of vertically resolved organic aerosols to provide insights into particle size evolution for future atmospheric models.

ArcticShark↗

A More Transparent Infrared Window

The infrared window region (780–1,250 cm –1 , 12.8 to 8.0 μm) is of great importance to Earth's climate due to its high transparency and thermal energy. We present here a new investigation of the transparency of this spectral region based on observations by interferometers of downwelling surface radiance at two DOE Atmospheric Radiation Measurement program sites. We focus on the dominant source of absorption in this region, the water vapor continuum, and derive updated values of spectral absorption coefficients for both the self and foreign continua. Our results show that the self continuum is too strong in the previous version of Mlawer-Tobin_Clough-Kneizys-Davies (MT_CKD) water vapor continuum model, a result that is consistent with other recent analyses, while the foreign continuum is too weak in MT_CKD. In general, the weaker self continuum derived in this study results in an overall increase in atmospheric transparency in the window, although in atmospheres with low amounts of water vapor the transparency may slightly decrease due to the increase in foreign continuum absorption. These continuum changes lead to a significant decrease in downwelling longwave flux at the surface for moist atmospheres and a modest increase in outgoing longwave radiation. The increased fraction of surface-leaving radiation that escapes to space leads to a notable increase (~5–10%) in climate feedback, implying that climate simulations that use the new infrared window continuum will show somewhat less warming than before. This study also points out the possibly important role that aerosol absorption may play in the longwave radiative budget.

58 GEOSCIENCES↗

Photo-induced structural dynamics of o -nitrophenol by ultrafast electron diffraction

The photo-induced dynamics of o-nitrophenol, particularly its photolysis, has garnered significant scientific interest as a potential source of nitrous acid in the atmosphere. Although the photolysis products and preceding photo-induced electronic structure dynamics have been investigated extensively, the nuclear dynamics accompanying the non-radiative relaxation of o-nitrophenol on the ultrafast timescale, which include an intramolecular proton transfer step, have not been experimentally resolved. Herein, we present a direct observation of the ultrafast nuclear motions mediating photo-relaxation using ultrafast electron diffraction. Here, this work spatiotemporally resolves the loss of planarity which enables access to a conical intersection between the first excited state and the ground state after the proton transfer step, on the femtosecond timescale and with sub-Angstrom resolution. Our observations, supported by ab initio multiple spawning simulations, provide new insights into the proton transfer mediated relaxation mechanism in o-nitrophenol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comparative study with two reduced test facilities for air-cooled RCCS scaling law

Korea Atomic Energy Research Institute (KAERI) and Argonne National Laboratory (ANL) conducted steady-state tests at constant heat flux conditions to investigate the scale effect on heat removal behavior of the air-cooled Reactor Cavity Cooling System (RCCS). Two differently scaled-down test facilities were used in this study, KAERI’s ¼ scale with 4.0-m cavity height, and ANL’s ½ scale with 6.8-m cavity height. Two scaling laws were proposed to simulate the radiation across the cavity and the buoyancy-driven duct flow in the riser, respectively. The test matrix for KAERI and ANL test facilities focused on the scaling effect of the heated riser length with constant heat flux at the PMR200 RCCS design. The test results showed that mixed convection in the riser duct is an important factor for accurately extrapolating the thermo-fluid behavior in the prototype from the test results in the scale-down facilities. The system analysis code, GAMMA+, with improved heat transfer models predicted fairly well the air-cooled RCCS test data from two facilities. GAMMA+ analysis results showed that the predicted radiation fractions on the heated plate in the scale-down test conditions were larger than those on the reactor vessel in the prototype. The scaling law for air-cooled RCCS was improved by considering the mixed convection in the riser duct and the same radiation fraction on the heated plate. The mixed convection effect was calculated by the height ratio. The same radiation fraction provided a conservative extrapolation from the scale-down test results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Extracting symplectic maps for space-charge dominated beams

Symplecticity of transfer maps is important for reliable evaluation of space-charge dominated beams in accelerators. Unfortunately, most simulation codes that include collective effects, such as space charge, do not use canonical phase-space variables and therefore are not symplectic in the presence of electromagnetic fields. In this paper, we present a numerical method to extract local linear symplectic transfer maps using particle tracking simulation code for space-charge dominated beams. We demonstrate this method for the photoinjector (113 MHz SRF gun) section of the Coherent electron Cooling (CeC) Proof of Principle (POP) experiment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Super-resolution Imaging of Plasmonic Near-Fields: Overcoming Emitter Mislocalizations

Plasmonic nano-objects have shown great potential in enhancing applications like biological/chemical sensing, light harvesting and energy transfer, and optical/quantum computing. Therefore, an extensive effort has been vested in optimizing plasmonic systems and exploiting their field enhancement properties. Super-resolution imaging with quantum dots (QDs) is a promising method to probe plasmonic near-fields but is hindered by the distortion of the QD radiation pattern. Here, we investigate the interaction between QDs and “L-shaped” gold nanoantennas and demonstrate both theoretically and experimentally that this strong interaction can induce polarization-dependent modifications to the apparent QD emission intensity, polarization, and localization. Based on FDTD simulations and polarization-modulated single-molecule microscopy, we show that the displacement of the emitter’s localization is due to the position-dependent interference between the emitter and the induced dipole, and can be up to 100 nm. Our results help pave a pathway for higher precision plasmonic near-field mapping and its underlying applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular dynamics study of grain boundaries as defect sinks under irradiation in LiAlO 2 and LiAl 5 O 8

Lithium aluminate ceramics, LiAlO 2 and LiAl 5 O 8 , show promise in nuclear environments due to their excellent radiation tolerance. Molecular dynamics simulations investigate grain boundaries (GB) and their role in defect evolution. Results reveal that GBs act as efficient defect sinks, with Li and Al atoms exhibiting distinct behaviors during displacement cascades. Tritium migration in LiAlO 2 is also studied, showing rapid diffusion and stable configurations with oxygen, corroborated by ab initio simulations from the literature. The calculated tritium diffusion coefficient of 1.33 × 10 - ¹⁴ m²/s aligns with the literature, validating the model. LiAl 5 O 8 demonstrates superior defect healing compared to LiAlO 2 , attributed to enhanced atomic transfer between grains and GBs. These findings reveal key insights into defect dynamics, providing essential insights for their application in tritium-producing burnable absorber rods (TPBARs).

36 MATERIALS SCIENCE↗

Preparation of high orbital angular momentum Rydberg states by optical-millimeter-wave STIRAP

Rydberg states of molecules are intrinsically challenging to study due to the presence of fast non-radiative decay pathways, such as predissociation. However, selectively exciting Rydberg states with values of the orbital angular momentum (ℓ) ℓ ≳ 3 is a productive strategy to minimize this rapid decay and to populate molecular Rydberg states with lifetimes that approach those of atoms. Here, in this proof-of-principle demonstration, we transfer population to an nf Rydberg state of the calcium atom by stimulated Raman adiabatic passage, in which an optical and a millimeter-wave field couple the initial and final states via an intermediate nd Rydberg state. Numerical simulations reproduce the observed time and frequency dependences of the population transfer and suggest the utility of this scheme to populate high-ℓ Rydberg states of molecules.

74 ATOMIC AND MOLECULAR PHYSICS↗

Cluster Dynamics Simulations of Intra-Granular Fission Gas Bubble Size and Pressure Evolution in UO 2

Fission gases such as xenon (Xe) play a critical role in determining the behavior and response of nuclear fuel. Given that Xe has little solubility in UO 2 , it accumulates and forms bubbles, which significantly impact fuel performance. Intra- and inter-granular bubble nucleation and growth can lead to fuel swelling, and once bubbles interconnect at grain boundaries, fission gas can be released into the plenum. At low temperatures, limited uranium vacancy mobility can restrict swelling, therefore causing the bubbles to become highly pressurized. Consequently, this can induce micro-cracking, promote fission gas release (increasing the likelihood of cladding failure), and even lead to fuel pulverization under accident conditions such as a loss of coolant accident. As bubble evolution is strongly influenced by local temperature and fission rate, markedly different behavior occurs across the radial profile of the fuel pellet. Capturing the mechanisms that underpin bubble evolution is therefore important to predict these behaviors in the fuel. Previous models describing important mechanisms informed by lower length scale simulations have been developed under the NEAMS program. These can describe the evolution of a single bubble type (i.e., single value for radius and pressure) at each position in the pellet, for instance using the Centipede cluster dynamic code. However, in reality, a full distribution in bubble sizes and pressures exists within the microstructure at a given position in the pellet. To address this the cluster dynamics code Xolotl, which can predict Xe and vacancy phase space (i.e., bubble distributions) for intra-granular bubbles, has been used before. Prior work benchmarked the Xolotl code against the Centipede cluster dynamics code to ensure compatibility and to verify that mobile defect properties are adequately transferred between the two codes, along with some physics improvements. In this work, we go further by introducing a physics-based set of improvements that will allow us to accurately predict bubble size distributions and internal bubble pressures under representative UO 2 irradiation conditions. The improvements include (i) coupling bubble-defect reaction energies to a virial equation of state (EOS), (ii) including a bubble surface tension contribution, (iii) incorporating radiation-induced re-solution of Xe and vacancies, (iv) enabling pressure-driven dislocation loop punching through an effective emission of interstitial clusters informed by interstitial loop energetics, (v) accounting for radiation induced athermal diffusion of Xe, and (vi) implementing a Booth-type grain boundary sink representation for all mobile defects and defect clusters. After these modifications, we observe good agreement of Xolotl fission gas bubble size and concentration predictions with legacy experimental measurements. Additionally, it allows the distribution of Xe bubble pressures and radius to also be predicted and compared to data produced through the Advanced Fuels Campaign (AFC) program. Here, we have done this by running simulations under conditions similar to the AFC post-irradiation examination (PIE) samples irradiated at North Anna 2 light water reactor (LWR). Our results shows excellent agreement with these experimental measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Efficient emulation of relativistic heavy ion collisions with transfer learning

Measurements from the Large Hadron Collider (LHC) and the Relativistic Heavy Ion Collider (RHIC) can be used to study the properties of quark-gluon plasma. Systematic constraints on these properties must combine measurements from different collision systems and methodically account for experimental and theoretical uncertainties. Such studies require a vast number of costly numerical simulations. While computationally inexpensive surrogate models (“emulators”) can be used to efficiently approximate the predictions of heavy ion simulations across a broad range of model parameters, training a reliable emulator remains a computationally expensive task. We use transfer learning to map the parameter dependencies of one model emulator onto another, leveraging similarities between different simulations of heavy ion collisions. By limiting the need for large numbers of simulations to only one of the emulators, this technique reduces the numerical cost of comprehensive uncertainty quantification when studying multiple collision systems and exploring different models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Investigation of the Failure of Marcus Theory for Hydrated Electron Reactions

Reactions of the hydrated electron with a wide variety of substrates have been found to exhibit unusually similar activation energies in a manner incompatible with Marcus electron transfer theory. Given the fundamental linear response assumption of Marcus theory, one possible explanation for this apparent failure is that the underlying free energy surfaces governing the reactions are not harmonic, i.e., hydrated electron structural fluctuations exhibit non-Gaussian behavior. In this work, we test this hypothesis by using simulations to calculate the hydrated electron vertical detachment energy distribution. Furthermore, we consider both cavity and non-cavity models for the hydrated electron, between which the actual hydrated electron behavior is expected to lie. Our results identify a possible origin for non-Gaussian behavior of the hydrated electron, but show that it is not of sufficient magnitude to explain the failure of Marcus theory to describe its reactions. Thus, other explanations must be sought.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗