Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large deviations”

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

Magnetic field-induced non-trivial electronic topology in Fe3− x GeTe2

The anomalous Hall, Nernst, and thermal Hall coefficients of the itinerant ferromagnet Fe3−xGeTe2 display anomalies upon cooling that are consistent with a topological transition that could induce deviations with respect to the Wiedemann–Franz (WF) law. This law has not yet been validated for the anomalous transport variables, with recent experimental studies yielding material-dependent results. Nevertheless, the anomalous Hall and thermal Hall coefficients of Fe3−xGeTe2 are found, within our experimental accuracy, to satisfy the WF law for magnetic fields μ0H applied along its c axis. Remarkably, large anomalous transport is also observed for μ0H||a axis with the field aligned along the gradient of the chemical potential generated by thermal gradients or electrical currents, a configuration that should not lead to their observation. These anomalous planar quantities are found to not scale with the component of the planar magnetization (M||), showing instead a sharp decrease beyond μ0H||= 4 T or the field required to align the magnetic moments along μ0H||. We argue that chiral spin structures associated with Bloch domain walls lead to a field-dependent spin chirality that produces a novel type of topological transport in the absence of interaction between the magnetic field and electrical or thermal currents. Locally chiral spin structures are captured by our Monte Carlo simulations incorporating small Dzyaloshinskii–Moriya and biquadratic exchange interactions. These observations reveal not only a new way to detect and expose topological excitations, but also a new configuration for heat conversion that expands the current technological horizon for thermoelectric energy applications.

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

Instabilities of a vortex-ring-bright soliton in trapped binary three-dimensional Bose-Einstein condensates

Instabilities of vortex-ring-bright coherent structures in harmonically trapped two-component three-dimensional Bose-Einstein condensates are studied numerically within the coupled Gross-Pitaevskii equations and interpreted analytically. Interestingly, the filled vortex core with a sufficiently large amount of the bright component is observed to reduce the parametric interval of stability of the vortex ring. Here, we have identified the mechanisms of several linear instabilities and one nonlinear parametric instability in this connection. Two of the linear instabilities are qualitatively different from ones reported earlier and are associated with azimuthal modes of m = 0 and m = 1, i.e., deviations of the vortex from the stationary ring shape. Our nonlinear parametric resonance instability occurs between the m = 0 and m = 2 modes and signals the exchange of energy between them.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Performance Characterization of FB-Line Neutron Multiplicity Counter and Large Neutron Multiplicity Counter

Savanah River National Laboratory’s (SRNL) Nuclear Measurements group was tasked with characterizing the performance of two neutron multiplicity counters located at SRNL. Characterization measurements were made to determine the gate width, pre-delay, deadtime parameters, triples and doubles gate fractions, detector efficiency, and operating high voltage for the Large Neutron Multiplicity Counter (LNMC) and the FB Line Neutron Multiplicity Counter (FBLNMC). The parameters were determined, shown below, and were, as to be expected, slightly different than the previous calibrations, which were performed over 20 years ago. Several Pu samples were measured to validate the characterizations of the FBLNMC and LNMC. The measurements determined the sample Pu-240 mass within <2% deviation for the pure plutonium samples and ~8% for the mixed oxide sample. The pure Pu samples had significantly better accuracy compared with the impure mixed oxide sample due to the lack of induced fission or alpha,n neutrons from impurities. Overall, the characterization of the neutron multiplicity counters, and the determination of their operability has been completed successfully.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1975

Abundances, charge states, and energy spectra of helium and heavy ions during solar particle events

Carbon and oxygen energy spectra observed during several solar events indicate a systematic deviation of these spectra from a simple power law: the spectra bend down below about 100 keV/nuc and the amount of this bending is highly correlated with the size of the flare, as measured by the 'event averaged' flux of 130-220 keV protons. The energy spectra of helium computed for the same time periods do not show a similar feature. A large variability of the alpha/CNO ratio is found from event to event (from 2 to about 20 at 40 keV/nuc), and in all cases examined the carbon and oxygen nuclei are nearly fully stripped. These results are interpreted as evidence for storage of energetic ions in hot coronal regions, followed by strong adiabatic deceleration.

