Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical modelling”

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 595 records · Page 33

Reversible structural transformations in supercooled liquid water from 135 to 245 K

Water has many anomalous properties compared to “simple” liquids, and these anomalies are typically enhanced in supercooled water.1-3 While numerous models have been proposed, including the liquid-liquid critical point,4,5 the singularity-free scenario,6 and the stability limit conjecture,1 a molecular-level understanding remains elusive. The main difficulty in determining which, if any, of these models is correct is the limited amount of data in the relevant temperature and pressure ranges. For water at ambient pressures, which is the focus of this work, data is largely missing from 160 – 232 K (“No Man’s Land”) due to rapid crystallization.2,3 Whether rapid crystallization is just an experimental obstacle, or a fundamental problem signaling the inability of water to thermally equilibrate prior to crystallization is also a major unanswered question.5,7 Here, we investigate the structural transformations of transiently-heated, supercooled water with nanosecond time resolution using infrared vibrational spectroscopy. The experiments demonstrate three key results. First, water’s structure relaxes from its initial configuration to a “steady-state” configuration prior to the onset of crystallization over a wide temperature range. Second, water’s steady-state structure can be reproduced by a linear combination of two, temperature-independent structures that correspond to a “high-temperature liquid” and a “low-temperature liquid.” Third, the observed structural changes are reversible over the full temperature range. Taken together, these results show that supercooled water can equilibrate prior to crystallization for temperatures from the homogeneous nucleation temperature, TH ~232 K,3 down to the glass transition temperature (Tg ~ 136 K). Second, the results provide support for the hypothesis that supercooled water can be described as a mixture of two, structurally-distinct, interconvertible liquids from 135 K to 245 K.5,8-18

Kringle, Loni M.↗

Volumetric Rendering on Wavelet-Based Adaptive Grid

Numerical modeling of physical phenomena frequently involves processes across a wide range of spatial and temporal scales. In the last two decades, the advancements in wavelet-based numerical methodologies to solve partial differential equations, combined with the unique properties of wavelet analysis to resolve localized structures of the solution on dynamically adaptive computational meshes, make it feasible to perform large-scale numerical simulations of a variety of physical systems on a dynamically adaptive computational mesh that changes both in space and time. Volumetric visualization of the solution is an essential part of scientific computing, yet the existing volumetric visualization techniques do not take full advantage of multi-resolution wavelet analysis and are not fully tailored for visualization of a compressed solution on the wavelet-based adaptive computational mesh. Our objective is to explore the alternatives for the visualization of time-dependent data on space-time varying adaptive mesh using volume rendering while capitalizing on the available sparse data representation. Two alternative formulations are explored. The first one is based on volumetric ray casting of multi-scale datasets in wavelet space. Rather than working with the wavelets at the finest possible resolution, a partial inverse wavelet transform is performed as a preprocessing step to obtain scaling functions on a uniform grid at a user-prescribed resolution. As a result, a solution in physical space is represented by a superposition of scaling functions on a coarse regular grid and wavelets on an adaptive mesh. An efficient and accurate ray casting algorithm is based just on these coarse scaling functions. Additional details are added during the ray tracing by taking an appropriate number of wavelets into account based on support overlap with the interpolation point, wavelet coefficient magnitude, and other characteristics, such as opacity accumulation (front to back ordering) and deviation from frontal viewing direction. The second approach is based on complementing of wavelet-based adaptive mesh to the traditional Adaptive Mesh Refinement (AMR) mesh. Both algorithms are illustrated and compared to the existing volume visualization software for Rayleigh-Benard thermal convection and electron density data sets in terms of rendering time and visual quality for different data compression of both wavelet-based and AMR adaptive meshes.

Vezolainen, Alexei V.↗

Automated power-following control for nuclear thermal propulsion startup and shutdown using MOOSE-based applications

