Search NASA⌕ Search

SEARCH · Search NASA

Results for “TRIANGULATION”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

BDEM (Discrete-element-simulator for high-solids granular flows) [SWR-22-72]

BDEM is a discrete element method based simulation tool developed specifically for modeling high-solids granular flows that include polydispersity, heat-transfer, moving boundaries and chemistry. Our solver provides facilities for simulating spherical/non-spherical particles with modified contact and friction models in complex dynamic geometries defined using level-sets or triangulated files. The solver is developed on top of NREL's open-source performance portable library, AMReX, providing parallel execution capabilities on current and upcoming high-performance-computing (HPC) architectures. Simulations at the scale of several millions to billion particles have been performed using this software on large scale computing resources. This software can be applied to non-reacting solids dominant flows in silos, hoppers and screw conveyors as well as in high temperature reacting systems such as screw kilns and auger reactors.

Sitaraman, Hariswaran↗
DOE OSTI2020

pnnl/pyraabe

Automated Raabe table generationPyRaabe is a Python package to automate Raabe table -- a graph representation of airway geometry -- generation from a triangulated surface of an airway geometry. With a single command, it extracts a centerline from an input 3D surface mesh and subsequently computes the relevant attributes for the Raabe table, including length, diameter, bifurcation angle, and gravity angle. These can be used to inform downstream models, such as the Multiple-Path Particle Dosimetry (MPPD) model.

Colby, Sean M.↗
DOE OSTI2021

Model, data, and code for paper "Modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x"

The data package includes data, model, and code that support the analyses and conclusions in the paper titled “modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x”. The primary goal of this paper is to demonstrate that key streamflow properties such as water depth, flow velocity, and dynamic pressure in a natural river at 30-kilometer scale over 5 years can be reliably and efficiently modeled using the computational framework presented in this paper. To support the paper, various data types from remote sensing, field observations, and computational models are used. Specific details are described as follows. Firstly, the river bathymetry data was obtained from a Light Detection and Ranging (LiDAR) survey. This data is then converted to a triangulated surface format, STL, for mesh generation in OpenFOAM. The STL data can be found in Model_Setups/BaseCase_2013To2015/constant/triSurface. The OpenFOAM mesh generated using this STL file can be found in constant/polyMesh. Other model setups, boundary and initial conditions can be found in /system and /0.org under folder BaseCase_2013To2015. A similar data structure can also be found in BaseCase_2018To2019 for the simulations during 2018 and 2019. Secondly, the OpenFOAM simulations need the upstream discharge and water depth information at the upstream boundary to drive the model. These data are generated from a one-dimensional hydraulic model and the data can be found under the folder Model_Setups /1D model Mass1 data. The mass1_65.csv and mass1_191.csv files include the results of the 1D model at the model inlet and outlet, respectively. The Matlab source code Mass1ToOFBC20182019.m is used to convert these data into OpenFOAM boundary condition setups.With the above OpenFOAM model, it can generate data for water surface elevation, flow velocity, and dynamic pressure. In this paper, the water surface elevation was measured at 7 locations during different periods between 2011 and 2019. The exact survey locations (see Fig1_SurveyLocations.txt) can be found in folder Fig_1. The variation of water stage over time at the 7 locations can be found in folder /Observation_WSE. The data type include .txt, .csv, .xlsx, and .mat. The .mat data can be loaded by Matlab.We also measured the flow velocities at 12 cross-sections along the river. At each cross-section, we recorded the x, y locations, depth, three velocity components u,v,w. These data are saved to a Matlab format which can be found under folder /Observation_Velocity and /Fig_1. The relative locations of velocity survey locations to the river bathymetry can be found in Figure 1c.The water stage data at the 7 locations from OpenFOAM, 1D, and 2D hydraulic models are also provided to evaluate the long-term performance of 3D models vs 1D/2D models. The water stage data for the 7 locations from OpenFOAM have been saved to .mat format and can be found in /OpenFOAM_WSE. The water stage data from the 1D model are saved in .csv format and can be found in /Mass1_WSE. The water stage from the 2D model is saved as .mat format and can be found in / Mass2_WSEIn addition, the OpenFOAM model outputs the information of hydrostatic and hydrodynamic pressure. They are saved as .mat format under folder /Fig_11/2013_1. As the files are too large, we only uploaded the data for January 2013. The area of different ratio of dynamic pressure to static pressure for all simulation range, i.e., 2013-2015, are saved to .mat format. They can be found in /Fig_11/PA. Further, the data of wall clock time versus the solution time of the OpenFOAM modeling are also saved to .mat format under folder /Fig_13/LogsMat. In summary, the data package contains seven data types, including .txt, .csv, .xlsx, .dat, .stl, .m, and .mat. The former 4 types can be directly open using a text editor or Microsoft Office. The .mat format needs to be read by Matlab. The Matlab source code .m files need to be run with Matlab. The OpenFOAM setups can be visualized in ParaView. The .stl file can be opened in ParaView or Blender. The data in subfolders Fig_1 to Fig_10 and Fig_12 are copied from the aforementioned data folders to generate specific figures for the paper. A readME.txt file is included in each subfolder to further describe how the data in each folder are generated and used to support the paper.Please use the data package's DOI to cite the data package. Please contact yunxiang.chen@pnnl.gov if you need more data related to the paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Stereo Cameras for Clouds (STEREOCAM) Instrument Handbook

