Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion couples”

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 991 records · Page 55

Self-organization of nanoparticles and molecules in periodic Liesegang-type structures

Chemical organization in reaction-diffusion systems offers a strategy for the generation of materials with ordered morphologies and structural hierarchy. Periodic structures are formed by either molecules or nanoparticles. On the premise of new directing factors and materials, an emerging frontier is the design of systems in which the precipitation partners are nanoparticles and molecules. We show that solvent evaporation from a suspension of cellulose nanocrystals (CNCs) and L -(+)-tartaric acid [ L -(+)-TA] causes phase separation and precipitation, which, being coupled with a reaction/diffusion, results in rhythmic alternation of CNC-rich and L -(+)-TA–rich rings. The CNC-rich regions have a cholesteric structure, while the L -(+)-TA–rich bands are formed by radially aligned elongated bundles. The moving edge of the pattern propagates with a finite constant velocity, which enables control of periodicity by varying film preparation conditions. This work expands knowledge about self-organizing reaction-diffusion systems and offers a strategy for the design of self-organizing materials.

36 MATERIALS SCIENCE↗

Dynamics of heavy quarks in strongly coupled $\mathcal{N}$ = 4 SYM plasma

We calculate the probability distribution P(k) for a heavy quark with velocity v propagating through strongly coupled N = 4 SYM plasma in the ’t Hooft limit (N c → ∞, λ = g 2 N c → ∞) at a temperature T to acquire a momentum k due to interactions with the plasma. This distribution encodes the well-known drag coefficient η D and the transverse and longitudinal momentum diffusion coefficients κ T and κ L . The jet quenching parameter $\hat{q}$ can be extracted from P(k) for v = 1. Going beyond these known Gaussian characteristics of P(k), our calculation determines all of the higher order and mixed moments to leading order in 1/$\sqrt{λ}$ for the first time. These non-Gaussian features of P(k) include qualitatively novel correlations between longitudinal energy loss and transverse momentum broadening at nonzero v. We show that all higher moments scale characteristically with an effective temperature of the boosted plasma in the heavy quark rest frame, and we demonstrate that these non-Gaussian characteristics can be sizable in magnitude and even dominant in physically relevant situations. We use these results to derive a Kolmogorov equation for the evolution of the probability distribution for the total momentum of a heavy quark that propagates through strongly coupled plasma. This evolution equation accounts for all higher order correlations between transverse momentum broadening and longitudinal energy loss, which we have calculated from first principles. It reduces to a Fokker-Planck equation when truncated to only include the effects of η D , κ T and κ L . Remarkably, while heavy quarks do not reach kinetic equilibrium with the plasma if evolved with this Fokker-Planck equation, by showing that the Boltzmann distribution is a static solution of the all-order Kolmogorov equation that we have derived we demonstrate that heavy quarks do reach kinetic equilibrium if evolved with this equation. Our results thus provide a dynamically complete framework for understanding the thermalization of a heavy quark that may be initially far from equilibrium in the strongly coupled N = 4 SYM plasma — as well as new insight into heavy quark transport and equilibration in quark-gluon plasma.

Holography and Hydrodynamics↗

Biobased Poly(dodecylene Furanoate) with Inherent Advantages in Performance and Circularity

Biobased polymers are gaining traction toward more sustainable flexible-film packaging, yet overcoming trade-offs between their performance properties and end-of-life (EoL) options still remains a challenge. Here, it is shown that biobased poly(dodecylene 2,5-furanoate) (PDDF), synthesized via both step-growth polycondensation and chain-growth ring-opening polymerization methods, exhibits advantages not only in gas barrier properties but also in EoL options due to its biodegradability and closed-loop chemical circularity. Specifically, PDDF displays significantly lower oxygen and carbon dioxide permeability than commercial poly(butylene adipate-co-terephthalate) (PBAT) and linear low-density polyethylene , alongside a markedly higher modulus (by ≈3 ×) and reduced water vapor transmission rate compared to PBAT. This superior performance is attributed to the inherently rigid, polar, H-bonding furan rings that enhance chain interaction, packing and crystallinity and thus reduce free volume impeding gas diffusion, while the long hydrophobic dodecylene segments inhibit water permeation. Furthermore, PDDF can be recycled back to its cyclic monomer by base-catalyzed depolymerization or diester and diol monomers by simple methanolysis. These superior barrier properties, coupled with biodegradation and closed-loop circularity, highlight the potential of the biobased PDDF as a more sustainable alternative for packaging.