This paper describes an investigation of various automated control strategies applied to a full-core multiphysics Griffin/Bison/RELAP-7 model of a prototypical nuclear thermal propulsion system. In all cases, control is achieved by actuating control drums based on the demanded power and predicted quantities from the numerical model. One key finding is that hybrid proportional integral derivative controllers – a novel type of controller that uses both power and reactivity predicted signals – can demonstrate a level of performance rivaling that of period-generated control. The former requires parameter tuning, while the latter mostly necessitates providing reactivity coefficients and temperature rates of change, which could be very challenging to accurately measure in real-time. In addition, decay heat plays an important role in determining cooling requirements during cooldown phases. A decay heat model, accounting for burn time and throttling, was derived, and is incorporated within the model to simulate the steady-state and shutdown phases and satisfyingly follow the power demand. However, temperature overshoots and non-constant specific impulse during throttling will necessitate further improvements.

33 ADVANCED PROPULSION SYSTEMS↗

Viscous Behavior of Clay‐Rich Rocks and Its Role in Focused Fluid Flow

Abstract Focused fluid flow is common in sedimentary basins worldwide, where flow structures often penetrate through sandy reservoir rocks, and clay‐rich caprocks. To better understand the mechanisms forming such structures, the impacts of the viscoelastic deformation and strongly nonlinear porosity‐dependent permeability of clay‐rich materials are assessed from an experimental and numerical modeling perspective. The experimental methods to measure the poroviscoelastic and transport properties of intact and remolded shale have been developed, and the experimental data is used to constrain the numerical simulations. It is demonstrated that viscoelastic deformation combined with nonlinear porosity‐dependent permeability triggers the development of localized flow channels, often imaged as seismic chimneys. The permeability inside a channel increases by several orders of magnitude compared to the background values. In addition, the propagation time scale and the channel size strongly depend on the material properties of the fluid and the rock. The time‐dependent behavior of the clay‐rich rock may play a key role in the long‐term integrity of the subsurface formations.

58 GEOSCIENCES↗

Progress Report on Model Development for Aerosol Transport through Divergent Cracks Paths

This report summarizes the progress in developing a phenomenological model of aerosol transport, deposition, and plugging through microchannels. The purpose of this effort is to introduce to a user community—involving researchers, regulators, and industry—a generic, reliable numerical model for the prediction of aerosol transport while accounting for potential deposition and plugging of the leak paths to model spent nuclear fuel (SNF) release from postulated Stress Corrosion Cracks in canisters under storage or transportation. Current work focuses on expanding the model to predict aerosol (and gas) flow through complex microchannel (nozzle) geometries in addition to the rectangular and cylindrical geometries, as we approach more realistic Stress Corrosion Crack conditions. In this regard, a divergent nozzle geometry, used at Sandia National Laboratories (SNL) for testing aerosol release and retention (Durbin et al. 2021), was added to the model. The model was then validated with blowdown data for the particular microchannel, from experiments conducted at SNL. The report also presents preliminary non-benchmarked aerosol penetration fraction and mass flow rate for a monodisperse 10-micron (AED) particle concentration of 1.7e-08 kg/m 3 being released through the divergent microchannel geometry from the SNL’s aerosol experimental tank setup.Future work will involve validating the aerosol model and updating the model’s Graphical User Interface (GUI) to include the divergent nozzle geometry. This is expected to help stakeholders perform quick and easy first principles calculations without needing to understand the underlying MATLAB script.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Summary Report of the FY24 DOE Contributions to the GIF VHTR CMVB

The Generation-IV Forum (GIF) Very-High-Temperature Reactor-Computational Methods Validation and Benchmark (VHTR-CMVB) initiative, involving organizations from Korea Atomic Energy Research Institute (KAERI) (South Korea), Institute of Nuclear and New Energy Technology of Tsinghua University (INET) (China), U.S. Department of Energy (DOE) (U.S.), Joint Research Centre (JRC) (Europe), and Japan Atomic Energy Agency (JAEA) (Japan), is dedicated to the verification and validation of tools for High-Temperature Gas-Cooled Reactors (HTGRs) analysis, using data shared by Computational Methods Validation and Benchmark (CMVB) signatories. For FY24, the US DOE CMVB has committed to several critical activities. Under WP1, led by the US, the integration of the High Temperature Gas Cooled Reactor - Pebble-Bed Module (HTR-PM) Phenomena Identification and Ranking Table (PIRT) into the comparative PIRT is progressing, with a new draft of the comparison tables issued earlier this year and currently being utilized by INET for their contribution. Neutronic validation efforts under WP3 include the preparation of the burnup analysis benchmark, preliminary calculations, and the development of reference models and results. In WP2, a validation exercise for hot gas mixing in the lower plenum of HTR-PM is in progress, using experimental data from INET (China) to validate modeling approaches. A model of the experimental facility has been developed using StarCCM+, with initial calculations slated for presentation at the GIF CMVB meeting this fall. Another WP2 activity focuses on validating numerical models for air-cooled Reactor Cavity Cooling System (RCCS) with experimental data from the Wisconsin Madison RCCS facility. A high-fidelity model, developed using NEK-RS, is currently being validated with available data from a low power forced convection test. These efforts are aimed at enhancing and confirming the accuracy of HTGR analysis tools, ensuring their alignment with experimental data and regulatory requirements.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Analysis of heat transfer and AuNPs-mediated photo-thermal inactivation of E. coli at varying laser powers using single-phase CFD modeling

