Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “scale coupling”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Air-sea interactions and precipitation over the tropical oceans

In this lecture, the author principally discusses air-sea exchanges that are relevant to climate and global problems. The processes of interest are those acting over time scales of months to decades, which in some instances are influenced by smaller-time-scale processes, down to the diurnal time scale. The repsective influence of these processes varies with regions, seasons and scales over which they occur and, because these processes are mostly nonlinear, scale interactions can be quite complex. Owing to the breadth of the topic addressed, the discussion is mostly focused on the tropical regions where air-sea interactions and precipitation processes eventually affect the entire globe. This allows a look in more detail at some air-sea processes, such as those associated with the El Nino southern oscillation (ENSO). This oscillation, which affects the climate of the entire globe, acts over periods of a year or longer and is caused, primarily, by sea surface temperature (SST) variations in the tropical Pacific. As a result, SST variability is often used as an indicator of coupled ocean-atmosphere low-frequency variability. Global or basin scale processes can uniquely be observed from space-born instruments with the coverage required. Space based techniques have been developed during the last decade which can now be used to illustrate the scientific issues presented and the presentation concludes with an overview of some Earth Observing System (EOS) capabilities for addressing air-sea interactions and hydrology issues.

Gautier, C.↗
DOE OSTI2021

Impacts of Anisotropic Porosity on Heat Transfer and Off-Gassing during Biomass Pyrolysis

The pore structure of biogenic materials imbues the ability to deliver water and nutrients through a plant from root to leaf. This anisotropic pore granularity can also play a significant role in processes such as biomass pyrolysis that are used to convert these materials into useful products like heat, fuel, and chemicals. Evolutions in modeling of biomass pyrolysis as well as imaging of pore structures allow for further insights into the concerted physics of phase change-induced off-gassing, heat transfer, and chemical reactions. In this work, we report a biomass single particle model which incorporates these physics to explore the impact of implementing anisotropic permeability and diffusivity on the conversion time and yields predicted for pyrolysis of oak and pine particles. Simulation results showed that anisotropic permeability impacts predicted conversion time more than 2 times when the Biot number is above 0.1 and pyrolysis numbers (Py 1 , Py 2 ) are less than 20. Pore structure significantly impacts predicted pyrolytic conversion time (>8 times) when the Biot number is above 1 and the pyrolysis number is below 1, i.e., the “conduction controlled” regime. Therefore, these nondimensional numbers reflect that when internal heat conduction limits pyrolysis performance, internal pyrolysis off-gassing further retards effective heat transfer rates as a closely coupled phenomenon. Overall, this study highlights physically meaningful opportunities to improve particle-scale pyrolysis modeling and experimental validation relevant to a variety of feedstock identities and preparations, guiding the future design of pyrolyzers for efficient biomass conversion.

09 BIOMASS FUELS↗
DOE OSTI2023

Automated Symbolic Upscaling: 1. Model Generation for Extended Applicability Regimes

Abstract In porous media theory, upscaling techniques are fundamental to deriving rigorous Darcy‐scale models for flow and reactive transport in subsurface systems. Due to limitations in classical techniques, a number of ad hoc approaches have been proposed to address physical regimes in which reactive time scales are similar to, or faster than, diffusive time scales. In Part 1 of this two part series, we present a strategy for expanding the applicability of classical homogenization theory by generalizing the assumed closure form. We detail the implementation of this strategy on two reactive mass transport problems with moderately reactive physics. The strategy produces nontrivial homogenized models with emergent terms and effective parameters that couple reactive, diffusive, and advective transport. The differences in equation forms between the macroscopic and pore‐scale descriptions advise caution to further studies where the forms of macroscopic equations are assumed, as opposed to rigorously derived. Numerical validation is provided for each problem to show that the error estimates of homogenization theory are satisfied, and to justify the implemented strategy. In Part 2, the presented strategy is automated using symbolic computing to expedite its implementation.

Pietrzyk, Kyle↗
DOE OSTI2026

A Scalable Reduced‐Order Model for the Steady Navier–Stokes Equations

