Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Cross Flow”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Apparatus for a laser welding system

A laser welding system for welding a component and reducing defects in the weld by ensuring uniform, laminar gas flow over a process area of the system. The laser welding system comprises a laser for welding the component, a platform for supporting the component, an enclosure surrounding the platform, a first actuatable barrier, a second actuatable barrier, an actuator, and a controller. The enclosure includes a plurality of walls, one of the walls having an inlet and another wall having an outlet. The inlet and outlet each having an opening having a cross-sectional area for letting gas flow through. The first and second barriers are configured to modify the cross-sectional areas of the openings when actuated. The actuator is configured to actuate the barriers, and the controller is configured to direct the actuator to actuate the barriers so that the cross-sectional area of the first opening is larger than the cross-sectional area of the second opening so that a pressure at the inlet is greater than a pressure at the outlet.

Sartin, Bryan↗
DOE OSTI2021

Turbulence-driven flow dynamics in general axisymmetric toroidal geometry

The present work gives the equations governing the generation of toroidally axisymmetric flows by turbulent Reynolds and Maxwell stresses in finite aspect ratio, general cross section tokamak plasmas. Inclusion of the divergence-free flow constraint in lowest order changes the time scale for evolution of the poloidal flows driven by turbulence by substantial factors. In the pedestal region for present day machines, comparing to earlier cylindrical models, the time scale evaluated using a large aspect ratio circular cross section model can be two orders of magnitude longer while the present, general geometry result can be about one order of magnitude longer. Inclusion of gyroviscosity in the calculation shows that the only lowest order radial velocity fluctuations that enter the problem are those due to fluctuating E×B flows. Toroidal and poloidal flow effects on the toroidally axisymmetric flows are inextricably coupled due to the neoclassical poloidal viscosity. Accordingly, the physics is inherently three dimensional and measurements of all three velocity components are required to obtain the information needed to quantitatively test the theory. The parallel and angular momentum equations for the lowest order, toroidally axisymmetric flows look like radial transport equations when the turbulence is included. The turbulence terms provide the radial transport fluxes. In standard neoclassical theory, the parallel flow equation is local on each flux surface; there is no radial derivative term. However, adding turbulence gives a way, in principle, for radial transport to lead to poloidal flows that deviate from the neoclassical prediction. As a result, the inclusion of the Maxwell stress provides a mechanism for MHD fluctuations to alter the toroidally axisymmetric flows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Evolution of in-medium baryon-baryon scattering cross sections and stiffness of dense nuclear matter from Bayesian analyses of FOPI proton-flow excitation functions

Within a Bayesian statistical framework using a Gaussian process (GP) emulator for an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model simulator of heavy-ion reactions with momentum-independent Skyrme interactions, we infer from the proton directed and elliptical flow in mid-central Au + Au reactions, measured at beam energies from 150 to 1200 MeV/nucleon by the FOPI Collaboration, the posterior probability distribution functions (PDFs) of the in-medium baryon-baryon scattering cross section (BBSCS) modification factor 𝑋 (with respect to their free-space values) and the stiffness parameter 𝐾 of dense nuclear matter. Here, we find that the most probable value of 𝑋 evolves from around 0.7 to 1.0 as the beam energy 𝐸 beam /𝐴 increases. On the other hand, the posterior PDF⁡(𝐾) may have dual peaks having roughly the same height or extended shoulders at high 𝐾 values. More quantitatively, the posterior PDF⁡(𝐾) changes from having a major peak around 220 MeV, characterizing a soft EOS in the reaction at 𝐸 beam /𝐴 = 150 MeV, to one that peaks around 320 MeV, indicating a stiff EOS in the reactions at 𝐸 beam /𝐴 higher than about 600 MeV. The transition from soft to stiff happens in mid-central Au + Au reactions at beam energies around 250 MeV/nucleon in which 𝐾 = 220 MeV and 𝐾 = 320 MeV are approximately equally probable. Altogether, the FOPI proton-flow excitation function data indicate a gradual hardening of hot and dense nuclear matter as its density and temperature increase in reactions with higher beam energies.

nuclear matter↗
DOE OSTI2021

