Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FRACTALS”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Imaging antiferromagnetic domain fluctuations and the effect of atomic scale disorder in a doped spin-orbit Mott insulator

Correlated oxides can exhibit complex magnetic patterns. Understanding how magnetic domains form in the presence of disorder and their robustness to temperature variations has been of particular interest, but atomic scale insight has been limited. We use spin-polarized scanning tunneling microscopy to image the evolution of spin-resolved modulations originating from antiferromagnetic (AF) ordering in a spin-orbit Mott insulator perovskite iridate Sr 3 Ir 2 O 7 as a function of chemical composition and temperature. We find that replacing only several percent of lanthanum for strontium leaves behind nanometer-scale AF puddles clustering away from lanthanum substitutions preferentially located in the middle strontium oxide layer. Thermal erasure and reentry into the low-temperature ground state leads to a spatial reorganization of the AF puddles, which nevertheless maintain scale-invariant fractal geometry in each configuration. Our experiments reveal multiple stable AF configurations at low temperature and shed light onto spatial fluctuations of the AF order around atomic scale disorder in electron-doped Sr 3 Ir 2 O 7 .

36 MATERIALS SCIENCE↗
DOE OSTI2025

Diffusion power spectra as a window into dynamic materials architecture

Chemical recycling of commodity and specialty polymers presents a multifaceted challenge for industrial societies. On one hand, macromolecular architectures must be engineered to yield durable products that, on the other hand, rapidly deconstruct to recyclable monomers under pre-determined conditions. Polymer deconstruction is a chemical process that requires deep understanding of molecular reactivity in heterogeneous media, where porous material architectures evolve in both space and time. To build this understanding, we develop herein experimental and analytical methods describing sets of diffusive eigenmodes that exist within time-varying, non-Euclidean boundary conditions, a situation commonly encountered in the reactive deconstruction of polymers where chain fragments splay, alter their local dynamics, and evolve in their confinement of reacting media. Diffusion power spectra, discerned experimentally by NMR, yield polymer and solvent frequency-domain velocity autocorrelation functions that are analyzed in the context of physical models for chemical reactions parameterized with fractal mathematics. The results connect local motion in polymers to chemical reactivity during acidolysis of circular elastomers.

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

Mapping structures and dynamics with frequency-correlated diffusion exchange

Understanding molecular motion in diffusion-driven complex environments is critical for designing sustainable materials and improving chemical processes. Here, we introduce a multidimensional nuclear magnetic resonance (NMR) method that captures how molecular populations exchange across different dynamic regimes. By extending the modulated gradient spin-echo technique to include frequency-frequency correlations, our approach reveals diffusion pathways that are otherwise obscured in heterogeneous systems. Implemented on a unilateral NMR magnet, the method eliminates gradient pulsing constraints and accesses dynamics in the kilohertz regime. We apply this technique to swelling and acid-catalyzed deconstruction of cross-linked and linear polymers to observe how structural heterogeneity evolves over time. By linking molecular motion to topology and chemical state, we extract physical metrics such as fractal surface dimensionality and reaction wavefront velocity, properties inaccessible with standard diffusion measurements. This work expands the capabilities of NMR for probing soft matter, with implications for polymer recycling and materials design.