Scaling up new scientific technologies from laboratory to industry often involves demonstrating performance on a larger scale. Computer simulations can accelerate design and predictions in the deployment process, though traditional numerical methods are computationally intractable even for intermediate pilot plant scales. Recently, the component reduced order modeling method has been developed to tackle this challenge by combining projection reduced order modeling and discontinuous Galerkin domain decomposition. However, while many scientific or engineering applications involve nonlinear physics, this method has only been demonstrated for various linear systems. In this work, the component reduced order modeling method is extended to steady Navier–Stokes flow, with application to general nonlinear physics in view. The large‐scale, global domain is decomposed into a combination of small‐scale unit component. Linear subspaces for flow velocity and pressure are identified via proper orthogonal decomposition over sample snapshots collected from each small‐scale unit component. Velocity bases are augmented with a pressure supremizer to satisfy the inf–sup condition for stable pressure prediction. Two different nonlinear reduced order modeling methods are employed and compared for efficient evaluation of nonlinear advection: A third‐order tensor projection operator and the empirical quadrature procedure. The proposed method is demonstrated on the flow over arrays of five different unit objects, achieving a 23‐fold speedup with less than 4% relative error in domains up to 256 times larger than the unit components. Furthermore, a numerical experiment with the pressure supremizer strongly indicates the need for a supremizer for stable pressure prediction. A comparison between the tensorial approach and the empirical quadrature procedure revealed a slight advantage of the empirical quadrature procedure. The framework is compared with an alternating Schwarz‐based reduced‐order approach, demonstrating improved efficiency and robustness for the DG‐based global solver while retaining flexibility for sub‐scale iterative solvers. The method is further extended to a coupled advection–diffusion and Navier–Stokes system, illustrating its applicability to multi‐physics problems and its potential for more general, inter‐coupled nonlinear systems.

42 ENGINEERING↗
DOE OSTI2020

Spontaneous breaking of U(1) symmetry in coupled complex SYK models

As shown in [1], two copies of the large N Majorana SYK model can produce spontaneous breaking of a Z 2 symmetry when they are coupled by appropriate quartic terms. In this paper we similarly study two copies of the complex SYK model coupled by a quartic term preserving the U(1) x U(1) symmetry. We also present a tensor counterpart of this coupled model. When the coefficient α of the quartic term lies in a certain range, the coupled large N theory is nearly conformal. We calculate the scaling dimensions of fermion bilinear operators as functions of α. We show that the operator c$^†_{1i}$c 2i , which is charged under the axial U(1) symmetry, acquires a complex dimension outside of the line of fixed points. We derive the large N Dyson-Schwinger equations and show that, outside the fixed line, this U(1) symmetry is spontaneously broken at low temperatures because this operator acquires an expectation value. We support these findings by exact diagonalizations extrapolated to large N.

1/N expansion↗
DOE OSTI2023

Design, fabrication, and testing of a water-cooled tungsten-foam high-heat flux plasma-facing module

Here, we present a comprehensive study detailing the design, fabrication and testing a new modular “finger” component, featuring a hexagonal tungsten foam armor tile bonded to the top surface of the water-cooled finger unit. The mechanical, thermal, and particle transport properties of the tungsten foam armor were determined for a wide range of porosity and cell size, and then combined in an equivalent solid material with “effective” properties suitable for large-scale simulations. A series of 3-D multiphysics simulations for coupled thermal and mechanical behavior of the “finger” module are used to optimize the design before fabrication. The simulation results show that, for the expected experimental conditions, a safe range of heat flux and foam porosity is possible for both water-cooling and helium-cooling. Based on these multiphysics simulations, an optimum finger cooling module has been fabricated at Ultramet, Inc. A series of experiments were designed to study the thermomechanical performance of the “finger” unit under steady-state and transient thermal loading at surface heat flux up to 16 MW/m 2 in the High Energy Flux Test facility (HEFTY) at UCLA. Testing conditions included low-power, steady-state plasma exposure at 6 MW/m 2 , followed by severe high heat flux cycles at 16 MW/m 2 . Results of these tests are explained through comparisons with multiphysics simulations. To eliminate the effects of thermomechanical damage carried out on a module on its thermo-mechanical response on a subsequent test, three different yet identically designed and fabricated finger modules were tested. The results show that the current design is capable of cooling the walls of fusion devices at surface heat flux up to 6 MW/m 2 without failure, and that the failure heat flux limit is below 16 MW/m 2 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Catalytic consequences of hydrogen addition events and solvent-adsorbate interactions during guaiacol-H 2 reactions at the H 2 O-Ru(0 0 0 1) interface