In the wake of the COVID-19 pandemics, the demand for innovative and effective methods of bacterial inactivation has become a critical area of research, providing the impetus for this study. The purpose of this research is to analyze the AuNPs-mediated photothermal inactivation of E. coli. Gold nanoparticles irradiated by laser represent a promising technique for combating bacterial infection that combines high-tech and scientific progress. The intermediate aim of the work was to present the calibration of the model with respect to the gold nanorods experiment. The purpose of this work is to study the effect of initial concentration of E. coli bacteria, the design of the chamber and the laser power on heat transfer and inactivation of E. coli bacteria. Using the CFD simulation, the work combines three main concepts. 1. The conversion of laser light to heat has been described by a combination of three distinctive approximations: a- Discrete particle integration to take into account every nanoparticle within the system, b- Rayleigh-Drude approximation to determine the scattering and extinction coefficients and c- Lambert–Beer–Bourger law to describe the decrease in laser intensity across the AuNPs. 2. The contribution of the presence of E. coli bacteria to the thermal and fluid-dynamic fields in the microdevice was modeled by single-phase approach by determining the effective thermophysical properties of the water-bacteria mixture. 3. An approach based on a temperature threshold attained at which bacteria will be inactivated, has been used to predict bacterial response to temperature increases. The comparison of the thermal fields and temporal temperature changes obtained by the CFD simulation with those obtained experimentally confirms the accuracy of the light-heat conversion model derived from the aforementioned approximations. The results show a linear relationship between maximum temperature and variation in laser power over the range studied, which is in line with previous experimental results. It was also found that the temperature inside the microchamber can exceed 55 °C only when a laser power higher than 0.8 W is used, so bacterial inactivation begins. The experimental data allows to determinate the concentration of nanoparticles. This parameter is introduced into the mathematical model obtaining the same number of AuNPs. However, this assumption introduces a certain simplification, as in the mathematical model the distribution of nanoparticles is uniform. This work is directly connected to the use of gold nanoparticles for energy conversion, as well as the field of bacterial inactivation in microfluidic systems such as lab-on-a-chip. Presented mathematical and numerical models can be extended to the entire spectrum of wavelengths with particular use of white light in the inactivation of bacteria. This work represents a significant advancement in the field, as to the best of the authors’ knowledge, it is the first to employ a single-phase computational fluid dynamics (CFD) approach specifically combined with the thermal inactivation of bacteria. Moreover, this research pioneers the use of a numerical simulation to analyze the temperature threshold of photothermal inactivation of E. coli mediated by gold nanorods (AuNRs). The integration of these methodologies offers a new perspective on optimizing bacterial inactivation techniques, making this study a valuable contribution to both computational modeling and biomedical applications.

36 MATERIALS SCIENCE↗

Preliminary Thermal-Hydraulic Analysis of the UTA–2 Subcritical Linear Accelerator–Driven System