Poly(dodecylene 2,5-furanoate) (PDDF)↗

A meshless stochastic method for Poisson–Nernst–Planck equations

A plethora of biological, physical, and chemical phenomena involve transport of charged particles (ions). Its continuum-scale description relies on the Poisson–Nernst–Planck (PNP) system, which encapsulates the conservation of mass and charge. The numerical solution of these coupled partial differential equations is challenging and suffers from both the curse of dimensionality and difficulty in efficiently parallelizing. We present a novel particle-based framework to solve the full PNP system by simulating a drift–diffusion process with time- and space-varying drift. We leverage Green’s functions, kernel-independent fast multipole methods, and kernel density estimation to solve the PNP system in a meshless manner, capable of handling discontinuous initial states. The method is embarrassingly parallel, and the computational cost scales linearly with the number of particles and dimension. We use a series of numerical experiments to demonstrate both the method’s convergence with respect to the number of particles and computational cost vis-à-vis a traditional partial differential equation solver.

Chemistry↗

Neutronics and thermal-hydraulics simulation of generic pebble-bed fluoride-salt-cooled high-temperature reactor (gFHR)

The fluoride-salt-cooled high-temperature reactor (FHR) is one type of the advanced reactors and has been attracting great interest from the research institutes and commercial companies in the recent years. However, currently the technology is relatively immature. To facilitate the design and safety analysis of FHRs, Kairos Power published the generic FHR (gFHR) benchmark. In this paper, a hybrid method combining the stochastic code and deterministic code is developed to simulate the gFHR benchmark. Serpent 2 is employed as a few-group cross section generator and the cross sections are applied to the finite difference neutron diffusion code AGREE. A good consistency with the gFHR benchmark is achieved. The agreement between the Serpent 2 results and the AGREE results shows that the hybrid method is applicable to FHRs. The thermal-hydraulics are coupled to neutronics for the steady-state calculations and good agreement between AGREE and SAM is achieved. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

PDF approach for compressible turbulent reacting flows

The objective of the present work is to develop a probability density function (pdf) turbulence model for compressible reacting flows for use with a CFD flow solver. The probability density function of the species mass fraction and enthalpy are obtained by solving a pdf evolution equation using a Monte Carlo scheme. The pdf solution procedure is coupled with a compressible CFD flow solver which provides the velocity and pressure fields. A modeled pdf equation for compressible flows, capable of capturing shock waves and suitable to the present coupling scheme, is proposed and tested. Convergence of the combined finite-volume Monte Carlo solution procedure is discussed, and an averaging procedure is developed to provide smooth Monte-Carlo solutions to ensure convergence. Two supersonic diffusion flames are studied using the proposed pdf model and the results are compared with experimental data; marked improvements over CFD solutions without pdf are observed. Preliminary applications of pdf to 3D flows are also reported.

Hsu, A. T.↗

Dissolution kinetics of a lunar glass simulant at 25 degrees C: the effect of pH and organic acids

The dissolution kinetics of a simulated lunar glass were examined at pH 3, 5, and 7. Additionally, the pH 7 experiments were conducted in the presence of citric and oxalic acid at concentrations of 2 and 20 mM. The organic acids were buffered at pH 7 to examine the effect of each molecule in their dissociated form. At pH 3, 5, and 7, the dissolution of the synthetic lunar glass was observed to proceed via a two-stage process. The first stage involved the parabolic release of Ca, Mg, Al, and Fe, and the linear release of Si. Dissolution was incongruent, creating a leached layer rich in Si and Ti which was verified by transmission electron microscopy (TEM). During the second stage the release of Ca, Mg, Al, and Fe was linear. A coupled diffusion/surface dissolution model was proposed for dissolution of the simulated lunar glass at pH 3, 5, and 7. During the first stage the initial release of mobile cations (i.e., Ca, Mg, Al, Fe) was limited by diffusion through the surface leached layer of the glass (parabolic release), while Si release was controlled by the hydrolysis of the Si-O-Al bonds at the glass surface (linear release). As dissolution continued, the mobile cations diffused from greater depths within the glass surface. A steady-state was then reached where the diffusion rate across the increased path lengths equalled the Si release rate from the surface. In the presence of the organic acids, the dissolution of the synthetic lunar glass proceeded by a one stage process. The release of Ca, Mg, Al, and Fe followed a parabolic relationship, while the release of Si was linear. The relative reactivity of the organic acids used in the experiments was citrate > oxalate. A thinner leached layer rich in Si/Ti, as compared to the pH experiments, was observed using TEM. Rate data suggest that the chemisorption of the organic anion to the surface silanol groups was responsible for enhanced dissolution in the presence of the organic acids. It is proposed that the increased rate of Si release is responsible for the one stage parabolic release of mobile cations and the relatively thin leached layer compared to experiments at pH 3 and 5.

