Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “smoothing”

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 1,441 records · Page 80

Machine Learning-Assisted Recovery of Delicate Kinetic Information from Transient Reactor Experiments

Identifying active sites and their roles in chemical reaction steps remains a vital challenge in heterogeneous catalysis. Transient experiments offer a unique way to probe active sites and distinguish subtle kinetic features. Although physics-based analysis methods may be well-developed, they can be highly susceptible to experimental noise, and smoothing methods may erase or even distort important features; a smooth curve is not always the best curve. We demonstrate a new workflow for the direct interpretation of intrinsic kinetic information from exit flux curves measured in transient reactor experiments. This workflow contains three artificial neural networks (ANNs), including a noise reducer, a concentration predictor, and a rate predictor to analyze experimental data, followed by the virtual TAP (VTAP) physics-based reactor model and density functional theory (DFT) calculations of adsorption energies on specific sites. We use this workflow to analyze the data from experiments titrating Pt/Al 2 O 3 and Pt/SiO 2 catalysts with carbon monoxide (CO) in the temporal analysis of products (TAP) reactor. Our workflow separates the time-evolving chemical reaction and mass transfer information contained in the TAP pulse response. The existence of strong- and weak-binding sites on the Pt/Al 2 O 3 catalyst is observed in the catalyst titration experiment in the transient reactor. The structures of the strong- and weak-binding sites are then identified by using DFT calculations. We find that the Pt/SiO 2 catalyst has only strong-binding sites, which aligns with the inactive support effect of SiO 2 . We demonstrate how machine learning methods provide unique insights with high-resolution data analysis that cannot be achieved by using state-of-the-art physics-based methods.

Adsorption↗

Superposition Principle in Auger Recombination of Charged and Neutral Multicarrier States in Semiconductor Quantum Dots

Application of colloidal semiconductor quantum dots (QDs) in optical and optoelectronic devices is often complicated by unintentional generation of extra charges, which opens fast nonradiative Auger recombination pathways whereby the recombination energy of an exciton is quickly transferred to the extra carrier(s) and ultimately dissipated as heat. Previous studies of Auger recombination have primarily focused on neutral and, more recently, negatively charged multicarrier states. Auger dynamics of positively charged species remains more poorly explored due to difficulties in creating, stabilizing, and detecting excess holes in the QDs. Here we apply photochemical doping to prepare both negatively and positively charged CdSe/CdS QDs with two distinct core/shell interfacial profiles (“sharp” versus “smooth”). Using neutral and charged QD samples we evaluate Auger lifetimes of biexcitons, negative and positive trions (an exciton with an extra electron or a hole, respectively), and multiply negatively charged excitons. Using these measurements, we demonstrate that Auger decay of both neutral and charged multicarrier states can be presented as a superposition of independent elementary three-particle Auger events. As one of the manifestations of the superposition principle, we observe that the biexciton Auger decay rate can be presented as a sum of the Auger rates for independent negative and positive trion pathways. Furthermore, by comparing the measurements on the QDs with the “sharp” versus “smooth” interfaces, we also find that while affecting the absolute values of Auger lifetimes, manipulation of the shape of the confinement potential does not lead to violation of the superposition principle, which still allows us to accurately predict the biexciton Auger lifetimes based on the measured negative and positive trion dynamics. Our findings indicate considerable robustness of the superposition principle as applied to Auger decay of charged and neutral multicarrier states, suggesting its generality to quantum-confined nanocrystals of arbitrary compositions and complexities.

36 MATERIALS SCIENCE↗

Understanding Multiscale Structural Changes During Dilute Acid Pretreatment of Switchgrass and Poplar

Biofuels produced from lignocellulosic biomass hold great promise as a renewable alternative energy and fuel source. To realize a cost and energy efficient approach, a fundamental understanding of the deconstruction process is critically necessary to reduce biomass recalcitrance. Here, the structural and morphological changes over multiple scales (5–6000 Å) in herbaceous (switchgrass) and woody (hybrid poplar) biomass during dilute sulfuric acid pretreatment were explored using neutron scattering and X-ray diffraction. Switchgrass undergoes a larger increase (20–84 Å) in the average diameter of the crystalline core of the elementary cellulose fibril than hybrid poplar (19–50 Å). Switchgrass initially forms lignin aggregates with an average size of 90 Å that coalesce to 200 Å, which is double that observed for hybrid poplar, 55–130 Å. Switchgrass shows a smooth-to-rough transition in the cell wall surface morphology unlike the diffuse-to-smooth transition of hybrid poplar. Yet, switchgrass and hybrid poplar pretreated under the same experimental conditions result in pretreated switchgrass producing higher glucose yields (~76 wt %) than pretreated hybrid poplar (~60 wt %). This observation shows that other aspects like cellulose allomorph transitions, cellulose accessibility, cellular biopolymer spatial distribution, and enzyme–substrate interactions may be more critical in governing the enzymatic hydrolysis efficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nematic director reorientation at solid and liquid interfaces under flow: SAXS studies in a microfluidic device

