Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Flux pinning”

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 181 records · Page 10

Noble gas constraints on spent fuel irradiation histories

Fission gas isotopic compositions are sensitive to a variety of reactor operating parameters that include the neutron flux, neutron fluence at different neutron energies, and operating temperature. Measurements of fission gas isotopic compositions thus have potential to constrain reactor simulations for nuclear forensics and safeguards applications. In this paper, we present Kr and Xe isotope measurements from a suite of samples obtained from locations that span the axial length a fuel pin with a well-characterized irradiation history and compare these data to spatially resolved reactor simulations. Here, we observed positive correlations between fluence sensitive isotopic ratios and burnup and between a flux sensitive ratio and power, although some discrepancies are observed between the measured data and model predictions. These differences may be due to simplifications in the model and/or inaccuracies in the cross sections. A much broader measurement to model comparison is required to better understand the discrepancies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

ML-Based Pebble Power Reconstruction for Pebble Bed Reactor Analysis

Pebble power reconstruction has been explored to complement the conventional homogenized modeling approach in pebble bed reactor (PBR) analysis, as detailed heterogeneous geometry calculations are computationally expensive. The random distribution of pebble fuels within the core challenges the application of conventional pin power reconstruction methods. To address this, we introduce a machine learning approach based on the transformer model, composed of encoder and decoder layers, to estimate the flux and power form functions for reconstructing individual pebble neutron fluxes and powers. The homogeneous neutron flux distribution within each spectral zone (SZ) is obtained from finite element solutions of global diffusion or transport calculations. Verification tests demonstrate that the trained transformer model accurately predicts power form functions over a range of conditions, including variations in pebble enrichment, location, type, SZ size, and burnup. In particular, verification using a three-dimensional PBR benchmark with burned pebbles shows good agreement in heterogeneous pebble power distributions between Griffin and Serpent. These results highlight the potential of applying conventional pin power reconstruction approaches to PBR cores with randomly distributed pebbles.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The importance of explicitly representing the streambed in watershed models

Abstract The streambed is the critical interface between the aquatic and terrestrial systems and hosts important biogeochemical hot spots within river corridors. Although the streambed characteristics are significantly different from those of its surrounding soil, the streambed itself has not been explicitly represented in watershed models. Here, we explicitly incorporated a streambed layer into an integrated hydrologic model through model parameterization and discretization. We examined the hydrological effects of streambed characteristics, including hydraulic conductivity ( K ), layer thickness, and resolution, on the exchange fluxes across the streambed as well as the streamflow at the watershed outlet. The numerical experiments were performed in the American River Watershed, a headwater, mountainous watershed within the Yakima River Basin in central Washington. Despite having a negligible effect on the watershed streamflow, an explicit representation of the streambed with distinctive properties dramatically changed the magnitude and variability of the exchange flux. In general, a larger streambed K along with a thicker streambed layer induced larger exchange fluxes. The exchange flux was most sensitive to the streambed resolution. A finer streambed resolution increased exchange fluxes per unit area while reducing the overall exchange volumes across the entire streambed. The amount of baseflow decreased by 6% as the streambed resolution increased from 250 to 50 m. This finding is important because these hydrological changes may, in turn, affect the exchange of nutrients and contaminants between surface water and groundwater and the associated biogeochemical processes. Our work demonstrated the importance of representing streambeds in fully distributed, process‐based watershed models to better capture the exchange flow dynamics in river corridors.

54 ENVIRONMENTAL SCIENCES↗

Investigation of hot-spots due to trapped flux in niobium superconducting radiofrequency cavities

One of the significant sources of residual losses in superconducting radio-frequency cavities is magnetic flux trapped during the cool-down due to the incomplete Messier effect. If the trapped vortices are non-uniformly distributed on the cavity surface, the temperature mapping revealed the “hotspots” at the location of high density of pinned vortices. Here, we performed a rf test on 1.3 GHz single cell cavity with the combination of the temperature mapping system. The temperature mapping reveled the development of the hot spots with the increase in rf field inside the cavity. When magnetic field is trapped locally on the surface of cavity, the hot-spots strength increase rapidly, showing the direct correlation of vortex induced hot spot and corresponding rf loss.

