Finite Angle Research Articles

Optimally Suppressed Phonon Tunneling in van der Waals Graphene-WS2 Heterostructure with Ultralow Thermal Conductivity.

Van der Waals heterostructures have great potential for realizing ultimately low thermal conductivity because defectless interfaces can be constructed at a length scale smaller than the phonon wavelength, allowing modulation of coherent phonon transport. In this Letter, we demonstrate the mechanism for thermal conductivity reduction at a mode-resolved level. The graphene-WS2 heterostructure with the lowest cross-plane thermal conductivity of 0.048 W/(m·K) is identified from 16,384 candidates by combining Bayesian optimization and molecular dynamics simulations. Then, the angle-resolved phonon transmission is calculated using the mode-resolved atomistic Green's function. The results reveal that the optimal heterostructure nearly completely terminates phonon transport with finite incident angles, owing to the reduced critical incident angle and suppression of phonon tunneling.

Influence of the acceptance angle on the evaluation of reflectance data of randomly rough surfaces using scalar diffraction theory

Apart from coherent reflectance, which corresponds to specular reflection, the values obtained by real spectrophotometers also include contribution from incoherent reflectance, which represents light scattered by the samples and registered by the detector due to its finite acceptance angle. This work aims to investigate the influence of this second part on reflectance spectra measured for samples with randomly rough surfaces. Three silicon samples with roughened surfaces are investigated. The reflectance is measured using a commercial spectrophotometer with acceptance angles restricted by apertures placed in the incident and reflected beam. The proposed method is based on the simultaneous processing of spectral dependencies of reflectance measured with differently-sized apertures. The utilized theoretical approach is based on the scalar diffraction theory. Because the dependencies on both wavelength and acceptance angle are considered, a model providing correct predictions for these dependencies should also correctly describe how is the total reflectance separated into its coherent and incoherent parts. It is shown that the theoretical predictions for incoherent reflectance are consistent with the changes in the diameter of the apertures. It was possible to determine the RMS value of the heights as well as the estimate for the autocorrelation length and additional parameter controlling the course of the autocorrelation function. A short discussion comparing our results with those achieved using methods employed in earlier works is also provided.

Approach and rotation of reconnecting topological defect lines in liquid crystal

Topological defects are a universal concept across many disciplines, such as crystallography, liquid-crystalline physics, low-temperature physics, cosmology, and even biology. In nematic liquid crystals, topological defects called disclinations have been widely studied. For their three-dimensional (3D) dynamics, however, only recently have theoretical approaches dealing with fully 3D configurations been reported. Further, recent experiments have observed 3D disclination line reconnections, a phenomenon characteristic of defect line dynamics, but detailed discussions were limited to the case of approximately parallel defects. In this paper, we focus on the case of two disclination lines that approach at finite angles and lie in separate planes, a more fundamentally 3D reconnection configuration. Observing and analyzing such reconnection events, we find the square-root law of the distance between disclinations and the decrease of the interdisclination angle over time. We compare the experimental results with theory and find qualitative agreement on the scaling of distance and angle with time, but quantitative disagreement on distance and angle relative mobilities. To probe this disagreement, we derive the equations of motion for systems with reduced twist constant and also carry out simulations for this case. These, together with the experimental results, suggest that deformations of disclinations may be responsible for the disagreement. Published by the American Physical Society 2024

Wetting and interfacial behavior of aqueous Deep eutectic solvents on Graphite, silica and calcite surfaces: Molecular dynamics study

Understanding the wetting behavior of liquids/fluids on solid surfaces at the nanoscale can help develop more effective and efficient technologies in various applications, including microfluidics, biomedical, enhanced oil recovery, and many others. Deep eutectic solvents (DESs) have recently gained attention due to their physical and chemical properties and tunability, making them useful in many applications. This article reports the investigation of the wetting behavior of aqueous deep eutectic solvent (DES) on graphite, silicon dioxide (SiO2), and calcite (CaCO3) using molecular dynamics (MD) simulation. The study's main objective is to understand the interfacial behavior of aqueous DES on these surfaces. In the current study, two different DESs, choline-chloride (ChCl) as hydrogen bond acceptor (HBA) and two different hydrogen bond donors (HBD), urea (U) and ethylene glycol (EG) with molar ratio 1:2are used. Aqueous DES has been studied for three different concentrations of10, 20, and 30 wt% of DES molecules to explore the effect of concentration.AqueousDES system on the graphite makes a finite contact angle for all concentrations of DES.In contrast, for SiO2 and CaCO3, droplets spread and attain a contact angle below ∼ 10° due to high hydrophilicity. We confirmed that our simulation had attained equilibrium by computing the interaction energy with the substrate and system molecules. The wetting behavior of aqueous DES thin film on the surfaceis characterizedby a density profile along the z-axis (normal to the surface). Furthermore, hydrogen bond (HB) analysisis performedto elucidate the structural relationship of these aqueous-DES systems. Overall, this study provides insights into the wetting behavior of nanodroplets/film of aqueous DES on different surfaces and illustrates the structural characteristics of these systems through various properties such as contact angle, density profiles, and hydrogen bonds.

