Thermal Conduction in Polymeric Nanofluids and Nanosolids Controlled by Interfacial Scattering:Solutions to Some Selected Problems

Dept.of Instrumentation,Cochin University of Science and Technology

The role of copper in the thermal conductivity of thermoelectric oxychalcogenides: do lone pairs matter?

Understanding the underlying mechanisms that suppress thermal conduction in solids is of paramount importance for the targeted design of materials for thermal management and thermoelectric energy conversion applications. Bismuth copper oxychalcogenides, BiOCuQ (Q = Se, Te), are highly crystalline thermoelectric materials with an unusually low lattice thermal conductivity of approx. 0.5 Wm-1K-1, a value normally found in amorphous materials. Here we unveil the origin of the unusual thermal transport properties of these phases. First principles calculations of the vibrational properties combined with analysis of in-situ neutron diffraction data, demonstrate that weak bonding of copper atoms within the structure leads to an unexpected vibrational mode at low frequencies, which is likely to be a major contributor to the low thermal conductivity of these materials. In addition, we show that anharmonicity and the large Grüneisen parameter in these oxychalcogenides are mainly related to the low frequency copper vibrations, rather than to the Bi3+ lone pairs.

Thermal, stress, and traps effects in AlGaN/GaN HEMTs

GaN y AlN son materiales semiconductores piezoeléctricos del grupo III-V. La heterounión AlGaN/GaN presenta una elevada carga de polarización tanto piezoeléctrica como espontánea en la intercara, lo que genera en su cercanía un 2DEG de grandes concentración y movilidad. Este 2DEG produce una muy alta potencia de salida, que a su vez genera una elevada temperatura de red. Las tensiones de puerta y drenador provocan un stress piezoeléctrico inverso, que puede afectar a la carga de polarización piezoeléctrica y así influir la densidad 2DEG y las características de salida. Por tanto, la física del dispositivo es relevante para todos sus aspectos eléctricos, térmicos y mecánicos. En esta tesis se utiliza el software comercial COMSOL, basado en el método de elementos finitos (FEM), para simular el comportamiento integral electro-térmico, electro-mecánico y electro-térmico-mecánico de los HEMTs de GaN. Las partes de acoplamiento incluyen el modelo de deriva y difusión para el transporte electrónico, la conducción térmica y el efecto piezoeléctrico. Mediante simulaciones y algunas caracterizaciones experimentales de los dispositivos, hemos analizado los efectos térmicos, de deformación y de trampas. Se ha estudiado el impacto de la geometría del dispositivo en su auto-calentamiento mediante simulaciones electro-térmicas y algunas caracterizaciones eléctricas. Entre los resultados más sobresalientes, encontramos que para la misma potencia de salida la distancia entre los contactos de puerta y drenador influye en generación de calor en el canal, y así en su temperatura. El diamante posee une elevada conductividad térmica. Integrando el diamante en el dispositivo se puede dispersar el calor producido y así reducir el auto-calentamiento, al respecto de lo cual se han realizado diversas simulaciones electro-térmicas. Si la integración del diamante es en la parte superior del transistor, los factores determinantes para la capacidad disipadora son el espesor de la capa de diamante, su conductividad térmica y su distancia a la fuente de calor. Este procedimiento de disipación superior también puede reducir el impacto de la barrera térmica de intercara entre la capa adaptadora (buffer) y el substrato. La muy reducida conductividad eléctrica del diamante permite que pueda contactar directamente el metal de puerta (muy cercano a la fuente de calor), lo que resulta muy conveniente para reducir el auto-calentamiento del dispositivo con polarización pulsada. Por otra parte se simuló el dispositivo con diamante depositado en surcos atacados sobre el sustrato como caminos de disipación de calor (disipador posterior). Aquí aparece una competencia de factores que influyen en la capacidad de disipación, a saber, el surco atacado contribuye a aumentar la temperatura del dispositivo debido al pequeño tamaño del disipador, mientras que el diamante disminuiría esa temperatura gracias a su elevada conductividad térmica. Por tanto, se precisan capas de diamante relativamente gruesas para reducer ele efecto de auto-calentamiento. Se comparó la simulación de la deformación local en el borde de la puerta del lado cercano al drenador con estructuras de puerta estándar y con field plate, que podrían ser muy relevantes respecto a fallos mecánicos del dispositivo. Otras simulaciones se enfocaron al efecto de la deformación intrínseca de la capa de diamante en el comportamiento eléctrico del dispositivo. Se han comparado los resultados de las simulaciones de la deformación y las características eléctricas de salida con datos experimentales obtenidos por espectroscopía micro-Raman y medidas eléctricas, respectivamente. Los resultados muestran el stress intrínseco en la capa producido por la distribución no uniforme del 2DEG en el canal y la región de acceso. Además de aumentar la potencia de salida del dispositivo, la deformación intrínseca en la capa de diamante podría mejorar la fiabilidad del dispositivo modulando la deformación local en el borde de la puerta del lado del drenador. Finalmente, también se han simulado en este trabajo los efectos de trampas localizados en la superficie, el buffer y la barrera. Las medidas pulsadas muestran que tanto las puertas largas como las grandes separaciones entre los contactos de puerta y drenador aumentan el cociente entre la corriente pulsada frente a la corriente continua (lag ratio), es decir, disminuir el colapse de corriente (current collapse). Este efecto ha sido explicado mediante las simulaciones de los efectos de trampa de superficie. Por su parte, las referidas a trampas en el buffer se enfocaron en los efectos de atrapamiento dinámico, y su impacto en el auto-calentamiento del dispositivo. Se presenta también un modelo que describe el atrapamiento y liberación de trampas en la barrera: mientras que el atrapamiento se debe a un túnel directo del electrón desde el metal de puerta, el desatrapamiento consiste en la emisión del electrón en la banda de conducción mediante túnel asistido por fonones. El modelo también simula la corriente de puerta, debida a la emisión electrónica dependiente de la temperatura y el campo eléctrico. Además, también se ilustra la corriente de drenador dependiente de la temperatura y el campo eléctrico. ABSTRACT GaN and AlN are group III-V piezoelectric semiconductor materials. The AlGaN/GaN heterojunction presents large piezoelectric and spontaneous polarization charge at the interface, leading to high 2DEG density close to the interface. A high power output would be obtained due to the high 2DEG density and mobility, which leads to elevated lattice temperature. The gate and drain biases induce converse piezoelectric stress that can influence the piezoelectric polarization charge and further influence the 2DEG density and output characteristics. Therefore, the device physics is relevant to all the electrical, thermal, and mechanical aspects. In this dissertation, by using the commercial finite-element-method (FEM) software COMSOL, we achieved the GaN HEMTs simulation with electro-thermal, electro-mechanical, and electro-thermo-mechanical full coupling. The coupling parts include the drift-diffusion model for the electron transport, the thermal conduction, and the piezoelectric effect. By simulations and some experimental characterizations, we have studied the device thermal, stress, and traps effects described in the following. The device geometry impact on the self-heating was studied by electro-thermal simulations and electrical characterizations. Among the obtained interesting results, we found that, for same power output, the distance between the gate and drain contact can influence distribution of the heat generation in the channel and thus influence the channel temperature. Diamond possesses high thermal conductivity. Integrated diamond with the device can spread the generated heat and thus potentially reduce the device self-heating effect. Electro-thermal simulations on this topic were performed. For the diamond integration on top of the device (top-side heat spreading), the determinant factors for the heat spreading ability are the diamond thickness, its thermal conductivity, and its distance to the heat source. The top-side heat spreading can also reduce the impact of thermal boundary resistance between the buffer and the substrate on the device thermal behavior. The very low electrical conductivity of diamond allows that it can directly contact the gate metal (which is very close to the heat source), being quite convenient to reduce the self-heating for the device under pulsed bias. Also, the diamond coated in vias etched in the substrate as heat spreading path (back-side heat spreading) was simulated. A competing mechanism influences the heat spreading ability, i.e., the etched vias would increase the device temperature due to the reduced heat sink while the coated diamond would decrease the device temperature due to its higher thermal conductivity. Therefore, relative thick coated diamond is needed in order to reduce the self-heating effect. The simulated local stress at the gate edge of the drain side for the device with standard and field plate gate structure were compared, which would be relevant to the device mechanical failure. Other stress simulations focused on the intrinsic stress in the diamond capping layer impact on the device electrical behaviors. The simulated stress and electrical output characteristics were compared to experimental data obtained by micro-Raman spectroscopy and electrical characterization, respectively. Results showed that the intrinsic stress in the capping layer caused the non-uniform distribution of 2DEG in the channel and the access region. Besides the enhancement of the device power output, intrinsic stress in the capping layer can potentially improve the device reliability by modulating the local stress at the gate edge of the drain side. Finally, the surface, buffer, and barrier traps effects were simulated in this work. Pulsed measurements showed that long gates and distances between gate and drain contact can increase the gate lag ratio (decrease the current collapse). This was explained by simulations on the surface traps effect. The simulations on buffer traps effects focused on illustrating the dynamic trapping/detrapping in the buffer and the self-heating impact on the device transient drain current. A model was presented to describe the trapping and detrapping in the barrier. The trapping was the electron direct tunneling from the gate metal while the detrapping was the electron emission into the conduction band described by phonon-assisted tunneling. The reverse gate current was simulated based on this model, whose mechanism can be attributed to the temperature and electric field dependent electron emission in the barrier. Furthermore, the mechanism of the device bias via the self-heating and electric field impact on the electron emission and the transient drain current were also illustrated.