In this work we investigate the interplay between flow and boundary condition effects on the orientation field of a thermotropic nematic liquid crystal under flow and confinement in a microfluidic device. Two types of experiments were performed using synchrotron small-angle X-ray-scattering (SAXS). In the first, a nematic liquid crystal flows through a square-channel cross section at varying flow rates, while the nematic director orientation projected onto the velocity/velocity gradient plane is measured using a 2D detector. At moderate-to-high flow rates, the nematic director is predominantly aligned in the flow direction, but with a small tilt angle of ~±11° in the velocity gradient direction. The director tilt angle is constant throughout most of the channel width but switches sign when crossing the center of the channel, in agreement with the Ericksen–Leslie–Parodi (ELP) theory. At low flow rates, boundary conditions begin to dominate, and a flow profile resembling the escaped radial director configuration is observed, where the director is seen to vary more smoothly from the edges (with homeotropic alignment) to the center of the channel. In the second experiment, hydrodynamic focusing is employed to confine the nematic phase into a sheet of liquid sandwiched between two layers of Triton X-100 aqueous solutions. The average nematic director orientation shifts to some extent from the flow direction toward the liquid boundaries, although it remains unclear if one tilt angle is dominant through most of the nematic sheet (with abrupt jumps near the boundaries) or if the tilt angle varies smoothly between two extreme values (~90 and 0°). Lastly, the technique presented here could be applied to perform high-throughput measurements for assessing the influence of different surfactants on the orientation of nematic phases and may lead to further improvements in areas such as boundary lubrication and clarifying the nature of defect structures in LC displays.

36 MATERIALS SCIENCE↗

Crust and Upper Mantle Structure Beneath the Eastern United States

Abstract The Eastern United States (EUS) has a complex geological history and hosts several seismic active regions. We investigate the subsurface structure beneath the broader EUS. To produce reliable images of the subsurface, we simultaneously invert smoothed P‐wave receiver functions, Rayleigh‐wave phase and group velocity measurements, and Bouguer gravity observations for the 3D shear‐wave speed. Using surface‐wave observations (3–250 s) and spatially smoothed receiver functions, our velocity models are robust, reliable, and rich in detail. The shear‐wave velocity models fit all three types of observations well. The resulting velocity model for the eastern U.S. shows thinner crust beneath New England, the east coast, and the Mississippi Embayment (ME). A relatively thicker crust was found beneath the stable North America craton. A relatively slower upper mantle was imaged beneath New England, the east coast, and western ME. A comparison of crust thickness derived from our model against four recent published models shows first‐order consistency. A relatively small upper mantle low‐speed region correlates with a published P‐wave analysis that has associated the anomaly with a 75 Ma kimberlite volcanic site in Kentucky. We also explored the relationship between the subsurface structure and seismicity in the eastern U.S. We found that earthquakes often locate near regions with seismic velocity variations, but not universally. Not all regions of significant subsurface wave speed changes are loci of seismicity. A weak correlation between upper mantle shear velocity and earthquake focal mechanism has been observed.

58 GEOSCIENCES↗

Spontaneous formation of nanoparticles on electrospun nanofibres

We report the spontaneous formation of nanoparticles on smooth nanofibres in a single-step electrospinning process, as an inexpensive and scalable method for producing high-surface-area composites. Layers of nanofibres, containing the proton conducting electrolyte, caesium dihydrogen phosphate, are deposited uniformly over large area substrates from clear solutions of the electrolyte mixed with polymers. Under certain conditions, the normally smooth nanofibres develop caesium dihydrogen phosphate nanoparticles in large numbers on their external surface. The nanoparticles appear to originate from the electrolyte within the fibres, which is transported to the outer surface after the fibres are deposited, as evidenced by cross-sectional imaging of the electrospun fibres. The presence of nanoparticles on the fibre surface yields composites with increased surface area of exposed electrolyte, which ultimately enhances electrocatalytic performance. Indeed, solid acid fuel cells fabricated with electrodes from processed nanofibre-nanoparticle composites, produced higher cell voltage as compared to fuel cells fabricated with state-of-the-art electrodes.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Mild-to-wild plastic transition is governed by athermal screw dislocation slip in bcc Nb