Geometric Optimization of an Electrochemical Purification Cell to Prevent Corrosion in CSP Plants During Operation

When exposed to moisture or oxygen, molten chloride salts produce corrosive impurities which degrade containment alloys. As a result, this can significantly decrease the lifetime and increase costs of molten-salt-based systems. To overcome this barrier, we designed and modeled an electrochemical purification cell to remove the corrosive impurity MgOH+. Various reactor architectures, including constant stirred tank reactors (CSTRs) and plug flow reactors (PFRs) were investigated. Steady-state thermoelectric properties were evaluated, allowing assessment of the effects of structure and design parameters such as flow rate, cell length, and cross-sectional area of molten salt. The results suggest that our design could most effectively increase reliability and decrease costs of molten-chloride-salt-based systems by protecting them during continuous operation using an annular plug flow reactor.

chloride salts↗
DOE OSTI2022

Proactive Posturing of Large Power Grid for Mitigating Hurricane Impacts

In the past decade, natural disasters such as hurricanes have challenged the operation and control of U.S. power grid more. It is crucial to develop proactive strategies to assist grid operators for better emergency response and minimized electricity service interruptions. In this paper, we propose a proactive posturing methodology of power system elements, and formulated a Security-Constrained Optimal Power Flow (SCOPF) that informed by cross-domain hurricane modeling and its potential impacts on the grid elements. Simulation results based on real-world power grid and historical hurricane event validated the applicability of the proposed optimization formulation, and show potential to enable grid operators and planners with interactive cross-domain data analytics in mitigating hurricane impacts.

power system emergency response, security constrai↗
DOE OSTI2022

Feynman rules for scalar conformal blocks

We complete the proof of “Feynman rules” for constructing M-point conformal blocks with external and internal scalars in any topology for arbitrary M in any spacetime dimension by combining the rules for the blocks (based on their Witten diagram interpretation) with the rules for the construction of conformal cross ratios (based on the OPE and “flow diagrams”). The full set of Feynman rules leads to blocks as power series of the hypergeometric type in the conformal cross ratios. We then provide a proof by recursion of the Feynman rules which relies heavily on the first Barnes lemma and the decomposition of the topology of interest in comb structures. Finally, we provide a nine-point example to illustrate the rules.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Stepped pressure equilibrium with relaxed flow and applications in reversed-field pinch plasmas

The multi-region relaxed magnetohydrodynamics (MRxMHD) has been successful in the construction of equilibria in three-dimensional (3D) configurations. In MRxMHD, the plasma is sliced into sub-volumes separated by ideal interfaces, each undergoing relaxation, allowing the formation of islands and chaos. The resulting equilibrium has a stepped pressure profile across sub-volumes. The stepped pressure equilibrium code (SPEC) was developed to calculate MRxMHD equilibria numerically. In this work, we have extended the SPEC code to compute MRxMHD equilibria with field-aligned flow and rotation, following the theoretical development to incorporate cross-helicity and angular momentum constraints. The code has been verified for convergence and compared to a Grad–Shafranov solver in 2D. We apply our new tool to study the flow profile change before and after the sawtooth crash of a reversed-field pinch discharge, in which data of the parallel flow is available. We find the promising result that under the constraints of cross-helicity and angular momentum, the parallel flow profile in post-crash SPEC equilibrium is flat in the plasma core and the amplitude of the flow matches experimental observations. Lastly, we provide an example equilibrium with a 3D helical field structure as the favored lower energy state. This will be the first 3D numerical equilibrium in which the flow effects are self-consistently calculated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Reduced Order Models for Liquid Hydrogen Pooling and Vaporization Supported by Experiments

In the event of a leak of liquid hydrogen, a pool can form that vaporizes, disperses, and eventually dilutes to a non-flammable mixture. In this work, we describe fast-running models for the pooling and vaporization of liquid hydrogen in a steady cross-wind. Several pooling models from the literature are compared to solve for the flow and extent of the pool. The size of the pool can serve as the source for a separate dispersion model, which builds upon the existing one-dimensional Gaussian plume model in HyRAM+. Additional terms for the effects of a cross-wind on momentum and entrainment were added so that the model could handle the effects of a cross-wind on a low-speed flow. The models are compared to experimental data on pooling extent and downwind dispersion for steady flow rates of liquid hydrogen in a steady cross-wind. In the two compared experiments, liquid flow rates of 15 and 45 g/s were spilled onto concrete in cross-winds of approximately 1.8 m/s. The rate of growth of the pool and the downwind concentration boundaries are compared to the models, showing good agreement, although additional tuning is needed. These models can contribute to the advancement of codes and standards for liquid hydrogen systems.

