SEARCH · Engineering Papers

Results for “Rotational dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Sliding dynamics of skyrmion molecular crystals

Using both atomistic and particle-based simulations, we investigate the current-driven dynamics of skyrmions on two-dimensional periodic substrates when there are multiple skyrmions per substrate minimum. At zero drive, the system forms pinned skyrmion molecular crystal states consisting of dimers, trimers, or dimer-trimer mixtures that have both positional and orientational order. On a square substrate lattice, the motion above depinning occurs via a running soliton that travels completely transverse to the applied current. This motion is generated by a torque from the Magnus force, which rotates the n-mer states perpendicular to the applied current. At higher drives, the flow becomes disordered while the Hall angle diminishes and gradually approaches the intrinsic value. In some cases, we also find directional locking where the Hall angle becomes locked to certain symmetry directions of the substrate over a range of currents. The transitions into and out of directionally locked states are accompanied by negative differential mobility in which the net velocity decreases as the drive increases. In conclusion, on a triangular substrate, we find no transverse mobility effects, but still observe directionally locked motion.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Understanding the lattice thermal conductivity of SrTiO 3 from an ab initio perspective

Here, we present a detailed analysis of the structure dependence of the lattice thermal conductivity of SrTiO 3 . We have used both ab initio molecular dynamic simulations and density functional theory calculations to decouple the effect of different structural distortions on the thermal conductivity. We have identified two main mechanisms for tuning the thermal conductivity when a distortion is applied. First, the modification of the acoustic-modes energy dispersion when a change in the lattice parameters is imposed and second, the low energy polar modes. In particular and counterintuitively, we have found that an increase in the angle of the oxygen octahedral rotations increases the thermal conductivity due to its coupling to these polar modes.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Numerical Modeling of a Two-Stage Ocean Current Turbine

The Equinox Ocean Turbines (EQOT) current energy converter has a unique design with power generation in two small-diameter turbines attached to the tips of a large-diameter passive rotor. This configuration offers some key advantages for capturing ocean currents. With no centrally placed generator, almost no reaction torque is required at the nacelle of the main large-diameter rotor, and the small-diameter tip turbine generators operate at a higher speed and lower torque. The physics that determine the performance and loads on the turbine are also unique. The interactions of the flow field between the two stages and the general architecture of the system cannot be captured with traditional mid-fidelity modeling tools. For design iterations and large sets of load cases, it is important to have mid-fidelity models that can capture the important phenomenon with enough accuracy to identify global trends. This work uses a limited set of high-fidelity computational fluid dynamics (CFD) simulations to help inform the selection of and construction of a custom mid-fidelity model. Mid-fidelity modeling approaches were verified by comparing key turbine performance quantities to those found with the CFD model. Hydrodynamic interactions of the two-stage rotor were identified through high-fidelity CFD modeling. This highlighted the impact of the main rotor tip vortex and wake on the secondary rotor apparent inflow. This results in a relative flow rotation and sharp deficit, that change the optimal secondary rotor rotation speed and adds unsteadiness to the blade loading respectively. Multiple mid-fidelity approaches were evaluated for their ability to capture these effects. A simple approximation of the combined-stage performance based on single-stage BEM provides a reasonable rough prediction, especially near the peak TSR values, with some larger discrepancy at higher TSRs. Predicting the combined-stage performance based on single-stage CFD data improves this prediction across the TSR range. Although the combined-stage modeling in OLAF was not successful in this stage of the project, it showed promise as a mid-fidelity method, assuming the parameters can be tuned to account for the significant differences in time and length scales between the main and secondary rotors. This may be addressed through code changes in future work. A significant finding from the OLAF work was the agreement between the vortex core radius values found independently via a parameter space search and via CFD. The technique of using single-stage secondary rotor BEM, with a custom inflow taken from single-stage main rotor CFD or OLAF, provides an efficient method to capture one-way coupled flow interactions. This method provided generally good predictions of the impact of the flow rotation on the secondary rotor but struggled to accurately predict the peaks of the unsteady load progression. Future work could include some superposition of a tuned main rotor trailing edge viscous wake into the custom inflow to better predict this interaction.

