SEARCH · Engineering Papers

Results for “Thermodynamic limit”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Accelerating Thermochemical Equilibrium Calculations for Nuclear Reactor Applications

Thermochemical properties play a key role in modeling and simulation of several key phenomena in nuclear reactors. There has been an increasing interest in incorporating CALPHAD-based formulations in multiphysics simulations including for Molten Salt Reactors where knowledge of phase evolution of the salt and the chemical potentials of various elements are of utmost importance in source term analyses and redox control. However, the size of such simulations is often limited by the high computational cost of full thermodynamic equilibrium calculations. This work discusses the current efforts aimed at accelerating thermochemical equilibrium calculations for multiphysics simulations performed using the open-source finite element / finite volume code Multiphysics Object Oriented Simulation Environment (MOOSE) [1]. While several methods have been proposed for accelerating phase equilibrium calculations [2], most focus on relatively small systems and often rely on a- priori knowledge of the state-space of the system. Nuclear materials, however, are often multi-component systems owing to the evolution of composition under irradiation and an approach based on a-priori mapping of phase diagram is often not enough. This work is aimed at demonstrating an on-the-fly surrogate modeling framework that uses active learning to reduce the number of full equilibrium calculations that must be performed. By combining with efficient coupling approaches, the surrogate framework helps in reducing the computational cost of thermodynamic equilibrium informed multiphysics simulations of nuclear materials. The performance is benchmarked against full coupling with the thermochemistry library Thermochimica [3]. This work uses a machine learning based approach for constructing surrogate models to predict the stable phases in a multicomponent system. The surrogates were constructed using neural networks and Gaussian process classification. In this work, we compare the relative performance of the two methods. We also demonstrate the use of caching previous calculations by interpolating the values from nearest neighbors. References [1] Lindsay, A.D., et al. "2.0 – MOOSE: Enabling massively parallel multiphysics simulation", SoftwareX, 20 (2022): 101202. [2] Roos, W.A. and Zietsman J.H. "Accelerating complex chemical equilibrium calculations – A Review", Calphad, 77 (2022): 102380. [3] Piro, M.H.A., et al. "The thermochemistry library Thermochimica", Computational Materials Science, 67 (2013): 266-272.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Thermodynamics of rotating quantum matter in the virial expansion

We characterize the high-temperature thermodynamics of rotating bosons and fermions in two-dimensional (2D) and three-dimensional (3D) isotropic harmonic trapping potentials. In this work, we begin by calculating analytically the conventional virial coefficients $b_n$ for all $\textit{n}$ in the noninteracting case, as functions of the trapping and rotational frequencies. We also report on the virial coefficients for the angular momentum and associated moment of inertia. Using the $b_n$ coefficients, we analyze the deconfined limit (in which the angular frequency matches the trapping frequency) and derive explicitly the limiting form of the partition function, showing from the thermodynamic standpoint how both the 2D and 3D cases become effectively homogeneous 2D systems. To tackle the virial coefficients in the presence of weak interactions, we implement a coarse temporal lattice approximation and obtain virial coefficients up to third order.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1994

Lunar base heat pump, phase 1

