Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Adiabatic approximation”

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

Vortex Structure of Head Bubble in Convective Cloud Starting Plume

Abstract A developing cumulus cloud (Cu) was modeled, and dynamic, thermodynamic, and microphysical properties of an ascending head bubble reproducing the upper part of a developing Cu were investigated. The data for analysis are taken from 10-m-resolution LES of trade wind Cu under BOMEX conditions. The detection of a rising bubble is carried out using wavelet filtering of the velocity fields and microphysical fields, while a low-frequency signal of the filtering is associated with the convective-scale structure of cloud. We substantiate and discuss the representation of the bubble as a vortex ring, and estimate the parameters of this vortex ring. The simplest Hill’s vortex was chosen as a model of a vortex ring inside cloud. Analytical approximations of the radial profiles of the vertical velocity and of conservative quantities (such as total water mixing ratio and liquid water potential temperature inside and outside the bubble) are obtained. The spatial structure of these quantities is investigated using analytical expressions. Analytical models for spatial distributions of liquid water content (LWC) and adiabatic fraction (AF) are also designed and analyzed. The results demonstrate the existence of a cloud core with high values of LWC and AF up to the height of 1800 m. The horizontally averaged value of the adiabatic fraction, calculated analytically using the Hill’s vortex concept, is evaluated as 0.39, which is the typical AF value in the upper parts of such Cu. The vertical profiles of different important quantities characterizing cloud structure are presented. The analysis performed in this study allows us to conclude that a rising vortex ring plays the dominating role in formation of the thermodynamic and microphysical structure of developing Cu. Significance Statement 1) Dynamic and thermodynamic fields of a developing cumulus cloud simulated by high-resolution LES with spectral bin microphysics are separated into convective and turbulent components by means of the wavelet technique. 2) The analysis of convective component of the cloud revealed the existence of a vortex ring at the developing stage of the cloud and evaluate its parameters. 3) The analysis performed in this study allows us to conclude that a rising vortex ring plays the important role in formation of the thermodynamic and microphysical structure of developing Cu. 4) The study provides a novel insight into the cloud–environment interaction. 5) The approximating equations describing the vortex ring can be usefully applied for developing new schemes of convective parameterization.

Meteorology & Atmospheric Sciences↗
NASA NTRS1993

Mass loading and velocity diffusion models for heavy pickup ions at comet Grigg-Skjellerup

We compare model predictions of cometary water group ion densities and the solar wind slow down with measurements made by the Giotto Johnstone plasma analyzer implanted ion sensor at the encounter with comet Grigg-Skjellerup (G-S) on July 10, 1992. The observed slope of the ion density profile on approach to the comet is unexpectedly steep. Possible explanations for this are discussed. We present also a preliminary investigation of the quasilinear velocity-space diffusion of the implanted heavy ion population at G-S using a transport equation including souce, convection, adiabatic compression, and velocity diffusion terms. Resulting distributions are anisotropic, in agreement with observations. We consider theoretically the waves that may be generated by the diffusion process for the observed solar wind conditions. At initial ion injections, waves are generated at omega approximately Omega(sub i) the ion gyrofrequency, and lower frequencies are predicted for diffusion toward a bispherical shell.

Huddleston, D. E.↗
NASA NTRS2003

On the Critical Behaviour, Crossover Point and Complexity of the Exact Cover Problem

Research into quantum algorithms for NP-complete problems has rekindled interest in the detailed study a broad class of combinatorial problems. A recent paper applied the quantum adiabatic evolution algorithm to the Exact Cover problem for 3-sets (EC3), and provided an empirical evidence that the algorithm was polynomial. In this paper we provide a detailed study of the characteristics of the exact cover problem. We present the annealing approximation applied to EC3, which gives an over-estimate of the phase transition point. We also identify empirically the phase transition point. We also study the complexity of two classical algorithms on this problem: Davis-Putnam and Simulated Annealing. For these algorithms, EC3 is significantly easier than 3-SAT.

Morris, Robin D.↗
DOE OSTI2022

Conformer-Specific Dissociation Dynamics in Dimethyl Methylphosphonate Radical Cation