Gloeckler, G.↗
NASA NTRS1994

R(sub nu)-dependent optical and near-ultraviolet extinction

We have derived extinctions A(lambda)/A(V) at the wavelengths of the uvby filters for 22 stars, with a range of values of R(sub nu), from the sample of Cardelli, Clayton, & Mathis (1989, hereafter CCM). We have fit these extinctions, and also UBVRIJHKL, IUE and ANS extinction measurements, with linear relations A(lambda)/A/(V) = a+b/R(sub nu) and fit a and b as a function of x(=1/lambda) with polynomials to obtain an R(sub nu)-dependent mean extinction law (A(x)/A(V) = a(x) + b(x)/R(sub nu))in the optical and near-ultraviolet (1.1/micrometer less than or equal to 3.3/micrometer). This law is virtually identical to the CCM extinction law for large values of R(sub nu)(R(sub nu) approximately 5) but is slightly lower in the near-ultraviolet for smaller R(sub nu) (R(sub nu) approximately 3). The extinction law presented here agrees much better with a high-resolution extinction curve for the diffuse interstellar medium (R(sub nu) approximately 3.1), presented by Bastiaansen (1992), than CCM. The deviations of individual extinction curves from the mean are dominated by observational errors. The wavelength resolution of this work is not high enough to show evidence for or against the existence of very broad structure in optical extinction curves.

O'Donnell, James E.↗
NASA NTRS2014

An Architecture to Enable Autonomous Control of Spacecraft

Autonomy is required for manned spacecraft missions distant enough that light-time communication delays make ground-based mission control infeasible. Presently, ground controllers develop a complete schedule of power modes for all spacecraft components based on a large number of factors. The proposed architecture is an early attempt to formalize and automate this process using on-vehicle computation resources. In order to demonstrate this architecture, an autonomous electrical power system controller and vehicle Mission Manager are constructed. These two components are designed to work together in order to plan upcoming load use as well as respond to unanticipated deviations from the plan. The communication protocol was developed using "paper" simulations prior to formally encoding the messages and developing software to implement the required functionality. These software routines exchange data via TCP/IP sockets with the Mission Manager operating at NASA Ames Research Center and the autonomous power controller running at NASA Glenn Research Center. The interconnected systems are tested and shown to be effective at planning the operation of a simulated quasi-steady state spacecraft power system and responding to unexpected disturbances.

space power systems↗
DOE OSTI2022

Numerical study of multi-component flow and mixing in a scaled fission product venting system

Numerical simulation of isothermal multi-component gas flow was conducted to study fluid flow and mixing in a scaled low-pressure, low-temperature test facility for the fission product venting system (FPVS) of a gas cooled fast reactor (GFR). The geometry of FPVS test facility was an open loop including gradual expansion coupling and two 90°pipe elbows. First, the validation of large eddy simulation (LES) wall-adapting local eddy-viscosity (WALE) turbulent model was performed for transitional flow in a circular pipe with gradual expansion. The reactingFOAM solver in OpenFOAM v8 was employed. By introducing appropriate turbulence disturbance at the flow inlet, such as turbulence intensity, the numerical results showed good agreement in (1) the velocity profile downstream of the pipe expansion measured by particle image velocimetry (PIV) in this study, and (2) the reattachment length reported in the literature. Using this validated model, the numerical simulation of flow and component mixing was performed for the FPVS with an inlet flowrate corresponding to a Reynolds number of 2,400 to investigate the flow behavior and component mixing. Standard deviation of the mass fraction of component, referred to as absolute mixing index (AMI), was calculated to quantify the mixing of components, showing that the component mixing length determined by AMI is related to the reattachment length downstream of the expansion. Additionally, this study was able to identify the best location in the FPVS test facility to measure the velocity and concentration profies where multicomponent flow is well mixed.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Control over Conformational Landscapes of Polypeptoids by Monomer Sequence Patterning