16 TIDAL AND WAVE POWER↗

DOE OSTI2025

Electronic Melting of Silicon in Nanostructures using X-ray Forbidden Bragg Reflections

Abstract We carried out a short beamtime at the Pohang Accelerator Laboratory x-ray Free Electron Laser to perform a pump-probe (PP) laser excitation diffraction experiment on the silicon (222) forbidden Bragg peak. To limit the x-ray penetration, we used a “device layer” silicon film wafer bonded to a silicon substrate. The sample, specially fabricated by MEMC Electronic Materials, had a Si(100) substrate bonded to a 170 nm Si(100) film rotated at 45° for crystallographic isolation. A second sample was reactive-ion-etched down to 52 nm thickness. In the silicon lattice, the covalent bonds are seen exclusively at the 222 reflection. Upon laser excitation, these electrons are expected to be excited to the valence band on femtosecond electronic time scales. The Si(222) reflection is therefore expected to be extinguished on this fast time scale, while the electron–phonon coupled acoustic response is determined by the lattice dynamics. The latter is determined by the speed of sound over the device thickness, which is in the mid-picosecond range.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2023

Hysteretic Responses of Skyrmion Lattices to Electric Fields in Magnetoelectric Cu 2 OSeO 3

Electric field control of topologically nontrivial magnetic textures, such as skyrmions, provides a paradigm shift for future spintronics beyond current silicon-based technology. While significant progress has been made by X-ray and neutron scattering studies, direct observation of such nanoscale spin structures and their dynamics driven by external electric fields remains a challenge to understand underlying mechanisms and harness functionalities. Here, using Lorentz transmission electron microscopy combined with in situ electric and magnetic fields at liquid helium temperatures, we report the crystallographic-orientation dependent skyrmion responses to electric fields in thin slabs of magnetoelectric Cu 2 OSeO 3 . We show electric fields not only stabilize the hexagonally packed skyrmion lattices in the entire sample in a hysteretic manner but also induce the rotation of their reciprocal vector discretely by 30°. Furthermore, the nonvolatile and energy-efficient skyrmion-lattice control by electric fields demonstrated in this work provides an important foundation for designing skyrmion-based qubits and memory devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

A simple catch: Fluctuations enable hydrodynamic trapping of microrollers by obstacles

It is known that obstacles can hydrodynamically trap bacteria and synthetic microswimmers in orbits, where the trapping time heavily depends on the swimmer flow field and noise is needed to escape the trap. Here, we use experiments and simulations to investigate the trapping of microrollers by obstacles. Microrollers are rotating particles close to a bottom surface, which have a prescribed propulsion direction imposed by an external rotating magnetic field. The flow field that drives their motion is quite different from previously studied swimmers. We found that the trapping time can be controlled by modifying the obstacle size or the colloid-obstacle repulsive potential. We detail the mechanisms of the trapping and find two remarkable features: The microroller is confined in the wake of the obstacle, and it can only enter the trap with Brownian motion. While noise is usually needed to escape traps in dynamical systems, here, we show that it is the only means to reach the hydrodynamic attractor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Unconventional colloidal aggregation in chiral bacterial baths

Abstract When in equilibrium, thermal forces agitate molecules, which then diffuse, collide and bind to form materials. However, the space of accessible structures in which micron-scale particles can be organized by thermal forces is limited, owing to the slow dynamics and metastable states. Active agents in a passive fluid generate forces and flows, forming a bath with active fluctuations. Two unanswered questions are whether those active agents can drive the assembly of passive components into unconventional states and which material properties they will exhibit. Here we show that passive, sticky beads immersed in a bath of swimmingEscherichia colibacteria aggregate into unconventional clusters and gels that are controlled by the activity of the bath. We observe a slow but persistent rotation of the aggregates that originates in the chirality of theE. coliflagella and directs aggregation into structures that are not accessible thermally. We elucidate the aggregation mechanism with a numerical model of spinning, sticky beads and reproduce quantitatively the experimental results. We show that internal activity controls the phase diagram and the structure of the aggregates. Overall, our results highlight the promising role of active baths in designing the structural and mechanical properties of materials with unconventional phases.