This report describes the Phase 1 process and analysis used to select a refrigerant and thermodynamic cycle as the basis of a vapor compression heat pump requiring a high temperature lift, then to perform a preliminary design to implement the selected concept, including major component selection. Use of a vapor compression heat pump versus other types was based on prior work performed for the Electric Power Research Institute. A high lift heat pump is needed to enable a thermal control system to remove heat down to 275 K from a habitable volume when the external thermal environment is severe. For example, a long-term lunar base habitat will reject heat from a space radiator to a 325 K environment. The first step in the selection process was to perform an optimization trade study, quantifying the effect of radiator operating temperature and heat pump efficiency on total system mass; then, select the radiator operating temperature corresponding to the lowest system mass. Total system mass included radiators, all heat pump components, and the power supply system. The study showed that lunar night operation, with no temperature lift, dictated the radiator size. To operate otherwise would require a high mass penalty to store power. With the defined radiation surface, and heat pump performances assumed to be from 40 percent to 60 percent of the Carnot ideal, the optimum heat rejection temperature ranged from 387 K to 377 K, as a function of heat pump performance. Refrigerant and thermodynamic cycles were then selected to best meet the previously determined design conditions. The system was then adapted as a ground-based prototype lifting temperature to 360 K (versus 385 K for flight unit) and using readily available commercial-grade components. Over 40 refrigerants, separated into wet and dry compression behavioral types, were considered in the selection process. Refrigerants were initially screened for acceptable critical temperature. The acceptable refrigerants were analyzed in ideal single and two-stage thermodynamic cycles. Top candidates were analyzed assuming realistic component limits and system pressure drops, and were evaluated for other considerations such as safety, environmental impact, and commercial availability. A maximum coefficient of performance (COP) of 56 percent of the Carnot ideal was achievable for a three-stage CFC-11 cycle operating under the flight conditions above. The program was completed by defining a control scheme and by researching and selecting the major components, compressor and heat exchangers, that could be used to implement the thermodynamic cycle selected. Special attention was paid to using similar technologies for the SIRF and flight heat pumps resulting in the commercially available equivalent of the flight unit. A package concept was generated for the components selected and mass and volume estimated.

Goldman, Jeffrey H.↗

NASA NTRS1993

Space Life Support Engineering Program

This report covers the second year of research relating to the development of closed-loop long-term life support systems. Emphasis was directed toward concentrating on the development of dynamic simulation techniques and software and on performing a thermodynamic systems analysis in an effort to begin optimizing the system needed for water purification. Four appendices are attached. The first covers the ASPEN modeling of the closed loop Environmental Control Life Support System (ECLSS) and its thermodynamic analysis. The second is a report on the dynamic model development for water regulation in humans. The third regards the development of an interactive computer-based model for determining exercise limitations. The fourth attachment is an estimate of the second law thermodynamic efficiency of the various units comprising an ECLSS.

Seagrave, Richard C.↗

DOE OSTI2024

Dual-module humidity pump for efficient air dehumidification: Demonstration and performance limitations

Condensation dehumidification in conventional air conditioning technologies is energy-intensive, accounting for up to 50% of building cooling energy used in some climates. Selective vacuum membrane dehumidification (VMD) is one of the leading alternative dehumidification technologies due to its potential for significant energy savings, and the “dual-module humidity pump” is one of the most promising VMD concepts. Here, this work is the first to provide experimental proof-of-concept for the dual-module humidity pump system and provides the first thermodynamic modeling framework that accounts for realistic steady-state operating limitations, both of which are lacking in the current literature. Additionally, this work is the first to provide a system design solution that overcomes practical challenges associated with air accumulation in the vacuum channels. The experimental results in this work show that the current prototype can remove up to 45% of the water vapor in the air stream, and the vapor pressure difference in the vapor rejection module needs to be approximately 2–4 times greater than that of the dehumidification module in order to maintain balanced mass transfer. The thermodynamic model applied to typical air conditioning conditions shows that the ideal dehumidification (latent) COPs can reach up to 40, but practical COPs are limited to approximately 10. Furthermore, the model shows that the overall energy efficiency increases as the membrane air selectivity increases, though this improvement gradually starts to diminish when the membrane selectivity is increased above 10,000.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRS1995

Interface Instabilities in the Interstellar Medium

In the present communication, we reexamine two limiting cases of star-forming mechanisms involving self-gravity, thermodynamics, and velocity fields, that we believe must be ubiquitous in the ISM -- the generally oblique collision of supersonic gas streams or turbulent eddies. The general case of oblique collisions has not yet been examined. However, two limiting cases have been studied in detail: (1) The head-on collision of two identical gas streams that form dense, cool accretion shocks that become unstable and may form Jeans mass clouds, which subsequently undergo collapse. (2) Linearly unstable tangential velocity discontinuities, which result in Kelvin-Helmholtz (K-H) instabilities and related phenomena. The compressible K-H instabilities exhibit rich and unexpected behaviors. Moreover a new thermal-dynamic (T-D) mode was discovered that arises from the coupling of the perturbed thermal behavior and the unperturbed flow. The T-D mode has the curious characteristic that it may be strongly unstable to interface modes when the global modes in either medium are absolutely thermally stable. In the present communication additional models of case 1 are described and discussed, and self-gravity is added in the linear theory of tangential discontinuities, case 2. We prove that self-gravity fundamentally changes the behavior of interfacial modes -- density discontinuities (or steps) are inherently unstable on roughly the free-fall timescale of the denser medium to perturbations of all wavelengths.