Design and construction of precision heat fluxmeters

The possibility of local elastic instabilities is considered in a first¿order structural phase transition, typically a thermoelastic martensitic transformation, with associated interfacial and volumic strain energy. They appear, for instance, as the result of shape change accommodation by simultaneous growth of different crystallographic variants. The treatment is phenomenological and deals with growth in both thermoelastic equilibrium and in nonequilibrium conditions produced by the elastic instability. Scaling of the transformed fraction curves against temperature is predicted only in the case of purely thermoelastic growth. The role of the transformation latent heat on the relaxation kinetics is also considered, and it is shown that it tends to increase the characteristic relaxation times as adiabatic conditions are approached, by keeping the system closer to a constant temperature. The analysis also reveals that the energy dissipated in the relaxation process has a double origin: release of elastic energy Wi and entropy production Si. The latter is shown to depend on both temperature rate and thermal conduction in the system.

Discretized integral hydrodynamics

Using an interpolant form for the gradient of a function of position, we write an integral version of the conservation equations for a fluid. In the appropriate limit, these become the usual conservation laws of mass, momentum, and energy. We also discuss the special cases of the Navier-Stokes equations for viscous flow and the Fourier law for thermal conduction in the presence of hydrodynamic fluctuations. By means of a discretization procedure, we show how the integral equations can give rise to the so-called particle dynamics of smoothed particle hydrodynamics and dissipative particle dynamics.