Atomic wetting of oil droplets into hexagons and stripes

HypothesisPrevious studies showed that liquid droplets can be spherical caps or polygonal shapes on patterned surfaces. We hypothesize that these phenomena can be narrowed down to nanoscale and affecting liquid molecular wetting. We anticipate observing preferential wetting on the underlying substrate with atomic crystalline lattice. ExperimentsThis study explores the wetting behavior of oil nanodroplets deposited on highly ordered pyrolytic graphite (HOPG) with atomic force microscope (AFM). Oil droplets were scanned with hydrophobic and hydrophilic AFM tips in water in peak force mode AFM, carefully controlling the scanning force. Simple modeling was carried out using Surface Evolver to simulate how oil droplets wet on HOPG with molecular lattice. FindingsThe results showed two distinct wetting behaviors: the hexagonal nanodroplets and stripe nanometric films. In both cases, the droplets appear to wet preferentially along preferred orientations which reflect the atomic periodicity of the underlying substrate. This preferential wetting occurs because the atomic roughness is significant at the length scale of the nanodroplets, while the eventual outcome of the wetting (i.e. whether the droplet settles to a finite contact angle, or wets completely into a nanometric film) is dependent on the spreading parameter of the droplet-substrate system, independent of scanning direction.

Binary discrimination through next-to-leading order

Binary discrimination between well-defined signal and background datasets is a problem of fundamental importance in particle physics. With detailed event simulation and the advent of extensive deep learning tools, identification of the likelihood ratio has typically been reserved as a computational problem. However, this approach can obscure overtraining or excessive sensitivity to tuned features of the simulation that may not be well-defined theoretically. Here, we present the first analysis of binary discrimination for signal and background distributions for which their likelihood ratio is infrared and collinear safe, and can therefore be calculated order-by-order in perturbation theory. We present explicit, general formulas for receiver operator characteristic curves and the area under it through next-to-leading order. These results can then establish absolute upper bounds on discrimination performance because any realistic implementation will have measurement errors, undetected particles, or restrictions on fiducial phase space. As a demonstration of this formalism, we apply it to discrimination of highly-boosted Higgs decays from gluon splitting to bottom quarks. Effects at next-to-leading order are first sensitive to the flow of color in the jet and significantly modify discrimination performance at leading-order. In the limit of infinite boost, these events can be perfectly discriminated because only the gluon will radiate at finite angles from the bottom quarks, and we find that large effects persist at energies accessible at the Large Hadron Collider. Next-to-leading order is therefore required to qualitatively understand results using machine-learning methods.

Mathematical fundamentals of spherical kinematics of plate tectonics in terms of quaternions

To be a quantitative and testable tectonic model, plate tectonics requires spherical geometry and spherical kinematics in terms of finite rotations conveniently parametrized by their angle and axis and described by unit quaternions. In treatises on “Plate Tectonics” infinitesimal, instantaneous, and finite rotations, absolute and relative rotations are said to be applied to model the motion of tectonic plates. Even though these terms are strictly defined in mathematics, they are often casually used in geosciences. Here, their definitions are recalled and clarified as well as the terms rotation, orientation, and location on the sphere. For instance, infinitesimal rotations refer to a mathematical limit, when the angle of rotation tends to zero. Their rules do not apply to finite rotations, no matter how small their finite angles of rotation are. Mathematical approaches applying appropriate and feasible assumptions to model spherical motion of tectonic plates over geological times of hundreds of millions of years are derived including (i) sequences of incremental finite rotations, (ii) sequences of accumulating successive concatenations of finite rotations, and (iii) continuous rotations in terms of fully transient quaternions. The incremental and the accumulating approaches provide complementary views. While the relative Euler pole appears to migrate in the latter, it appears fixed in the former. Path, mean, and instantaneous velocity of the migrating Euler pole are derived as well as the angular and trajectoral velocity of the rotational motion about it. The approaches are illustrated by a geological example with actual data and a numerical yet geologically inspired example with artificial data. The former revisits the three‐plate scenario with stationary axes of two “absolute” rotations implying transient “relative” rotations about a migrating Euler pole and employs a proper plate‐circuit argument to determine them numerically without resorting to approximations. The latter applies an involved interplay of incremental and accumulating modeling inducing split–join cycles to approximate sinusoidal trajectories as reported to record plates' motion during the Gondwana breakup.