Hunter, J. H., Jr.↗

DOE OSTI2022

Desalinating a Real Hyper-Saline Pre-Treated Produced Water via Direct-Heat Vacuum Membrane Distillation

Membrane distillation (MD) is an emerging thermal desalination technology capable of desalinating waters of any salinity. During typical MD processes, the saline feedwater is heated and acts as the thermal energy carrier; however, temperature polarization (as well as thermal energy loss) contributes to low distillate fluxes, low single-pass water recovery and poor thermal efficiency. An alternative approach is to integrate an extra thermal energy carrier as part of the membrane and/or module assembly, which can channel externally provided heat directly to the membrane-feedwater interface and/or along the feed channel length. This direct-heat delivery has been demonstrated to increase single-pass water recovery and enhance the overall thermal efficiency. We developed a bench-scale direct-heated vacuum MD (DHVMD) process to desalinate pre-treated oil and gas "produced water" with an initial total dissolved solids of 115,500 ppm at a feed temperature ranging between 24 and 32 degrees C. We evaluated both water flux and specific energy consumption (SEC) as a function of water recovery. The system achieved a 50% water recovery without significant scaling, with an average flux >6 kg m-2 hr-1 and a SEC as low as 2,530 kJ kg-1. The major species of mineral scales (i.e., NaCl, CaSO4, and SrSO4) that limited the water recovery to 68% were modeled in terms of thermodynamics and identified by scanning electron microscopy and energy-dispersive X-ray spectroscopy. In addition, we further developed and employed a physics-based process model to estimate temperature, salinity, water transport and energy flows for full-scale vacuum MD and DHVMD modules. Model results show that a direct-heat input rate of 3,600 W can increase single-pass water recovery from 2.1% to 3.1% while lowering the thermal SEC from 7,800 kJ kg-1 to 6,517 kJ kg-1 in an unoptimized module. Finally, the scaling up potential of DHVMD process is briefly discussed.

brine↗

DOE OSTI2020

Decoupling Kinetics and Thermodynamics of Interfacial Catalysis at a Chemically Modified Black Silicon Semiconductor Photoelectrode

Understanding the interplay between the kinetics of interfacial catalytic reactions and the thermodynamics of an underlying semiconductor electrode is imperative for rational construction of efficient photoelectrocatalytic systems. Current understanding of the thermodynamic effects of molecular catalyst attachment to semiconductor electrodes is limited. We report the immobilization of a molecular cobalt bis(benzenedithiolate) proton reduction catalyst onto nanoporous black silicon (b-Si) electrodes through p–p interactions with a series of aromatic molecules covalently attached to the surface. Furthermore, intensity-modulated high-frequency resistivity and linear sweep voltammetry measurements are used to show that the kinetics of proton reduction are decoupled from the thermodynamic properties of the underlying b-Si photoelectrode.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

From electronic structure to model application of key reactions for gasoline/alcohol combustion: Hydrogen-atom abstraction by $CH_3O\dot{O}$ radicals