The three pairs of stereo camera setups aim to provide synchronized and stereo calibrated time series of images that can be used for 3D cloud mask reconstruction. Each camera pair is positioned at approximately 120 degrees from the other pair, with a 17o-19o pitch angle from the ground, and at 5-6 km distance from the U.S. Department of Energy (DOE) Central Facility at the Atmospheric Radiation Measurement (ARM) Climate Research Facility Southern Great Plains (SGP) observatory to cover the region from northeast, northwest, and southern views. Images from both cameras of the same stereo setup can be paired together to obtain 3D reconstruction by triangulation. 3D reconstructions from the ring of three stereo pairs can be combined together to generate a 3D mask from surrounding views. This handbook delivers all stereo reconstruction parameters of the cameras necessary to make 3D reconstructions from the stereo camera images.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Point Cloud of Cloud Points (PCCP) Value-Added Product Report

In this report, we refer to a pair of cameras that capture synchronized pictures with overlapping fields of view (FOV) as a stereo pair. Each stereo pair independently performs a stereo reconstruction of cloud points. Currently, there are three U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility stereo pairs (six cameras in total) positioned around the Southern Great Plains (SGP) observatory’s Central Facility (CF; Romps and Oktem 2017). Time-synchronized pictures from the two cameras in a stereo pair can be paired together to obtain a three-dimensional (3D) reconstruction of feature points by triangulation. This document explains how we use ARM stereo cameras and stereophotogrammetric principles to generate the Point Cloud of Cloud Points (PCCP) Value-Added Product (VAP). The PCCP VAP is essentially a set of 3D positions representing the locations of cloud features in the sky. Thousands of cloud features can be reconstructed instantly in each stereo pair's FOV, which covers an area of tens of square kilometers. Cloud base and cloud top heights can be extracted from the PCCP product.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Retrieving Point Cloud of Cloud Points (PCCP) Value-Added Product from Stereo Cameras

In this report, we refer to a pair of cameras that capture synchronized pictures with overlapping fields of view (FOV) as a stereo pair. Each stereo pair independently performs a stereo reconstruction of cloud points. Currently, there are three U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility stereo pairs (six cameras in total) positioned around the Southern Great Plains (SGP) observatory’s Central Facility (CF; Romps and Oktem 2017). Time-synchronized pictures from the two cameras in a stereo pair can be paired together to obtain a three-dimensional (3D) reconstruction of feature points by triangulation. This document explains how we use ARM stereo cameras and stereophotogrammetric principles to generate the Point Cloud of Cloud Points (PCCP) Value-Added Product (VAP). The PCCP VAP is essentially a set of 3D positions representing the locations of cloud features in the sky. Thousands of cloud features can be reconstructed instantly in each stereo pair's FOV, which covers an area of tens of square kilometers. Cloud base and cloud top heights can be extracted from the PCCP product.