Accelerator Physics↗

Preliminary Evaluation of Driver Fuel Plates and Reduced Sodium Temperature to Increase Fast Test Reactor Experimental Neutron Irradiation Capabilities

All sodium-cooled fast test reactors operated so far have used essentially the same fuel geometrical ar-rangement, namely cylindrical fuel pins assembled into a triangular array in hexagonal fuel assemblies. Fuel compositions and dimensions may vary but the overall geometrical pattern remains identical. Fur-thermore, most sodium-cooled reactors have operated, and still operate, with core-average inlet and outlet sodium temperatures of, respectively, 350-400°C and 500-550°C. The present preliminary study revisits the standard fuel assembly geometry as well as operating tempera-tures to evaluate potential impact on fast test reactor neutron irradiation capabilities. The reference fuel is U-20Pu-10Zr. Pin (reference) and plate geometries are considered together with two sets of core-average sodium inlet/outlet temperatures: 350/500°C (reference) and 200/400°C. More prototypic environments could be obtained inside ad-hoc irradiation testing vehicles that are independent from the reactor primary coolant. Everything else being the same, fuel temperature is significantly lower in plates than in pins. Therefore, plates could in principle accommodate higher plutonium content, thus, reducing the need for uranium enrichment or allowing the use of lower quality plutonium (U-xPu-10Zr thermal conductivity and solidus temperature decrease as x increases, hence, necessitating additional thermal margins). For the reference pin configuration, lowering the inlet/outlet sodium temperatures from the reference 350/500°C down to 200/400°C provides additional thermal margins that can be used to increase the peak fast flux from about 4.5 × 10 15 n/cm 2 -s to 6 ×10 15 n/cm 2 -s for the same core power of 300 MW. Assuming 300 Equivalent Fuel Power Days (EFPD) of operation per calendar year, a (steel) test article could accu-mulate up to 75 dpa/year. The use of fuel plates provides even more thermal margins which may allow the peak fast flux to reach values as high as 8 × 10 15 n/cm 2 -s for both sets of inlet/outlet temperatures. In this environment, a test article could accumulate up to 95 dpa/year assuming 300 EFPD/year, hence greatly accelerating irradiation testing. Allowing a fast test reactor to operate over a wide range of inlet/outlet temperatures could add significant flexibility to its neutron irradiation capabilities. Steel present in the driver fuel assemblies is not expected to accumulate more than 90 displacements per atom (dpa) over its lifetime. Hence, steel alloys that are not suited for very high dpa but can operate over a wide range of temperatures (from low to high) would be appropriate for this application. Finally, because pin and plate assemblies fit on the same grid plate, a fast test reactor could in principle start-up with standard fuel pin assemblies and, later on, move to a plate-type Mark-II fuel to further in-crease its neutron irradiation capabilities while potentially reducing the level of uranium enrichment needed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

In Situ Visualization of the Pinning Effect of Planar Defects on Li Ion Insertion

Longevity of Li ion batteries strongly depends on the interaction of transporting Li ions in electrode crystals with defects. However, detailed interactions between the Li ion flux and structural defects in the host crystal remain obscure due to the transient nature of such interactions. Here, by in situ transmission electron microscopy and density function theory calculations, we reveal how the diffusion pathways and transport kinetics of a Li ion can be affected by planar defects in a tungsten trioxide lattice. We uncover that changes in charge distribution and lattice spacing along the planar defects disrupt the continuity of ion conduction channels and dramatically increase the energy barrier of Li diffusion, thus, arresting Li ions at the defect sites and twisting the lithiation front. In conclusion, the atomic-scale understanding holds critical implications for rational interface design in solid-state batteries and solid oxide fuel cells.

25 ENERGY STORAGE↗

Sr(Ag 1− x Li x ) 2 Se 2 and [Sr 3 Se 2 ][(Ag 1− x Li x ) 2 Se 2 ] Tunable Direct Band Gap Semiconductors