No Major Differences Found between the Effects of Microwave-Based and Conventional Heat Treatment Methods on Two Different Liquid Foods

Extension of shelf life and preservation of products are both very important for the food industry. However, just as with other processes, speed and higher manufacturing performance are also beneficial. Although microwave heating is utilized in a number of industrial processes, there are many unanswered questions about its effects on foods. Here we analyze whether the effects of microwave heating with continuous flow are equivalent to those of traditional heat transfer methods. In our study, the effects of heating of liquid foods by conventional and continuous flow microwave heating were studied. Among other properties, we compared the stability of the liquid foods between the two heat treatments. Our goal was to determine whether the continuous flow microwave heating and the conventional heating methods have the same effects on the liquid foods, and, therefore, whether microwave heat treatment can effectively replace conventional heat treatments. We have compared the colour, separation phenomena of the samples treated by different methods. For milk, we also monitored the total viable cell count, for orange juice, vitamin C contents in addition to the taste of the product by sensory analysis. The majority of the results indicate that the circulating coil microwave method used here is equivalent to the conventional heating method based on thermal conduction and convection. However, some results in the analysis of the milk samples show clear differences between heat transfer methods. According to our results, the colour parameters (lightness, red-green and blue-yellow values) of the microwave treated samples differed not only from the untreated control, but also from the traditional heat treated samples. The differences are visually undetectable, however, they become evident through analytical measurement with spectrophotometer. This finding suggests that besides thermal effects, microwave-based food treatment can alter product properties in other ways as well.

Aspetti innovativi per la valutazione e la misura della porosità nei materiali dell'edilizia antica e moderna. La geometria frattale della porosità