Catalytic reactions of biomass-derived phenolics and H 2 occur on transition metal surfaces via competitive C–O cleavage and ring saturation pathways, with both requiring multiple hydrogen addition events before forming their respective rate limiting transition states. These events are markedly affected by solvent chemical identity, with polar protic solvents ionizing hydrogen adatoms (H*) to interfacial protons (H + ) and opening up new catalytic routes. Here, we establish the reaction coordinate space for guaiacol-H 2 reactions on Ru(0 0 0 1) using density functional theory and describe the atomic-scale effect of a polar protic solvent, H 2 O. Coupled H + and H* attack leads to quasi-equilibrated enol and keto intermediates as the precursors for C–O cleavage and ring saturation, respectively. For C–O cleavage, H 2 O solvent enables a lower energy pathway via concomitant transfer of the hydroxyl H + to the methoxy oxygen during C–OCH 3 cleavage, forming a charge separated [Ru(s)–(C 6 H 5 O – )…(H + )…OCH 3 ] transition state and reducing the barrier by up to 0.8 eV as compared to unassisted C–OCH 3 cleavage. For ring saturation, H* attack onto an unsaturated meta carbon is rate limiting with no direct solvent participation, suggesting that protic polar solvents selectively promote the C–O cleavage pathway. Taken together, we show that activating guaiacol for either C–O bond cleavage or ring saturation product formation depends on the reactive hydrogen identity (H* or H + ), enol/keto isomerization equilibrium, and accessibility of the proton assisted Car–OCH 3 cleavage transition state. Here, all such factors are tunable via changes to the solvent or metal identity.

09 BIOMASS FUELS↗
DOE OSTI2020

Vertical Strain-Driven Antiferromagnetic to Ferromagnetic Phase Transition in EuTiO 3 Nanocomposite Thin Films

Three-dimensional (3D) strain induced in self-assembled vertically aligned nanocomposite (VAN) epitaxial films provides an unrivaled method to induce very large strains in thin films. Here, by growing VAN films of EuTiO 3 (ETO)–Eu 2 O 3 (EO) with different EO fractions, the vertical strain was systematically increased in ETO, up to 3.15%, and the Eu–Ti–Eu bond angle along $\langle$111$\rangle$ decreased by up to 1°, leading to a weakening of the antiferromagnetic interactions and switching from antiferromagnetic to ferromagnetic behavior. Our work has shown for the first time that Eu–Ti–Eu superexchange interactions play a key role in determining the magnetic ground state of ETO. Broadly, our work serves as an exemplar to show that multifunctionalities in strong spin–lattice coupling perovskite oxides can be uniquely tuned at the atomic scale using simple VAN structures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Visualizing Light-Induced Microstrain and Phase Transition in Lead-Free Perovskites Using Time-Resolved X-Ray Diffraction

Metal halide perovskites have emerged as promising materials for optoelectronic applications in the last decade. A large amount of effort has been made to investigate the interplay between the crystalline lattice and photoexcited charge carriers as it is vital to their optoelectronic performance. Among them, ultrafast laser spectroscopy has been intensively utilized to explore the charge carrier dynamics of perovskites, from which the local structural information can only be extracted indirectly. Here, we have applied a time-resolved X-ray diffraction technique to investigate the structural dynamics of prototypical two-dimensional lead-free halide perovskite Cs 3 Bi 2 Br 9 nanoparticles across temporal scales from 80 ps to microseconds. We observed a quick recoverable (a few ns) photoinduced microstrain up to 0.15% and a long existing lattice expansion (~a few hundred nanoseconds) at mild laser fluence. Once the laser flux exceeds 1.4 mJ/cm 2 , the microstrain saturates and the crystalline phase partially transfers into a disordered phase. This photoinduced transient structural change can recover within the nanosecond time scale. Furthermore, these results indicate that photoexcitation of charge carriers couples with lattice distortion, which fundamentally affects the dielectric environment and charge carrier transport.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Multiple Metrics Informed Projections of Future Precipitation in China