Abstract Synthesizing solids in molten fluxes enables the rapid diffusion of soluble species at temperatures lower than in solid‐state reactions, leading to crystal formation of kinetically stable compounds. In this study, we demonstrate the effectiveness of mixed hydroxide and halide fluxes in synthesizing complex Sr/Ag/Se in mixed LiOH/LiCl. We have accessed a series of two‐dimensional Sr(Ag 1− x Li x ) 2 Se 2 layered phases. With increased LiOH/LiCl ratio or reaction temperature, Li partially substituted Ag to form solid solutions of Sr(Ag 1− x Li x ) 2 Se 2 with x up to 0.45. In addition, a new type of intergrowth compound [Sr 3 Se 2 ][(Ag 1− x Li x ) 2 Se 2 ] was synthesized upon further reaction of Sr(Ag 1− x Li x ) 2 Se 2 with SrSe. Both Sr(Ag 1− x Li x ) 2 Se 2 and [Sr 3 Se 2 ][(Ag 1− x Li x ) 2 Se 2 ] exhibit a direct band gap, which increases with increasing Li substitution ( x ). Therefore, the band gap of Sr(Ag 1− x Li x ) 2 Se 2 can be precisely tuned via fine‐tuning x that is controlled by only the flux ratio and temperature.

Zhou, Xiuquan↗

Midwest US Croplands Determine Model Divergence in North American Carbon Fluxes

Large uncertainties in North American terrestrial carbon fluxes hinder regional climate projections. Terrestrial biosphere models (TBMs), the essential tools for understanding continental-scale carbon cycle, diverge on whether temperate forests or croplands dominate carbon uptake in North America. Evidence from novel photosynthetic proxies, such as those based on chlorophyll fluorescence, has cast doubt on the “weak cropland, strong forest” carbon uptake patterns simulated by most TBMs. However, no systematic evaluation of TBMs has yet been attempted to pin down space-time patterns that are most consistent with regional CO 2 observational constraints. Here, we leverage atmospheric CO 2 observations and satellite-observed photosynthetic proxies to understand emergent space-time patterns in North American carbon fluxes from a large suite of TBMs and data-driven models. To do so, we evaluate how well the atmospheric signals resulting from carbon flux estimates reproduce the space-time variability in atmospheric CO 2 , as is observed by a network of continuous-monitoring towers over North America. Models with gross or net carbon fluxes that are consistent with the observed CO 2 variability share a salient feature of growing-season carbon uptake in Midwest US croplands. Conversely, the remaining models place most growing-season uptake in boreal or temperate forests. Differences in model explanatory power depend mainly on the simulated annual cycles of cropland uptake—especially, the timing of peak uptake—rather than the distribution of annual mean fluxes across biomes. Our results suggest that improved model representation of cropland phenology is crucial to robust, policy-relevant estimation of North American carbon exchange.

54 ENVIRONMENTAL SCIENCES↗

Numerical study of heat transfer enhancement in laminar and turbulent flow in enhanced tubes with uniform and non-uniform thermal boundary conditions

Various types of tube enhancements have been proposed in order to increase heat transfer while minimizing any associated increases in pressure drop. This paper presents a numerical study of two different types of enhanced tubes: the relatively common twisted tape insert tube, and a novel tube with airfoil shaped pins on the inner surface. Each tube is studied in both the laminar and turbulent flow regime. In addition to the typical constant temperature and constant flux boundary conditions, a non-constant heat flux boundary condition is applied on one half of the tube, a situation which may be encountered when absorbing solar radiation. In general, each type of enhanced tube offers a thermal performance factor greater than unity at certain Re. In addition, the non-uniform heat flux boundary condition appears to function similarly as the uniform heat flux boundary condition.

Sandlin, Matthew J.↗

Numerical study of heat transfer enhancement in laminar and turbulent flow in enhanced tubes with uniform and non-uniform thermal boundary conditions