New concepts on porosity appraisal in ancient and modern construction materials. The role of Fractal Geometry on porosity characterization and transport phenomena. This work studied the potential of Fractal Geometry to the characterization of porous materials. Besides the descriptive aspects of the pore size distribution, the fractal dimensions have led to the development of rational relations for the prediction of permeability coefficients to fluid and heat transfer. The research considered natural materials used in historical buildings (rock and earth) as well as currently employed materials as hydraulic cement and technologically advanced materials such as silicon carbide or YSZ ceramics. The experimental results of porosity derived from the techniques of mercury intrusion and from the image analysis. Data elaboration was carried out according to established procedures of Fractal Geometry. It was found that certain classes of materials are clearly fractal and respond to simple patterns such as Sierpinski and Menger models. In several cases, however, the fractal character is not recognised because the microstructure of the material is based on different phases at different dimensional scales, and in consequence the “fractal dimensions” calculated from porosimetric data do not come within the standard range (less than 3). Using different type and numbers of fractal units is possible, however, to obtain “virtual” microstructures that have the fraction of voids and pore size distribution equivalent with the experimental ones for almost any material. Thus it was possible to take the expressions for the permeability and the thermal conduction which does not require empirical “constants”, these expressions have also provided values that are generally in agreement with the experimental available data. More problematic has been the fractal discussion of the geometry of the rupture of the material subjected to mechanical stress both external and internal applied. The results achieved on these issues are qualitative and prone to future studies. Keywords: Materials, Microstructure, Porosity, Fractal Geometry, Permeability, Thermal conduction, Mechanical strength.

High-Speed Laser Processing of Thin Single and Multi-Layer Films

Theoretical models are developed for the continuous-wave and pulsed laser incision and cut of thin single and multi-layer films. A one-dimensional steady-state model establishes the theoretical foundations of the problem by combining a power-balance integral with heat flow in the direction of laser motion. In this approach, classical modelling methods for laser processing are extended by introducing multi-layer optical absorption and thermal properties. The calculation domain is consequently divided in correspondence with the progressive removal of individual layers. A second, time-domain numerical model for the short-pulse laser ablation of metals accounts for changes in optical and thermal properties during a single laser pulse. With sufficient fluence, the target surface is heated towards its critical temperature and homogeneous boiling or "phase explosion" takes place. Improvements are seen over previous works with the more accurate calculation of optical absorption and shielding of the incident beam by the ablation products. A third, general time-domain numerical laser processing model combines ablation depth and energy absorption data from the short-pulse model with two-dimensional heat flow in an arbitrary multi-layer structure. Layer removal is the result of both progressive short-pulse ablation and classical vaporisation due to long-term heating of the sample. At low velocity, pulsed laser exposure of multi-layer films comprising aluminium-plastic and aluminium-paper are found to be characterised by short-pulse ablation of the metallic layer and vaporisation or degradation of the others due to thermal conduction from the former. At high velocity, all layers of the two films are ultimately removed by vaporisation or degradation as the average beam power is increased to achieve a complete cut. The transition velocity between the two characteristic removal types is shown to be a function of the pulse repetition rate. An experimental investigation validates the simulation results and provides new laser processing data for some typical packaging materials.

Hydrodynamic stability theory of double ablation front structures in inertial confinement fusion.

The aim of inertial confinement fusion is the production of energy by the fusion of thermonuclear fuel (deuterium-tritium) enclosed in a spherical target due to its implosion. In the direct-drive approach, the energy needed to spark fusion reactions is delivered by the irradiation of laser beams that leads to the ablation of the outer shell of the target (the so-called ablator). As a reaction to this ablation process, the target is accelerated inwards, and, provided that this implosion is sufficiently strong a symmetric, the requirements of temperature and pressure in the center of the target are achieved leading to the ignition of the target (fusion). One of the obstacles capable to prevent appropriate target implosions takes place in the ablation region where any perturbation can grow even causing the ablator shell break, due to the ablative Rayleigh-Taylor instability. The ablative Rayleigh-Taylor instability has been extensively studied throughout the last 40 years in the case where the density/temperature profiles in the ablation region present a single front (the ablation front). Single ablation fronts appear when the ablator material has a low atomic number (deuterium/tritium ice, plastic). In this case, the main mechanism of energy transport from the laser energy absorption region (low density plasma) to the ablation region is the electron thermal conduction. However, recently, the use of materials with a moderate atomic number (silica, doped plastic) as ablators, with the aim of reducing the target pre-heating caused by suprathermal electrons generated by the laser-plasma interaction, has demonstrated an ablation region composed of two ablation fronts. This fact appears due to increasing importance of radiative effects in the energy transport. The linear theory describing the Rayleigh-Taylor instability for single ablation fronts cannot be applied for the stability analysis of double ablation front structures. Therefore, the aim of this thesis is to develop, for the first time, a linear stability theory for this type of hydrodynamic structures.