Predicting how regional precipitation will respond to future warming is among the most challenging undertaking in climate change projection. Despite sustained efforts to improve modeling and understanding of precipitation, the overall uncertainty in projecting regional precipitation has not been reduced substantially. In this work, the potential for more robust precipitation projections is demonstrated through the use of discriminating metrics to subsample the multimodel ensemble. Using a two-dimensional metric of precipitation and its relationship with large-scale circulation indices in East Asia, 31 models in the Coupled Model Intercomparison Project Phase 5 (CMIP5) are classified into three groups. Models in the top performing group projected statistically significant increasing trends in precipitation and the regional precipitation patterns are more similar to each other than to the patterns in the bottom performing group. In contrast, models in the bottom performing group projected diverse responses, with overall small drying or no significant trends in precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Deep learning ferroelectric polarization distributions from STEM data via with and without atom finding

Over the last decade, scanning transmission electron microscopy (STEM) has emerged as a powerful tool for probing atomic structures of complex materials with picometer precision, opening the pathway toward exploring ferroelectric, ferroelastic, and chemical phenomena on the atomic scale. Analyses to date extracting a polarization signal from lattice coupled distortions in STEM imaging rely on discovery of atomic positions from intensity maxima/minima and subsequent calculation of polarization and other order parameter fields from the atomic displacements. Here, we explore the feasibility of polarization mapping directly from the analysis of STEM images using deep convolutional neural networks (DCNNs). In this approach, the DCNN is trained on the labeled part of the image (i.e., for human labelling), and the trained network is subsequently applied to other images. We explore the effects of the choice of the descriptors (centered on atomic columns and grid-based), the effects of observational bias, and whether the network trained on one composition can be applied to a different one. This analysis demonstrates the tremendous potential of the DCNN for the analysis of high-resolution STEM imaging and spectral data and highlights the associated limitations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Operator learning for predicting multiscale bubble growth dynamics

We report simulating and predicting multiscale problems that couple multiple physics and dynamics across many orders of spatiotemporal scales is a great challenge that has not been investigated systematically by deep neural networks (DNNs). Herein, we develop a framework based on operator regression, the so-called deep operator network (DeepONet), with the long-term objective to simplify multiscale modeling by avoiding the fragile and time-consuming “hand-shaking” interface algorithms for stitching together heterogeneous descriptions of multiscale phenomena. To this end, as a first step, we investigate if a DeepONet can learn the dynamics of different scale regimes, one at the deterministic macroscale and the other at the stochastic microscale regime with inherent thermal fluctuations. Specifically, we test the effectiveness and accuracy of the DeepONet in predicting multirate bubble growth dynamics, which is described by a Rayleigh–Plesset (R–P) equation at the macroscale and modeled as a stochastic nucleation and cavitation process at the microscale by dissipative particle dynamics (DPD). First, we generate data using the R–P equation for multirate bubble growth dynamics caused by randomly time-varying liquid pressures drawn from Gaussian random fields (GRFs). Our results show that properly trained DeepONets can accurately predict the macroscale bubble growth dynamics and can outperform long short-term memory networks. We also demonstrate that the DeepONet can extrapolate accurately outside the input distribution using only very few new measurements. Subsequently, we train the DeepONet with DPD data corresponding to stochastic bubble growth dynamics. Although the DPD data are noisy and we only collect sparse data points on the trajectories, the trained DeepONet model is able to predict accurately the mean bubble dynamics for time-varying GRF pressures. Taken together, our findings demonstrate that DeepONets can be employed to unify the macroscale and microscale models of the multirate bubble growth problem, hence providing new insight into the role of operator regression via DNNs in tackling realistic multiscale problems and in simplifying modeling with heterogeneous descriptions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Interpretation of electromagnetic modes in the sub-TAE frequency range in JET plasmas with elevated monotonic q-profiles