Various types of tube enhancements have been proposed in order to increase heat transfer while minimizing any associated increases in pressure drop. This paper presents a numerical study of two different types of enhanced tubes: the relatively common twisted tape insert tube, and a novel tube with airfoil shaped pins on the inner surface. Each tube is studied in both the laminar and turbulent flow regime. In addition to the typical constant temperature and constant flux boundary conditions, a non-constant heat flux boundary condition is applied on one half of the tube, a situation which may be encountered when absorbing solar radiation. In general, each type of enhanced tube offers a thermal performance factor greater than unity at certain Re. In addition, the non-uniform heat flux boundary condition appears to function similarly as the uniform heat flux boundary condition.

Sandlin, M↗

Electron Anisotropies in Magnetotail Dipolarization Events

Anisotropic electron distributions can be a major source of free energy for the generation of microinstabilities and waves. Here we investigate specifically anisotropies of the suprathermal electrons associated with dipolarization events in the geomagnetic tail. The investigation is based on an MHD simulation of magnetotail reconnection, flow bursts and dipolarization. This simulation is used to trace test particles in the electromagnetic fields of the MHD simulation and investigate their acceleration and contributions to flux increases. The simulated velocity distributions yield anisotropies which can be dominantly parallel (“cigar”-shaped), perpendicular (“pancake”-shaped) to the magnetic field, or a combination thereof (e.g., “rolling-pin”-shaped), depending on location, relative timing, and energy.

electron anisotropies↗

Progress on Demonstration of a MOOSE-Based Coupled Capability for Hot Channel Factors in Fast Reactors

Hot channel factors (HCFs) are computed values that account for the impact on predicted peak fuel, cladding, and coolant temperatures due to uncertainties in the as-built reactor’s material properties and geometry as well as uncertainties due to modeling approximations. Reduction in computed HCF values via reduction or elimination of modeling approximations may translate to significant economic savings if the reactor power can be raised due to the extra temperature margin gained. While limited historical datasets exist for sodium-cooled fast reactors (SFRs), there are no available HCF data for lead-cooled fast reactors (LFRs) outside of work generated previously within NEAMS. The computation of HCFs involves insights from reactor physics, thermal fluids and heat conduction calculations to determine how the peak temperatures respond to various uncertainties in the design. Due to the significant advantages for multi-physics coupling offered by the MOOSE framework, Griffin (MOOSE-based reactor physics code), MOOSE Heat Conduction Module, and Cardinal (MOOSE-wrapped multi-physics application which includes the NekRS thermal fluids code) are being coupled together using the MOOSE MultiApp System to develop a highfidelity multi-physics modeling capability for HCF simulations. This high-fidelity coupling workflow may also be beneficial for other fast reactor applications in the future. In previous work, Griffin and NekRS were individually assessed to ensure the necessary capabilities were in place. This work describes initial efforts to couple the codes (including folding in the MOOSE Heat Conduction Module) and determining the workflow for the perturbed calculations which will leverage the Stochastic Tools Module (STM). To our knowledge, this is the first coupling of Griffin and NekRS as well as the first exploratory use of Stochastic Tools Module for Cardinal. In this report, the neutronics code Griffin, the heat conduction solver in MOOSE, and the MOOSE-wrapped application containing NekRS (Cardinal) are linked together to demonstrate the coupled capability. Griffin and Cardinal are linked dynamically by specifying shared libraries. Different coupling hierarchies are tested for selecting the most appropriate coupling strategy. A coupling scheme is selected based on the efficiency of calculation and ease of data communication. Multiple tests are performed to choose suitable mesh structure, model configurations, scheme setup and boundary conditions to avoid loss of energy due to data interpolation between different modules or weak imposition of fluxes in finite element codes. Computational experiments are performed to study the tolerance control of each type of iteration to avoid false convergence. The coupled capability is demonstrated in both single pin and 7-pin models based on LFR materials and geometry. The study finds that the use of too large a time step size in the heat conduction module can lead to temperature oscillation even though the heat conduction equation does not have a time-derivative kernel, but only the time-dependent boundary condition. A 7-pin model without duct region achieved good convergence in the coupled calculation while a 7 pin model with duct region experienced data communication issues which need to be resolved.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Examination of stiff ion temperature gradient mode physics in simulations of DIII-D H-mode transport