The ability to program chain conformation and structure through control over the monomer sequence of synthetic polymers has broad implications for next-generation material design. While related problems of protein-folding and de novo design have generated accurate predictions of 3D folded chain structures, generalization to synthetic polymers remains intractable due to the requirement of large structural databases and the intrinsically disordered nature of polymer building blocks. In this work, polypeptoids, a class of peptidomimetic synthetic polymers, are utilized to build a general workflow for the study of relationships between monomer sequence and dynamic 3D chain structure in solution. Furthermore, this work demonstrates how control over the monomer sequence can alter the conformational landscape of synthetic polymers to deviate dramatically from classical chain statistics. Specifically, the distribution of end-to-end distances, as measured by double electron-electron resonance spectroscopy in dilute solvent, is systematically skewed towards shorter distances with an increasing number of hydrophobes and further refined by hydrophobe arrangement in amphiphilic polypeptoid chains.

36 MATERIALS SCIENCE↗
DOE OSTI2022

‘Flux+Mutability’: a conditional generative approach to one-class classification and anomaly detection

Abstract Anomaly Detection is becoming increasingly popular within the experimental physics community. At experiments such as the Large Hadron Collider, anomaly detection is growing in interest for finding new physics beyond the Standard Model. This paper details the implementation of a novel Machine Learning architecture, called Flux+Mutability, which combines cutting-edge conditional generative models with clustering algorithms. In the ‘flux’ stage we learn the distribution of a reference class. The ‘mutability’ stage at inference addresses if data significantly deviates from the reference class. We demonstrate the validity of our approach and its connection to multiple problems spanning from one-class classification to anomaly detection. In particular, we apply our method to the isolation of neutral showers in an electromagnetic calorimeter and show its performance in detecting anomalous dijets events from standard QCD background. This approach limits assumptions on the reference sample and remains agnostic to the complementary class of objects of a given problem. We describe the possibility of dynamically generating a reference population and defining selection criteria via quantile cuts. Remarkably this flexible architecture can be deployed for a wide range of problems, and applications like multi-class classification or data quality control are left for further exploration.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Beehive haloscope for high-mass axion dark matter

We propose a new haloscope geometry that can arbitrarily increase the resonator volume for a given target axion mass. This geometry consists of closely packed, overlapping coaxial cavities operating as a single resonator. While the resonant frequency is still determined by the dimensions of the individual “cells,” the strong interactions between the cells encourage the entire “beehive” to oscillate in phase, a phenomenon expected of tightly coupled harmonic oscillators. This synchronization behavior allows the construction of a singly connected large volume resonator at high frequency by simply increasing the number of the cells. Using direct numerical simulations, we verify the existence of a global eigenmode that has a high (40%) form factor in a 169-element beehive resonator. The resonant frequency of the eigenmode is tunable by moving the center rods laterally in unison. The form factor is very tolerant to dimensional deviations and misalignment as a result of mode hybridization due to strong coupling. The beehive haloscope inherits many appealing properties from the conventional coaxial cavity: a high quality factor, compatibility with a solenoid magnet, and ease of fabrication, tuning, and coupling. We argue that this geometry is an excellent candidate for high-mass axion searches covering the post-inflationary parameter space (>5 GHz).

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Interacting dark sector within ETHOS: Cosmological constraints from SPT cluster abundance with DES and HST weak lensing data

We use galaxy cluster abundance measurements from the South Pole Telescope enhanced by multicomponent matched filter confirmation and complemented with mass information obtained using weak-lensing data from Dark Energy Survey Year 3 (DES Y3) and targeted Hubble Space Telescope observations for probing deviations from the cold dark matter paradigm. Concretely, we consider a class of dark sector models featuring interactions between dark matter (DM) and a dark radiation (DR) component within the framework of the effective theory of structure formation (ETHOS). We focus on scenarios that lead to power suppression over a wide range of scales, and thus can be tested with data sensitive to large scales, as realized, for example, for DM–DR interactions following from an unbroken non-Abelian 𝑆⁢𝑈⁡(𝑁) gauge theory (interaction rate with power-law index 𝑛 = 0 within the ETHOS parametrization). Cluster abundance measurements are mostly sensitive to the amount of DR interacting with DM, parametrized by the ratio of DR temperature to the cosmic microwave background (CMB) temperature, 𝜉 DR = 𝑇 DR /𝑇 CMB . We find an upper limit 𝜉 DR < 17% at 95% credibility. When the cluster data are combined with Planck 2018 CMB data along with baryon acoustic oscillation (BAO) measurements we find 𝜉 DR < 10%, corresponding to a limit on the abundance of interacting DR that is around 3 times tighter than that from CMB + BAO data alone. We also discuss the complementarity of weak lensing informed cluster abundance studies with probes sensitive to smaller scales, explore the impact on our analysis of massive neutrinos, and comment on a slight preference for the presence of a nonzero interacting DR abundance, which enables a physical solution to the 𝑆 8 tension.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1992