Self-similar motion of laser fusion plasmas. Absorption in an unbounded plasma

The one-dimensional motion generated in a cold, infinite, uniform plasma of density na by the absorption, in a certain plane, of a linear pulse of energy per unit time and area = 4>0t/r, 0< t< r, is considered, the analysis allows for thermal conduction and viscosity of ions and electrons, their energy exchange, and an electron heat flux limiter The resulting motion is self-similar and governed by a single nondimensional parameter a«(n0 2T/0)2/3 Detailed asymptotic results are obtained for both a < l and a > l , the general behavior of the solution for arbitrary a is discussed The analysis can be extended to the case of a plasma initially occupying a half-space, and throws light on how to optimize the hydrodynamics of laser fusion plasmas Known approximate results corresponding to motion of a plasma submitted to constant irradiation (<()) are recovered in the present work under appropriate limiting processes

Nonlocal electron heat-flux

Electron thermal conduction in a not quite collisional unmagnetlzed plasma is analysed. The failure of classical results for temperature scale-length up to 100 times larger than thermal mean-free-path for electron scattering, and large ion-charge number Z , is discussed. Recent results from a nonlocal model of conduction at large Z are reviewed. Closed form expressions for Braginskii's coefficients a ,/3 , y for Z =0(1) are derived. An extension of the nonlocal model for Z =0(1) is discussed.

Joule heating in ferromagnetic nanostripes with a notch

The temperature in a ferromagnetic nanostripe with a notch subject to Joule heating has been studied in detail. We first performed an experimental real-time calibration of the temperature versus time as a 100 ns current pulse was injected into a Permalloy nanostripe. This calibration was repeated for different pulse amplitudes and stripe dimensions and the set of experimental curves were fitted with a computer simulation using the Fourier thermal conduction equation. The best fit of these experimental curves was obtained by including the temperature-dependent behavior of the electrical resistivity of the Permalloy and of the thermal conductivity of thesubstrate(SiO2). Notably, a nonzero interface thermal resistance between the metallic nanostripe and thesubstrate was also necessary to fit the experimental curves. We found this parameter pivotal to understand ourresults and the results from previous works. The higher current density in the notch, together with the interface thermal resistance, allows a considerable increase of the temperature in the notch, creating a large horizontal thermal gradient. This gradient, together with the high temperature in the notch and the larger current density close to the edges of the notch, can be very influential in experiments studying the current assisted domain wall motion.

Contêineres metálicos para canteiros de obras: análise experimental de desempenho térmico e melhorias na transferência de calor pela envoltória.

Os contêineres metálicos foram desenvolvidos para a utilização no setor de logística e transporte, mas por sua escala adaptável à das edificações e pela mobilidade e praticidade de instalação, tiveram sua utilização apropriada também pelo setor da construção civil. Essas instalações possuem diversas qualidades ambientalmente amigáveis, mas seu aspecto térmico é extremamente insuficiente: sem isolamento térmico, demandam alta carga térmica de refrigeração e aquecimento, no verão e inverno, respectivamente e, consequentemente, um alto consumo energético. Tal característica foi crucial para que se determinassem como objetivos da presente pesquisa investigar o comportamento térmico dessas construções metálicas, avaliar seus parâmetros de desempenho, conforto e estresse térmicos, por meio de uma ampla coleta de dados experimentais. O experimento com duração de um ano - contou com três tipologias de contêiner em escala real, sendo o primeiro em aço Tipo X sem isolamento térmico, o segundo com um isolamento térmico para o fenômeno da condução e o terceiro com isolamento térmico para o fenômeno da radiação. Os diferentes tipos de tratamentos térmicos proporcionaram melhorias à envoltória dos contêineres, chegando a uma diferença nas temperaturas internas de até 9 °C. Constatou-se a extrema necessidade de adequação do tipo de isolamento térmico dos contêineres ao uso a que tais instalações se destinam escritório ou alojamento, no caso dos canteiros de obras para que as características da envoltória minimizem de fato a demanda ou mesmo atinjam a eliminação da necessidade de condicionamento artificial.