NASA Discipline Life Support Systems↗

Rain Reevaporation, Boundary Layer Convection Interactions, and Pacific Rainfall Patterns in an AGCM

Sensitivity experiments with an atmospheric general circulation model (AGCM) show that parameterized rain re-evaporation has a large impact on simulated precipitation patterns in the tropical Pacific, especially on the configuration of the model s intertropical convergence zone (ITCZ). Weak re-evaporation leads t o the formation of a "double ITCZ" during the northern warm season. The double ITCZ is accompanied by strong coupling between precipitation and high-frequency vertical motion in the planetary boundary layer (PBL). Strong reevaporation leads to a better overall agreement of simulated precipitation with observations. The model s double ITCZ bias is reduced. At the same time, correlation between high-frequency vertical motion in the PBL and precipitation is reduced. Experiments with modified physics suggest that evaporative cooling by rain near the PBL top weakens the coupling between precipitation and vertical motion. This may reduce the model s tendency to form double ITCZs. The strength of high-frequency vertical motions in the PBL was also reduced directly through the introduction of a diffusive cumulus momentum transport (DCMT) parameterization. The DCMT had a visible impact on simulated precipitation in the tropics, but did not reduce the model s double bias in all cases.

Bacmeister, Julio T.↗

Coupled Gel/Seeded Growth of Cs2AgBiBr6 Single Crystals for Radiation Detection

Abstract— The all-inorganic lead free Cs2AgBiBr6 perovskite represents a highly promising environment-friendly radiation detector material. The first attempt to grow Cs2AgBiBr6 single crystals using the unique gel growth method is reported in this study. This technique enables controlled reagent diffusion in a gel medium at room temperature, fostering the growth of visible single crystals. It offers advantages over conventional methods, leading to fewer equilibrium and non-equilibrium defects and simplified execution. Interestingly, the as-grown crystals exhibited variations in size and colors within the gel matrix, indicating the influence of Bi solutes in Bi-poor and Bi-rich Cs2AgBiBr6 crystals. Once the crystals were confirmed by powder xray diffraction, a series of comprehensive material characterizations was performed. The resistivity of Cs2AgBiBr6 single crystals was measured at 1.12x107 Ω•cm. Trap density analysis yielded a value of 2.82 x1010 cm-3. The charge carrier mobility was determined to be 32.9 cm2V-1S-1, suggesting the suitability for radiation detector applications. Additional growth optimization is going on and further characterization will be reported during the presentation. Our work helps accelerate the development of lead-free perovskite materials toward radiation detector applications.

James, Ralph B.↗

Enhanced diffusion bonding of alloy 617 using electric field-assisted sintering