42 ENGINEERING↗
DOE OSTI2021

Frequency Comb-Based Remote Methane Observation Network (Final Scientific/Technical Report)

ARPA-E investment in long-range, frequency comb-based methane detection has resulted in the first scalable, cost-effective continuous methane emissions system. At the onset of this ARPA-E MONITOR award, the dual-frequency comb spectrometer was a sprawling, expensive, laboratory-confined device which had only recently been demonstrated capable of measuring atmospheric trace gases over open-paths. Bringing the technology to the point of being able to monitor for and characterize methane emissions at oil and natural gas production facilities seemed an almost impossible task. Over the subsequent months and years, our team transitioned the Nobel-prize-winning frequency-comb-laser spectrometer technology to a field-deployed regional methane leak detection system. The spectrometer enables highly sensitive near-infrared absorption measurements of methane along kilometer-scale laser beam paths. The measurements are coupled with an atmospheric modeling and inversion framework to triangulate the location of emission sources and quantify the emission rates.

03 NATURAL GAS↗
DOE OSTI2022

A Mountaintop View Requires Minimal Sorting: A Faster Contour Tree Algorithm

Consider a scalar field f : M → R, where M is a triangulated simplicial mesh in R d . A level set, or contour, at value v is a connected component of f –1 (v). As v is changed, these contours change topology, merge into each other, or split. Contour trees are concise representations of f that track this contour behavior. The vertices of these trees are the critical points of f, where the gradient is zero. The edges represent changes in the topology of contours. It is a fundamental data structure in data analysis and visualization, and there is significant previous work (both theoretical and practical) on algorithms for constructing contour trees. Suppose M has n vertices, N facets, and t critical points. A classic result of Carr, Snoeyink, and Axen (2000) gives an algorithm that takes O(n log n+Nα(N)) time (where α(·) is the inverse Ackermann function). A further improvement to O(t log t + N) time was given by Chiang et al. All these algorithms involve a global sort of the critical points, a significant computational bottleneck. Unfortunately, lower bounds of Ω(t log t) also exist. We present the first algorithm that can avoid the global sort and has a refined time complexity that depends on the contour tree structure. Intuitively, if the tree is short and fat, we get significant improvements in running time. For a partition of the contour tree into a set of descending paths, P, our algorithm runs in O($\Sigma$ pϵP |p| log |p| + tα(t) + N). This is at most O(t log D + N), where D is the diameter of the contour tree. Moreover, it is O(tα(t) + N) for balanced trees, a significant improvement over the previous complexity. Our algorithm requires numerous ideas: partitioning the contour tree into join and split trees, a local growing procedure to iteratively build contour trees, and the use of heavy path decompositions for the time complexity analysis. There is a crucial use of a family of binomial heaps to maintain priorities, ensuring that any comparison made is between comparable nodes of the contour tree. We also prove lower bounds showing that the $\Sigma$ pϵP |p| log |p| complexity is inherent to computing contour trees.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Rapidly Deployable Acoustic Monitoring and Localization System Based on a Low-Cost Wave Buoy Platform