Hydrogen atom abstraction by methyl peroxy ($CH_3O\dot{O}$) radicals can play an important role in gasoline/ethanol interacting chemistry for fuels that produce high concentrations of methyl radicals. Detailed kinetic reactions for hydrogen atom abstraction by $CH_3O\dot{O}$ radicals from the components of FGF-LLNL (a gasoline surrogate) including cyclopentane, toluene, 1-hexene, n -heptane, and isooctane have been systematically studied in this work. Here, the M06-2X/6-311 ++ G(d,p) level of theory was used to obtain the optimized structure and vibrational frequency for all stationary points and the low-frequency torsional modes. The 1-D hindered rotor treatment for low-frequency torsional modes was treated at M06-2X/6-31G level of theory. The UCCSD(T)-F12a/cc-pVDZ-F12 and QCISD(T)/CBS level of theory were used to calculate single point energies for all species. High pressure limiting rate constants for all hydrogen atom abstraction channels were performed using conventional transition state theory with unsymmetric tunneling corrections. Individual rate constants are reported in the temperature range from 298.15 to 2000 K. Our computed results show that the abstraction of allylic hydrogen atoms from 1-hexene is the fastest at low temperatures. When the temperature increases, the hydrogen atom abstraction reaction channel at the primary alkyl site gradually becomes dominant. Thermodynamics properties for all stable species and high-pressure limiting rate constants for each reaction pathway obtained in this work were incorporated into the latest gasoline surrogate/ethanol model to investigate the influence of the rate constants calculated here on model predicted ignition delay times.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2021

Thermodynamics of boron incorporation in BGaN

We study the thermodynamics of boron (B) incorporation into gallium nitride (GaN) using first-principles calculations. In the dilute limit, we have calculated the formation energies of different configurations of the B impurity in GaN and found that substitution on the cation site is favored over substitution on the anion site. Furthermore, under p-type conditions, interstitial boron can become the more favorable configuration and will ultimately limit the p-type conductivity. At higher B concentrations we use the generalized quasi-chemical approximation to elucidate the thermodynamic stability of boron gallium nitride (BGaN) alloys. We also investigate the effects of strain, which will be present if BGaN alloys are grown pseudomorphically on a GaN substrate.Without strain, B incorporation at typical growth conditions is limited to about 1.4% at 800 °C. Pseudomorphic strain raises the limit to 3.0% at the same temperature, close to experimentally observed levels of B incorporation.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Thermal and electric multidomain dynamic model for integration of power grid distribution with behind-the-meter devices

As renewable energy sources like solar and wind power become more integrated into the grid, coordinated control of behind-the-meter devices is crucial for enhancing grid flexibility and reliability and for meeting cost targets, with standardized models being developed to support this transition. The increasing flexibility and uncertainty of integrated renewable energy grids, along with interactions between various subsystems, make traditional steady-state modeling insufficient to capture transient and dynamic behaviors. Current models (e.g., composite load and battery equivalent models) focus on thermodynamic or electrical characteristics but overlook critical electromechanical interactions. This limits the ability to share performance information for grid services and hampers fast dynamic simulations. In addition, motor stalling is usually triggered by a fault event and attributed to the characteristics of the mechanical torque of the motor, resulting in absorption of a large amount of reactive power during the stalling period. Further, this significant withdrawal of reactive power will deteriorate the dynamic voltage stability of power grids and cause delayed voltage recovery. Therefore, an in-depth modeling of the thermodynamics or mechanical torque is essential to study the impacts of the realistic torque characteristics of those behind-the-meter devices on power system voltage stability. This study developed a dynamic multidomain model for building HVAC systems, such as air-source heat pumps, to simulate their thermal and electrical responses to grid transients. The model can accurately predict power metrics with a mean absolute percentage error of 10%, by validating against with power system computer-aided design performance data. Case studies demonstrate the model capability of capturing the transient response to sudden voltage changes, rapid load fluctuations, and system shutdowns respectively. During a sudden voltage drop (30% for 0.1s), a fully loaded heat pump’s motor speed dropped, continued declining, and shut down after 3.6s, with severe power oscillations and a torque spike. A partially loaded unit experienced temporary oscillations but stabilized. Under higher building loads, compressor speed increased from 64% to 100%, with power and torque rising before stabilizing. In safety-triggered shutdowns, power decreased after minor fluctuations, and torque briefly spiked before dropping to zero.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2020

The Role of Demixing and Crystallization Kinetics on the Stability of Non‐Fullerene Organic Solar Cells