Physics↗

DOE OSTI2023

Local structure, bonding, and asymmetry of $\mathrm{((NH_2)_2CH)PbBr_3}$, $\mathrm{CsPbBr_3}$, and $\mathrm{(CH_3MH_3)PbBr_3}$

Here we report local structure measurements for CsPbBr 3 and ((NH 2 ) 2 CH)PbBr 3 (FAPbBr 3 ) and compare them with recent results for (CH 3 NH 3 )PbBr 3 (MAPbBr 3 ). The Pb-Br bonding is similar for all three systems; the effective spring constants, κ, are comparable (ranging from 1.20 to 1.95 eV/Å 2 ), but small in magnitude indicating very soft materials. However, there are also important differences between the three systems. Static disorder is very small for CsPbBr3 but increases somewhat with the size of the organic molecular ions MA + and FA + . At room temperature, dynamic disorder dominates in all compounds. The thermal disorder of the Pb-Br pair distribution function (PDF), i.e., the Debye-Waller factor σ 2 follows a correlated Debye or Einstein model up to 300 K in CsPbBr 3 (orthorhombic phase), but for FAPbBr 3 and MAPbBr 3 , there is a break in the σ 2 (T) curve at the orthorhombic-tetragonal transition (o-t) near 150 K, indicating a small change in the spring constant κ. κ increases for MAPbBr 3 but decreases in FAPbBr 3 at this transition. These changes are attributed to changes in the H-bonding between Br - and MA + or FA + at this transition, as a result of librations or rotations of these molecular cations. In addition, the Pb-Br PDF becomes asymmetric at a relatively low temperature for FAPbBr 3 and MAPbBr 3 , while this effect is significantly smaller for CsPbBr 3 . Finally, we address the question of a model to explain the asymmetric PDF. Two main models are discussed in the literature, an anharmonic pair potential and a split-pair distribution, possibly driven by the presence of a lone pair on the Pb ion. We show that the fourth cumulant C 4 can differentiate between these two models and other possible models. Experimentally C 4 is positive at 250 K and above, for all three systems and that is inconsistent with a split-peak model, for which C 4 is negative for splittings larger than 0.12 Å.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Diversifying and perennializing plants in agroecosystems alters retention of new C and N from crop residues

Abstract Managing soils to retain new plant inputs is key to moving toward a sustainable and regenerative agriculture. Management practices, like diversifying and perennializing agroecosystems, may affect the decomposer organisms that regulate how new residue is converted to persistent soil organic matter. Here we tested whether 12 years of diversifying/perennializing plants in agroecosystems through extended rotations or grassland restoration would decrease losses of new plant residue inputs and, thus, increase retention of carbon (C) and nitrogen (N) in soil. We tracked dual‐labeled ( 13 C and 15 N), isotopically enriched wheat ( Triticum aestivum ) residue in situ for 2 years as it decomposed in three agroecosystems: maize–soybean (CS) rotation, maize–soybean–wheat plus red clover and cereal rye cover crops (CSW2), and spring fallow management with regeneration of natural grassland species (seven to 10 species; SF). We measured losses of wheat residue (C wheat and N wheat ) in leached soil solution and greenhouse gas fluxes, as well as how much was recovered in microbial biomass and bulk soil at 5‐cm increments down to 20 cm. CSW2 and SF both had unique, significant effects on residue decomposition and retention dynamics that were clear only when using nuanced metrics that able to tease apart subtle differences. For example, SF retained a greater portion of C wheat in 0–5 cm surface soils (155%, p = 0.035) and narrowed the C wheat to N wheat ratio ( p < 0.030) compared to CS. CSW2 increased an index of carbon‐retention efficiency, C wheat retained in the mesocosm divided by total measured, from 0.18 to 0.27 (49%, p = 0.001), compared to CS. Overall, we found that diversifying and extending the duration of living plants in agroecosystems can lead to greater retention of new residue inputs in subtle ways that require further investigation to fully understand.

McDaniel, Marshall D.↗