The dynamics of the dimethyl methylphosphonate (DMMP) radical cation after production by strong field adiabatic ionization have been investigated. Pump-probe experiments using strong field 1300 nm pulses to adiabatically ionize DMMP and a 800 nm non-ionizing probe induce coherent oscillations of the parent ion yield with a period of about 45 fs. The yields of two fragments, PO 2 C 2 H 7 + and PO 2 CH 4 + , oscillate approximately out of phase with the parent ion, but with a slight phase shift relative to each other. We use electronic structure theory and nonadiabatic surface hopping dynamics to understand the underlying dynamics. The results show that while the cation oscillates on the ground state along the P=O bond stretch coordinate, the probe excites population to higher electronic states that can lead to fragments PO 2 C 2 H 7 + and PO 2 CH 4 + . The computational results combined with the experimental observations indicate that the two conformers of DMMP that are populated under experimental conditions exhibit different dynamics after being excited to the higher electronic states of the cation leading to different dissociation products. These results highlight the potential usefulness of these pump-probe measurements as a tool to study conformer-specific dynamics in molecules of biological interest.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1972

Diffusion models for Jupiter's radiation belt

Solutions are given for the diffusion of trapped particles in a planetary magnetic field in which the first and second adiabatic invariants are preserved but the third is not, using as boundary conditions a fixed density at the outer boundary (the magnetopause) and a zero density at an inner boundary (the planetary surface). Losses to an orbiting natural satellite are included and an approximate evaluation is made of the effects of the synchrotron radiation on the energy of relativistic electrons. Choosing parameters appropriate to Jupiter, the electrons required to produce the observed synchrotron radiation are explained. If a speculative mechanism in which the diffusion is driven by ionospheric wind is the true explanation of the electrons producing the synchrotron emission it can be concluded that Jupiter's inner magnetosphere is occupied by an energetic proton flux that would be a serious hazard to spacecraft.

Jacques, S. A.↗
NASA NTRS2011

Evolution of the Equatorial Oscillation in Saturn's Stratosphere Between 2005 and 2010 from Cassini/CIRS Limb Data Analysis

We present an analysis of thermal infrared spectra acquired in limb viewing geometry by Cassini/CIRS in February 2010. We retrieve vertical profiles of Saturn's stratospheric temperature from 20 hPa to 10 (exp -2) hPa, at 9 latitudes between 20 deg N and 20 deg S. Using the gradient thermal wind equation, we derive a map of the zonal wind field. Both the temperature and the zonal wind vertical profiles exhibit an oscillation in the equatorial region. These results are compared to the temperature and zonal wind maps obtained from 2005-2006 CIRS limb data, when this oscillation was first reported. In both epochs, strong temperature anomalies at the equator (up to 20K) are consistent with adiabatic heating (cooling) due to a sinking (rising) motion at a speed of 0.1 - 0.2 mm/s. Finally, we show that the altitude of the maximum eastward wind has moved downwards by 1.3 scale heights in 4.2 years, hence with a 'phase' speed of approximately 0.5 mm/s. This rate is consistent with the estimated period of 14.7 years for the equatorial oscillation, and requires a local zonal acceleration of 1.1 x 10(exp -6) m.s(exp -2) at the 2.5 hPa pressure level. This downward propagation of the oscillation is consistent with it being driven by absorption of upwardly propagating waves.

Guerlet, S.↗
DOE OSTI2023

Symmetric Kondo Lattice States in Doped Strained Twisted Bilayer Graphene