Antipodal hotspot pairs on the earth

The results of statistical analyses performed on three published hotspot distributions suggest that significantly more hotspots occur as nearly antipodal pairs than is anticipated from a random distribution, or from their association with geoid highs and divergent plate margins. The observed number of antipodal hotspot pairs depends on the maximum allowable deviation from exact antipodality. At a maximum deviation of not greater than 700 km, 26 to 37 percent of hotspots form antipodal pairs in the published lists examined here, significantly more than would be expected from the general hotspot distribution. Two possible mechanisms that might create such a distribution include: (1) symmetry in the generation of mantle plumes, and (2) melting related to antipodal focusing of seismic energy from large-body impacts.

Rampino, Michael R.↗
NASA NTRS2007

Multi-Grid-Cell Validation of Satellite Aerosol Property Retrievals in INTEX/ITCT/ICARTT 2004

Aerosol transport off the US Northeast coast during the Summer 2004 International Consortium for Atmospheric Research on Transport and Transformation (ICARTT) Intercontinental Chemical Transport Experiment (INTEX) and Intercontinental Transport and Chemical Transformation (ITCT) experiments produced a wide range of aerosol types and aerosol optical depth (AOD) values, often with strong horizontal AOD gradients. In these conditions we flew the 14-channel NASA Ames Airborne Tracking Sun photometer (AATS) on a Jetstream 31 (J31) aircraft. Legs flown at low altitude (usually less than 100 m ASL) provided comparisons of AATS AOD spectra to retrievals for 90 grid cells of the satellite radiometers MODIS-Terra, MODIS-Aqua, and MISR, all over the ocean. Characterization of the retrieval environment was aided by using vertical profiles by the J31 (showing aerosol vertical structure) and, on occasion, shipboard measurements of light scattering and absorption. AATS provides AOD at 13 wavelengths lambda from 354 to 2138 nm, spanning the range of aerosol retrieval wavelengths for MODIS over ocean (466-2119 nm) and MISR (446-866 nm). Midvisible AOD on low-altitude J31 legs in satellite grid cells ranged from 0.05 to 0.9, with horizontal gradients often in the range 0.05 to 0.13 per 10 km. When possible, we used ship measurements of humidified aerosol scattering and absorption to estimate AOD below the J31. In these cases, which had J31 altitudes 60-110 m ASL (typical of J31 low-altitude transects), estimated midvisible AOD below the J31 ranged from 0.003 to 0.013, with mean 0.009 and standard deviation 0.003. These values averaged 6 percent of AOD above the 53 1. MISR-AATS comparisons on 29 July 2004 in 8 grid cells (each -17.6 km x 17.6 km) show that MISR versions 15 and 16 captured the AATS-measured AOD gradient (correlation coefficient R2 = 0.87 to 0.92), but the MISR gradient was somewhat weaker than the AATS gradient. The large AOD (midvisible values up to -0.9) and differing gradients in this case produced root-mean-square (RMS) MISR-AATS AOD differences of 0.03 to 0.21 (9 to 31%). MISR V15 Angstrom exponent alpha (= -dlnAOD/dln lambda) was closer to AATS than was MISR V16. MODIS-AATS AOD comparisons on 8 overpasses using 61 grid cells (each nominally 10 km x 10 km) had R2 approximately 0.97, with RMS AOD difference approximately 0.03 (approximately 20%). About 87% of the MODIS AOD retrievals differed from AATS values by less than the predicted MODIS over-ocean uncertainty, Delta tau = plus/minus 0.03 plus/minus 0.05 tau. In contrast to the small MODIS-AATS differences in AOD, MODIS-AATS differences in Angstrom exponent alpha were large: RMS differences for alpha (553, 855 nm) were 0.28 for MODIS-Terra and 0.64 for MODIS-Aqua; RMS differences for alpha (855, 2119 nm) were larger still, 0.61 for MODIS-Terra and 1.14 for MODIS-Aqua. The largest MODIS-AATS Angstrom exponent differences were associated with small AOD values, for which MODIS AOD relative uncertainty is large. Excluding cases with AOD(855 nm) less than 0.1 reduced MODIS-AATS alpha differences substantially. In one grid cell on 21 July 2004, smoke over cloud appeared to impair the MODIS-Aqua cloud mask, resulting in retrieved AODs that significantly exceeded AATS values. Experiments with extending MODIS retrievals into the glint mask yielded MODIS AODs consistently less than AATS AODs, especially at long wavelength, indicating that the current MODIS glint mask limits should not be reduced to the extent tried here. The sign of the AOD differences within the glint mask (MODIS AOD less than AATS AOD) is consistent with ship-measured wind speeds there.