DOE OSTI2024

Analysis of single-mode Richtmyer–Meshkov instability using high-order incompressible vorticity-streamfunction and shock-capturing simulations

In this study, two- and three-dimensional simulation results obtained using a new high-order incompressible, variable-density vorticity–streamfunction (VS) method and data from previous ninth-order weighted essentially nonoscillatory (WENO) shock-capturing simulations are used to investigate the nonlinear dynamics of single-mode Richtmyer–Meshkov instability using a model of a Mach 1.3 air(acetone)/SF6 shock tube experiment. A comparison of the density fields from both simulations with the experimental images demonstrates very good agreement in the large-scale structure with both methods but differences in the small-scale structure. The WENO method captures the small-scale disordered structure observed in the experiment, while the VS method partially captures such structure and yields a strong rotating core. The perturbation amplitude growth from the simulations generally agrees well with the experiment. The simulation bubble and spike amplitudes agree well at early times. At later times, the WENO bubble amplitude is smaller than the VS amplitude and vice versa for the spike amplitude. The predictions of nonlinear single-mode instability growth models are shown to agree with the simulation amplitudes at early-to-intermediate times but underpredict the amplitudes at later times in the nonlinear regime. Visualizations of the mass fraction and enstrophy isosurfaces, velocity and vorticity fields, and baroclinic vorticity production and vortex stretching terms from the three-dimensional simulations indicate that, with the exception of the small-scale structure within the rollups, the VS and WENO results are in good agreement.

42 ENGINEERING↗

DOE OSTI2021

Analyzing Gas Turbine-Generator Performance of the Hybrid Power System

We report kinetic energy available from the turbomachinery may boost short-term ramp rates during load shifting between generator and fuel cell that may be exploited to optimize the performance of turbine-fuel cell hybrid systems during load following. The paper starts with modeling a gas turbine, gearbox, and generator sitting on the rigid shafts as a two-mass model primarily for estimating the rotational inertia available from the turbogenerator. The gas turbine and its environment impose constraints of available space allowances, measurements from high-speed rotating parts, lack of accessibility to the shaft of measure, making replacements and/or required installation adjustments of the torque transducer difficult. It is demonstrated that in the absence of torque measurement, the data containing only a slow excursion of the rotor motion in response to a slight load change is sufficient to determine inertia and damping coefficients from the model. Furthermore, the generator’s reactive power capability was investigated because it serves as an asset in integrated hybrid systems to meet dynamic reactive power demands, reducing the burden on the fuel cell’s power electronics control in generating reactive power.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2023

Effects of Magnetic Fields on Gas Dynamics and Star Formation in Nuclear Rings

Nuclear rings at the centers of barred galaxies are known to be strongly magnetized. To explore the effects of magnetic fields on star formation in these rings and nuclear gas flows, we run magnetohydrodynamic simulations in which there is a temporally constant magnetized inflow to the ring, representing a bar-driven inflow. The mass inflow rate is 1M ⊙ yr –1 , and we explore models with a range of field strength in the inflow. We adopt the TIGRESS framework developed by Kim & Ostriker to handle radiative heating and cooling, star formation, and resulting supernova (SN) feedback. We find that magnetic fields are efficiently amplified in the ring due to rotational shear and SN feedback. Within a few 100 Myr, the turbulent component Btrb in the ring saturates at ~35 μG (in rough equipartition with the turbulent kinetic energy density), while the regular component Breg exceeds 50 μG. Expanding superbubbles created by clustered SN explosions vertically drag predominantly toroidal fields from near the midplane to produce poloidal fields in high-altitude regions. The growth of magnetic fields greatly suppresses star formation at late times. Simultaneously, strong magnetic tension in the ring drives radially inward accretion flows from the ring to form a circumnuclear disk in the central region; this feature is absent in the unmagnetized model.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Understanding the phase transformation mechanisms that affect the dynamic response of Fe-based microstructures at the atomic scales