Plastic deformation in crystals is mediated by the motion of line defects known as dislocations. For decades, dislocation activity has been treated as a homogeneous, smooth continuous process. However, it is now recognized that plasticity can be determined by long-range correlated and intermittent collective dislocation processes, known as avalanches. Here we demonstrate in body-centered cubic Nb how the long-range and scale-free dynamics at room temperature are progressively quenched out with decreasing temperature, eventually revealing intermittency with a characteristic length scale that approaches the Burgers vector itself. Plasticity is shown to be bimodal across the studied temperature regime, with conventional thermally-activated smooth plastic flow (‘mild’) coexisting with sporadic bursts (‘wild’) controlled by athermal screw dislocation activity, thereby violating the classical notion of temperature-dependent screw dislocation motion at low temperatures. An abrupt increase of the athermal avalanche component is identified at the critical temperature of the material. Our results indicate that plasticity at any scale can be understood in terms of the coexistence of these mild and wild modes of deformation, which could help design better alloys by suppressing one of the two modes in desired temperature windows.

36 MATERIALS SCIENCE↗

Electrochemical properties and structural evolution of O3-type layered sodium mixed transition metal oxides with trivalent nickel

The electrochemical properties of NaNi 0.5 Co 0.5 O 2 and NaNi 0.5 Fe 0.5 O 2 and their structural transitions as a function of Na extraction associated with redox reactions are investigated. Synthesized in the O3-type layered structure, both materials show reasonable electrochemical activities at room temperature, delivering approximately 0.5 Na per formula unit at C/10 discharge. More Na can be reversibly cycled in NaNi 0.5 Co 0.5 O 2 at elevated temperature and/or in an extended voltage window, while NaNi 0.5 Fe 0.5 O 2 shows significant capacity fading at a high voltage cutoff which is likely due to Fe 4+ migration. In situ X-ray diffraction shows that the structural changes in the two materials upon desodiation are very different. NaNi 0.5 Co 0.5 O 2 goes through many different two-phase reactions including three different O3-type and three different P3-type structures during cycling, producing a voltage profile with multiple plateau-like features. In contrast, NaNi 0.5 Fe 0.5 O 2 has a smooth voltage profile and shows the typical O3–P3 phase transition without lattice distortion seen in other materials. Finally, this different structural evolution upon desodiation and re-sodiation can be explained by the electronic structure of the mixed transition metals and how it perturbs the ordering between Na ions differently.

36 MATERIALS SCIENCE↗

Local spring effect in titanium-based layered oxides

Titanium-based layered oxides (TLOs) are one of the most promising electrode material families for sodium-ion batteries (NIBs) due to their smooth charge/discharge profiles and excellent cycle performance. However, the reaction mechanism of these materials, especially the reason for the disappearance of multiple voltage plateaus, is still not clear. Herein, two representative TLOs (Na 2/3 Ni 1/3 Ti 2/3 O 2 and Na 2/3 Co 1/3 Ti 2/3 O 2 ) with the same P2 crystal structure have been studied to scrutinize those unexplained issues. In situ synchrotron high-energy X-ray diffraction revealed a solid solution reaction mechanism for both, suggesting the absence of rigid phase transitions upon electrochemical cycling. An interesting “spring effect” of the TiO 6 octahedron, i.e., the reversible vibration of the central Ti atom inside the local octahedron upon electrochemical redox, was demonstrated by advanced X-ray absorption spectroscopy and theoretical calculations. Such an effect could suppress the rigid phase transitions, and result in smooth charge/discharge profiles and enhanced cycle stability. Finally, this work not only accounts for the disappearance of multiple voltage plateaus of TLOs for NIBs, but also provides an effective local-structure viewpoint to increase the cycle stability of electrode materials for other advanced battery systems.

25 ENERGY STORAGE↗

Computational general relativistic force-free electrodynamics: II. Characterization of numerical diffusivity