dispersion↗
DOE OSTI2025

Belowground cross-trophic networks impact CH4 and CO2 emissions in degraded alpine peatlands

Belowground organisms forming complex cross-trophic ecological networks are essential for maintaining peatland carbon stability and energy flow. However, how peatland degradation affects the biodiversity and cross-trophic ecological networks of soil communities remains poorly understood. Here, we examined the degradation effects on soil prokaryotes (i.e., bacteria, archaea), fungi and nematodes in alpine peatlands on the eastern Tibetan Plateau, characterized by varying water table depths (indicating degradation levels). We found that peatland degradation, accompanied by significant shifts in soil moisture and pH (P < 0.05), reduced the taxonomic richness and phylogenetic diversity of prokaryotes, fungi, and nematodes, particularly in deeper soil layers (20–50 cm). Crucially, peatland degradation weakened potential cross-trophic interactions within bipartite networks of prokaryotes-nematodes and fungi-nematodes, resulting in less than 6.5 %–28.8 % of unchanged modules. Degradation-induced changes in soil moisture and pH were identified as primary drivers of biodiversity loss and network restructuring. Furthermore, such changes of belowground cross-trophic networks (particularly prokaryote-nematode) were significantly correlated with greenhouse gas emissions, such as decreased CO2 emissions, maintained CH4 emissions (leading to a higher CH4/CO2 ratio in deep layers), and reduced temperature sensitivity (Q10) of soil respiration. These findings underscore the critical need to protect soil biodiversity and cross-trophic networks in peatlands, particularly under the threat of climate change, to preserve peatland carbon stocks and maintain ecosystem stability. Our findings highlight that belowground cross-trophic networks are pivotal to decipher soil carbon dynamics of degraded peatlands and project the fate of peatland carbon stocks under future climate change scenarios.

Cross-trophic networks↗
DOE OSTI2023

Fourier coefficients of noninterdependent collective motions in heavy-ion collisions

Here, we present a scenario in heavy-ion collisions where different modes of collective motions are noninterdependent, driven by factorized actions in the created nuclear medium. Such physics mechanisms could each dominate at a distinct evolution stage, or coexist simultaneously. If the probability of particle emission is modulated by each nondependent collective motion with a single-harmonic Fourier expansion, the particle azimuthal distribution should be the product of all these expansions. Consequently, nonleading cross terms between collectivity modes appear, and their contributions to experimental observables could be significant. In particular, we argue that the chiral magnetic effect (CME) and elliptic flow can develop separately, with their convolution affecting the observable that is sensitive to the shear-induced CME. We will use the event-by-event anomalous-viscous fluid dynamics model to illustrate the effects of this scenario. Besides giving insights into searches for the CME, we also propose feasible experimental tests based on conventional flow harmonics, and demonstrate the emergence of nonleading cross terms with a multiphase transport model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Elucidating Norrish type I reactive pathways by ultrafast X-ray absorption spectroscopy

Norrish type I reactions selectively cleave carbon–carbon bonds directly adjacent to carbonyl groups. Despite their broad use in combination with aromatic carbonyls for additive manufacturing and dental UV curing applications, the nature of the photochemically active state and its population mechanism remain insufficiently understood. Detailed mechanistic insight requires mapping of the photoexcited population flow involving internal conversion and intersystem crossing. We present a time-domain study of gas phase acetophenone as a prototypical aromatic carbonyl combining soft X-ray time-resolved near-edge X-ray absorption fine structure (TR-NEXAFS) spectroscopy at the oxygen K-edge with ab initio multiple spawning (AIMS) simulations. Exploiting the specific sensitivity of TR-NEXAFS spectroscopy to states with nπ* character, we observe population transfer from the initially excited 1 ππ* state to the 1 nπ* state with a time constant of (0.13 ± 0.02) ps after an initial induction period of (0.12 ± 0.02) ps without population transfer, in quantitative agreement with the AIMS simulations. The population in the 1 nπ* state subsequently decays via intersystem crossing, likely mediated by a 3 ππ* state, within (3.17 ± 0.66) ps to a long-lived 3 nπ* state, which is presumed to be active towards Norrish type I chemistry.