The development of compact heat exchangers (CHXs) has gained increasing interest in many industries owing to their high thermal efficiency and reduced size. Diffusion bonding (DB) is an advantageous technique for fabricating CHXs. Alloy 617 is a candidate for manufacturing CHXs for high-temperature advanced nuclear reactors due to its elevated-temperature properties. Previous endeavors in DB of Alloy 617 were conducted by hot pressing (HP), which reported precipitates at the diffusion-bond interface, limited grain boundary (GB) migration, and significantly reduced high-temperature mechanical properties. To overcome these challenges, this study investigated DB of Alloy 617 using electric field-assisted sintering (EFAS). Stacks composed of three sheets were bonded with EFAS using different temperatures, pressures, and hold times. DB using HP as the zero-current analog of EFAS was also performed for comparison. The result shows that Cr- and Mo-rich precipitates were formed at the interface of the hot-pressed samples. The electric current and temperature in EFAS play a significant role in precipitation and GB migration. The electric current coupled with correct temperatures can effectively prevent precipitate formation at the interface and achieve excellent GB migration. Nanoscale Al-rich oxide was formed at the interface of the samples made by both HP and EFAS, but grain boundaries can ignore the nanoscale Al-oxide and migrate across the interface. The temperature, pressure, and hold time also affected diffusion. The temperature is a prerequisite for a successful GB migration, and GB migration can be enhanced by increasing pressure and hold time.

36 MATERIALS SCIENCE↗

A non‐intrusive domain‐decomposition model reduction method for linear steady‐state partial differential equations with random coefficients

Abstract Domain decomposition methods have been proved to be an effective strategy to reduce the dimension of parametric partial differential equations (PDEs). However, existing domain decomposition methods for parametric PDEs are usually intrusive, which means domain decomposition based solvers need to be implemented from scratch for each target parametric PDE. To address this issue, we develop a new non‐intrusive domain‐decomposition model reduction method for linear steady‐state PDEs with random‐field coefficients. As a variant of our previous work by Mu and Zhang, the new method only needs access to the final linear system, that is, the global stiffness matrix and the right hand side, of a deterministic PDE solver, in order to build a domain‐decomposition‐based reduced model without intrusive implementation from scratch. The key idea is to remove the interface condition between sub‐domains and rely on the correlation between columns of the linear system to couple the sub‐domains. The non‐intrusive feature enables the applicability of the proposed method to a broader class of uncertainty quantification problems, where many legacy codes/solvers can be fully reused by our method. Two numerical examples including diffusion equations with random diffusivity and convection‐dominated transport with random velocity, are provided to demonstrate the effectiveness and efficiency of our method.

Zhang, Guannan↗

Pt Nanoparticle Disintegration at Oxide Interfaces Enhances CO Oxidation Catalysis

Understanding how supported metal nanoparticles dynamically evolve under reaction conditions is critical for controlling their catalytic function. Here, the mechanism behind the dynamic disintegration of Pt nanoparticles (NPs) supported on CeOx-TiO2 (CT) during CO oxidation is elucidated, leading to the formation of single atoms (SAs) and/or sub-nanometer clusters. Density functional theory (DFT) calculations reveal that strong Pt-CO interactions weaken Pt─Pt cohesion, while electronic coupling between Pt and Ce ions stabilizes Pt-CO* intermediates at the oxide interface. Surface oxygen vacancies kinetically trap Pt-CO*, but the vacancies are replenished under oxygen-rich conditions, enabling Pt-CO* surface diffusion and subsequent structural reorganization. In situ spectroscopic analyses confirm the oxygen-driven transformation of Pt NPs, correlating with a threefold increase in mass-specific activity at 150 °C. These findings demonstrate that interfacial oxygen dynamics and metal-support interactions can be leveraged to induce nanoparticle disintegration and optimize catalytic performance, highlighting the catalytic potential of interface-engineered Pt nanostructures.

CO oxidation↗

Scaling microstructural processes in the sintering of ionic ceramics