Fricke, Sophia N. [University of California, Berke↗
DOE OSTI2023

Topological aspects of matters and Langlands program

The Langlands program is a vast mathematical projection linking number theory and geometry. In high-energy physics, a connection with mirror symmetry has been suggested in string theory, but it has been little studied in low-energy physics. In the framework of the Langlands program, we present a unified description of the integer and fractional quantum Hall effect and the duality found in the fractal nature of the energy spectrum of two-dimensional Bloch electrons, statistical physics, and quantum computation. The new unified view of existing dualism presented in this paper raises the entirely new question of how each theory of physics is connected as a piece of the Langlands program.

Physics↗
DOE OSTI2020

Multiple Environmental Influences on the Lightning of Cold-Based Continental Cumulonimbus Clouds. Part I: Description and Validation of Model

In this two-part paper, influences from environmental factors on lightning in a convective storm are assessed with a model. In Part I, an electrical component is described and applied in the Aerosol–Cloud model (AC). AC treats many types of secondary (e.g., breakup in ice–ice collisions, raindrop-freezing fragmentation, rime splintering) and primary (heterogeneous, homogeneous freezing) ice initiation. AC represents lightning flashes with a statistical treatment of branching from a fractal law constrained by video imagery. The storm simulated is from the Severe Thunderstorm Electrification and Precipitation Study (STEPS; 19/20 June 2000). The simulation was validated microphysically [e.g., ice/droplet concentrations and mean sizes, liquid water content (LWC), reflectivity, surface precipitation] and dynamically (e.g., ascent) in our 2017 paper. Predicted ice concentrations (~10 L -1 ) agreed—to within a factor of about 2—with aircraft data at flight levels (-10° to -15°C). Here, electrical statistics of the same simulation are compared with observations. Flash rates (to within a factor of 2), triggering altitudes and polarity of flashes, and electric fields, all agree with the coincident STEPS observations. The “normal” tripole of charge structure observed during an electrical balloon sounding is reproduced by AC. It is related to reversal of polarity of noninductive charging in ice–ice collisions seen in laboratory experiments when temperature or LWC are varied. Positively charged graupel and negatively charged snow at most midlevels, charged away from the fastest updrafts, is predicted to cause the normal tripole. Total charge separated in the simulated storm is dominated by collisions involving secondary ice from fragmentation in graupel–snow collisions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Chemical and architectural intricacy from nanoscale tetrahedra and their analogues

The tetrahedron, as the simplest platonic shape, is a profound building block with the potential to create intricate superstructures. Noteworthy designs utilizing tetrahedral building blocks include the Sierpiński tetrahedron (the most fundamental three-dimensional fractal), a one-dimensional helical structure known as the tetrahelix, and various crystalline and quasicrystalline packings. Historically, the practicality of tetrahedral superstructures has been evident, providing stable, well-defined frameworks for various constructions, including truss bridges, tower cranes, and electricity transmission line pylons. In the field of self-assembled nanocrystal superlattices, tetrahedral nanocrystals, as building blocks, occupy a unique place among all the possible nanoscale particles. Mathematical models, simulation work, and experimental studies using nanocrystals in the laboratory have suggested that self-assembled structures derived from nanoscale tetrahedral building blocks are notably intricate, giving rise to new horizons of high-entropy nanocrystal superlattices. An important implication from previous works is that such tetrahedral nanocrystal superlattices form through highly delicate interparticle interactions, emphasizing the importance of the fine features of these nanocrystals. In this article, we summarize the advances in superlattices assembled from tetrahedral nanocrystals. Here, we first define the tetrahedron and tetrahedron analogues based on Conway’s transformation and graph theory, underscoring their relevance to the crystallization process producing tetrahedral nanocrystals. Then, we showcase previous reports on the synthesis of tetrahedral nanocrystals and the resulting nanocrystal superstructures. Finally, we conclude by offering insights and perspective into the chemical and architectural intricacy that could emerge from tetrahedral nanocrystals.

Nanostructure↗
DOE OSTI2020

Cellular Response to Proton Irradiation: A Simulation Study with TOPAS-nBio

The cellular response to ionizing radiation continues to be of significant research interest in cancer radiotherapy, and DNA is recognized as the critical target for most of the biologic effects of radiation. Incident particles can cause initial DNA damages through physical and chemical interactions within a short time scale. Initial DNA damages can undergo repair via different pathways available at different stages of the cell cycle. The misrepair of DNA damage results in genomic rearrangement and causes mutations and chromosome aberrations, which are drivers of cell death. This work presents an integrated study of simulating cell response after proton irradiation with energies of 0.5–500 MeV (LET of 60–0.2 keV/µm). A model of a whole nucleus with fractal DNA geometry was implemented in TOPAS-nBio for initial DNA damage simulations. The default physics and chemistry models in TOPAS-nBio were used to describe interactions of primary particles, secondary particles, and radiolysis products within the nucleus. The initial DNA double-strand break (DSB) yield was found to increase from 6.5 DSB/Gy/Gbp at low-linear energy transfer (LET) of 0.2 keV/µm to 21.2 DSB/Gy/Gbp at high LET of 60 keV/µm. A mechanistic repair model was applied to predict the characteristics of DNA damage repair and dose response of chromosome aberrations. It was found that more than 95% of the DSBs are repaired within the first 24 h and the misrepaired DSB fraction increases rapidly with LET and reaches 15.8% at 60 keV/µm with an estimated chromosome aberration detection threshold of 3 Mbp. The dicentric and acentric fragment yields and the dose response of micronuclei formation after proton irradiation were calculated and compared with experimental results.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2021

Diving into a holographic superconductor

Charged black holes in anti-de Sitter space become unstable to forming charged scalar hair at low temperatures T < T_\text{c} T < T c . This phenomenon is a holographic realization of superconductivity. We look inside the horizon of these holographic superconductors and find intricate dynamical behavior. The spacetime ends at a spacelike Kasner singularity, and there is no Cauchy horizon. Before reaching the singularity, there are several intermediate regimes which we study both analytically and numerically. These include strong Josephson oscillations in the condensate and possible `Kasner inversions’ in which after many e-folds of expansion, the Einstein-Rosen bridge contracts towards the singularity. Due to the Josephson oscillations, the number of Kasner inversions depends very sensitively on T T , and diverges at a discrete set of temperatures \{T_n\} { T n } that accumulate at T_c T c . Near these T_n T n , the final Kasner exponent exhibits fractal-like behavior.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Evaluating one-loop string amplitudes

We evaluate one-loop open-string amplitudes at finite \alpha' α ′ for the first time. Our method involves a deformation of the integration contour over the modular parameter \tau τ to a fractal contour introduced by Rademacher in the context of analytic number theory. This procedure leads to explicit and practical formulas for the one-loop four-point amplitudes in type-I superstring theory, amenable to numerical evaluation. We plot the amplitudes as a function of the Mandelstam invariants s s and t t and directly verify long-standing conjectures about their behaviour at high energies.

Eberhardt, Lorenz↗
DOE OSTI2023

Assessing the link between aerosol mixing state, structure and composition and their optical properties: Ascension Island as a testbed for the South‐East Atlantic aerosol regimes

Forest and savanna fires are a major source of pollution over the globe. Such fire events emit to the atmosphere a considerable amount of “fluffy” fractal black looking particles that are often referred to as black carbon (BC) or soot aerosols. In fact, BC particles (aerosols) emitted from savanna and agricultural fires over the African continent play an important role in the region’s and the global climate. These fires occur each year between July and October, which is referred to as the biomass burning (BB) season and are considered among the globe's largest man-made emission sources of BC particles. However, during these fire events, not only BC particles are emitted, but also other types of particles and gases that can interact with each other to create new particles, generating a slew of particles often generally referred to as BB aerosols. Although the majority of the fire emissions occur inland in the sub-Saharan part of Africa, their effect reaches far beyond the continent due to the prevailing easterly winds, which transport those aerosols over the South-East Atlantic (SEA) ocean, off the west-coast of Africa. The reason why these particles are so important to the regional and global climate is that they absorb sunlight radiation, which overall warms the climate. However, the magnitude of this warming effect is highly uncertain and depends upon a multitude of aspects. Our overarching goal for this project was to better understand the changes in the aerosol properties as they transport from fire sources downstream toward the ocean and how this can eventually affect the earth's radiative budget in terms of warming or cooling. Following the above, our work here was focused on the investigation of the African BB aerosol plumes mixing state (mixing of BB aerosols and other aerosol types), composition and size distribution as they transport from the African continent towards the SEA ocean, and their link with some of the commonly measurable bulk optical properties such as mass absorption coefficient (MAC) and single scattering albedo (SSA). We utilized three field missions that were conducted over the SEA ocean between 2016 and 2018. Specifically, we used (1) ground-based measurements from the DOE ARM Mobile Facility (AMF1), which was deployed at Ascension Island (ASI) for the LASIC (Layered Atlantic Smoke Interactions with Clouds) campaign between June 2016 and October 2017, (2) airborne measurements from the UK CLARIFY (Clouds and Aerosol Radiative impacts and Forcing) campaign during Aug-Sep-2017, and (3) airborne measurements from the ORACLES (Observations of Aerosols above Clouds and their interactions) campaign during Aug-2017, and Sep-Oct-2018. The three campaigns cover a relatively large region, from the western African coast on the east towards Ascension Island in the middle of the SEA to the west. We used the data gathered by these campaigns to find how BB properties and composition change from near emission sources (the ORACLES campaign flights) downstream (the CLARIFY flights and the ground based LASIC measurements). We used particle trajectory following methodologies to connect the measurements from the different campaigns, which allowed us to follow a certain aerosol plume (airmass) from near-source towards Ascension Island. We found out that atmospheric aging processes governed by the UV light from the sun (photochemistry) and cloud processes (aqueous chemistry) govern the changes seen in the particle composition during the transport were more dominant in dictating aerosol composition than specific emission composition or source types. We assessed how aerosol optical properties change during the BB season and what are the main drivers of this change. We found that the enhanced ratio of BC to CO is well correlated with single scattering albedo (SSA) and mass absorption coefficient (MACBC), providing a simple way to estimate the aerosol optical characteristics in the south-eastern Atlantic Ocean. From the analysis of the location of BB, the primary source fuel, the water content in the fuel, combined with the mean cloud cover and precipitation in the transport areas of the BB plume, we conclude that the increase in BC/CO from June to August is likely to be caused by burning becoming more flaming (hot fires), and the decrease in BC/CO in September and October may be caused by smoldering (colder) fires. We found that aerosol hygroscopicity increases with decreasing altitude below 2km, and that enhanced BB hygroscopicity at lower altitudes is mainly due to a lower organic aerosol (OA) fraction, increased sulfate fraction, and greater hygroscopicity parameter of OA at lower altitudes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Thermally Insulating Transparent Barrier (THINNER) coatings on single pane windows

Conventional silica aerogel monoliths can provide remarkable thermal insulation but the presence of large pores (> 30 nm) tends to scatter visible light and render the material opaque or translucent instead of transparent. In addition, they are prone to cracking during synthesis and handling which makes them difficult to integrate in products and in particular in window solutions. This project developed two new solgel synthesis methods using silica precursors or preformed silica nanoparticles and ambient drying to produce mesoporous organo-silica monoliths. The monoliths were (i) thermally insulating, (ii) optically transparent, (iii) flexible, and (iv) hydrophobic. They feature porosity ranging from 50% to 90% with narrow pore size distribution with pore less than 20 nm resulting in excellent optical clarity (haze < 2%) and very low thermal conductivity (< 30 mW/mK). Interestingly, not only porosity but also pore size and mass fractal dimension were found to affect the thermal conductivity of the mesoporous silica. The superior transparency of the monoliths was shown to be attributed to dependent scattering among silica nanoparticles. Flexibility was achieved through trimethylchlorosilane surface modification. Furthermore, process scale-up and integration of the ambigel monoliths into window solutions using optically clear adhesives were also demonstrated in 6”x6” double-pane windows. The aerogel and the window solution were shown to be durable to accelerated aging under UV, moisture, and/or temperature gradient and thermal cycling. Overall, the window solution achieved the technical performance and the cost target of $10/sqft set for the SHIELD program. However, scaling to industry relevant scale (> 10’x10’) remains a challenge due to requirements on the drying process and fume hood size and to the propensity of large aerogel slab to crack during drying.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

NLR CSP Optical Facilities: Illuminating the Path Forward Through Innovation and Impact: Agreement 38490

This initiative is a multi-faceted project at the National Laboratory of the Rockies (NLR) aimed at strengthening its Concentrating Solar Power (CSP) Optical Facilities to advance the development of low-cost, high-performance materials for solar and other applications. The project's strategy is built on three pillars: strategic stakeholder engagement, diligent facility maintenance and utilization, and the development of new research capabilities. The overarching goal is to ensure the facilities remain state-of-the-art resources for industry and academia, thereby accelerating the conversion of concentrated sunlight into energy. A key driver of the project is an international Advisory Board, which provides critical guidance on research priorities and industry needs, leading to new collaborations and secured funding. This external engagement, combined with proactive outreach to industry partners, ensures the lab's work remains aligned with real-world challenges, including materials durability and performance certification. Significant efforts in facility maintenance have addressed challenges with aging infrastructure. Notable achievements include the complete refurbishment of the hail-damaged Ultra-Accelerated Weathering System (UAWS) and the successful replacement of a failing 15-year-old Lambda 1050 spectrophotometer with a new-generation model, substantially upgrading material characterization capabilities. These maintenance activities were complemented by achieving a prestigious ISO 9001:2015 certification for the Advanced Optical Materials Labs, formally recognizing the quality and reliability of NLR's measurement capabilities. Despite these successes, challenges remain, including high demand for the High Flux Solar Furnace (HFSF) and intermittent failures of other key instruments. The project has delivered major advancements in research techniques and capabilities. At the Flatirons campus, a new indoor laboratory, was established to house advanced deflectometry and photogrammetry systems for heliostat characterization. For on-sun testing, a novel, actively cooled turning mirror was developed for the HFSF, enabling more realistic testing of particle receivers and components. A collaboration with Virginia Tech successfully demonstrated the high-temperature durability of a new solar absorber coating through extensive cyclic testing. Concurrently, new modeling took place to better predict material degradation on rough, fractal surfaces. In summary, this project has systematically enhanced NLR's CSP Optical Facilities through strategic upgrades, rigorous maintenance, and stakeholder-guided research. By overcoming equipment failures, budgetary constraints, and logistical hurdles, the project has reinforced NLR's role as a central hub for CSP innovation and materials testing. Future work will focus on securing diverse funding, expanding collaborations, and continuing to provide the critical infrastructure needed to accelerate the development and deployment of next-generation technologies.

14 SOLAR ENERGY↗
DOE OSTI2026

Analytical Dynamic Aperture Calculation Compared to Tracking Results

In studying dynamic aperture in the storage rings, it is becoming increasing beneficial to understand the growth of the fractal boundary of the dynamic aperture without tracking for millions of turns in a hadron ring. This note describes the J. Gao theory applied to the Electron Ion Collider (EIC) Hadron Storage Ring (HSR) and compares simulation results to theoretical predictions to show that the analytical method may prove useful in quick determinations of improved dynamic aperture based on sextupole and octupole configurations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Large-scale homogeneity in the distribution of quasars in the Hercules-Corona Borealis Great Wall region

In light of recent debates on the existence of a gigaparsec-scale structure traced by gamma-ray bursts, namely the Hercules-Corona Borealis Great Wall (HCBGW), we revisit large-scale homogeneity in the spatial distribution of quasars. Our volume-limited sample of quasars in the redshift range 1:6 < z ≤ 2:1, which is constructed from the data release 7 of the Sloan Digital Sky Survey quasar catalogue, covers about half of the suspected HCBGW region. We analyze the sample in two complementary ways: fractal analysis of determining the average scale of homogeneity and friends- of-friends analysis of identifying specific large-scale structures. The quasar distribution on average reaches homogeneity at rh = 136 ± 38h-1 Mpc and the richness and comoving size frequencies of large (> ~ 150h-1 Mpc) quasar groups are consistent with the prediction of homogeneous distribution. These results put constraints on the spatial extent of the HCBGW but do not contradict its existence since our quasar sample does not cover the entire HCBGW region.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

DEEP CELLULAR RECURRENT NEURAL ARCHITECTURE FOR EFFICIENT MULTIDIMENSIONAL TIME-SERIES DATA PROCESSING

Efficient processing of time series data is a fundamental yet challenging problem in pattern recognition. Though recent developments in machine learning and deep learning have enabled remarkable improvements in processing large scale datasets in many application domains, most are designed and regulated to handle inputs that are static in time. Many real-world data, such as in biomedical, surveillance and security, financial, manufacturing and engineering applications, are rarely static in time, and demand models able to recognize patterns in both space and time. Current machine learning (ML) and deep learning (DL) models adapted for time series processing tend to grow in complexity and size to accommodate the additional dimensionality of time. Specifically, the biologically inspired learning based models known as artificial neural networks that have shown extraordinary success in pattern recognition, tend to grow prohibitively large and cumbersome in the presence of large scale multi-dimensional time series biomedical data such as EEG. Consequently, this work aims to develop representative ML and DL models for robust and efficient large scale time series processing. First, we design a novel ML pipeline with efficient feature engineering to process a large scale multi-channel scalp EEG dataset for automated detection of epileptic seizures. With the use of a sophisticated yet computationally efficient time-frequency analysis technique known as harmonic wavelet packet transform and an efficient self-similarity computation based on fractal dimension, we achieve state-of-the-art performance for automated seizure detection in EEG data. Subsequently, we investigate the development of a novel efficient deep recurrent learning model for large scale time series processing. For this, we first study the functionality and training of a biologically inspired neural network architecture known as cellular simultaneous recurrent neural network (CSRN). We obtain a generalization of this network for multiple topological image processing tasks and investigate the learning efficacy of the complex cellular architecture using several state-of-the?art training methods. Finally, we develop a novel deep cellular recurrent neural network (CDRNN) architecture based on the biologically inspired distributed processing used in CSRN for processing time series data. The proposed DCRNN leverages the cellular recurrent architecture to promote extensive weight sharing and efficient, individualized, synchronous processing of multi-source time series data. Experiments on a large scale multi-channel scalp EEG, and a machine fault detection dataset show that the proposed DCRNN offers state-of-the-art recognition performance while using substantially fewer trainable recurrent units.

Vidyaratne, Lasitha S.↗
NASA NTRS1984

Nonlinear problems in flight dynamics

A comprehensive framework is proposed for the description and analysis of nonlinear problems in flight dynamics. Emphasis is placed on the aerodynamic component as the major source of nonlinearities in the flight dynamic system. Four aerodynamic flows are examined to illustrate the richness and regularity of the flow structures and the nature of the flow structures and the nature of the resulting nonlinear aerodynamic forces and moments. A framework to facilitate the study of the aerodynamic system is proposed having parallel observational and mathematical components. The observational component, structure is described in the language of topology. Changes in flow structure are described via bifurcation theory. Chaos or turbulence is related to the analogous chaotic behavior of nonlinear dynamical systems characterized by the existence of strange attractors having fractal dimensionality. Scales of the flow are considered in the light of ideas from group theory. Several one and two degree of freedom dynamical systems with various mathematical models of the nonlinear aerodynamic forces and moments are examined to illustrate the resulting types of dynamical behavior. The mathematical ideas that proved useful in the description of fluid flows are shown to be similarly useful in the description of flight dynamic behavior.

Chapman, G. T.↗
NASA NTRS1984

The distribution of catchment coverage by stationary rainstorms

The occurrence of wetted rainstorm area within a catchment is modeled as a Poisson arrival process in which each storm is composed of stationary, nonoverlapping, independent random cell clusters whose centers are Poisson-distributed in space and whose areas are fractals. The two Poisson parameters and hence the first two moments of the wetted fraction are derived in terms of catchment average characteristics of the (observable) station precipitation. The model is used to estimate spatial properties of tropical air mass thunderstorms on six tropical catchments in the Sudan.

Eagleson, P. S.↗
NASA NTRS1985

Observations, theoretical ideas and modeling of turbulent flows: Past, present and future

Turbulence was analyzed in a historical context featuring the interactions between observations, theoretical ideas, and modeling within three successive movements. These are identified as predominantly statistical, structural and deterministic. The statistical movement is criticized for its failure to deal with the structural elements observed in turbulent flows. The structural movement is criticized for its failure to embody observed structural elements within a formal theory. The deterministic movement is described as having the potential of overcoming these deficiencies by allowing structural elements to exhibit chaotic behavior that is nevertheless embodied within a theory. Four major ideas of this movement are described: bifurcation theory, strange attractors, fractals, and the renormalization group. A framework for the future study of turbulent flows is proposed, based on the premises of the deterministic movement.

Chapman, G. T.↗