Scientific codes are an indispensable link between theory and experiment; in (astro-)plasma physics, such numerical tools are one window into the universe’s most extreme flows of energy. The discretization of Maxwell’s equations – needed to make highly magnetized (astro)physical plasma amenable to its numerical modeling – introduces numerical diffusion. It acts as a source of dissipation independent of the system’s physical constituents. Understanding the numerical diffusion of scientific codes is the key to classifying their reliability. It gives specific limits in which the results of numerical experiments are physical. We aim at quantifying and characterizing the numerical diffusion properties of our recently developed numerical tool for the simulation of general relativistic force-free electrodynamics by calibrating and comparing it with other strategies found in the literature. Our code correctly models smooth waves of highly magnetized plasma. We evaluate the limits of general relativistic force-free electrodynamics in the context of current sheets and tearing mode instabilities. We identify that the current parallel to the magnetic field (j ∥ ), in combination with the breakdown of general relativistic force-free electrodynamics across current sheets, impairs the physical modeling of resistive instabilities. We find that at least eight numerical cells per characteristic size of interest (e.g., the wavelength in plasma waves or the transverse width of a current sheet) are needed to find consistency between resistivity of numerical and of physical origins. High-order discretization of the force-free current allows us to provide almost ideal orders of convergence for (smooth) plasma wave dynamics. The physical modeling of resistive layers requires suitable current prescriptions or a sub-grid modeling for the evolution of j ∥ .

79 ASTRONOMY AND ASTROPHYSICS↗

Anomalies in the topology of the temperature fluctuations in the cosmic microwave background: An analysis of the NPIPE and FFP10 data releases

We present a topological analysis of the temperature fluctuation maps from the Planck 2020 Data Release 4 NPIPE dataset and the Planck 2018 Data Release 3 FFP10 dataset. We performed a multiscale analysis in terms of the homology characteristics of the maps, invoking relative homology to account for the analysis in the presence of masks. We performed our analysis for a range of smoothing scales spanning sub- and super-horizon scales corresponding to a full width at half maximum (FWHM) of 5',10',20',40',80',160',320', and 640', and employed simulations based on the standard model for comparison, which assumes the initial fluctuation field to be an isotropic and homogeneous Gaussian random field. Examining the behavior of topological components, represented by the 0D homology group, we find the observations to be approximately 2σ or less deviant from the simulations for all resolutions and scales for the NPIPE dataset. For the FFP10 dataset, we detect a 2.96σ deviation between the observations and simulations at N = 128, FWHM = 80'. For the topological loops, represented by the first homology group, the simulations and observations are consistent within 2σ for most resolutions and scales for both the datasets. However, for the NPIPE dataset, we observe a high deviation between the observation and simulations in the number of loops at FWHM = 320', but at a low dimensionless threshold ν = –2.5. Under a Gaussian assumption, this would amount to a deviation of ~4σ. However, the distribution in this bin is manifestly non-Gaussian and does not obey Poisson statistics either. In the absence of a true theoretical understanding, we simply note that the significance is higher than what may be resolved by 600 simulations, yielding an empirical p-value of at most 0.0016. Specifically in this case, our tests indicate that the numbers arise from a statistically stable regime, despite being based on small numbers. For the FFP10 dataset, the differences are not as strong as for the NPIPE dataset, indicating a 2.77σ deviation at this resolution and threshold. The Euler characteristic, which is the alternating sum of the ranks of relative homology groups, reflects the deviations in the components and loops. To assess the significance of combined levels for a given scale, we employed the empirical and theoretical versions of the χ 2 test as well as the nonparametric Tukey depth test. Although all statistics exhibit a stable distribution, we favor the empirical version of the χ 2 test in the final interpretation, as it indicates the most conservative differences. For the NPIPE dataset, we find that the components and loops differ at more than 95%, but agree within the 99% confidence level with respect to the base model at N = 32, FWHM = 320'. The Euler characteristic at this resolution displays a per mil deviation. In contrast, the FFP10 dataset shows that the observations are consistent with the base model within the 95% confidence level, at this and smaller scales. This is consistent with the observations of the Planck analysis pipeline via Minkowski functionals. For the largest smoothing scale, N = 16, FWHM = 640', both datasets exhibit an anomalous behavior of the loops, where FFP10 data exhibit a deviation that is larger by an order of magnitude than that of the NPIPE dataset. In contrast, the values for the topological components and the Euler characteristic agree between observations and model to within a confidence level of 99%. However, for the largest scales, the statistics are based on low numbers and may have to be regarded with caution. Even though both datasets exhibit mild to significant discrepancies, they also exhibit contrasting behaviors at various instances. Therefore, we do not find it feasible to convincingly accept or reject the null hypothesis. Disregarding the large-scale anomalies that persist at similar scales in WMAP and Planck, observations of the cosmic microwave background are largely consistent with the standard cosmological model within 2σ.