A multi-scale framework, combining a multiphase field formulation and large deformation mechanics, was developed as a stepping stone to perform the data analytics of the microstructural level kinetics of a sintering solid. Relevant microstructural information from this framework, such as grain, stress, and porosity statistics, was scaled up to describe the macroscopic level sintering kinetics. Here, the developed formulation was applied to describe the electric field assisted sintering of Y 2 O 3 . Microstructural inhomogeneities in a multi-granular solid result in the formation of a field of compressive stress networks, which interleave with low compression and weakly tensile regions, defining a scaffolding for sintering concentration regions to develop. A Poisson effect-induced lateral stress network is also naturally self-induced as a result of the mechanical constraints imposed by the sintering apparatus. For long sintering times, localized shear stresses enhancing mass flow along grain boundaries and internal surfaces develop. Three-sided pores are removed by either vacancy transport to the surrounding pores, or move towards the external surfaces through grain boundary diffusion. Four- and higher order-sided pores stabilize because an equal amount of vacancies are gained and lost through the connecting grain boundaries. Grain dewetting contributes to pore coalescence, suggesting that pore kinetics and grain growth are coupled and should be analyzed in concert. The combined sintering and grain growth kinetics define six regimes of sintering behavior: (1) T, the transient regime; (2) E$_Υ$, the surface energy dominated, early sintering regime, where the grain growth exponent, p = 1, and the stress concentration factor, $f$ ~ $1/\hat{ρ}^{4.6}$; (3) E S , the stress dominated, early sintering regime, where p = 1 and $f$ ~ $1/\hat{ρ}^{4}$; (4) I$_Υ$, the surface energy dominated, intermediate sintering regime, where p = 2 and $f$ ~ $1/\hat{ρ}^{4.6}$; (5) I S , the stress dominated, intermediate sintering regime, where p = 2 and $f$ ~ $1/\hat{ρ}^{4}$; and (6) L, the late sintering regime, where p = 3 and $f$ ~ 1. At the macroscopic level, the rapid densification and suppression of grain growth observed in the electric field assisted sintering process is a consequence of the compounding effects of the underlying stress-, transport-, and interfacial-energy-induced energy minimization kinetics, as predicted by the multi-scale framework.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Probing the role of local tunnel variations in early-stage lithiation of α-MnO₂ nanowires via in situ TEM

Understanding lithium-ion transport in tunnel-structured manganese oxides is essential for designing high-performance lithium-ion battery electrode materials. Here, we elucidate the early-stage lithiation mechanism of potassium-stabilized α-MnO 2 nanowires using in situ transmission electron microscopy (TEM) coupled with electron energy-loss spectroscopy (EELS), high-resolution TEM (HRTEM), and geometric phase analysis (GPA). Real-time TEM imaging reveals clear volume expansion at the reaction front, while EELS analysis uncovers lithium-ion diffusion far beyond this region, where no visible expansion is observed, indicating fast, defect-assisted transport. GPA and HRTEM analyses show that localized tensile and compressive strain fields, originating from pre-existing local tunnel structural variations, persist after lithiation. The tensile-strained regions enable lithium-ion insertion with minimal lattice distortion, offering additional free volume that facilitates rapid lithium-ion accommodation ahead of the structural transformation. Our results demonstrate a local tunnel variation-mediated fast diffusion pathway that precedes bulk reaction, underscoring the critical role of local strain in enabling early-stage lithium transport. Given the structural versatility of MnO 2 and its ability to accommodate diverse atomic arrangements beyond the well-known tunnel phases (β-, γ-, δ-, λ-, R-phases), our findings highlight the importance of understanding and engineering local structural environments. This work provides fundamental insights into the interplay between defects, strain, and ion dynamics, and presents defect engineering as a promising approach to enhance both rate performance and structural stability in manganese-based cathodes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A mechanistic interpretation of Nelson curves for PVP failures under high temperature hydrogen attack

As an empirically established design criterion, Nelson curves that relate the service temperature and the allowable hydrogen partial pressure have been developed and utilized for more than sixty years in pressure vessels and piping (PVP) safety design. Despite a relatively clear thermodynamic understanding of the high-temperature-hydrogen-attack (HTHA) problem, the detailed fracture process on the microstructural length scales, however, remains elusive, and a quantitative assessment of the PVP lifetime under HTHA from the available creep fracture dataset is still not possible. This work develops a microstructure-informed and micromechanics-based model by incorporating a synergy between hydrogen transport and intergranular-cavity-based fracture process. Based on the available creep lifetime data of C-0.5Mo steels, we are able to calibrate material constitutive parameters, and then conduct nonlinear finite element simulations that reveal a real-time stress-induced hydrogen diffusional transport along grain boundaries, coupled with a microstructure-explicit failure process, from which Nelson curves can be computed. Such failure analyses allow us to delineate two distinct regimes on the Nelson curves, i.e., dislocation-creep-controlled or grain boundary diffusion-assisted cavity growth. More importantly, we found that a small change of the pipe thickness and applied stresses can significantly shift these lifetime curves. However, these two parameters are usually not provided in Nelson curves, thus limiting their usage in material selection and safety design. In conclusion, this discrepancy can clearly be mitigated by extensive parametric studies from our micromechanical modeling/simulation framework.