The National Nuclear Security Administration’s mission includes establishing a reliable supply of 99 Mo without highly enriched uranium. Oak Ridge National Laboratory (ORNL) supports this objective through collaborative research and development with industrial partners. Niowave Inc., a current partner, is currently designing a subcritical linear accelerator-driven system (ADS) and preparing for the US Nuclear Regulatory Commission’s licensing process. Niowave’s system has the potential to efficiently supply medical radioisotopes. The technology includes a superconducting electron accelerator and a pile of both natural and low-enriched uranium (LEU) targets. The process fissions uranium and many valuable isotopes can then be extracted from the targets. Niowave is currently iterating through conceptual and detailed design processes for several system sizes. This report discusses UTA–2, which is at the demonstration stage. UTA–2 will validate numerical modeling results with experimental measurements before progressing to the detailed design of UTA–3, the commercial-sized ADS. The thermal-hydraulic behavior of the UTA–2 core design was numerically investigated using STAR-CCM+, as described in this report. STAR-CCM+, a state-of-the-art computational fluid dynamics (CFD) software that was commercially developed by Siemens, has an extensive user base and a set of validation studies. It is also compliant with the American Society of Mechanical Engineers’ Nuclear Quality Assurance 1 standard. Three cases are investigated in this report. Case 1 quantified the temperature field within the UTA–2 assembly and water tank using only conduction as the method of thermal energy transport. This simplified approach was overly conservative and yielded wetted cladding temperatures above the coolant saturation temperature. In this case, the maximum temperature of the wetted cladding surface of the highest power rod exceeded the saturation temperature of water by 63.8°C. Because of the overly conservative approach taken in Case 1 and its negative subcooled margin, buoyancy-driven natural circulation flow physics were implemented in Case 2. Adding coolant motion significantly distributed the thermal energy of the system through convective heat transfer. This relatively small amount of convective heat transfer significantly reduced system temperatures and increased the subcooled margin from –63.8 to 61.3°C. This margin confirmed that no boiling was expected during normal operation of UTA–2 at 230 W. Case 3 had no additional physics models but considered an overpower event during which the power of each LEU and natural uranium rod was at its respective peak values. This resulted in a study with the total assembly power equal to 176% of the nominal power of 230 W considered in Cases 1 and 2. The resulting natural circulation flows were slightly enhanced. The subcooled margin decreased slightly to 46.3°C, which is still a significant margin to local boiling of the water within the tank. This margin confirmed that no boiling was expected during an abnormal operation of UTA–2 at 406 W.

07 ISOTOPE AND RADIATION SOURCES↗

Confinement/deconfinement transition in the D0-brane matrix model — A signature of M-theory?

We study the confinement/deconfinement transition in the D0-brane matrix model (often called the BFSS matrix model) and its one-parameter deformation (the BMN matrix model) numerically by lattice Monte Carlo simulations. Our results confirm general expectations from the dual string/M-theory picture for strong coupling. In particular, we observe the confined phase in the BFSS matrix model, which is a nontrivial consequence of the M-theory picture. We suggest that these models provide us with an ideal framework to study the Schwarzschild black hole, M-theory, and furthermore, the parameter region of the phase transition between type IIA superstring theory and M-theory. A detailed study of M-theory via lattice Monte Carlo simulations of the D0-brane matrix model might be doable with much smaller computational resources than previously expected.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Modelling of the laser amplification process with allowance for the effect of the temperature distribution in an Yb : YAG gain element on the thermophysical and lasing characteristics of the medium

A time-dependent three-dimensional model for the laser amplification process has been constructed with allowance for the effect of the temperature distribution on the thermophysical and lasing characteristics of gain media. We have performed numerical modelling of the laser amplification process in the gain elements of a two-stage subjoule-level cryogenic laser amplifier operating at a pulse repetition rate of up to 1 kHz. It has been shown that taking into account the temperature distribution is of critical importance in calculation of cryogenically cooled laser amplifiers pumped with high-power diodes. We have found optimal diode pump parameters at which the maximum achievable pulse energy at the amplifier output can reach 300 and 570 mJ at pulse repetition rates of 1000 and 500 Hz, respectively. (extreme light fields and their interaction with matter)

42 ENGINEERING↗

Closing the loop between in-situ stress complexity and EGS fracture complexity

The goal of the project is to employ a combination of high-fidelity simulations and true-triaxial block fracturing tests at high temperature to explore the intricate relationship between in-situ stress and hydraulic fracture patterns and better characterize the in-situ stress at Utah FORGE. Laboratory experiments are used to validate numerical models and better understand the impact of various parameters (e.g., temperature, well orientation, stress, etc) on the hydraulic fracture nucleation process in EGS. The project will have a significant impact as it will (1) improve the characterization of the in-situ stress field at FORGE; (2) demonstrate how high-fidelity modeling tools can be employed to analyze data for EGS; (3) provide a unique set of high-temperature hydraulic fracturing experimental results for EGS to help identifying the most important components for in-situ stress estimation and to validate the current in-situ stress estimation theory and establish a strong foundation for proposing a new one; (4) provide a validated set of numerical tools in an open source simulation framework (GEOS) that can be used to model EGS.