Recent JET deuterium experiments with an advanced tokamak scenario using an internal transport barrier (ITB) have been performed to clearly observe destabilized toroidicity-induced Alfvén eigenmodes (TAEs) by fast ions; interestingly, these also exhibit unstable electromagnetic (EM) perturbations in the sub-TAE frequency range. We identify such EM perturbations to be beta-induced ion temperature gradient (BTG) eigenmodes and not beta-induced Alfvén eigenmodes nor beta-induced Alfvén acoustic eigenmodes, which are often unstable in such high-beta plasmas with high-power neutral beam injection (NBI). The BTG modes are the most unstable modes due to the high thermal ion temperature gradient related to the ITB, high thermal ion temperature compared to thermal electron temperature (high T i /T⁠ e ), and a high ion beta regime. BTG mode experimental characteristics match the analytical theory, i.e., location in the vicinity of a rational magnetic surface with a low magnetic shear, mode frequency scaling with the ion diamagnetic frequency (⁠ω$^*_i$⁠), and a coupling among Alfvén and drift waves. We also perform linear gyrokinetic simulations with validated plasma profiles and equilibrium, and find a mode kinetically driven by thermal ions with similar characteristics as the experimental BTG modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Spheroidization of Nd–Fe–B particles

Spherical powders are required for many advanced manufacturing techniques due to their inherent requirement of flowability, either within feed tubes or during powder spreading. As advanced manufacturing of magnets continues to develop, new production methods for feedstocks are also sought. Plasma spheroidization is a high-yield method to produce spherical Nd–Fe–B powders from irregularly shaped particles, with advantages including high throughput and a well-controlled size distribution. Highly spherical Nd–Fe–B powders with large scale production (i.e., kg) have been demonstrated using an inductively coupled thermal plasma system; however, the magnetic properties of the output powder display significant degradation. The coercivity was decreased from the initial 8 kOe (636 kA/m) of the as-received to 0.7 kOe (55 kA/m) for spheroidized powders. Microstructural investigation reveals 6% Nd depletion caused by the extreme temperatures of the plasma, leading to the formation of low-coercivity α-Fe and a subsequent decrease in energy product. In conclusion, post-spheroidization heat treatments with Nd can partially mitigate the coercivity degradation, increasing to 1.7 kOe (135 kA/m), potentially offering a pathway toward spherical powders for a range of applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

An Initial Microstructurally Informed Model of High Burnup Structure Formation in UO 2 Fuel

The microstructure of a UO 2 fuel pellet changes as burnup increases, impacting fuel performance. Predicting and characterizing high burnup structure (HBS) and dark zone formation is a key part of supporting burnup limit extensions for light water reactors. This paper describes a model developed through fitting radially resolved pellet data obtained from recently published microstructural characterization data. The model predicts grain size and grain character, in addition to pore density and size, with fitting dependencies on power history variables. Separately fitting power history variables to microstructural parameters allows for insight into the underlying physical phenomena for future model development. Additionally, experimental data have been correlated to an HBS fraction to facilitate the development of a model capable of predicting a total fuel restructured fraction at the engineering scale. In conclusion, this two-step approach provides a coupling from reactor power history to microstructural data to fractional HBS and creates a basis to model HBS-dependent parameters in a fuel performance code.

High burnup structure↗
DOE OSTI2025

On the role of ion temperature gradient turbulence in driving ion thermal transport in neutral beam injection-heated L-mode plasmas in a superconducting tokamak