Industrial computed tomography for three-dimensional cultural relic model reconstruction based on L1-αL2+TV norm minimization

Industrial computed tomography (ICT) has emerged as a valuable tool for the conservation of cultural relics since it enables the extraction of information about the internal structure of artifacts without causing any damage. However, sparse and finite angle projection computed tomography (CT) often generates images with a large number of artifacts caused by insufficient data. This paper introduces a new cone-beam CT hybrid model based on compressed sensing (CS) theory and total variance (TV) denoising using the Split Bregman algorithm to address this problem, for the first time using compressed sensing theory for CT reconstruction of cultural relics. The model combines the L1 and L2 norms and optimizes the total variance of the image as the objective function. To enhance the reconstruction of L1-L2, the paper employs the continuation strategy of α and proposes a new combination function to update α. Consequently, the reconstruction results are compared with those of traditional common methods, and the effect of different α update functions is observed on the reconstruction results. The results from the model and real CT images demonstrate that the proposed L1-αL2 + TV method yields superior reconstruction results compared to ART + TV, SART, and L1/L2 algorithms based on objective evaluation indexes such as RMSE, SSIM, and PSNR, the proposed L1-αL2 + TV method achieves SSIM improvement of at least 15.2 %. Moreover, the newly proposed combined function is better at preserving the grain of the image, resulting in a much-improved reconstruction of the contours, making it highly beneficial for subsequent conservation and restoration of heritage. Overall, this proposed L1-αL2 + TV mixture model provides low-noise and high-quality reconstruction results, while maintaining the grain of the original image, making it a great boon for heritage conservation and restoration.

Development and evaluation of a standard for absorbed dose to water from alpha-particle-emitting radionuclides

Targeted radionuclide therapy and brachytherapy with alpha particles has gained significant clinical relevance recently. Absorbed dose traceability to a standard is currently lacking in the dosimetry chain. The short range of alpha particles in water of <100 μm complicates the absorbed dose measurements in the form of significant attenuation and perturbation effects. The aim of this work was to develop and validate a standard for absorbed dose to water from alpha-particle-emitting radionuclides. A dosimetric formalism to realize surface absorbed dose to water per unit activity using a windowless cylindrical parallel-plate ion chamber (IC) was introduced. IC-based and Monte Carlo (MC)-based correction factors were calculated for a planar circular 210Po alpha-particle emitter. The measured absorbed dose to air was compared to the MC-calculated absorbed dose. A parallel-plate IC with a nominal 4 mm collector diameter composed of polystyrene-equivalent material was utilized as a standard. MC simulations were performed using the TOPAS MC code and finite source size, backscatter, and emission angle divergence correction factors were calculated by modeling the IC and the source. Multiple measurement trials were performed to measure ionization current at air gaps in the 0.3 mm to 0.525 mm range. The proposed dosimetric formalism was employed to calculate the surface absorbed dose to water from a point-like 210Po source. The recombination and polarity correction factors were measured to be <0.50% when a 150 V mm−1 electric field strength was applied. The MC-calculated and measured absorbed dose to air agreed within 2.05%. The finite source size and backscatter corrections were <10% with the emission angle divergence correction being in the 93%–239% range. The surface absorbed dose to water was measured to be 2.8913 × 10−6 Gy s−1 Bq−1 with a combined uncertainty of 4.23% at k = 1. This work demonstrated the ability of a windowless parallel-plate IC as a standard for absorbed dose from alpha-particle-emitting radionuclides.

The signed helical angle: A technique for characterizing midfoot motion during gait