Large-scale molecular dynamics (MD) simulations were carried out to investigate the shock-induced evolution of microstructure in Fe-based systems comprising single-crystal and layered Cu/Fe alloys with a distribution of interfaces. The shock compression of pure single-crystal Fe oriented along [110] above a threshold pressure results in a BCC (α)→HCP (ε) phase transformation behavior that generates a distribution of ε phase variants in the phase transformed region of the microstructure behind the shock front. The propagation of the release wave through a phase transformed ε phase causes a reverse ε→α phase transformation and renders a distribution of twins for the [110] oriented Fe that serve as void nucleation sites during spall failure. The simulations reveal that the α→ε→α transformation-induced twinning for shock loading along the [110] direction is due to a dominant ε phase variant formed during compression that rotates on the arrival of the release wave followed by a reverse phase transformation to twins in the α phase. The modifications in the evolution of the ε phase variants and twins in Fe behavior are also studied for Cu–Fe layered microstructures due to the shock wave interactions with the Cu/Fe interfaces using a newly constructed Cu–Fe alloy potential. Here, the MD simulations suggest that interfaces affect the observed variants during shock compression and, hence, distributions of twins during shock release that affects the void nucleation stresses in the Fe phase of Cu/Fe microstructures.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Control Design for a Marine Hydrokinetic Cycloturbine Vehicle

Abstract Marine Hydrokinetic (MHK) cycloturbines generate sustainable power by exploiting tidal currents. By powering the turbines and using pitching foils for control, a vehicle comprised of MHK cycloturbines also has the ability to station keep and maneuver. The vehicle consists of four counter-rotating cycloturbines, with hydrofoils oriented perpendicular to the flow in a paddlewheel configuration. Lift and drag generated from these foils sum together to produce thrust. An experimentally tuned simulation model that solves the six-degrees-of-freedom rigid body equations of motion for the MHK vehicle subject to hydrodynamic, hydrostatic, and propulsive forces is used to aid the design of vehicle controllers. Global feedback controllers are initially designed by applying classical control methods to an approximate linear model of the system dynamics. A higher performing nonlinear controller is designed using the nonlinear dynamic inversion (NDI) method. NDI accounts for the nonlinearities of the MHK system and therefore is suitable for a wide range of operating conditions. The response of the classical and NDI controllers to speed, depth, roll, pitch, and yaw commands are evaluated and compared in simulation. The classical controller outperforms the NDI controller for small amplitude maneuvers, although the degradation with NDI is minor. However, in the nonlinear operating regime the NDI controller outperforms the classical controller and the classical controller exhibits instability.

Automation & Control Systems↗

DOE OSTI2021

Distributionally Safe Path Planning: Wasserstein Safe RRT

In this paper, we propose a Wasserstein metric-based random path planning algorithm. Wasserstein Safe RRT (W-Safe RRT) provides finite-sample probabilistic guarantees on the safety of a returned path in an uncertain obstacle environment. Vehicle and obstacle states are modeled as distributions based upon state and model observations. Additionally, we define limits on distributional sampling error so the Wasserstein distance between a vehicle state distribution and obstacle distributions can be bounded. This enables the algorithm to return safe paths with a confidence bound through combining finite sampling error bounds with calculations of the Wasserstein distance between discrete distributions. W-Safe RRT is compared against a baseline minimum encompassing ball algorithm, which ensures balls that minimally encompass discrete state and obstacle distributions do not overlap. The improved performance is verified in a 3D environment using single, multi, and rotating non-convex obstacle cases, with and without forced obstacle error in adversarial directions, showing that W-Safe RRT can handle poorly modeled complex environments.

42 ENGINEERING↗

DOE OSTI2020

Creep Deformation and Dynamic Grain Growth in an Interstitial-Free Steel

