Search NASA⌕ Search

SEARCH · Search NASA

Results for “Linear Models”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

ACDC (Automated Campbell Diagram Code) [SWR-26-042]

This application provides a web-based graphical user interface to generating Campbell Diagrams and visualizing mode shapes for OpenFAST turbine models. Determining the aeroelastic stability and dynamic characteristics of wind turbines is a critical step in turbine design and analysis. Historically, extracting natural frequencies and mode shapes from OpenFAST—the industry-standard whole-turbine simulation code—has been a fragmented and tedious process. It required manual model configuration, command-line linearization execution, and complex post-processing via proprietary scripts to handle rotating-frame dynamics. To address these workflow bottlenecks, we present the Automated Campbell Diagram Code (ACDC), an open-source graphical software tool developed by the National Laboratory of the Rockies (NLR) under the DOE-funded Distributed Wind Aeroelastic Modeling (dWAM) project. ACDC streamlines the end-to-end linearization and stability analysis workflow into a single, intuitive cross-platform application. The software guides users through OpenFAST model configuration, definition of operating points, and the automated execution of steady-state trim and linearization simulations. Under the hood, ACDC automates the complex mathematical post-processing steps required for rotating systems, including Multi-Blade Coordinate (MBC) transformations, eigenanalysis, and advanced modal tracking utilizing the Modal Assurance Criterion (MAC) and spectral clustering. Finally, ACDC processes these results to automatically generate Campbell diagrams and features a robust 3D visualization engine to animate full-system mode shapes. By eliminating the reliance on external post-processing environments and manual data manipulation, ACDC significantly accelerates dynamic analysis and lowers the barrier to entry for wind energy researchers and engineers.