Russell, P. B.↗
DOE OSTI2025

Probing the Higgs boson CP properties in vector-boson fusion production in the H → τ + τ − channel with the ATLAS detector

The CP properties of the Higgs boson are studied in the vector-boson fusion production mode. The analysis exploits the decay mode of the Higgs boson into two τ-leptons using 140 fb −1 of proton-proton collision data at $\sqrt{s}$ = 13 TeV collected by the ATLAS experiment at the Large Hadron Collider. Results are obtained using the Optimal Observable method. CP-violating interactions between the Higgs boson and electroweak gauge bosons are considered in the effective field theory framework, with the interaction strength described in the HISZ basis by $\tilde{d}$, and in the Warsaw basis by c H$\tilde{W}$ , c H$\tilde{B}$ , and c H$\tilde{W}$B . No deviations relative to the Standard Model are observed, and limits are obtained on the strength parameters. The $\tilde{d}$ parameter is constrained to the interval [−0.012, 0.044] at the 95% confidence level while c H$\tilde{W}$ is constrained to [−0.24, 0.83], when considering both linear and quadratic effects of physics beyond the Standard Model.

Hadron-Hadron Scattering↗
DOE OSTI2026

Direct Ab Initio Simulation of the Synthesis of BaZrO 3 and the Microstructure Impacts on Proton Transport

Controlling and predicting the processing-structure-performance relationship in functional materials is a grand challenge in materials science, with important implications for a wide range of emerging applications; a high fidelity understanding of the performance impact of microstructures formed under synthesis conditions is required to develop advanced materials, such as solid-state fuel cells and electrolyzers. Using the ceramic BaZrO 3 as a case study, we directly simulate the synthesis and investigate how proton transport is dictated by microstructures. We develop a framework that couples density functional theory (DFT), machine-learning interatomic potential (MLIP) driven molecular dynamics, and grand canonical Monte Carlo to perform large-scale, microstructure-resolved, atomistic simulations of proton transport in experimentally representative polycrystalline structures. Our fully ab initio approach, using a MLIP as a proxy for DFT, allows us to quantify the competition between two distinct diffusion mechanisms: one associated with grain-boundary regions and another within grains. When the impacts of grain boundaries are taken into account, proton transport exhibits substantial deviation from the bulk oxide limit. This addresses long-standing discrepancies between theory and experiments. Our integrated approach provides atomistic insight into microstructure-dependent proton pathways in BaZrO 3 and establishes a general protocol for predicting processing-structure-performance relationships.

organic↗
DOE OSTI2023

WPEC SG50: Developing an Automatically Readable, Comprehensive and Curated Experimental Nuclear Reaction Database