15 GEOTHERMAL ENERGY↗

Methods and devices for preventing computationally explosive calculations in a computer for model parameters distributed on a hierarchy of geometric simplices

A computer-implemented method of preventing computationally explosive calculations. The method includes obtaining, by a processor of the computer, measured data of one of a physical process or a physical object; performing hierarchical numerical modeling of a physical process inclusive of an Earth model containing at least one of (a) infrastructure in the ground and (b) a formation feature in the ground, wherein predicted data is generated; comparing the measured data to the predicted data to calculate an estimated error; analyzing the estimated error via an inversion process to update the at least one of the Earth model and infrastructure model so as to reduce the estimated error and to determine a final composite Earth model of at least one of the infrastructure and the feature; and using the final composite Earth model to characterize at least one of the process and the physical object.

Source record↗

A two-fluid single-column model of turbulent shallow convection. Part II: Single-column model formulation and numerics

The two-fluid single-column model of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) is extended to include moisture and horizontal wind shear. Turbulent kinetic energy is introduced as a prognostic variable, dependence on a diagnosed boundary-layer height is removed, and subfilter fluxes are approximated using a two-fluid version of a Mellor–Yamada scheme. Three mechanisms for entrainment and detrainment processes are introduced, which represent entrainment of unstable air at the surface, forced detrainment of air at the top of the boundary/cloud layers, and turbulent mixing that relaxes the convective fluid to a reference profile. A semi-implicit Eulerian discretization replaces the semi-implicit semi-Lagrangian implementation of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) to improve numerical stability and conservation. The equations for the implicit time step are solved using a quasi-Newton method, which is shown to perform well in numerical tests for conservation and convergence. The two-fluid single-column model presented in this article will be applied to simulations of shallow cumulus convection in Part III.

54 ENVIRONMENTAL SCIENCES↗

Evaluating Climate Models with the CLIVAR 2020 ENSO Metrics Package

El Niño–Southern Oscillation (ENSO) is the dominant mode of interannual climate variability on the planet, with far-reaching global impacts. It is therefore key to evaluate ENSO simulations in state-of-the-art numerical models used to study past, present, and future climate. Recently, the Pacific Region Panel of the International Climate and Ocean: Variability, Predictability and Change (CLIVAR) Project, as a part of the World Climate Research Programme (WCRP), led a community-wide effort to evaluate the simulation of ENSO variability, teleconnections, and processes in climate models. The new CLIVAR 2020 ENSO metrics package enables model diagnosis, comparison, and evaluation to 1) highlight aspects that need improvement; 2) monitor progress across model generations; 3) help in selecting models that are well suited for particular analyses; 4) reveal links between various model biases, illuminating the impacts of those biases on ENSO and its sensitivity to climate change; and to 5) advance ENSO literacy. By interfacing with existing model evaluation tools, the ENSO metrics package enables rapid analysis of multipetabyte databases of simulations, such as those generated by the Coupled Model Intercomparison Project phases 5 (CMIP5) and 6 (CMIP6). The CMIP6 models are found to significantly outperform those from CMIP5 for 8 out of 24 ENSO-relevant metrics, with most CMIP6 models showing improved tropical Pacific seasonality and ENSO teleconnections. Only one ENSO metric is significantly degraded in CMIP6, namely, the coupling between the ocean surface and subsurface temperature anomalies, while the majority of metrics remain unchanged.

54 ENVIRONMENTAL SCIENCES↗

On measuring stimulated photon–photon scattering using multiple ultraintense lasers