We use the topological heavy fermion (THF) model and its Kondo lattice (KL) formulation to study the possibility of a symmetric Kondo (SK) state in twisted bilayer graphene. Via a large-N approximation, we find a SK state in the KL model at fillings v = 0,±1,±2 where a KL model can be constructed. In the SK state, all symmetries are preserved and the local moments are Kondo screened by the conduction electrons. At the mean-field level of the THF model at v = 0,±1,±2,±3 we also find a similar symmetric state that is adiabatically connected to the symmetric Kondo state. We study the stability of the symmetric state by comparing its energy with the ordered (symmetry-breaking) states found in [H. Hu et al., Phys. Rev. Lett. 131, 026502 (2023)., Z.-D. Song and B. A. Bernevig, Phys. Rev. Lett. 129, 047601 (2022).] and find the ordered states to have lower energy at v = 0,±1,±2. However, moving away from integer fillings by doping the light bands, our mean-field calculations find the energy difference between the ordered state and the symmetric state to be reduced, which suggests the loss of ordering and a tendency toward Kondo screening. In order to include many-body effects beyond the mean-field approximation, we also performed dynamical mean-field theory calculations on the THF model in the nonordered phase. The spin susceptibility follows a Curie behavior at v = 0,±1,±2 down to ~2 K where the onset of screening of the local moment becomes visible. This hints to very low Kondo temperatures at these fillings, in agreement with the outcome of our mean-field calculations. At noninteger filling v = ±0.5,±0.8,±1.2 dynamical mean-field theory shows deviations from a 1/T susceptibility at much higher temperatures, suggesting a more effective screening of local moments with doping. Lastly, we study the effect of a C 3⁢z -rotational-symmetry-breaking strain via mean-field approaches and find that a symmetric phase (that only breaks C 3⁢z symmetry) can be stabilized at sufficiently large strain at v = 0,±1,±2. Our results suggest that a symmetric Kondo phase is strongly suppressed at integer fillings, but could be stabilized either at noninteger fillings or by applying strain.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Nonadiabatic transitions during a passage near a critical point

The passage through a critical point of a many-body quantum system leads to abundant nonadiabatic excitations. Here, we explore a regime, in which the critical point is not crossed although the system is passing slowly very close to it. Here, we show that the leading exponent for the excitation probability can then be obtained by standard arguments of the Dykhne formula, but the exponential prefactor is no longer simple and behaves as a power law on the characteristic transition rate. We derive this prefactor for the nonlinear Landau–Zener model by adjusting Dykhne’s approach. Then, we introduce an exactly solvable model of the transition near a critical point in the Stark ladder. We derive the number of excitations for it without approximations and find qualitatively similar results for the excitation scaling.

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

Toward efficient elastocaloric systems: Predicting material thermal properties with high fidelity

A critical need to accurately model thermal behaviors of materials that exhibit strong elastocaloric effects, including predicting the effects themselves at varying stresses and temperatures, has been addressed using a simple and versatile physics-based approach. The key factor leading to the high precision is approximating the underlying elastic phase transition as a smooth modification of lattice entropy between coexisting phases. Once the phase transformation entropy is modeled to match experimentally measured strain as a function of temperature and applied stress, estimating the heat capacity, entropy, and isothermal entropy and adiabatic temperature changes in temperature-stress coordinates becomes straightforward. This approach provides insight into how thermal properties of elastocaloric materials vary through the transition based on strain measurements that are simple to perform and interpret. In addition to aiding in the rapid evaluation of new and existing elastocaloric materials, this advance is expected to prove invaluable for accurate heat transfer modeling aimed at designing efficient regenerative elastocaloric cooling devices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Diabatic quantum annealing for the frustrated ring model

Abstract Quantum annealing (QA) is a continuous-time heuristic quantum algorithm for solving or approximately solving classical optimization problems. The algorithm uses a schedule to interpolate between a driver Hamiltonian with an easy-to-prepare ground state and a problem Hamiltonian whose ground state encodes solutions to an optimization problem. The standard implementation relies on the evolution being adiabatic: keeping the system in the instantaneous ground state with high probability and requiring a time scale inversely related to the minimum energy gap between the instantaneous ground and excited states. However, adiabatic evolution can lead to evolution times that scale exponentially with the system size, even for computationally simple problems. Here, we study whether non-adiabatic evolutions with optimized annealing schedules can bypass this exponential slowdown for one such class of problems called the frustrated ring model. For sufficiently optimized annealing schedules and system sizes of up to 39 qubits, we provide numerical evidence that we can avoid the exponential slowdown. Our work highlights the potential of highly-controllable QA to circumvent bottlenecks associated with the standard implementation of QA.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Lattice holography on a quantum computer