The Organisation for Economic Co-operation and Development (OECD) Nuclear Energy Agency (NEA) Working Party on International Nuclear Data Evaluation Co-operation (WPEC) subgroup (SG) 50 was formed in 2020 to develop an automatically readable, comprehensive and curated experimental nuclear reaction database. This database is called MEDUSAL (Machine-readable Experimental Data User Application & Library), and will draw from EXFOR. The EXFOR database preserves experimental nuclear reaction data true to its original documentation and information from the authors of the data. MEDUSAL will deviate from EXFOR by storing additional information from users of the data for their fields of work (evaluation, model development, validation, etc.). This includes expert judgment on the data sets, identification of data points as outliers, renormalization of the data to the newest monitor reactions, and estimations of missing uncertainty sources. The format for MEDUSAL is being developed to enable easy automatic parsing of large amounts of data. Here, we will summarize the use cases, high-level requirements and first steps towards developing the database MEDUSAL and the API to access it.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Validation of FAST.Farm Against Full-Scale Turbine SCADA Data for a Small Wind Farm

FAST.Farm is a new midfidelity engineering tool developed by the National Renewable Energy Laboratory targeted at accurately and efficiently predicting wind turbine power production and structural loading in wind farm settings, including wake interactions between turbines. FAST.Farm is based on some of the principles of the Dynamic Wake Meandering model—including passive tracer modeling of wake meandering—but addresses many of the limitations of previous Dynamic Wake Meandering (DWM) implementations. Previous FAST.Farm verification studies show the similarities and differences between FAST.Farm and large-eddy simulations for rigid and flexible turbines. In this validation study, FAST.Farm turbine responses are compared to multiturbine measurements from a subset of a full-scale wind farm. FAST.Farm predictions of turbine generator power, rotor speed, and blade pitch for five-turbine simulations are compared to supervisory control and data acquisition results. Results reveal that FAST.Farm generator power mean and standard deviation results reasonably match measured data for upstream and downstream turbines, as well as the mean rotor speed and blade pitch above rated wind speeds. However, FAST.Farm generally underpredicts the mean rotor speed and overpredicts the mean blade pitch below rated operation. These errors are likely related to inaccuracies in the generic controller simulated. Despite controller differences, FAST.Farm predicts the same overall relative rotor power trends for all waked turbines at all wind speeds.

17 WIND ENERGY↗
DOE OSTI2021

A new method for measuring the 3D turbulent velocity dispersion of molecular clouds

ABSTRACT The structure and star formation activity of a molecular cloud are fundamentally linked to its internal turbulence. However, accurately measuring the turbulent velocity dispersion is challenging due to projection effects and observational limitations, such as telescope resolution, particularly for clouds that include non-turbulent motions, such as large-scale rotation. Here, we develop a new method to recover the 3D turbulent velocity dispersion (σv,3D) from position–position–velocity (PPV) data. We simulate a rotating, turbulent, collapsing molecular cloud, and compare its intrinsic σv,3D with three different measures of the velocity dispersion accessible in PPV space: (1) the spatial mean of the 2nd-moment map, σi, (2) the standard deviation of the gradient/rotation-corrected 1st-moment map, σ(c − grad), and (3) a combination of (1) and (2), called the ‘gradient-corrected parent velocity dispersion’, $\sigma _{\mathrm{(p}-\mathrm{grad)}}=(\sigma _{\mathrm{i}}^2+\sigma _{(\mathrm{c}-\mathrm{grad)}}^2)^{1/2}$. We show that the gradient correction is crucial in order to recover purely turbulent motions of the cloud, independent of the orientation of the cloud with respect to the line of sight. We find that with a suitable correction factor and appropriate filters applied to the moment maps, all three statistics can be used to recover σv,3D, with method 3 being the most robust and reliable. We determine the correction factor as a function of the telescope beam size for different levels of cloud rotation, and find that for a beam full width at half-maximum f and cloud radius R, the 3D turbulent velocity dispersion can best be recovered from the gradient-corrected parent velocity dispersion via $\sigma _{v,\mathrm{3D}}= \left[(-0.29\pm 0.26)\, f/R + 1.93 \pm 0.15\right] \sigma _{\mathrm{(p}-\mathrm{grad)}}$ for f/R < 1, independent of the level of cloud rotation or LOS orientation.

Stewart, Madeleine↗