Quantifying motion in the midfoot during gait and other movements is important for a variety of applications, but challenging due to the complexity of the multiple small articulations involved. The most common motion capture based techniques are limited in their ability to characterize the non-planar nature of the midfoot joint axes. In this study we developed a novel Signed Helical Angle (SHA) to quantify midfoot angular displacement. Motion capture data from 40 healthy subjects walking at a controlled speed were used to calculate finite helical axes and angles from a two-segment foot model. Axes were classified as either pronation or supination based on their orientation, and given a sign, thus either adding to or subtracting from the angular displacement. Analysis focused on insights from axis orientation and comparisons to other techniques. Results showed that when transitions were excluded, pronation and supination axes were fairly well clustered in the transverse plane. The resulting SHA midfoot angle waveform was comparable to sagittal plane Euler and helical component waveforms, but with 39% (approximately 3°) greater range of motion in pronation and 25% (approximately 4°) greater in supination, due to the direct measurement of the motion path and the influence of the other planes. The proposed SHA method may provide an intuitive and useful method to analyze midfoot motion for a variety of applications, particularly when interventions cause subtle changes that may be diluted in planar analyses.

The magnetic field clock angle departure in the Venusian magnetosheath and its response to IMF rotation

We investigate the characteristics of interplanetary magnetic field (IMF) draping in the Venusian magnetosheath using both Venus Express (VEX) observations and magnetohydrodynamics simulations. The distributions of magnetosheath field clock angle illustrate the nearly symmetric morphology of draped magnetic field with respect to the solar wind electric field, and the departure of the IMF clock angle is larger at closer distances. Based on VEX data, the sheath field clock angle departures are found to be &lt;45 degrees for 90% of the instances under steady IMF and this parameter can respond almost immediately to the unsteady IMF. We suggest the magnetosheath field just slips around the planet without significant pileup or bending. Our time-dependent simulations indicate that the response time of sheath field to IMF variation is not more than 1 min and it depends on the involved regions of magnetosheath: the timescale in the inner part of magnetosheath adjacent to the induced magnetosphere is longer than that in the outer part. We find this timescale is controlled by the convection velocity in the magnetosheath, emphasizing the magnetohydrodynamic characteristics of the behavior of the sheath field. The finite magnetosheath field clock angle departure and its quick response to IMF variation suggest that the magnetic field clock angle measured within the Venusian magnetosheath can be used as a reasonable proxy for the upstream IMF clock angle.

A Virtual-Flux State Observer-Based Inductance Identification Method for Model Predictive Control of Grid-Tied Inverters With a Finite Phase Angle Set-Based PLL

To reduce the calculation burden of the conventional finite phase angle set-based phase-locked loop (PLL) and enhance the parameter robustness of the model predictive control for grid-tied inverters (GTIs), an improved finite phase angle set-based PLL and a new inductanceidentification method basedonvirtual-flux state observer (VFSO) are proposed in this paper. First, three improved PLLs are proposed based on the principles of the conventional finite phase angle set-based PLL. Second, the proposed PLLs are further improved by redesigning the cost function, which can search the optimal phase angle within the region of [0, π] instead of [0, 2π], reducing the iteration times. Third, a new VFSO with only one observer parameter is proposed, based on which a new inductance identification method is presented. Fourth, by analyzing the characteristics of the proposed VFSO, it is found that a larger observer gain can improve the dynamic convergence speed of the proposed VFSO, while a smaller one can enhance its robustness against grid frequency deviation. So, an observer gain adaptive method is proposed, aiming to improve both the dynamic and steady state control performance of the proposed inductance identification method. Finally, comparative experimental studies verify the effectiveness of the proposed strategies.

Effect of horizontal magnetic field on Küppers–Lortz instability

We investigate the effect of an external horizontal magnetic field on the Küppers–Lortz instability (KLI) in rotating Rayleigh–Bénard convection of Boussinesq fluids using weakly nonlinear theory along with linear theory. By the KLI, we mean the instability where the two-dimensional roll solutions of the system occurring at the onset of convection become unstable against the perturbations by rolls oriented at different angles with the previous one as the rotation rate exceeds a critical value. The governing parameters, namely, the Prandtl number (Pr), the Taylor number (Ta), and the Chandrasekhar number (Q), are varied in the ranges 0.8≤Pr&amp;lt;∞, 0&amp;lt;Ta≤104, and 0≤Q≤104, respectively, by considering the vanishingly small magnetic Prandtl number limit. In the Pr→∞ limit, magnetic field is found to inhibit the KLI by enhancing the critical Taylor number (Tac) for its onset. On the other hand, for finite Prandtl number fluids, the KLI is favored for lower Q, and it is inhibited for higher Q. Interestingly, in the finite Prandtl number range, both KLI and small angle instability are manifested depending on the Prandtl number. No small-angle instability is observed for Pr≥50, and the rotation-induced KLI is inhibited predominantly by the magnetic field, while, for Pr&amp;lt;50, along with the Küppers–Lortz instability, small-angle instability is also observed. However, in this case, the KLI is favored for lower Q, while it is inhibited for higher Q.