Summerville, Brent [National Laboratory of the Roc↗
DOE OSTI2024

Optimal Electrification Using Renewable Energies: Microgrid Installation Model with Combined Mixture k-Means Clustering Algorithm, Mixed Integer Linear Programming, and Onsset Method

Optimal planning and design of microgrids are priorities in the electrification of off-grid areas. Indeed, in one of the Sustainable Development Goals (SDG 7), the UN recommends universal access to electricity for all at the lowest cost. Several optimization methods with different strategies have been proposed in the literature as ways to achieve this goal. This paper proposes a microgrid installation and planning model based on a combination of several techniques. The programming language Python 3.10 was used in conjunction with machine learning techniques such as unsupervised learning based on K-means clustering and deterministic optimization methods based on mixed linear programming. These methods were complemented by the open-source spatial method for optimal electrification planning: onsset. Four levels of study were carried out. The first level consisted of simulating the model obtained with a cluster, which is considered based on the elbow and k-means clustering method as a case study. The second level involved sizing the microgrid with a capacity of 40 kW and optimizing all the resources available on site. The example of the different resources in the Togo case was considered. At the third level, the work consisted of proposing an optimal connection model for the microgrid based on voltage stability constraints and considering, above all, the capacity limit of the source substation. Finally, the fourth level involved a planning study of electrification strategies based mainly on microgrids according to the study scenario. The results of the first level of study enabled us to obtain an optimal location for the centroid of the cluster under consideration, according to the different load positions of this cluster. Then, the results of the second level of study were used to highlight the optimal resources obtained and proposed by the optimization model formulated based on the various technology costs, such as investment, maintenance, and operating costs, which were based on the technical limits of the various technologies. In these results, solar systems account for 80% of the maximum load considered, compared to 7.5% for wind systems and 12.5% for battery systems. Next, an optimal microgrid connection model was proposed based on the constraints of a voltage stability limit estimated to be 10% of the maximum voltage drop. The results obtained for the third level of study enabled us to present selective results for load nodes in relation to the source station node. Finally, the last results made it possible to plan electrification using different network technologies and systems in the short and long term. The case study of Togo was taken into account. The various results obtained from the different techniques provide the necessary leads for a feasibility study for optimal electrification of off-grid areas using microgrid systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Machine Learning-Based Process Control for Injection Molding of Recycled Polypropylene

The increased interest in artificial intelligence in manufacturing has driven the adoption of machine learning to optimize processes and improve efficiency. A key challenge in injection molding is the variability of recycled materials, which affects part quality and processing stability. This study presents a novel closed-loop process control approach for injection molding, leveraging machine learning to adaptively predict processing inputs and quality outcomes. The methodology was tested on five blends of recycled polypropylene (rPP), using artificial neural networks (ANNs), linear regression, and polynomial regression to model the relationships between material properties and process parameters. The dataset was split 80/20 into training and testing sets. The ANN model was implemented using TensorFlow and Keras, with six hidden layers of 32 neurons per layer, ReLU activation, and an Adam optimizer. Empirical tuning and early stopping were used to optimize performance and prevent overfitting. Predictions were evaluated based on mean absolute error (MAE), mean squared error (MSE), and percentage error. The results showed that yield stress, ultimate elongation, and part weight were accurately predicted within a 5% error for linear and polynomial regression models and within a 10% error for the ANN. However, modulus predictions were less reliable, with errors of ~11% for ANN and linear regression and ~40% for polynomial regression, reflecting the inherent variability of this property in rPP blends. Predictions of processing inputs had errors ranging from 3% to 25%, depending on the model and response variable. No single modeling approach was consistently superior across all responses, highlighting the complexity of the relationship between material properties, process parameters, and quality metrics. Overall, the work demonstrates that closed-loop process control, powered by machine learning, can effectively predict key quality parameters in injection molding of recycled materials. The proposed approach can improve process stability and material utilization, facilitating increased adoption of sustainable materials.

Krantz, Joshua↗
DOE OSTI2025

Sparse measurement medical CT reconstruction using multi-fused block matching denoising priors

A major challenge for medical X-ray CT imaging is reducing the number of X-ray projections to lower radiation dosage and reduce scan times without compromising image quality. However these under-determined inverse imaging problems rely on the formulation of an expressive prior model to constrain the solution space while remaining computationally tractable. Traditional analytical reconstruction methods like Filtered Back Projection (FBP) often fail with sparse measurements, producing artifacts due to their reliance on the Shannon-Nyquist Sampling Theorem. Consensus Equilibrium, which is a generalization of Plug and Play, is a recent advancement in Model-Based Iterative Reconstruction (MBIR), has facilitated the use of multiple denoisers are prior models in an optimization free framework to capture complex, non-linear prior information. However, 3D prior modelling in a Plug and Play approach for volumetric image reconstruction requires long processing time due to high computing requirement. Instead of directly using a 3D prior, this work proposes a BM3D Multi Slice Fusion (BM3D-MSF) prior that uses multiple 2D image denoisers fused to act as a fully 3D prior model in Plug and Play reconstruction approach. Our approach does not require training and are thus able to circumvent ethical issues related with patient training data and are readily deployable in varying noise and measurement sparsity levels. In addition, reconstruction with the BM3D-MSF prior achieves similar reconstruction image quality as fully 3D image priors, but with significantly reduced computational complexity. We test our method on clinical CT data and demonstrate that our approach improves reconstructed image quality.

Hossain, Maliha [ORNL]↗
DOE OSTI2025

The SRG/eROSITA All-Sky Survey: Constraints on the structure growth from cluster number counts

Recent advancements in methods used in wide-area surveys have demonstrated the reliability of the number density of galaxy clusters as a viable tool for precision cosmology. Beyond testing the current cosmological paradigm, cluster number counts can also be used to investigate the discrepancies currently affecting cosmological measurements. In particular, cosmological studies based on cosmic shear and other large-scale structure probes routinely find a value for the amplitude of the fluctuations in the universe S​ 8 = σ 8 (Ω m /0.3) 0.5 smaller than the one inferred from the primary cosmic microwave background. In this work, we investigate this tension by measuring structure evolution across cosmic time as probed by the number counts of massive halos with the first SRG/eROSITA All-Sky Survey cluster catalog in the western Galactic hemisphere, complemented with the overlapping Dark Energy Survey Year-3, Kilo-Degree Survey, and Hyper Suprime-Cam data for weak lensing mass calibration, by implementing two different parameterizations and a model-agnostic method. In the first model, we measured the cosmic linear growth index as γ = 1.19 ± 0.21, which is in tension with the standard value of γ = 0.55 but in good statistical agreement with other large-scale structure probes. The second model is a phenomenological scenario in which we rescale the linear matter power spectrum at low redshift to investigate a potential reduction of structure formation, and it provided similar results. Finally, in a third strategy, we considered a standard ΛCDM cosmology, but we separated the cluster catalog into five redshift bins, measuring the cosmological parameters in each and inferring the evolution of the structure formation, finding hints of a reduction. Interestingly, the S ​8 value inferred from the number counts of the cluster eRASS1 when we add a degree of freedom to the matter power spectrum recovers the value inferred by cosmic shear studies. The observed reduction in the growth rate or systematic uncertainties associated with various measurements may account for the discrepancy in the S 8 ​ values suggested between cosmic shear probes and eROSITA cluster number counts and Planck CMB measurements.

cosmological parameters↗
DOE OSTI2024

Implementation of the Glued Sphere Discrete Element Model for Non-Spherical Particles in MFiX Software

To enhance solver capabilities, simulation flexibility and model validation within the MFiX software, the U.S. Department of Energy (DOE) is funding efforts to develop and integrate the glued-sphere discrete element method into the latest version of MFiX as a dedicated computational module. The glued-sphere discrete element method is a numerical technique to depict the behavior of non-spherical particles in granular flows or particulate systems by representing them as a collection of component spheres. These spheres are bonded together to approximate the shape and mechanical/chemical properties of a more complex particle. The method effectively reuses the existing sphere-sphere collision algorithm, interphase momentum and heat transfer calculations utilized in the traditional discrete element method, extending these capabilities to non-spherical particles. Additionally, this method explicitly resolves intra-particle temperature and species distributions. The MFiX glued-sphere computational module includes tools for generating glued sphere configurations, a dedicated solver, and visualization capabilities in post-processing. More specifically within the computational module, collision detection and calculations were first performed on component spheres and then mapped onto non-spherical particles. The linear spring-dashpot model was utilized to simulate the sphere-sphere interactions.

Ke, Renjie↗
DOE OSTI2026

Dihadron azimuthal asymmetry and light-quark dipole moments at the Electron-Ion Collider

Here, we propose a novel method to probe light-quark dipole moments by examining the azimuthal asymmetries between a collinear pair of hadrons in semi-inclusive deep inelastic lepton scattering off an unpolarized proton target at the Electron-Ion Collider. These asymmetries provide a means to observe transversely polarized quarks, which arise exclusively from the interference between the dipole and the Standard Model interactions, thereby depending linearly on the dipole couplings. We demonstrate that this novel approach can enhance current constraints on light-quark dipole operators by an order of magnitude, free from contamination of other new physics effects. Furthermore, it allows for a simultaneous determination of both the real and imaginary parts of the dipole couplings, offering a new avenue for investigating potential CP-violating effects at high energies.

dihadron fragmentation↗
DOE OSTI2026

Flow dynamics and heat transfer in simplified battery energy storage systems with heated battery modules

Large-scale energy storage systems (ESSs) composed of batteries show promise in addressing current energy challenges, but dissipation of generated heat is important. Here, this paper focuses on buoyant convective flows in simplified ESS battery racks. Natural convection is not generally the primary cooling strategy but can be important in abnormal scenarios where there is module overheat or potentially thermal runaway. We use computational fluid dynamics to investigate the flow dynamics and heat transfer mechanisms in a simplified parameterized rack design. Despite its simplicity, this configuration produces many of the relevant features expected in real ESSs without details of module geometry or hardware, allowing broad conclusions independent of manufacture-specific designs. We start by providing visualizations of the flowfield and measurements of entrainment, heat flux, and pressure. To characterize the dependence on the system parameters, we develop an integral-scale analysis of the average temperature equation to highlight the dominant source terms. We use results from this analysis to derive a steady network model composed of simple algebraic expressions to provide first-order predictions of entrainment through the rack. The network model leads to a linear scaling of the Reynolds number based on convective mass flux with respect to the Grashof number based on the heat source. We deduce empirical relationships that relate the heat exchanged between modules using a surface-averaged Nusselt number as a function of the local Reynolds and Rayleigh numbers. Lastly, we investigate how space between the modules and rack in the spanwise direction creates flow bypass, resulting in different flow pathways.

Battery thermal management↗
DOE OSTI2025

Coupled geomechanical investigation of depletion-induced fault reactivation

Fault reactivation during subsurface fluid production pose significant challenges to safe and sustainable resource extraction. Here, this study presents a three-dimensional coupled geomechanical framework to investigate the processes driving fault reactivation, capturing the interactions between reservoir dynamics and geomechanical responses. Verification against theoretical estimations based on linear poroelasticity confirms the model's capacity in representing reservoir background stress responses. However, the study reveals that relying solely on background stress states can underestimate or overestimate fault reactivation potential, emphasizing the importance of including localized stress perturbations such as differential compaction and stress redistribution. Applied to a fault (M1) inspired by the geological characteristics of the Groningen field, the model shows slip initiation at 2965 m depth with 16.0 MPa depletion, aligning with field observations where seismicity occurred at approximately 3 km depth after 15.8 MPa depletion. Parametric studies reveal: (1) inelastic reservoir compaction delays fault reactivation and mitigates fault slip by reducing stress concentration, (2) higher intermediate in-situ stress magnitudes decrease the Coulomb Failure Stress (CFS) increase rate and reduce fault slip, (3) larger fault offsets amplify shear stress near the offset zone, promoting earlier reactivation and longer rupture propagation, and (4) fault permeability significantly influences pressure diffusion, with low-permeability faults leading to sharper stress changes and earlier fault destabilization. These insights highlight the critical role of geological and mechanical parameters in fault reactivation and provide a predictive framework for mitigating induced seismicity risks.

Coupled geomechanical modeling↗
DOE OSTI2025

Causal Directions Matter: How Environmental Factors Drive Convective Cloud Detrainment Heights

This study investigates how environmental factors influence the level of maximum detrainment (LMD) in deep convective clouds. Through a novel application of the Linear Non‐Gaussian Acyclic Model (LiNGAM), we discover causal structures between environmental variables and LMD, observed at six tropical sites operated by the Atmospheric Radiation Measurement (ARM) user facility. LiNGAM effectively identifies causal directions among variables of interest, revealing robust relationships such as those among the lifting condensation level (LCL), level of free convection (LFC), and convective inhibition (CIN), aligning with prior knowledge. Relative humidity is shown to directly influence LMD; however, this relationship exhibits strong nonlinearity and becomes difficult to detect when the contrast between oceanic and continental environments is excluded from the analysis. This study highlights the importance of establishing causal relationships before performing statistical inference.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Deciphering yield modification of hadron-triggered semi-inclusive recoil jets in heavy-ion collisions

Jet quenching is recognized as critical evidence for the existence of the quark-gluon plasma (QGP) and serves as an essential probe to study its transport properties. Measurements of hadron-triggered semi-inclusive recoil jets have gained popularity due to its capability to probe jets over an extended phase space at low transverse momenta (p T ) and large radii. Recent ALICE measurements showed that the I AA , yield ratio of recoil jets between heavy-ion and p+p collisions, rises with jet p T and exceeds unity at high p T , contradicting conventional expectations that jet quenching should result in I AA values less than one. In this contribution, we re-examine the surface bias and study the effects of energy losses for both trigger hadrons and recoil jets on I AA , employing the Linear Boltzmann Transport (LBT) model to simulate jet-medium interactions. Our findings suggest that a large portion of hadrons used for the triggers undergoes substantial energy loss, despite surface bias. In particular, the energy loss of the trigger hadrons elevates the I AA baseline, corresponding to the case of no energy loss for recoil jets, to be greatly larger than unity. This enhancement of the baseline implies that the measured I AA values being larger than unity could still signal jet quenching.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Global stellarator coil optimization with quadratic constraints and objectives

Most present stellarator designs are produced by costly two-stage optimization: the first for an optimized equilibrium, and the second for a coil design reproducing its magnetic configuration. Few proxies for coil complexity and forces exist at the equilibrium stage. Rapid initial state finding for both stages is a topic of active research. Most present convex coil optimization codes use the least square winding surface method by Merkel (NESCOIL), with recent improvements in conditioning, regularization, sparsity, and physics objectives. While elegant, the method is limited to modeling the norms of linear functions in coil current. We present QUADCOIL, a global coil optimization method that targets combinations of linear and quadratic functions of the current. It can directly constrain and/or minimize a wide range of physics objectives unavailable in NESCOIL and REGCOIL, including the Lorentz force, magnetic energy, curvature, field-current alignment, and the maximum density of a dipole array. QUADCOIL requires no initial guess and runs nearly $10$ 2 x faster than filament optimization. Integrating it in the equilibrium optimization stage can potentially exclude equilibria with difficult-to-design coils, without significantly increasing the computation time per iteration. QUADCOIL finds the exact, global minimum in a large parameter space when possible, and otherwise finds a well-performing approximate global minimum. It supports most regularization techniques developed for NESCOIL and REGCOIL. We demonstrate QUADCOIL’s effectiveness in coil topology control, minimizing non-convex penalties, and predicting filament coil complexity with three numerical examples.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Asymmetric jet shapes with two-dimensional jet tomography

Two-dimensional (2D) jet tomography is a promising tool to study jet medium modification in high-energy heavy-ion collisions. It combines gradient (transverse) and longitudinal jet tomography for selection of events with localized initial jet production positions. It exploits the transverse asymmetry and energy loss that depend, respectively, on the transverse gradient and jet path length inside the quark-gluon plasma (QGP). In this study, we employ 2D jet tomography to study medium modification of the jet shape of γ -triggered jets within the linear Boltzmann transport (LBT) model for jet propagation in heavy-ion collisions. Our results show that jets with small transverse asymmetry ( A N n ⃗ ) or small γ -jet asymmetry ( x J γ = p T jet / p T γ ) exhibit a broader jet shape than those with larger A N n ⃗ or x J γ , since the former are produced at the center and go through longer path lengths while the later are off center and close to the surface of the QGP fireball. In events with finite values of A N n ⃗ , jet shapes are asymmetric with respect to the event plane. Hard partons at the core of the jet are deflected away from the denser region while soft partons from the medium response at large angles flow toward the denser part of QGP. Future experimental measurements of these asymmetric features of the jet shape can be used to study the transport properties of jets and medium responses. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Spin dynamics in natural multiferroic pyroxene NaFeSi 2 ⁢O 6

Spin dynamics in the natural mineral aegirine, NaFeSi 2⁢ O 6 , a member of the pyroxene family, was studied by elastic and inelastic neutron scattering. Magnetization and specific-heat measurements as well as single-crystal neutron diffraction maps, taken in the temperature range 2–20 K, confirm two successive magnetic transitions at 8.8 and 5.8 K, consistent with previous studies. The observed spin-wave excitations emerge from the incommensurate magnetic Bragg peaks corresponding to the propagation vector 𝑘 ICM = (0,0.77,0), and extend up to energies of about 1.5 meV. In the low-temperature helical phase, the spin dynamics of the Fe3+ ions is well described by a simple linear spin-wave model. The observed excitations can be modeled using a spin Hamiltonian that includes three primary exchange interactions—intrachain coupling 𝐽 = 0.142⁢(2) ⁢meV, and interchain couplings 𝐽 1 = 0.083⁢(1) ⁢meV and 𝐽 2 = 0.186⁢(1) meV—and an easy-plane anisotropy 𝐷 = 0.020⁢(6)⁢ meV. Furthermore, our results show that no single exchange interaction dominates the spin dynamics. The similar strengths of the intrachain and interchain couplings point to the fact that the magnetic interactions in aegirine are three dimensional rather than confined along one direction. As a result, the system cannot be considered quasi one dimensional, as previously suggested, and calls for further investigations.

Magnetic coupling↗
DOE OSTI2025

Optimal Economic Dispatch and Load-Following Strategies for Nuclear Integrated Energy Systems

The need for distributed and adaptable energy resources that can handle the growing unpredictability in both supply and demand is rising as the power system continues to modernize. In order to satisfy those needs and maintain grid resilience, nuclear power plants can dynamically control their output, despite typically being used as baseload generators. By incorporating energy storage and renewable energy sources, nuclear integrated energy systems are designed to satisfy the electrical and thermal demands of different end-user applications while ensuring flexible power operation. These systems generate revenue by participating in both wholesale and ancillary services electricity markets, as well as commodity markets for various byproducts generated from coupled industrial processes. This study addresses the economic dispatch efficiency of a tightly coupled nuclear integrated energy system comprising a gigawatt-scale light water reactor, commercialized in the U.S., a high-temperature steam electrolysis unit, a district heating network, and specified electrical loads. To demonstrate the nuclear power plant’s flexibility within the day-ahead unit commitment and economic dispatch framework, while maintaining equilibrium even during periods of refueling outages, this paper develops a mixed-integer linear programming framework that models the subsystems and components of its nuclear steam supply system. A systematic comparative analysis of flexible versus baseload nuclear power plant operation under varying levels of renewable energy integration indicates that flexible operation enhances system profitability by more than 18% while also increasing energy storage utilization, improving reactor responsiveness to load fluctuations, and allowing for greater participation across numerous electricity markets.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2026

Comparative Analysis of Pendulum and Flywheel Power Take-Off Mechanisms for Wave Energy Conversion

Here, this study compares three power-take-off (PTO) mechanisms for wave energy converters—a flywheel -type pitch resonator and two pendulum variants (normal and inverted)—all designed to fit within a fixed envelope. Performance is assessed via linear frequency-domain models and bi-conjugate impedance matching under a representative sea state. The inverted-pendulum PTO achieves 27.9% of the incident wave energy, versus 16.5% for the normal pendulum and 4.2% for the flywheel—that is, 6.6 × and 3.9 × improvements, respectively. These gains arise because pendulum-based systems generate higher torque for a given buoy motion and maintain impedance matching across a wider frequency band. Parametric studies of mass, pendulum arm length, spring stiffness, generator friction, and motor constant quantify the impact of key design choices on power output. By locating the pendulum pivot near the buoy's center of gravity and positioning the center of rotation below the COG, surge–pitch coupling is improved and energy transfer is further boosted.

Keow, Alicia Li Jen [Sandia National Laboratories ↗
DOE OSTI2024

Assessing the Reliability Benefits of Energy Storage as a Transmission Asset

Utilizing energy storage solutions to reduce the need for traditional transmission investments has been recognized by system planners and supported by federal policies in recent years. This work demonstrates the need for detailed reliability assessment for quantitative comparison of the reliability benefits of energy storage and traditional transmission investments. First, a mixed-integer linear programming expansion planning model considering candidate transmission lines and storage technologies is solved to find the least-cost investment decisions. Next, operations under the resulting system configuration are simulated in a probabilistic reliability assessment which accounts for weather-dependent forced outages. The outcome of this work, when applied to TPPs, is to further equalize the consideration of energy storage compared to traditional transmission assets by capturing the value of storage for system reliability.

co-optimization↗
DOE OSTI2026

Altered morphology and diffusivity of water confined in MXenes: Machine learning–accelerated computations combined with experiments

Nanoconfined water exhibits unique properties compared to bulk water due to limited quantities, frustrated hydrogen bonding, and surface interactions, which are fundamental for energy storage and transport applications. We integrate machine learning–accelerated ab initio molecular dynamics with x-ray diffraction (XRD) and inelastic neutron scattering (INS) to systematically analyze the thermodynamic and dynamic behavior of water confined between functionalized (-F, -O, and -OH) two-dimensional (2D) Ti 3 C 2 T x MXene layers. As water intercalates between layers, the interlayer spacing exhibits layer-dependent staging characteristics. The water polarization can be flipped by the count and morphology of intercalated molecules interacting with MXene surface groups, resulting in varying electrostatic potential profiles. On the basis of interfacial electrostatic potential, hydrogen bond lifetime, and molecular orientation, we establish a linear combination of exponential model describing water diffusivity. These computational insights align well with experimental x-ray and neutron measurements, suggesting strategies for tuning water morphology and transport by tailoring MXene surface chemistry and water content for electrochemical energy storage and nanofluidic applications.

Tang, Jiawei [Southeast Univ., Nanjing (China)] (O↗