In this paper, we report a detailed experimental study of the role of ion-temperature-gradient driven turbulence in driving ion thermal transport in neutral beam injection (NBI)-heated L-mode plasmas in Experimental Advanced Superconducting Tokamak (EAST) (Wan et al 2000 Nucl. Fusion 40 1057) coupled with linear and nonlinear gyrokinetic simulations. Significant ion-scale turbulence k ⟂ < cm −1 (k ⟂ < 1.5, where k ⟂ where is the perpendicular wavenumber and ρ s is the ion gyroradius calculated using local electron temperature T e ), measured with a microwave reflectometer, is observed in the plasma core at for r/a ≈ 0.265 and 0.5 (where r is half the diameter of the closed flux surface at a given radial position, and a is half the diameter of the last closed flux surface). Local linear stability analysis with experimental equilibrium quantities at these two radial locations using the GS2 gyrokinetic code shows that the most unstable ion-scale micro-instability is the ion temperature gradient (ITG) mode. Since the computed maximum ion-scale linear growth rates are larger than the local Waltz-Miller E × B shearing rate (Waltz and Miller 1999 Phys. Plasmas 6 4265), the E × B shear is unable to suppress ion-scale turbulence, consistent with the experimental observation of ion-scale turbulence. Experimental ion and electron thermal transports, calculated with power balance analysis, are both anomalous with ion thermal transport being the dominant channel. Nonlinear gyrokinetic simulations using the electrostatic global particle-in-cell gyrokinetic tokamak simulation (GTS) code (Weixing Wang et al 2010 Phys. Plasmas 17 072511) show good quantitative agreement between predicted ion thermal transport and that from the power balance analysis. These results demonstrate for the first time the important role of ITG turbulence in driving thermal transport in NBI-heated L-mode plasmas in EAST, showing the effectiveness of GTS code in simulating electrostatic turbulence in EAST and the need of suppressing ITG turbulence to improve energy confinement in EAST.

Sun, P. J. [Chinese Academy of Sciences (CAS), Hef↗
DOE OSTI2025

Unprecedented Beaufort Sea ice loss in late summer 2021 and its relationship to an extended period of unusually stormy weather

Previous case studies have linked cyclone-induced atmospheric forcing and/or upper-ocean processes to notable Arctic sea ice loss events in the summers of 2012 and 2016. This study examines a more recent and noteworthy case in late summer 2021 in which substantial sea ice loss followed a period of surface meteorological extremes in the Beaufort Sea region of the Arctic. We focus on the period from mid-August to mid-September 2021 that coincided with the Office of Naval Research THINICE Pilot Field Campaign and investigate stormy and windy conditions with respect to air-sea processes impacting sea ice conditions. We find that during the stormy first half of the campaign, cyclone-induced energy fluxes into the marginal ice zone and surrounding waters preconditioned the ice pack for more rapid melt later in the campaign. The second half of the campaign, in contrast, was marked by non-cyclone wind events that enhanced turbulent (namely sensible) heat fluxes into the ice and upper ocean that increased melt. Moreover, this latter period had enhanced advection of the Beaufort Sea ice pack into above-freezing waters, increasing bottom melt to >1 cm d −1 over the remainder of the campaign. While findings are shown to vary by surface type and at relatively small (i.e. ice-floe) scales, insights are offered on the roles of late summer coupled processes on rapid ice loss events in today’s Arctic environment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Mitigation of birefringence in cavity-based quantum networks using frequency-encoded photons

Atom-cavity systems offer unique advantages for building large-scale distributed quantum computers by providing strong atom-photon coupling while allowing for high-fidelity local operations of atomic qubits. However, in prevalent schemes where the photonic state is encoded in polarization, cavity birefringence introduces an energy splitting of the cavity eigenmodes and alters the polarization states, thus limiting the fidelity of remote entanglement generation. To address this challenge, we propose a scheme that encodes the photonic qubit in the frequency degree-of-freedom. The scheme relies on resonant coupling of multiple transverse cavity modes to different atomic transitions that are well-separated in frequency. We numerically investigate the temporal properties of the photonic wavepacket, two-photon interference visibility, and atom-atom entanglement fidelity under various cavity polarization-mode splittings and find that our scheme is less affected by cavity birefringence. Finally, we propose practical implementations in two trapped ion systems, using the fine structure splitting in the metastable D state of 40 Ca + , and the hyperfine splitting in the ground state of 225 Ra + . Furthermore, our study presents an alternative approach for cavity-based quantum networks that is less sensitive to birefringent effects, and is applicable to a variety of atomic and solid-state emitter-cavity interfaces.

Cavity quantum electrodynamics↗