A systematic evaluation of gyrokinetic and gyrofluid model predictions of ion temperature gradient (ITG) stability and transport using parameters from DIII-D high confinement mode (H-mode) plasmas has been performed. The nonlinear CGYRO code is used to make the gyrokinetic predictions, and the quasilinear TGLF model for the corresponding gyrofluid predictions. The assessments are made at three radii (normalized toroidal flux ρtor = 0.4, 0.55, and 0.7) in three different plasma scenarios with varying levels of neutral beam heating and torque. For each of the nine cases (3 radii × 3 scenarios) considered, ITG turbulence is found to be the dominant long-wavelength instability and transport mechanism. The inclusions of both transverse magnetic fluctuations and dynamic fast beam ions are stabilizing for all cases considered, with strongest effects seen at ρor = 0.4 where the fast ion population and normalized plasma pressure β = 2μ0nT/B2 are highest. The further inclusion of parallel magnetic fluctuations does not have a meaningful impact on the ITG turbulence in these scenarios, but does destabilize (in combination with fast ions) new high-frequency instabilities at ρtor = 0.4 in the high power scenarios. In each case the linear and nonlinear ITG critical gradients are predicted to be lower than the measured ITG scale lengths and their associated uncertainties. Inclusion of equilibrium flow shear in the transport predictions generally leads to an upshift in effective critical gradient rather than a qualitative change in the predicted stiffness, with stronger responses typically seen in the gyrokinetic predictions than in the gyrofluid results. However, in most cases these upshifted gradients still remain below the measured values and their uncertainties. Although the predicted critical gradients are below the measured gradients, both models predicted flux-matching gradients consistent with measured values in six of the nine cases considered, with no clear systematic over- or underprediction. Thus, while the experimental ion temperature profiles do not appear to be closely pinned to the ITG critical gradient, both gyrokinetic and gyrofluid models are able to accurately match the measured gradients reasonably well in most cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sensor Reduction for Diversion Detection in a Realistic Heat Pipe Microreactor Using Supervised Machine Learning

Microreactors are designed as a smaller, cheaper, and safer alternative to traditional nuclear power plants. Their non-traditional characteristics and prospect of mass production and deployment will likely require new approaches to nuclear safeguards. The primary proliferation concern with microreactors is the diversion of fuel material. Such diversion may produce measurable defects in key physical attributes like neutron flux, which may in turn be detectable using machine learning models. Preliminary work has demonstrated this ability for modeled nominal and diversion scenarios using large quantities of energy integrated neutron flux data. In practice, the number of available sensors for such measurements will be limited and energy integrated flux information will not be available. This work explores the ability of tree-based gradient boosted ensemble models to classify a given microreactor core is nominal or diversion, and determine the number of fuel pins diverted in the case of diversion with reduced numbers of sensors and more realistic detector responses. Classification accuracy of greater than 98% and regression errors as low as 5% of the total number of fuel pins were achieved with as few as 15 sensors, compared to 99% and 4.1% with a maximum of 240 sensors.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

FY26 Progress on Demonstration of a Multiphysics Steady State Capability for Modeling Core Radial Expansion in SFRs

Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism in liquid metal-cooled fast reactors and involves multiphysics effects including radiation transport, fluid flow, heat transfer, and mechanical response to temperature and flux gradients. This report summarizes recent progress on developing a multiphysics, MOOSE-based workflow to predict core bowing and associated reactivity feedback. Significant new capabilities in the reactor physics code Griffin - sodium backfill and pin power reconstruction for deformed geometries - were applied in this effort. This year’s work included verification, code comparisons, sensitivity studies, and coupled demonstrations that advance the state of MOOSE-based core bowing workflow. Griffin’s sodium backfill capability was verified by demonstrating that its automated treatment of geometry expansion and material-density updates reproduces manual calculations exactly, confirming solid mass conservation and proper coolant backfilling in expanded geometries. Reconstructed pin powers were compared for Griffin’s ductheterogeneous and ring-heterogeneous treatments in single-, seven-, and nineteen-assembly cases, with best agreement observed in lower-leakage configurations and the duct-heterogeneous approach offering substantially lower computational cost. Thermal-hydraulic sensitivity sensitivities showed that MOOSE SCM, SAM, and CFD are expected to produce similar deformation predictions despite variances in their temperature predictions, and that explicit treatment of inter-assembly flow becomes increasingly important as gap flow rate increases. Finally, coupled demonstrations on small multi-assembly configurations using Griffin, MOOSE Solid Mechanics, MOOSE SCM, and Heat Conduction produced physically consistent reactivity feedback from thermal expansion and bowing. The coupled demonstrations simulated grid plate expansion as well as resultant core bowing at full power conditions. Simplifications were made in current workflow, namely the assumption of instantaneous full power conditions following hot zero power, and pre-expanding the Griffin geometry axially due to lack of an axial fuel pin expansion model and temperature feedback to Griffin.