36 MATERIALS SCIENCE↗

Concentration Field of Reactants and Products Species in a Reacting Vortex Ring

The proposed paper will present experimental and numerical results on the concentration fields of both reactants and products species in a reacting vortex ring that is generated from the interaction between a diffusion flame and a laminar vortex ring. Flame-vortex interactions are canonical configurations used to study the underlying processes occurring in complicated turbulent reacting flows. This type of configuration contains many of the fundamental aspects of the coupling between fluid dynamics and combustion that could be investigated with more controllable conditions than are possible under direct investigations of turbulent flames. The current configuration has been studied experimentally by Chen and Dahm and Chen et al. under microgravity conditions, and by Park and Shin, and You et al. under normal gravity conditions. This configuration is similar to that used in the analyses of Karagozian and Manda of their 2-D vortex pair in which both fuel and entrained oxidizer are present. The vortex ring used in this study is generated by issuing methane into an air environment through the exit of an axisymmetric nozzle. The experiments were conducted under microgravity conditions in order to remove the undesirable effects of buoyancy that can affect both the flame structure and ring dynamics resulting in possibly asymmetric and nonrepeatable interactions. The experimental technique of diode laser wavelength modulation spectroscopy (WMS) is used to measure concentration fields of reactants, CH4 and O2, products, H2O, CO2, OH, and temperature fields which can be inferred from either line pairs of O2 or OH lines. This technique has been investigated previously by Silver and Bomse et al. This is the first time that the technique has been applied to reacting vortex rings under microgravity conditions. The effect of ring circulation and fuel volume on the species concentration fields will be investigated. The experimental results will be compared to the current numerical results, and used to validate the numerical studies. In addition, the existence of burned cores during the interactions will be determined, and the increase in reactant consumption with increased ring circulation will be examined. Numerical studies were also conducted by solving the Navier-Stokes and mixture fraction equations with the assumptions of unity Lewis and Schmidt numbers. Equilibrium chemistry and flamelet libraries were used to obtain the temperature and species mass fraction fields. The numerical results will serve as guidelines in conducting the experimental studies. Ring circulation and fuel volume effects on the interactions and species concentration fields will be investigated and compared to experimental results.

Chen, Shin-Juh↗

High density protein crystal growth

A protein crystal growth assembly including a crystal growth cell and further including a cell body having a top side and a bottom side and a first aperture defined therethrough, the cell body having opposing first and second sides and a second aperture defined therethrough. A cell barrel is disposed within the cell body, the cell barrel defining a cavity alignable with the first aperture of the cell body, the cell barrel being rotatable within the second aperture. A reservoir is coupled to the bottom side of the cell body and a cap having a top side is disposed on the top side of the cell body. The protein crystal growth assembly may be employed in methods including vapor diffusion crystallization, liquid to liquid crystallization, batch crystallization, and temperature induction batch mode crystallization.

Rouleau, Robyn↗

A Powerful Toolkit for Synthetic Biology: Over 3.8 Billion Years of Evolution

The combination of evolutionary with engineering principles will enhance synthetic biology. Conversely, synthetic biology has the potential to enrich evolutionary biology by explaining why some adaptive space is empty, on Earth or elsewhere. Synthetic biology, the design and construction of artificial biological systems, substitutes bio-engineering for evolution, which is seen as an obstacle. But because evolution has produced the complexity and diversity of life, it provides a proven toolkit of genetic materials and principles available to synthetic biology. Evolution operates on the population level, with the populations composed of unique individuals that are historical entities. The source of genetic novelty includes mutation, gene regulation, sex, symbiosis, and interspecies gene transfer. At a phenotypic level, variation derives from regulatory control, replication and diversification of components, compartmentalization, sexual selection and speciation, among others. Variation is limited by physical constraints such as diffusion, and chemical constraints such as reaction rates and membrane fluidity. While some of these tools of evolution are currently in use in synthetic biology, all ought to be examined for utility. A hybrid approach of synthetic biology coupled with fine-tuning through evolution is suggested

Rothschild, Lynn J.↗