The primary objective of this project is to develop a cost-effective, fit-for-purpose environmental monitoring system, “NoiseSpotter®,” that characterizes, classifies, and provides accurate location information for anthropogenic and natural sounds in near real-time. NoiseSpotter was developed to support the evaluation of potential acoustic effects of marine energy (ME) projects. By utilizing a compact array of three acoustic particle motion sensors, NoiseSpotter triangulates individual bearings to provide sound source localization to within 5% accuracy, allowing the ability to discern ME device sounds relative to other confounding sounds in the environment, while providing location estimates to nearby marine mammals for environmental mitigation purposes. The ME industry needs proven solutions to meet environmental impact assessment needs. The NoiseSpotter® includes off-the-shelf, modular components that are easy to assemble and disassemble. Its acoustic particle motion sensors are commercially available and the data logger and real-time telemetry system is designed to be plug-and-play. The entire system is relatively compact and can be deployed from small vessels. NoiseSpotter’s near real-time capability enables operational monitoring of ME sounds, particularly during early stages of technology adoption to facilitate mitigation of potential noise effects. Widespread adoption of the technology for acoustic monitoring of ME devices requires that it be cost effective; hence the anticipated commercial cost of system hardware is $35,000. This project contributes to reducing barriers to ME testing through support of scientific research focused on reducing or mitigating environmental risks and lowering costs and complexity of environmental monitoring. This project has developed an acoustic monitoring system, NoiseSpotter® (U.S. Patent No. 11,156,734 and U.S. Registered Trademark No. 6,442,313), to detect and characterize baseline noise and sounds from ME operations and support geolocation of detected sounds. The intended outcome of the project is to mitigate concerns about the potential for ME device noise to alter marine mammal or fish behavior. NoiseSpotter® enables cost-effective, near real-time acoustic monitoring of an operational ME device relative to ambient environmental noise and provides a technical basis for ME developers seeking to navigate the permitting process in an efficient manner.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Influence of Vegetation Canopy Architecture on Dynamics of Spatiotemporal Light Distribution and its Effect on Land-Biosphere-Atmosphere Fluxes of Energy, Carbon, and Water

The dynamic spatiotemporal distribution of light and micro climatic variables (wind speed, atmospheric CO 2 concentration, humidity, pressure, and air temperature) within the plant canopy as determined by the canopy architecture play a critical role in land-biosphere-atmosphere (LBA) photosynthetic carbon uptake, transpiratory water loss and energy partitioning. Vegetation models that include effects of sun and shade flecks, leaf temperature and energy balance partitioning between sensible and latent heat fluxes, and reduced turbulent mixing in dense vegetation canopies can accurately predict canopy performance, particularly its shaded leaves, which can constitute 80% of leaf area and contribute up to 50% of carbon update and 65% of water loss. However, current vegetation models with their bulk characterization of shaded leaves poorly represent the complex three dimensional (3D) canopy architecture are insufficient to capture the spatiotemporal dynamics of shaded leaves. Here, we propose to employ an x-ray computed tomography (XCT) imaging technique in combination with an automated, high throughput, image processing pipeline (3DArch, to be developed) to obtain an explicit 3D triangulated surface mesh of a plant canopy (Fig 1). 3DArch would provide a much needed technological jump in 3D LBA modeling where computer generated artificial 3D models suffer from lack of high resolution detailed canopy architecture parameters. The 3D canopy surface mesh output from 3DArch is coupled with an explicit 3D forward ray tracing model (fastTracer) to obtain a high resolution representation of the dynamic spatiotemporal light environment within the canopy over the course of a day (Fig 1). The results from the light attenuation model are then coupled with a process based multi-scale bio-physical canopy and leaf model of photosynthesis (PlantGro) to get fluxes of carbon, water and energy at the leaf and canopy scale, in addition to estimating whole plant carbon uptake, and water loss. Our proposed model framework would be the first of its kind to provide a breakthrough capability in the field of process-based plant modeling that bridges the key technological gap in shaded leaf modeling through explicit representation of 3D canopy architecture. Computer generated artificial 3D model canopy simulations on grasses show that the canopy light distribution does not follow an exponential behavior as assumed in most 1D vegetation models, but rather exhibits a heterogeneous behavior, with a small but significant portion of interior canopy leaves experiencing peak light levels due to sun flecks while the rest of the canopy receives substantially lower light levels. These results indicate a bias in traditional models that overestimates shaded leaf contribution towards whole canopy photosynthesis by up to 15% and underestimate transpiration by up to 12% with an associated increase in sensible heat fluxes since photosynthesis, transpiration and energy partitioning are tightly coupled processes. However, these test results need to be validated using the 3D architecture obtained from a real plant imaged in 3D using 3DArch, rather than an artificial model plant. Our proposed model framework would allow us to identify inefficiencies in canopy light distribution and prescribe optimal canopy architectures (leaf optical properties, leaf area, leaf angles, row spacing, planting density, etc.) to traditional and genetic engineering based crop breeders, for maximizing target functions such as light use efficiency, water use efficiency, drought tolerance, crop yield, crop resilience etc., with important applications in biofuel crops such as sorghum, and switchgrass. Our model can also be applied for designing optimal canopy architecture that increase surface albedo as a mitigation strategy for global warming and projections of crop yield for future food and water security under climate change.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Automatic Image Point Matching