Dynamic grain growth is demonstrated to be much faster than static grain growth in a body-centered-cubic, interstitial-free steel sheet material at 850°C. Dynamic grain growth occurs during concurrent plastic deformation at elevated temperature, whereas static grain growth occurs during static annealing. Grain growth during steady-state plastic flow in tension at 850°C to a true strain of 0.2 at a true-strain rate of 10 -4 s -1 doubled grain size, while static annealing for the same time produced no increase in grain size. This is described as dynamic normal grain growth (DNGG) because no abnormally large grains were observed. The recrystallized microstructure of the steel demonstrated a log-normal distribution of grain sizes. DNGG produced bimodal grain size distributions that deviate from the theoretical expectation of a simple shift to larger sizes during normal growth. The bimodal distributions contained a remnant of small grains that were not consumed during grain growth. DNGG produced a crystallographic texture that is unique from both the recrystallized material and that produced by lattice rotation alone. DNGG strengthened the {111}&lang;110&rang; and {111}&lang;112&rang; components of the strong γ-fiber component in the original recrystallization texture. Lattice rotation from tensile deformation, by contrast, strengthened the α-fiber components that intersect the original γ-fiber.

36 MATERIALS SCIENCE↗

DOE OSTI2023

$μ$ $\mathrm{SR}$ study of the dipole-octupole quantum spin ice candidate $\mathrm{Ce_2Zr_2O_7}$

The Ce 3+ pseudospin-1/2 degrees of freedom in Ce 2 Zr 2 O 7 possess both dipolar and octupolar character which enables the possibility of novel quantum spin liquid ground states in this material. In this study, we report muon spin relaxation and rotation (μSR) measurements on single-crystal samples of Ce 2 Zr 2 O 7 in zero magnetic field and in magnetic fields directed along the [1, $\overline{1}$, 0] and [1, 1, 1] crystallographic directions, and for magnetic fields directed both longitudinal and transverse to the direction of muon polarization. Our zero-field results show no signs of magnetic ordering or spin freezing, consistent with earlier zero-field μSR measurements on a powder sample of Ce 2 Zr 2 O 7 , and also with the expectations for a quantum spin ice. However, we measure a more gentle relaxation rate for Ce 2 Zr 2 O 7 in zero field at low temperatures than was previously reported. This difference in relaxation rate is likely due to the low oxidation and, correspondingly, the high stoichiometry of our singlecrystal samples. Longitudinal field measurements confirm that the magnetic dipole moments in Ce 2 Zr 2 O 7 remain dynamic at T = 0.1 K on the microsecond timescale. For both [1, $\overline{1}$, 0] and [1, 1, 1] magnetic fields, our μSR Knight shift measurements show a field-induced leveling off of the magnetic susceptibility at low temperature which is qualitatively consistent with corresponding calculations using the numerical-linked-cluster method in combination with recent estimates for the nearest-neighbor exchange parameters of Ce 2 Zr 2 O 7 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

Online task-space motion control for positioner-coordinated multi-robot manufacturing systems

Incorporating multiple robotic manipulators into large-scale manufacturing systems enhances production efficiency and expands manufacturing capabilities beyond those of single-robot systems. Workpiece positioners in robotic manufacturing have demonstrated significant benefits for process optimization, but coordination strategies for multi-robot systems with shared positioners have received limited attention. This work presents a task-space coordinated trajectory-tracking control framework for multi-robot manufacturing systems, in which robots coordinate their motions within a shared, dynamic workpiece positioning frame. A workpiece positioner actively adjusts the pose of the manufactured component to enable greater operational concurrency and improve overall production efficiency. The proposed motion-coordination scheme employs a distributed and scalable architecture, supporting coordination across heterogeneous multi-robot systems. Two optimization methodologies are introduced to manage kinematic redundancies and maintain continuous, near-optimal operation throughout the manufacturing process. The first strategy exploits a task-space dimensionality reduction to achieve locally optimal configurations by leveraging symmetry-axis rotations of the tool. The second strategy utilizes the workpiece positioner to drive the coordinated robots toward stable and kinematically favorable configurations. For both optimization strategies, multiple objectives are defined to improve key performance metrics, including manipulability, configuration consistency, proximity to mechanical limits, and motion efficiency. Addressing a key limitation of existing coordination approaches, the framework is designed around online setpoint modification, allowing coordinated robots to respond effectively to in-situ process feedback. The proposed control framework is validated using the Robot Operating System (ROS) middleware on a combination of physical and simulated multi-robot system hardware.

Arbogast, Alex [ORNL] (ORCID:0000000154740723)↗