Abstract With power conversion efficiency now over 17%, a long operational lifetime is essential for the successful application of organic solar cells. However, most non‐fullerene acceptors can crystallize and destroy devices, yet the fundamental underlying thermodynamic and kinetic aspects of acceptor crystallization have received limited attention. Here, room‐temperature (RT) diffusion coefficients of 3.4 × 10 −23 and 2.0 × 10 −22 are measured for ITIC‐2Cl and ITIC‐2F, two state‐of‐the‐art non‐fullerene acceptors. The low coefficients are enough to provide for kinetic stabilization of the morphology against demixing at RT. Additionally profound differences in crystallization characteristics are discovered between ITIC‐2F and ITIC‐2Cl. The differences as observed by secondary‐ion mass spectrometry, differential scanning calorimetry (DSC), grazing‐incidence wide‐angle X‐ray scattering, and microscopy can be related directly to device degradation and are attributed to the significantly different nucleation and growth rates, with a difference in the growth rate of a factor of 12 at RT. ITIC‐4F and ITIC‐4Cl exhibit similar characteristics. The results reveal the importance of diffusion coefficients and melting enthalpies in controlling the growth rates, and that differences in halogenation can drastically change crystallization kinetics and device stability. It is furthermore delineated how low nucleation density and large growth rates can be inferred from DSC and microscopy experiments which could be used to guide molecular design for stability.

Hu, Huawei↗

DOE OSTI2022

Dissociating the phononic, magnetic and electronic contributions to thermal conductivity: a computational study in alpha-iron

Computational tools to study thermodynamic properties of magnetic materials have, until recently, been limited to phenomenological modeling or to small domain sizes limiting our mechanistic understanding of thermal transport in ferromagnets. Herein, we study the interplay of phonon and magnetic spin contributions to the thermal conductivity in a-iron utilizing non-equilibrium molecular dynamics simulations. It was observed that the magnetic spin contribution to the total thermal conductivity exceeds lattice transport for temperatures up to two-thirds of the Curie temperature after which only strongly coupled magnon-phonon modes become active heat carriers. Characterizations of the phonon and magnon spectra give a detailed insight into the coupling between these heat carriers, and the temperature sensitivity of these coupled systems. Comparisons to both experiments and ab initio data support our inferred electronic thermal conductivity, supporting the coupled molecular dynamics/spin dynamics framework as a viable method to extend the predictive capability for magnetic material properties.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Pressure effects on reactivity and extinction of n -dodecane diffusion cool flame

Low-temperature combustion at high pressure has gained increasing interest owing to the growing demand for clean and efficient propulsion and energy systems. In particular, cool flame structure, extinction limits, and fuel reactivities at high pressure are critical in affecting the engine performance at near-limit conditions. Here, this paper investigates the effects of pressure on the cool flame extinction limit, structure, radical index, reactivity, and oxygen concentration dependence by experiments, analysis, and modeling. A large n-alkane, namely n-dodecane (n C 12 H 26 ), is selected to study its diffusion cool flame dynamics and reactivity in a high-pressure counterflow burner up to 10 atm. The experimental results show that higher pressure increases the cool flame extinction strain rates, and that the pressure-weighted extinction strain rate (aP) is proportional to the square of pressure, P 2 . A scaling analysis explains the relationship between the dependence of flame structure, heat release rate, and pressure-weighted cool flame strain rate on pressure. Furthermore, radical indexes at different pressures are measured by isolating the thermal and transport effects from the chemical contribution to diffusion cool flame extinction. The radical index clearly shows that the low-temperature reactivity increases with pressure. In addition, due to the critical role of multiple oxygen addition reactions in low-temperature chemistry, the relationship between the cool flame extinction limit and the oxygen concentration is explored. It is found that the cool flame extinction limits are proportional to the n th power of the oxygen concentration, [O 2 ] n , and increasing pressure leads to stronger extinction limit dependence (larger n) on the oxygen concentration. The present experiment and detailed kinetic analysis show clearly that increasing pressure promotes the low-temperature chemistry including the oxygen addition reactions, while the scaling analysis explains well the experimental results.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2020

Liquid–Solid Mixtures of Ga Metal Infused with Cu Microparticles and Nanoparticles for Microscale and Nanoscale Patterning of Solid Metals at Room Temperature