79 ASTRONOMY AND ASTROPHYSICS↗

Image segmentation for analyzing galaxy-galaxy strong lensing systems

The goal of this Letter is to develop a machine learning model to analyze the main gravitational lens and detect dark substructure (subhalos) within simulated images of strongly lensed galaxies. Using the technique of image segmentation, we turn the task of identifying subhalos into a classification problem, where we label each pixel in an image as coming from the main lens, a subhalo within a binned mass range, or neither. Our network is only trained on images with a single smooth lens and either zero or one subhalo near the Einstein ring. On an independent test set with lenses with large ellipticities, quadrupole and octopole moments, and for source apparent magnitudes between 17–25, the area of the main lens is recovered accurately. On average, only 1.3% of the true area is missed and 1.2% of the true area is added to another part of the lens. In addition, subhalos as light as 10 8.5 M ⊙ can be detected if they lie in bright pixels along the Einstein ring. Furthermore, the model is able to generalize to new contexts it has not been trained on, such as locating multiple subhalos with varying masses or more than one large smooth lens.

79 ASTRONOMY AND ASTROPHYSICS↗

BeyondPlanck: VII. Bayesian estimation of gain and absolute calibration for cosmic microwave background experiments

We present a Bayesian calibration algorithm for cosmic microwave background (CMB) observations as implemented within the global end-to-end BEYONDPLANCK framework and applied to the Planck Low Frequency Instrument (LFI) data. Following the most recent Planck analysis, we decomposed the full time-dependent gain into a sum of three nearly orthogonal components: one absolute calibration term, common to all detectors, one time-independent term that can vary between detectors, and one time-dependent component that was allowed to vary between one-hour pointing periods. Each term was then sampled conditionally on all other parameters in the global signal model through Gibbs sampling. The absolute calibration is sampled using only the orbital dipole as a reference source, while the two relative gain components were sampled using the full sky signal, including the orbital and Solar CMB dipoles, CMB fluctuations, and foreground contributions. We discuss various aspects of the data that influence gain estimation, including the dipole-polarization quadrupole degeneracy and processing masks. Comparing our solution to previous pipelines, we find good agreement in general, with relative deviations of -0.67% (-0.84%) for 30 GHz, 0.12% (-0.04%) for 44 GHz and -0.03% (-0.64%) for 70 GHz, compared to Planck PR4 and Planck 2018, respectively. We note that the BEYONDPLANCK calibration was performed globally, which results in better inter-frequency consistency than previous estimates. Additionally, WMAP observations were used actively in the BEYONDPLANCK analysis, which both breaks internal degeneracies in the Planck data set and results in an overall better agreement with WMAP. Finally, we used a Wiener filtering approach to smoothing the gain estimates. We show that this method avoids artifacts in the correlated noise maps as a result of oversmoothing the gain solution, which is difficult to avoid with methods like boxcar smoothing, as Wiener filtering by construction maintains a balance between data fidelity and prior knowledge. Although our presentation and algorithm are currently oriented toward LFI processing, the general procedure is fully generalizable to other experiments, as long as the Solar dipole signal is available to be used for calibration.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluation of Light-Element Reactions in the Resolved Resonance Region

Light-element reactions at low energies in the resolved resonance region are important for a range of applications in basic and applied sciences including nuclear reactors, nonproliferation, cultural heritage, forensics and environmental control, rare event investigations and nuclear astrophysics. In this paper, we report on an effort to evaluate charged-particle cross sections in the resolved resonance region and produce evaluated nuclear data files for further processing and inclusion in evaluated data libraries. We discuss the open issues in R-matrix calculations as we extend to higher energies, such as dealing with the rapidly growing number of open channels and merging with the regime of smooth cross sections described by the statistical model, and present attempts to address these issues in neutron-induced reactions relevant to nuclear reactor applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Ferromagnetic resonance of sputtered yttrium iron garnet nanometer films