We explore the potential application of quantum computers to the examination of lattice holography, which extends to the strongly coupled bulk theory regime. With adiabatic evolution, we compute the ground state of a spin system on a ( 2 + 1 )-dimensional hyperbolic lattice, and measure the spin-spin correlation function on the boundary. Notably, we observe that with achievable resources for coming quantum devices, the correlation function demonstrates an approximate scale-invariant behavior, aligning with the pivotal theoretical predictions of the anti–de Sitter/conformal field theory correspondence. Published by the American Physical Society 2024

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Observation of subcritical shocks in a collisional laboratory plasma: scale dependence near the resistive length

We present a study of subcritical shocks in a highly collisional laboratory plasma with a dynamically significant magnetic field. Shocks were produced by placing cylindrical obstacles into the supermagnetosonic ( $M_{{\rm MS}} \sim 1.9$ ) outflow from an inverse wire array z-pinch at the MAGPIE pulsed power facility ( $n_{e} \sim 8.5 \times 10^{17}\,{\rm cm}^{-3}$ , $v \sim 45\,{\rm km}\,{\rm s}^{-1}$ ). We demonstrate the existence of subcritical shocks in this regime and find that secondary stagnation shocks form in the downstream which we infer from interferometry and optical Thomson scattering measurements are hydrodynamic in nature. The subcritical shock width is found to be approximately equal to the resistive diffusion length and we demonstrate the absence of a jump in hydrodynamic parameters. Temperature measurements by collective optical Thomson scattering showed little temperature change across the subcritical shock ( ${<}10\,\%$ of the ion kinetic energy) which is consistent with a balance between adiabatic and Ohmic heating and radiative cooling. We demonstrate the absence of subcritical shocks when the obstacle diameter is less than the resistive diffusion length due to decoupling of the magnetic field from the plasma. These findings are supported by magnetohydrodynamic simulations using the Gorgon and AstroBEAR codes and discrepancies between the simulations and experiment are discussed.

Physics↗
NASA NTRS1993

Stormtime transport of ring current and radiation belt ions

This is an investigation of stormtime particle transport that leads to formation of the ring current. Our method is to trace the guiding-center motion of representative ions (having selected first adiabatic invariants mu) in response to model substorm-associated impulses in the convection electric field. We compare our simulation results qualitatively with existing analytically tractable idealizations of particle transport (direct convective access and radial diffusion) in order to assess the limits of validity of these approximations. For mu approximately less than 10 MeV/G (E approximately less than 10 keV at L equivalent to 3) the ion drift period on the final (ring-current) drift shell of interest (L equivalent to 3) exceeds the duration of the main phase of our model storm, and we find that the transport of ions to this drift shell is appropriately idealized as direct convective access, typically from open drift paths. Ion transport to a final closed drift path from an open (plasma-sheet) drift trajectory is possible for those portions of that drift path that lie outside the mean stormtime separatrix between closed and open drift trajectories, For mu approximately 10-25 MeV/G (110 keV approximately less than E approximately less than 280 keV at L equivalent to 3) the drift period at L equivalent to 3 is comparable to the postulated 3-hr duration of the storm, and the mode of transport is transitional between direct convective access and transport that resembles radial diffusion. (This particle population is transitional between the ring current and radiation belt). For mu approximately greater than 25 MeV/G (radiation-belt ions having E approximately greater than 280 keV at L equivalent to 3) the ion drift period is considerably shorter than the main phase of a typical storm, and ions gain access to the ring-current region essentially via radial diffusion. By computing the mean and mean-square cumulative changes in 1/L among (in this case) 12 representative ions equally spaced in drift time around the steady-state drift shell of interest (L equivalent to 3), we have estimated (from both our forward and our time-reversed simulations) the time-integrated radial-diffusion coefficients D(sup sim)(sub LL) for particles having selected values of mu approximately greater than 15 MeV/G. The results agree surprisingly well with the predictions (D(sup ql)(sub LL)) of quasilinear radial diffusion theory, despite the rather brief duration (approximately 3 hrs) of our model storm and despite the extreme variability (with frequency) of the spectral-density function that characterizes the applied electric field during our model storm. As expected, the values of D(sup sim)(sub LL) deduced (respectively) from our forward and time-reversed simulations agree even better with each other and with D(sup sim)(sub LL) when the impulse amplitudes which characterize the individual substorms of our model storm are systematically reduced.