This paper demonstrates and characterizes a simple ink for nanopatterning of solid metallic structures under ambient conditions by taking advantage of the low melting point of gallium and its affinity to form intermetallics with other metals. Bare copper micro- and nanoparticles readily mix with liquid gallium near room temperature to form a paste that ultimately amalgamates into a dense solid when mixed at the appropriate concentration ratios. The paste has sufficient shelf life that can be extruded from a nozzle at modest temperatures to form solid three-dimensional (3D) shapes. Additionally, the paste can be molded at room temperature to replicate feature sizes ranging from a few millimeters down to hundreds of nanometers. In situ X-ray diffraction (XRD) and thermo-mechanical analysis (TMA) data show gallium and copper readily interdiffuse to form the thermodynamically expected intermetallic phase. We describe the capabilities and limitations of a simple way to pattern solid metals in an additive fashion (syringe-based extrusion) and with high resolution (molding) at or near room temperature. Here, the use of a paste that solidifies provides a novel route for 3D printing of solid metals at ambient temperatures as well as the creation of micro- and nanostructured metallic surfaces that may be useful for optics, non-wetting surfaces, or electronic microcomponents.

36 MATERIALS SCIENCE↗

DOE OSTI2023

High-throughput calculations of charged point defect properties with semi-local density functional theory—performance benchmarks for materials screening applications

Abstract Calculations of point defect energetics with Density Functional Theory (DFT) can provide valuable insight into several optoelectronic, thermodynamic, and kinetic properties. These calculations commonly use methods ranging from semi-local functionals with a-posteriori corrections to more computationally intensive hybrid functional approaches. For applications of DFT-based high-throughput computation for data-driven materials discovery, point defect properties are of interest, yet are currently excluded from available materials databases. This work presents a benchmark analysis of automated, semi-local point defect calculations with a-posteriori corrections, compared to 245 “gold standard” hybrid calculations previously published. We consider three different a-posteriori correction sets implemented in an automated workflow, and evaluate the qualitative and quantitative differences among four different categories of defect information: thermodynamic transition levels, formation energies, Fermi levels, and dopability limits. We highlight qualitative information that can be extracted from high-throughput calculations based on semi-local DFT methods, while also demonstrating the limits of quantitative accuracy.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Overcoming immiscibility toward bimetallic catalyst library

Bimetallics are emerging as important materials that often exhibit distinct chemical properties from monometallics. However, there is limited access to homogeneously alloyed bimetallics because of the thermodynamic immiscibility of the constituent elements. Overcoming the inherent immiscibility in bimetallic systems would create a bimetallic library with unique properties. Here, we present a nonequilibrium synthesis strategy to address the immiscibility challenge in bimetallics. As a proof of concept, we synthesize a broad range of homogeneously alloyed Cu-based bimetallic nanoparticles regardless of the thermodynamic immiscibility. The nonequilibrated bimetallic nanoparticles are further investigated as electrocatalysts for carbon monoxide reduction at commercially relevant current densities (>100 mA cm -2 ), in which Cu 0.9 Ni 0.1 shows the highest multicarbon product Faradaic efficiency of ~76% with a current density of ~93 mA cm -2 . The ability to overcome thermodynamic immiscibility in multimetallic synthesis offers freedom to design and synthesize new functional nanomaterials with desired chemical compositions and catalytic properties.

25 ENERGY STORAGE↗

DOE OSTI2023

Charged point defect benchmark of Hybrid and GGA-PBE

Data used by the publication "High-throughput calculations of charged point defect properties with semi-local density functional theory - performance benchmarks for materials screening applications." This work presented an in-depth benchmark analysis of automated, semi-local point defect calculations with a-posteriori corrections, compared to 245 “gold standard” hybrid calculations previously published. We considered three different a-posteriori correction sets for semi-local calculations, implemented in a fully automated workflow, and consider the qualitative and quantitative differences for four different categories of defect information: thermodynamic transition levels, formation energies, fermi levels, and dopability limits. We highlighted the type of qualitative information about point defect properties that can be extracted from high-throughput calculations based on semi-local DFT methods, while also demonstrating the limits of quantitative accuracy that can be achieved by these approaches.

Broberg, Danny↗