Graßl, Martin [SLAC National Accelerator Laborator↗
DOE OSTI2020

Experimental investigation of internal two-phase flow structures and dynamics of quasi-stable sheet cavitation by fast synchrotron x-ray imaging

The quasi-stable sheet cavitation produced in a small Venturi channel is investigated using a fast synchrotron X-ray imaging technique aided with conventional high speed photography. The use of X-rays instead of visible light solves cavitation opacity related issues, and X-ray phase contrast-based edge enhancement enables high-definition visualization of the internal two-phase morphology. The simultaneous acquisition of time-resolved velocity and void fraction fields through post-processing the recorded X-ray images reveals, for the first time, the complex diphasic flow structures inside the sheet cavity, which is essentially divided into 6 characteristic parts. Distinct from the current mainstream view, the globallysteady sheet cavitation is found to be characterized by a weak but constantly-existing re-entrant flow that can penetrate the entire cavity. The turbulent velocity fluctuations inside the sheet cavity are also investigated. The turbulence level in the reverse flow region is observed to be as low as in the outer main flow demonstrating the relatively steady status of the reentrant flow. Unlike the streamwise and cross-stream fluctuations, the shear stress appears to be weakly correlated with the velocity gradient. The collapse of vapor phase and the vaporization at the upstream cavity interface are found to be the primary causes of shear stress intensification.

42 ENGINEERING↗
DOE OSTI2023

Single-Cell Discovery and Multiomic Characterization of Therapeutic Targets in Multiple Myeloma

Multiple myeloma (MM) is a highly refractory hematologic cancer. Targeted immunotherapy has shown promise in MM but remains hindered by the challenge of identifying specific yet broadly representative tumor markers. We analyzed 53 bone marrow (BM) aspirates from 41 MM patients using an unbiased, high-throughput pipeline for therapeutic target discovery via single-cell transcriptomic profiling, yielding 38 MM marker genes encoding cell-surface proteins and 15 encoding intracellular proteins. Of these, 20 candidate genes were highlighted that are not yet under clinical study, 11 of which were previously uncharacterized as therapeutic targets. The findings were cross-validated using bulk RNA sequencing, flow cytometry, and proteomic mass spectrometry of MM cell lines and patient BM, demonstrating high overall concordance across data types. Independent discovery using bulk RNA sequencing reiterated top candidates, further affirming the ability of single-cell transcriptomics to accurately capture marker expression despite limitations in sample size or sequencing depth. Target dynamics and heterogeneity were further examined using both transcriptomic and immuno-imaging methods. In summary, this study presents a robust and broadly applicable strategy for identifying tumor markers to better inform the development of targeted cancer therapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Experimental Measurements of Soot and Soot Precursors of Sustainable Aviation Fuels

Sustainable Aviation Fuels (SAFs) are promising for decarbonizing the challenging-to-electrify aerospace sector and reducing yearly greenhouse gas emissions by thousands of tons. Developing predictive detailed chemical kinetic models for SAF (and traditional) fuel mixtures necessitates a robust experimental database to elucidate the intricate pyrolysis and oxidation chemistry of SAF fuel blend components. Two critical components of SAF surrogate mixtures are iso-dodecane isomers and trans-decalin. This report characterizes the structure of three slightly sooting, laminar, non-premixed, nitrogen-diluted Planar Mixing Layer Flames (PMLFs) fueled by either pure ethylene, ethylene doped with 1500ppm of either 2,2,4,6,6-pentamethyl-heptane or trans-decalin. The flames have the same stoichiometric mixture fraction and total mole fraction of fuels in the fuel stream (X F,F = X C2H4,F + X Dopant,F = 0.26), resulting in approximately the same maximum temperature (T max ≈ 1800 K) to highlight the effect of doping on soot and its gaseous precursors. The PMLF is unconfined and established between adjacent planar jets of the fuel and oxidizer streams surrounded by annular shielding nitrogen. The flow is stabilized by exhausting the hot buoyant flame products through a slot in a downstream plate onto which the cold flow impinges. Importantly, any horizontal cross-section of the PMLF has a self-similar thermochemical structure of a reactive boundary layer that grows at increasing Height Above the Burner (HAB). Therefore, the horizontal cross-sections of a PMLF can be modeled as a One-Dimensional Counterflow Flame (1D-CF) with very low strain rates which are unachievable in traditional CFs and yield a structure that is several millimeters thick. In brief, the PMLF is ideal for performing spatially resolved measurements with minimal effects of the sampling-induced perturbations, to develop and validate detailed chemical kinetic models on time scales of a few tens of milliseconds that are relevant in aviation applications. The horizontal profiles of C0-C18 gaseous species at HAB = 50 mm and soot volume fraction at 25 and 50 mm are measured using capillary sampling followed by GC-MS analyses and Laser Induced Emission Spectroscopy (LIES), respectively. Additionally, the Elastic Laser Light Scattering (E-LLS) coefficient is measured at HAB = 50 mm and 25 mm to determine the strain rates of the 1D-CFs equivalent to the characterized PMLF horizontal cross sections and the profile of the E-LLS equivalent diameter of soot. The comparison of the results in the base PMLF fueled only by ethylene with those in either doped PMLF in which 1500 ppm of ethylene are replaced with 2,2,4,6,6-pentamethyl-heptane or trans-decalin, reveals that the doping increases the mole fraction of several Polycyclic Aromatic Hydrocarbons (PAHs) at both HABs by a factor of approximately 2.5, and the soot volume fraction at both HABs by factors of ~1.5 for 2,2,4,6,6-pentamethyl-heptane and ~1.25 for trans-decalin doping, respectively. Conversely, doping causes a near doubling of the soot equivalent diameter only at HAB = 50 mm, without significant effects at 25 mm. The experimental results partially validate the tested state-of-the-art detailed chemical kinetic model and also point toward the further improvement of its predictive capabilities.

10 SYNTHETIC FUELS↗
DOE OSTI2022

Header for high-pressure heat exchanger

A header for a high-pressure heat exchanger includes a first high-pressure transition section with inlets for multiple first high-pressure flow channels that are spaced from one another in a radial direction and collectively arranged in a substantially circular shape. The inlets for the multiple first high-pressure flow channels on a radially outer edge of the first high-pressure transition section are spaced further apart in a circumferential direction from adjacent inlets of the multiple first high-pressure flow channels than radially inward inlets are spaced from adjacent radially inward inlets of the multiple first high-pressure flow channels. The header also includes multiple first high-pressure flow channels extending from the first high-pressure transition section to a second-high pressure transition section that is configured to divide each of the multiple first high-pressure flow channels into at least two first high-pressure sub-flow channels. Each flow channel of the multiple first high-pressure flow channels and each sub-flow channel of the first high-pressure sub-flow channels have a round cross-sectional shape.

Kirsch, Kathryn L.↗
DOE OSTI2024

Curves, Coriolis, and Cross-Channel Circulation in the Hudson River Estuary

Abstract Despite its relatively small magnitude, cross-channel circulation in estuaries can influence the along-channel momentum balance, dispersion, and transport. We investigate spatial and temporal variation in cross-channel circulation at two contrasting sites in the Hudson River estuary. The two sites differ in the relative strength and direction of Coriolis and curvature forcing. We contrast the patterns and magnitudes of flow at the two sites during varying conditions in stratification driven by tidal amplitude and river discharge. We found well-defined flows during flood tides at both sites, characterized by mainly two-layer structures when the water column was more homogeneous and structures with three or more layers when the water column was more stratified. Ebb tides had generally weaker and less definite flows, except at one site where curvature and Coriolis reinforced each other during spring tide ebbs. Cross-channel currents had similar patterns, but were oppositely directed at the two sites, demonstrating the importance of curvature even in channels with relatively gradual curves. Coriolis and curvature dominated the measured terms in the cross-channel momentum balance. Their combination was generally consistent with driving the observed patterns and directions of flow, but local acceleration and cross-channel advection made some notable contributions. A large residual in the momentum balance indicates that some combination of vertical stress divergence, baroclinic pressure gradients, and along-channel and vertical advection must play an essential role, but data limitations prevented an accurate estimation of these terms. Cross-channel advection affected the along-channel momentum balance at times, with implications for the exchange flow’s strength. Significance Statement Currents that flow across the channel in an estuary move slower than those flowing along the channel, but they can transport materials and change water properties in important ways, affecting human uses of estuaries such as shipping, aquaculture, and recreation. We wanted to better understand cross-channel currents in the Hudson River estuary. We found that larger tides produced the strongest cross-channel currents with a two-layer pattern, compared to weaker currents with three layers during smaller tides. Higher or lower river flow also affected current strength. Comparing two locations, we saw cross-channel currents moving in opposite directions because of differences in the curvature of the river channel. Our results show how channel curvature and Earth’s rotation combine to produce cross-channel currents.

Conley, Margaret M.↗
DOE OSTI2026

CFD Modeling of an 81" Single Slope Barrier Wall Drainage Window for on Grade and Sag Conditions

The main goal of this study is to design a drainage window in an 81" single slope barrier wall for on-grade and sag conditions that could replace the standard ODOT I-3 barrier catch basin. The ODOT design guidelines are as follows: a) Examine configurations of the opening to the flume behind the wall under the wall section to maximize conveyance of flow to a roadside ditch section. b) Examine window lengths and configurations to maximize hydraulic efficiency on a grade up to an absolute maximum window length of 20 ft. c) Use a minimum opening height of 4 in. and consider tapering to an increased height at the outlet to the flume to reduce the potential for debris clogging beneath the barrier section. On-grade locations are analyzed for longitudinal slope up to 5% and spread up to 10 ft. Sag locations, represented by pavement at zero longitudinal slope, are analyzed for water depths up to 18 in. The roadway cross section consists of a 10-ft shoulder at a 4.0% cross slope and pavement at a 1.6% cross slope. All surfaces are represented with Manning’s n = 0.015. Computational fluid dynamics (CFD) is used to perform the analysis. The modeling approach follows previous studies by the authors. The new design is similar to a slotted inlet: a rectangular opening cut into the bottom of a barrier (e.g., a temporary barrier or single slope barrier wall), which typically has a uniform cross section. Findings from laboratory tests of flow under Florida DOT barrier walls were used for additional validation of the modeling approach used in this study.

42 ENGINEERING↗
DOE OSTI2021

Acceleration of Oxygen Ions in Dipolarization Events: 1. CPS Distributions

Abstract Using an MHD simulation of near tail reconnection associated with a flow burst and the collapse (dipolarization) of the inner tail in combination with test particle tracing we study the acceleration and flux increases of energetic oxygen ions (O + ). The characteristic orbits, distributions, and acceleration mechanisms are governed by the dimensionless parameter σ = ω ci t n , where ω ci is the ion gyro frequency and t n a characteristic Alfvén time of the MHD simulation. For σ < 1, central plasma sheet (CPS) populations after the passage of the dipolarization front are found to resemble half‐shells in velocity space oriented toward dusk. They originate from within the CPS and are energized typically by a single encounter of the region of enhanced cross‐tail electric field associated with the flow burst. For larger σ values ( σ > 1) the O + distributions resemble more closely those of protons, consisting of two counter‐streaming field‐aligned beams and an, albeit more tenuous and irregular, ring population perpendicular to the magnetic field. The existence of the beams, however, depends on suitable earthward moving source populations in the plasma sheet boundary layer or the adjacent lobes. The acceleration to higher energies is found to indicate a charge dependence, consistent with a dominance of more highly charged ions at energies of a few hundred keV. As in earlier simulations, the simulated fluxes show large anisotropies and nongyrotropic effects, phase bunching, and spatially and temporally localized beams.

Birn, J.↗