Chen, Margaret W.↗
NASA NTRS1994

Explosions of infalling comets in Jupiter's atmosphere

In view of the expected collision of comet Shoemaker-Levy 9 (1993e) with Jupiter in 1994 July, we calculate basic properties of the initial interaction for a simplified Jovian atmosphere. The comet is expected to impact Jupiter at 60 km/sec and at an angle of 45 deg to the zenith. The shock wave generated by the bolide should be optically thick once it has penetrated to an atmospheric density approximately 10(exp -6) gr/cm(exp 3), and we calculate the post-shock conditions assuming local thermodynamic equilibrium (LTE) for shock velocities v(sub sh) in the range 10 to 60 km/sec and preshock densities rho(sub a) = 10(exp -6) to 10(exp -2) gr/cm(exp 3). Our shock calculations include molecular hydrogen, atomic hydrogen, ionized hydrogen, neutral helium, and singly ionized helium. Even at the highest shock velocity, the gas is only partially ionized and the postshock temperature rises with preshock density in order to maintain the ionization. The value of the effective shock adiabatic index gamma(sub sh) varies from 1.17 (at low v(sub sh) and rho(sub a) to 1.40 (at high v(sub sh) and rho(sub a). The ablation rate is limited by the radiative flux that reaches the bolide surface. We argue that the ablated gas does not efficiently transfer its kinetic energy to the atmosphere, and it ultimately slows in a similar fashion to the comet material. As the bolide initially falls through the atmosphere, the character of the shock emission changes. At rho(sub a) approximately 10(exp -8) gr/cm (exp 3), the gas is optically thin and we expect line emission; in the optical spectrum, Balmer emission is expected from the shocked atmosphere and low-ionization metal lines from ablated cometary material. At rho(sub a) approximately 10(exp -6) gr/m(exp 3), the shocked gas is optically thick and the shock front near the bolide produces a blackbody spectrum. The temperature is favorable for ultraviolet (1000 to 3000 A) emission and the luminosity may be approximately 5 x 10(exp 23) ergs/sec for approximately 0.6 sec for a bolide 1 km in radius. At rho(sub a) approximately 10(exp -4) gr/cm(exp 3), the bolide has passed below the ultraviolet photosphere. The shock front emits considerable ionizing radiation, but it is absorbed in a narrow preshock region. The bolometric correction for the optical luminosity is large and we expect a 3000 to 8000 A luminosity of approximately 3 x 10(exp 23) ergs/sec for approximately 1 sec. The optical emission is strongly peaked in the vicinity of the bolide. The bolide does have a somewhat less luminous, optically thick trail extending greater than or equal to 10 km, but the radiation is characterized by a temperature of 4000 to 5000 K. From the fragmentation model of Chyba, Thomas, & Zahnle (1993), the bolide deposits most of its kinetic energy at rho(sub a) approximately 10(exp -3) gr/cm(exp 3) and this is the effective explosion site. The shock wave from such an explosion can move up about one density scale height. We examine the breakout of the shock front from the Jovian atmosphere and find that the shock acceleration in the decreasing density region is slow, so that the energy flux in the shock front is small. Higher velocities might be generated by shock acceleration along the channel left by the bolide if the shock motion can occur before the channel closes off as a result of radiative cooling. Hot gas created by the explosion ultimately rises due to buoyancy on a timescale of a minute. The luminosity is highest when the bubble first rises into the optically thin part of the atmosphere and may be approximately 1 x 10(exp 25) ergs/sec in the near-infrared. Roughly 1% of the initial bolide energy may be radiated in this way; the rest of the energy is lost to sound waves from the initial explosion and to work done by the bubble on the surrounding atmosphere.

Chevalier, Roger A.↗
DOE OSTI2023

Multimode theory of electron hole transverse instability

We present Vlasov–Poisson three-dimensional linear stability analysis of an initially planar electron hole structure, solving for the distribution function by integration along unperturbed orbits. The non-sinusoidal potential perturbation shape (parallel to $B$ ) is expanded in eigenfunctions of the adiabatic Poisson operator, generalizing the prior assumption of a rigid shift of the equilibrium. We show that the shiftmode is then modified by a second discrete mode plus an integral over a continuum of wave-like modes. A rigorous treatment shows that the continuum can be approximated effectively by a single mode that satisfies the external wave dispersion relation, thus making the perturbation a weighted sum of three modes. We find numerically the solution for the complex instability frequency, and the corresponding three mode amplitudes determining the perturbation eigenmode. This multimode analysis refines the accuracy of the prior single-mode results, giving slightly higher growth rates at most parameters, as expected from the extra mode shape freedom. Oscillating modes near stability boundaries have larger mode distortions which help explain particle-in-cell simulations that observe instability up to ${\sim }20$ % beyond the prior shiftmode thresholds, and narrowing of the perturbation. At high magnetic field, the multimode analysis predicts a reduction of the already small growth rate.

Physics↗
NASA NTRS1991

The Uranian system - An overview

Results of many years of observations from the ground and from the Voyager encounter with Venus are presented which clarify a number of issues that were subjects of speculation or partial information. The new issues that the encounter brought to light are outlined. The encounter revealed additional rings, extensive dust associated with them, and complexities, including azimuthal asymmetries. Uranus possesses a strong (about 1 gauss) magnetic field with a bizarre geometry that can be described as a dipole offset 0.3 RU from the center of the planet and tilted approximately 60 deg with respect to the axis of rotation. It is argued that it is produced by a magnetohydrodynamic dynamo at about 0.4 RU from the planet's center. It is inferred from this that at least part of the interior is a convective fluid, indicating an adiabatic temperature profile and a warm interior. The magnetic field produces a magnetosphere with some peculiar characteristics, owing to the unique orientation of the dipole axis with respect to the direction of the solar wind flow.

Bergstralh, Jay T.↗
NASA NTRS1992

Quasi-linear regime of gravitational instability as a clue to understanding the large-scale structure in the Universe

In the late seventies, an image of the large-scale structure in the Universe began to emerge as a result of the accumulation of the galaxy redshifts. Most of the galaxies are found to concentrate in large filaments and perhaps sheets leaving most of the volume empty. Similar structures were predicted theoretically in the frame of the adiabatic theory of galaxy formation (Zeldovich) and later in the hot dark matter cosmology. However, both scenarios have been ruled out by the observations. With these scenarios the dynamical part of the scenario was also erroneously rejected by many as well. In this talk, I derive the Zeldovich approximation from the exact dynamic equations and show that it is always better than the standard linear approximation. The advantage of the Zeldovich approximation is the greatest in the quasi-linear regime when delta(sub rms) is less than 1 (delta identical to delta(rho)/rho), but the displacement of the matter is essential. The range of scales in the quasi-linear regime depends upon the slope of the initial spectrum and increases with decreasing n, where n is the exponent, if the initial spectrum is approximated by a simple power law P(k) varies as k(exp n).

Shandarin, Sergei F.↗
NASA NTRS2005

An Expanding Radio Nebula Produced by a Giant Flare from the Magnetar SGR 1806-20

Soft gamma repeaters (SGRs) are "magnetars", a small class of slowly spinning neutron stars with extreme surface magnetic fields, B approx. 10(sup 15) gauss. On 2004 December 27, a giant flare was detected from the magnetar SGR 1806-20, the third such event ever recorded. This burst of energy, which resulted in the highest flux of gamma-rays ever measured from a celestial object, was detected by a variety of instruments and even caused an ionospheric disturbance in the Earth's upper atmosphere recorded around the globe. Here we report the detection of a very bright but rapidly-fading radio source at the position of SGR 1806-20 following this outburst. From day 6 to day 19 after the flare, we see a resolved, linearly polarized, radio nebula, expanding at a velocity of approximately 0.3c. To create this nebula, at least 4 x 10(exp 43) ergs of energy must have been emitted by the giant flare in the form of magnetic fields and relativistic particles. The steep decline of the radio flux may indicate episodic particle acceleration followed by adiabatic expansion, as could occur if the ejecta have energized a thin shell surrounding a preexisting cavity.

Gaensler, B. M.↗