Sparse Image Point Matching (SIPM) is a foundational technology for photo triangulation, structure from motion (SfM), Simultaneous Location and Mapping (SLAM), and data fusion. The goal of the matching is to automatically generate sets of image coordinates that identify the same feature across images. Ideally, the process should be robust to lighting, scale, perspective, and modality changes. The scope of the image matching topic in the field of remote sensing (RS) is enormous because of the variety of collection platforms, modalities, sensor types, applications, and subjects. In this work, we report the history of and assess the state of the art of visible-spectrum (panchromatic and color) image matching of the Earth’s surface. Work specific to large-format images (LFI) (e.g., metric aerial cameras and Earth-observing satellites) will be highlighted. However, the state of the art in this century will mostly be traced through machine vision research and benchmarks because research specific to LFI is rare.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Developing New Beam Constraints to Improve Sterile Neutrino Searches at NOvA

NOvA is a long-baseline neutrino oscillation experiment utilizing two functionally identical liquid scintillator detectors located 14 mrad off-axis from the NuMI beam at Fermilab. The position of the Near Detector (ND), 1 km from the neutrino production target, means it sees a range of off-axis angles (12.6-17.2 mrad), resulting in a noticeable energy gradient across the detector face. In addition, the neutrino production points vary widely along the decay pipe, which leads to a distribution of true neutrino baselines (L), and limits the sensitivity of our current sterile neutrino searches. This work presents a new "baseline reconstruction" technique designed to exploit our off-axis position to triangulate the precise origin of neutrinos within the NuMI beamline. By projecting the total momentum vector of Charged Current (CC) and Neutral Current (NC) interactions back to the central beam axis, we estimate the individual baseline for each event. We describe the technique, evaluate its performance, and demonstrate how this can improve the sensitivity of future searches for sterile neutrinos through reduced flux systematic uncertainties and improved L/E resolution.

Burns, Jessica [Cincinnati U.]↗
DOE OSTI2020

Lattice Gauge Theory and Two Dimensional Quantum Gravity

We construct a two dimensional SU(2) lattice gauge theory containing a single adjoint scalar field which, when tuned to a critical point associated to spontaneous breaking of center symmetry, possesses a long distance effective theory which describes two dimensional quantum gravity. The metric is emergent in the theory and depends on the gauge fields with the diffeomorphism symmetry of the continuum metric theory arising from lattice gauge symmetry. Apart from the intrinsic interest in formulating gravity as a lattice gauge theory the construction holds some advantages over traditional numerical approaches eg. dynamical triangulations since it uses only conventional cubic lattices and Wilson gauge links.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Visualization techniques for the gyrokinetic tokamak simulation code

Gyrokinetic simulations of plasma microturbulence in tokamaks are challenging to visualize because the compute grid follows the magnetic field lines that spiral around the torus. We have overcome this challenge by developing three new approaches that improve visualization of gyrokinetics. Our techniques work directly with the topology of magnetic flux surfaces where the simulation stores variables in concentric rings on poloidal planes (vertical cross sections of the torus). Our visualization preview step triangulates each consecutive pair of rings to display the data on a poloidal plane. The second visualization technique follows spiral field lines around the torus and constructs polygons to visualize a flux surface. Third, the poloidal triangles are connected between planes to form prisms that compose a 3-D model of the entire torus. The visualization workflow produces detailed geometry that matches the high resolution, irregular compute grid for every time step. The surface and solid models are displayed in scientific visualization programs to effectively explore and communicate the results, including fluctuation of electron density, ion temperature, and electrostatic potential. Highly detailed renderings verify plasma behavior along magnetic field lines over time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