Controlling Nanostructures for in-situ TEM Characterization

Low dimensional nanostructures, such as nanotubes and 2D sheets, have unique and promising material properties both from a fundamental science and an application standpoint. Theoretical modelling and calculations predict previously unobserved phenomena that experimental scientists often struggle to reproduce because of the difficulty in controlling and characterizing the small structures under real-world constraints. The goal of this dissertation is to controlling these structures so that nanostructures can be characterized in-situ in transmission electron microscopes (TEM) allowing for direct observation of the actual physical responses of the materials to different stimuli. Of most interest to this work are the thermal and electrical properties of carbon nanotubes, boron nitride nanotubes, and graphene. The first topic of the dissertation is using surfactants for aqueous processing to fabricate, store, and deposit the nanostructures. More specifically, thorough characterization of a new surfactant, ammonium laurate (AL), is provided and shows that this new surfactant outperforms the standard surfactant for these materials, sodium dodecyl sulfate (SDS), in almost all tested metrics. New experimental set-ups have been developed by combining specialized in-situ TEM holders with innovative device fabrication. For example, electrical characterization of graphene was performed by using an STM-TEM holder and depositing graphene from aqueous solutions onto lithographically patterned, electron transparent silicon nitride membranes. These experiments produce exciting information about the interaction between graphene and metal probes and the substrate that it rests on. Then, by adding indium to the backside of the membrane and employing the electron thermal microscopy (EThM) technique, the same type of graphene samples could be characterized for thermal transport with high spatial resolution. It is found that reduced graphene oxide sheets deposited onto a silicon nitride membrane and displaying high levels of wrinkling have higher than expected electrical and thermal conduction properties. We are clearly able to visualize the ability of graphene to spread heat away from an electronic hot spot and into the substrate.

Examining the growth and evolution of magnetic fields in clusters of galaxies: a numerical perspective

Magnetic fields are ubiquitous in galaxy cluster atmospheres and have a variety of astrophysical and cosmological consequences. Magnetic fields can contribute to the pressure support of clusters, affect thermal conduction, and modify the evolution of bubbles driven by active galactic nuclei. However, we currently do not fully understand the origin and evolution of these fields throughout cosmic time. Furthermore, we do not have a general understanding of the relationship between magnetic field strength and topology and other cluster properties, such as mass and X-ray luminosity. We can now begin to answer some of these questions using large-scale cosmological magnetohydrodynamic (MHD) simulations of the formation of galaxy clusters including the seeding and growth of magnetic fields. Using large-scale cosmological simulations with the FLASH code combined with a simplified model of the acceleration of cosmic rays responsible for the generation of radio halos, we find that the galaxy cluster frequency distribution and expected number counts of radio halos from upcoming low-frequency sur- veys are strongly dependent on the strength of magnetic fields. Thus, a more complete understanding of the origin and evolution of magnetic fields is necessary to understand and constrain models of diffuse synchrotron emission from clusters. One favored model for generating magnetic fields is through the amplification of weak seed fields in active galactic nuclei (AGN) accretion disks and their subsequent injection into cluster atmospheres via AGN-driven jets and bubbles. However, current large-scale cosmological simulations cannot directly include the physical processes associated with the accretion and feedback processes of AGN or the seeding and merging of the associated SMBHs. Thus, we must include these effects as subgrid models. In order to carefully study the growth of magnetic fields in clusters via AGN-driven outflows, we present a systematic study of SMBH and AGN subgrid models. Using dark-matter only cosmological simulations, we find that many important quantities, such as the relationship between SMBH mass and galactic bulge velocity dispersion and the merger rate of black holes, are highly sensitive to the subgrid model assumptions of SMBHs. In addition, using MHD calculations of an isolated cluster, we find that magnetic field strengths, extent, topology, and relationship to other gas quantities such as temperature and density are also highly dependent on the chosen model of accretion and feedback. We use these systematic studies of SMBHs and AGN inform and constrain our choice of subgrid models, and we use those results to outline a fully cosmological MHD simulation to study the injection and growth of magnetic fields in clusters of galaxies. This simulation will be the first to study the birth and evolution of magnetic fields using a fully closed accretion-feedback cycle, with as few assumptions as possible and a clearer understanding of the effects of the various parameter choices.