Weighted topology optimization of back-up supporting structures of antennas

In a working environment, the supporting structure usually suffers from asymmetric exterior loads with some finite elevation angles, leading to an asymmetric and irregular structure when applying traditional topology methods. To solve this problem, in this article, a weighted topology optimization method based on a shape feature function is presented, with an innovative concept of shape feature control function, which is employed to reshape the stiffness distribution and make the optimized structure symmetric for the manufacturing process. Topology optimization is also performed under fixed connection constraints. By introducing both the shape control and electromagnetic characteristics functions into the topological design domain, a relatively symmetric and regular topology structure is obtained through the weighted topology optimization method under the asymmetric exterior loads, which benefits the following detailed design and manufacturing process for reflector antennas. Numerical simulation demonstrates the effectiveness and advantages of the proposed method.

Reorientation processes of tilted skyrmion and spiral states in a bulk cubic helimagnet Cu2OSeO3

We present a systematic study of tilted spiral states obtained theoretically within the classical Dzyaloshinskii model for magnetic states in cubic non-centrosymmetric ferromagnets. Such tilted spirals are shown to stabilize under the competing effect of cubic and exchange anisotropies inherent to cubic helimagnets. By focusing on the internal structure of these spirals and their field-driven behaviour for different aspect ratios of the anisotropy coefficients, we are able to capture the main features of the experimental findings in a bulk cubic helimagnet Cu2OSeO3and to make a step further towards a complete quantitative model of this chiral magnet. In particular, we show that for strong anisotropy values (which experimentally correspond to low temperatures near zero) there exist an angular separation between the conical and tilted spirals, i.e., the conical spiral flips into a tilted state and immediately composes some finite angle with respect to the field direction. As the anisotropy ratio decreases, such a transition between two spiral states becomes almost continuous and corresponds to higher temperatures at the experiments. In addition, we investigate the field-driven reorientation of metastable skyrmion lattices induced by the competing anisotropies, which may be responsible for some peculiarities at the experimental phase diagrams of Cu2OSeO3.

Unveiling the Rolling to Kayak Transition in Propelling Nanorods with Cargo Trapping and Pumping

Magnetic nanorods driven by rotating fields in water can be rapidly steered along any direction while generating strong and localized hydrodynamic flow fields. Here we show that, when raising the frequency of the rotating field, these nanopropellers undergo a dynamic transition from a rolling to a kayak-like motion due to the increase in viscous drag and acquire a finite inclination angle with respect to the plane perpendicular to the bottom surface. We explain these experimental observations with a theoretical model which considers the nanorod as a pair of ferromagnetic particles hydrodynamically interacting with a close stationary surface. Further, we quantify how efficiently microscopic cargoes can be trapped or expelled from the moving nanorod and use numerical simulations to unveil the generated hydrodynamic flow field. These propulsion regimes can be implemented in microfluidic devices to perform precise operations based on the selective sorting of microscopic cargoes.

A HYBRID SOLVER FOR THE RADIATIVE TRANSFER EQUATION IN NONGRAY COMBUSTION GASES

The discrete ordinates method (DOM) or its variant, the finite angle method (FAM), is a popular solution method for the radiative transfer equation (RTE). Accurate solution to the RTE using the DOM or FAM requires many solid angles (directions) in multidimensional geometry. In combustion gases, where the absorption coefficient oscillates wildly and the RTE must be solved repeatedly, this method becomes computationally intractable. Here, the FAM is hybridized with the P&lt;sub&gt;1&lt;/sub&gt; approximation, which is efficient since it requires solution to a single partial differential equation as opposed to a set of directional RTEs in the FAM. The P&lt;sub&gt;1&lt;/sub&gt; approximation is accurate when the intensity field is fairly isotropic, as evidenced in optically thick media. Hence, the hybridization employs the FAM for optically thin and intermediate spectral intervals and employs the P&lt;sub&gt;1&lt;/sub&gt; approximation for optically thick spectral intervals. The objective of the present study is to determine optimal parameters for hybridization that can provide the best compromise between accuracy and efficiency. Using a narrowband-based box model for carbon dioxide and water vapor, the nongray radiative transfer equation is solved in media with nonuniform properties enclosed in multidimensional enclosures. Two different approaches-cutoff and filter optical thickness-are investigated for hybridizing. Several problems, both two- and three-dimensional, and with and without coupling to other modes of heat transfer are considered. The filter approach was found to be the best choice for prediction of the radiative source and temperature (in the case of a coupled mode), while the cutoff approach was found to be the best for prediction of wall radiative heat fluxes.