MOCMC: Method of Characteristics Moment Closure, a Numerical Method for Covariant Radiation Magnetohydrodynamics

In this work, we present a conservative numerical method for radiation magnetohydrodynamics with frequency-dependent full transport in stationary spacetimes. This method is stable and accurate for both large and small optical depths and radiation pressures. The radiation stress–energy tensor is evolved in flux-conservative form, and closed with a swarm of samples that each transport a multigroup representation of the invariant specific intensity along a null geodesic. In each zone, the enclosed samples are used to efficiently construct a Delaunay triangulation of the unit sphere in the comoving frame, which in turn is used to calculate the Eddington tensor, average source terms, and adaptively refine the sample swarm. Furthermore, radiation four-forces are evaluated in the moment sector in a semi-implicit fashion. The radiative transfer equation is solved in invariant form deterministically for each sample. Since each sample carries a discrete representation of the full spectrum, the cost of evaluating the transport operator is independent of the number of frequency groups, representing a significant reduction of algorithmic complexity for transport in frequency-dependent problems. The major approximation we make in this work is performing scattering in an angle-averaged way. Local adaptivity in samples also makes this scheme more amenable to nonuniform meshes than a traditional Monte Carlo method. We describe the method and present results on a suite of test problems. We find that Method of Characteristics Moment Closure converges at least as ~N -1 , rather than the canonical Monte Carlo N -1/2 , where N is the number of samples per zone.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Multi-sensor Observations of a High-velocity Fireball over the South Atlantic on 2026 April 1

A high-velocity fireball was detected over the South Atlantic (41.9°S, 54.7°W) on 2026 April 1 at 02:13:14 UTC by U.S. Government (USG) sensors, with peak brightness at 90.5 km altitude. The event was well-observed from geostationary orbit by two civilian lightning imagers with near-orthogonal geometry, the Geostationary Operational Environmental Satellite-East Geostationary Lightning Mapper (GLM) and the Exploitation of Meteorological Satellites Meteosat Third Generation Imager 1 Lightning Imager, and it produced low-frequency acoustic signatures. The GLM measured a total radiated energy of 2.4 × 10 10 J, corresponding to a calculated impact energy of 0.086 kt TNT equivalent. Stereoscopic triangulation of the imager tracks yields an independent pre-atmospheric velocity of ~57 km s −1 , some 18% below the USG-reported value. Because orbital provenance is acutely sensitive to the entry velocity, whose reported uncertainty could be substantial, this discrepancy is notable. We document the multi-sensor record and identify the analysis required to assess provenance, which a subsequent study will present.

Silber, Elizabeth Allaryce [Sandia National Lab. (↗
DOE OSTI2021

Geographic Mapping Tool for Power Systems

Power systems models are continuously updated to reflect changing infrastructure and modeling techniques. Open source geographic locations of power systems infrastructure do not include attributes that directly tie to electrical models. Moreover, relationships established between electric models and geospatial data sets do not apply to a future model because of nomenclature changes and architecture changes. The software titled “Geographic Mapping Tool for Power Systems” starts with a known set of locations, then runs multiple iterations to identify missing locations with a triangulation algorithm based on the branch connections and distances to the known locations.”

Snyder, Isabelle B↗
DOE OSTI2021

System and method of mapping topology

A system for determining a fingerprint of a structure is provided. The system includes a plurality of granules inserted in a structure having a plurality of fissures, fractures, and cracks (collectively apertures), each granule comprising a membrane, and at least one bubble of compressed gas formed in the membrane, the membrane selectively dissolving in presence of a predetermined fluid and thereby selectively bursting the at least one bubble, thereby generating a concussing vibration, at least at least i) three detection devices for two-dimensional mapping or ii) four detection devices for three-dimensional mapping placed proximate to the structure according to a predetermined placement schedule, and a computing device comprising a processor configured to receive data from the at least three or four detection devices and to determine location of the at least one bubble of each of the plurality of the granules at the time of bursting by triangulating the concussive vibration in order to determine location of the at least one bubble.

Pyrak-Nolte, Laura↗