Stimulated photon–photon scattering is a predicted consequence of quantum electrodynamics that has yet to be measured directly. Measuring the cross section for stimulated photon–photon scattering is the aim of a flagship experiment for NSF OPAL, a proposed laser user facility with two, 25-PW beamlines. We present optimized experimental designs for achieving this challenging and canonical measurement. A family of experimental geometries is identified that satisfies the momentum- and energy-matching conditions for two selected laser frequency options. Numerical models predict a maximum signal exceeding 1000 scattered photons per shot at the experimental conditions envisaged at NSF OPAL. Experimental requirements on collision geometry, polarization, cotiming and copointing, background suppression, and diagnostic technologies are investigated numerically. These results confirm that a beam cotiming shorter than the pulse duration and control of the copointing on a scale smaller than the shortest laser wavelength are needed to robustly scatter photons on a per-shot basis. Finally, we assess the bounds that a successful execution of this experiment may place on the mass scale of Born–Infeld nonlinear electrodynamics beyond the standard model of physics.

Beyond the Standard Model↗

Geomechanical properties of the Meletta sandstone - the high-temperature heat storage reservoir rock of DeepStor

The DeepStor project aims at storing excess heat at temperatures up to 140 °C in the depleted Leopoldshafen oil field at a depth of about 1300 m. In order to gain knowledge on the target horizons, the different layers of the Meletta sandstone, samples cored in a block retrieved in a quarry near Nussloch where the Meletta sandstone outcrops were studied. Several petrophysical properties were investigated including mineralogy, porosity, permeability, thermal conductivity, P- and S-wave velocities and quality factor. A mechanical study focused on compressive strength under uniaxial and triaxial stress conditions, tensile strength, critical pressure and stress-dependence of physical properties. Our results show that the Meletta sandstone is heterogeneous, anisotropic, mechanically weak, stress-sensitive and prone to water weakening. Petrophysical measurements on few Meletta sandstone cores retrieved in boreholes at about 1250 m depth showed that the outcrop samples are significantly less cohesive and weaker. Therefore, one should consider the Nussloch outcrop as a good reservoir analog only with great caution. Several models were applied to account for our experimental results on anisotropy, mechanical behavior and stress-dependence of the outcrop samples. Our study provides a valuable data set which can feed numerical models to simulate the behavior of the host formation during heat exchange operations in the DeepStor project.

Geomechanics↗

Hyper‐Local Temperature Prediction Using Detailed Urban Climate Informatics

The accurate modeling of urban microclimate is a challenging task given the high surface heterogeneity of urban land cover and the vertical structure of street morphology. Recent years have witnessed significant efforts in numerical modeling and data collection of the urban environment. Nonetheless, it is difficult for the physical‐based models to fully utilize the high‐resolution data under the constraints of computing resources. The advancement in machine learning (ML) techniques offers the computational strength to handle the massive volume of data. In this study, we proposed a modeling framework that uses ML approach to estimate point‐scale street‐level air temperature from the urban‐resolving meso‐scale climate model and a suite of hyper‐resolution urban geospatial data sets, including three‐dimensional urban morphology, parcel‐level land use inventory, and weather observations from a sensor network. We implemented this approach in the City of Chicago as a case study to demonstrate the capability of the framework. The proposed approach vastly improves the resolution of temperature predictions in cities, which will help the city with walkability, drivability, and heat‐related behavioral studies. Moreover, we tested the model's reliability on out‐of‐sample locations to investigate the modeling uncertainties and the application potentials to the other areas. This study aims to gain insights into next‐gen urban climate modeling and guide the observation efforts in cities to build the strength for the holistic understanding of urban microclimate dynamics.

54 ENVIRONMENTAL SCIENCES↗

Topological solitonic macromolecules

Being ubiquitous, solitons have particle-like properties, exhibiting behaviour often associated with atoms. Bound solitons emulate dynamics of molecules, though solitonic analogues of polymeric materials have not been considered yet. Here we experimentally create and model soliton polymers, which we call “polyskyrmionomers”, built of atom-like individual solitons characterized by the topological invariant representing the skyrmion number. With the help of nonlinear optical imaging and numerical modelling based on minimizing the free energy, we reveal how topological point defects bind the solitonic quasi-atoms into polyskyrmionomers, featuring linear, branched, and other macromolecule-resembling architectures, as well as allowing for encoding data by spatial distributions of the skyrmion number. Application of oscillating electric fields activates diverse modes of locomotion and internal vibrations of these self-assembled soliton structures, which depend on symmetry of the solitonic macromolecules. Our findings suggest new designs of soliton meta matter, with a potential for the use in fundamental research and technology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