Wozniak, Nicholas↗

Magnetic flux trapping in hydrogen-rich high-temperature superconductors

Recent discoveries of superconductivity in various hydrides at high pressures have shown that a critical temperature of superconductivity can reach near-room-temperature values. However, experimental studies are limited by high-pressure conditions, and electrical transport measurements have been the primary technique for detecting superconductivity in hydrides. Here we implement a non-conventional protocol for the magnetic measurements of superconductors in a SQUID magnetometer and probe the trapped magnetic flux in two near-room-temperature superconductors H 3 S and LaH 10 at high pressures. Contrary to traditional magnetic susceptibility measurements, the magnetic response from the trapped flux is almost unaffected by the background signal of the diamond anvil cell due to the absence of external magnetic fields. The behaviour of the trapped flux generated under zero-field-cooled and field-cooled conditions proves the existence of superconductivity in these materials. We reveal that the absence of a pronounced Meissner effect is associated with the very strong pinning of vortices inside the samples. This approach can also be a tool for studying multiphase samples or samples that have a low superconducting fraction at ambient pressure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Material-dependent photon ionizing radiation effects in Si and GaAs PIN diodes: A numerical investigation

We present a finite-element drift-diffusion-Poisson model in the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework to compare the radiation response of silicon (Si) and gallium arsenide (GaAs) PIN diodes under high-energy photon irradiation. The model solves coupled carrier continuity and Poisson’s equations with Shockley-Read-Hall recombination, and is verified against standard analytical J-V behavior. Using a simplified 1D geometry with ideal Ohmic contacts, we quantify device response under forward and reverse bias with a 100 MeV photon flux. Under forward bias, Si exhibits markedly greater radiation sensitivity than GaAs, including larger increases in current density, stronger local field and carrier-product perturbations, and higher recombination. Under reverse bias, GaAs shows larger radiation-induced photocurrent and broader current-density peaks near junctions, indicating an advantage for photodetection. Integrated steady-state recombination is consistently higher in Si across voltages. Under periodic photon pulses, GaAs produces higher-amplitude photoresponse and settles more rapidly than Si. These results highlight material-dependent trade-offs for radiation-tolerant, high-speed optoelectronics and provide guidance for selecting PIN architectures in aerospace, nuclear, and high-energy physics environments.

36 MATERIALS SCIENCE↗

Calibration of the NEAR XRS Solar Monitor

The NEAR-Shoemaker X-ray spectrometer measured elemental abundance ratios for six elements (Mg, Al, Si, S, Ca, Fe) on the surface of 433 Eros by remote-sensing x-ray fluorescence spectroscopy (1,2). The fluorescent x-rays are excited by x-rays from the solar corona. Therefore, accurate knowledge of the incident solar spectrum was critically important to interpreting the fluorescence spectra received by the NEAR detectors in terms of elemental ratios on the asteroid's surface. NEAR, therefore, carried two additional X-ray spectrometers dedicated to monitoring the solar spectrum between 1 and 10 keV. One of these, a Si PIN photodiode, failed prior to orbit insertion. The second consisted of a gas-filled proportional counter similar to the asteroid-pointing detectors, but covered by a graded filter (3) intended to attenuate the incoming solar flux and increase the dynamic range of the detector. Unfortunately, the response of this filter as a function of solar incidence angle was not adequately determined in the laboratory prior to launch. The response function had to be calculated from the assumed properties of the graded filter.

Lim, L. F.↗