Growth of nm-thick yttrium iron garnet (YIG) films by sputtering and ferromagnetic resonance (FMR) properties in the films were studied. The FMR linewidth of the YIG film decreased as the film thickness was increased from several nanometers to about 100 nm. For films with very smooth surfaces, the linewidth increased linearly with frequency. In contrast, for films with big grains on the surface, the linewidth-frequency response was strongly nonlinear. Films in the 7–26 nm thickness range showed a surface roughness between 0.1 nm and 0.4 nm, a 9.48-GHz FMR linewidth in the 6–10 Oe range, and a damping constant of about 0.001.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Short initial length quench on CICC of ITER TF coils

Previous quench studies performed for the International Thermonuclear Experimental Reactor (ITER) Toroidal Field (TF) Coils have led to identify two extreme families of quench: first 'severe' quenches over long initial lengths in high magnetic field, and second smooth quenches over short initial lengths in low field region. Detailed analyses and results on smooth quench propagation and detectability on one TF Cable In Conduit Conductor (CICC) with a lower propagation velocity are presented here. The influence of the initial quench energy is shown and results of computations with either a Fast Discharge (FD) of the magnet or without (failure of the voltage quench detection system) are reported. The influence of the central spiral of the conductor on the propagation velocity is also detailed. In the cases of a regularly triggered FD, the hot spot temperature criterion of 150 K (with helium and jacket) is fulfilled for an initial quench length of 1 m, whereas this criterion is exceed (Tmax ≈ 200 K) for an extremely short length of 5 cm. These analyses were carried out using both the Supermagnet(trade mark, serif) and Venecia codes and the comparisons of the results are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Experimental validation of advanced regulations for superconducting magnet cooling undergoing periodic heat loads

Superconducting magnets used in tokamaks undergo periodic heat load caused by cycling plasma operations inducing AC losses, neutrons fluxes and eddy currents in magnet structures. In the cryogenic system of JT60-SA tokamak, the Auxiliary Cold Box (ACB) distributes helium from the refrigerator to the cryogenic users and in particular to the superconducting magnets. ACB comprises a saturated helium bath with immersed heat exchangers, extracting heat from independent cooling loops. The supercritical helium flow in each cooling loop is driven by a cold circulator. In order to safely operate the refrigerator during plasma pulses, the interface between the ACB and the refrigerator shall be as stable as possible, with well-balanced bath inlet and outlet mass flows during cycling operation. The solution presented in this paper relies on a combination of regulations to smooth pulsed heat loads and to keep a constant refrigeration power during all the cycle. Two smoothing strategies are presented, both regulating the outlet mass flow of the bath: the first one using the bath as a thermal buffer and the second one storing energy in the loop by varying the cold circulator speed. The bath outlet mass flow is also controlled by an immersed resistive heater which enables a constant evaporation rate in the bath when power coming from the loops is decreasing. The refrigeration power is controlled so that the compensating power remains within an acceptable margin. Experimental validation is achieved using the HELIOS facility. This facility running at CEA Grenoble since 2010 is a scaled down model of the ACB bath and Central Solenoid magnet cooling loop of the JT60-SA tokamak. Test results show performances and robustness of the regulations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Degenerate Ising model for atomistic simulation of crystal-melt interfaces

One of the simplest microscopic models for a thermally driven first-order phase transition is an Ising-type lattice system with nearest-neighbour interactions, an external field, and a degeneracy parameter. The underlying lattice and the interaction coupling constant control the anisotropic energy of the phase boundary, the field strength represents the bulk latent heat, and the degeneracy quantifies the difference in communal entropy between the two phases. We simulate the (stochastic) evolution of this minimal model by applying rejection-free canonical and microcanonical Monte Carlo algorithms, and we obtain caloric curves and heat capacity plots for square (2D) and face-centred cubic (3D) lattices with periodic boundary conditions. Since the model admits precise adjustment of bulk latent heat and communal entropy, neither of which affect the interface properties, we are able to tune the crystal nucleation barriers at a fixed degree of undercooling and verify a dimension-dependent scaling expected from classical nucleation theory. We also analyse the equilibrium crystal-melt coexistence in the microcanonical ensemble, where we detect negative heat capacities and find that this phenomenon is more pronounced when the interface is the dominant contributor to the total entropy. The negative branch of the heat capacity appears smooth only when the equilibrium interface-area-to-volume ratio is not constant but varies smoothly with the excitation energy. Finally, we simulate microcanonical crystal nucleation and subsequent relaxation to an equilibrium Wulff shape, demonstrating the model's utility in tracking crystal-melt interfaces at the atomistic level.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