Radiative heat transfer calculations using full spectrum k-distribution method for benchmark test cases

In the present work, the full spectrum k-distribution method (FSK) has been adopted to calculate the radiative transfer in the presence of participating gaseous medium within an enclosure. The radiative properties of the medium is obtained from the HITEMP-2010 database. Further, the radiative properties have been assembled into a monotonically increasing function using the FSK. Moreover, a look-up table has been developed for these properties at different thermodynamic states and a multi-dimensional linear interpolation technique for unavailable thermodynamic states of gases. Furthermore, the FSK method is extended for mixture of gases using different mixing models such as superposition, multiplication and hybrid models. The radiation transfer equation (RTE) is solved by finite angle method to calculate the wall heat fluxes and the divergence of radiative heat flux for various test cases in the categories of homogeneous and non-homogeneous medium having different conditions of temperature and mole-fraction. The multiplication mixing model produces most accurate results among the mixture models used here. The results obtained from FSK has been validated against benchmark solution of line by line method (LBL). The FSK method has been also successfully applied to non-homogeneous gaseous medium for single gas or mixture of gases with almost LBL accuracy at extremely less computational cost.

Comparison of spherical harmonics method and discrete ordinates method for radiative transfer in a turbulent jet flame

In this study, we systematically compared the accuracy and computational cost of two popular solution methods for the radiative transfer equation (RTE): the spherical harmonics method (PN) and the discrete ordinates method (DOM). We first investigated convergence characteristics of different orders of PN and DOM in a series of 1D homogeneous configurations with varying optical thicknesses. Both solvers perform better for optically thicker cases. The accuracy of PN methods increases with its order, N, but the gain in accuracy reduces with the increase in N, i.e., improvement of P7 over P5 is less than that of P3 over P1. This decreasing trend becomes more prominent as the optical thickness decreases. On the other hand, DOM’s accuracy increases almost linearly with the increase in the number of ordinates (or polar angles in this study) in all cases. While comparing the directional profile of radiative intensity, both solvers perform better when the radiative intensity is more isotropic. These solvers were then connected with a full spectrum k-distribution (FSK) spectral model and used to perform radiation-coupled simulations of a turbulent jet flame in an axi-symmetric cylindrical domain. Results obtained from P1 to P7 approximations for PN, and 2 × 4, 4 × 4, 4 × 8, 8 × 8 finite angles for DOM are compared with that from an optically thin model, and a reference solution from line-by-line (LBL) photon Monte Carlo (PMC) method. The choice of radiation solver shows a noticeable impact on the temperature distribution of the flame. The PN solvers lead to slightly higher radiant fractions and the DOM solvers lead to slightly lower radiant fractions than the PMC benchmark solution. Finally, the computational costs of each of these solvers are also reported and an intermittent evaluation / time blending scheme to improve the computational efficiency of radiation solvers in radiation-coupled simulations are also demonstrated.

Magnus induced diode effect for skyrmions in channels with periodic potentials

Using a particle based model, we investigate the skyrmion dynamical behavior in a channel where the upper wall contains divots of one depth and the lower wall contains divots of a different depth. Under an applied driving force, skyrmions in the channels move with a finite skyrmion Hall angle that deflects them toward the upper wall for −x direction driving and the lower wall for +x direction driving. When the upper divots have zero height, the skyrmions are deflected against the flat upper wall for −x direction driving and the skyrmion velocity depends linearly on the drive. For +x direction driving, the skyrmions are pushed against the lower divots and become trapped, giving reduced velocities and a nonlinear velocity-force response. When there are shallow divots on the upper wall and deep divots on the lower wall, skyrmions get trapped for both driving directions; however, due to the divot depth difference, skyrmions move more easily under −x direction driving, and become strongly trapped for +x direction driving. The preferred −x direction motion produces what we call a Magnus diode effect since it vanishes in the limit of zero Magnus force, unlike the diode effects observed for asymmetric sawtooth potentials. We show that the transport curves can exhibit a series of jumps or dips, negative differential conductivity, and reentrant pinning due to collective trapping events. We also discuss how our results relate to recent continuum modeling on a similar skyrmion diode system.