Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “embedded systems”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Hierarchical Model-Free Transactional Control of Building Loads to Support Grid Services

A transition from generation on demand to consumption on demand is one of the solutions to overcome the many limitations associated with the higher penetration of renewable energy sources. Such a transition, however, requires a considerable amount of load flexibility in the demand side. Demand response (DR) programs can reveal and utilize this demand flexibility by enabling the participation of a large number of grid-interactive efficient buildings (GEB). Existing approaches on DR require significant modelling or training efforts, are computationally expensive, and do not guarantee the satisfaction of end users. To address the aforementioned limitations, this software code considers a scalable hierarchical model-free transactional control approach that incorporates elements of virtual battery, game theory, and model-free control (MFC) mechanisms. The developed approach separates the control mechanism into upper and lower levels. The MFC modulates the flexible GEB in the lower level with guaranteed thermal comfort of end users in response to the optimal pricing and power signals determined in the upper layer using a Stackelberg game integrated with aggregate virtual battery constraints. Additionally, the usage of MFC necessitates less burdensome computational and communication requirements, thus, it is easily deployable even on small, embedded devices. This software code enables the use of residential and small-size commercial buildings to offer a potentially substantial source of ancillary grid services that are currently underutilized. A hierarchical, model-free transactive building control provides a smooth interface between the grid service requests of utilities and the reliable control required by participating buildings. Model-free control, which supports distributed control architecture, will be a more scalable solution that can be deployed in neighborhood-size systems as well as individual buildings. The proposed approach contributes to the body of knowledge on two main aspects. First, it couples the MFC with the game-theoretic control and proposes a scalable model-free transactive control approach. MFC does not require any modeling effort or model training for the various building loads. This is very beneficial since deriving an accurate model for every single unit participating in DR programs and obtaining all the parameters about the units (e.g., thermal coefficients, standby losses) are infeasible. Also, MFC is computationally efficient, easily deployable even on small, embedded devices, and can be implemented in real time. Second, it integrates the concept of virtual battery into DR via the Stackelberg game. The concept of virtual battery enables efficient coordination and aggregation of a large number of flexible GEB with guaranteed thermal comfort of end users.

Olama, Mohammed [Oak Ridge National Lab. (ORNL), O↗
DOE OSTI2021

Nonlinear convergence in contact mechanics: Immersed boundary finite volume

In this report we present an immersed boundary finite volume (IBM) method for simulating quasistatic contact mechanics of linearly elastic domains at small strains. In IBM, all external boundaries and internal contacts of an object are represented by embedded surfaces inside a Cartesian mesh, which need not conform to the grid lines. The contact constraints consist of the non-penetrability condition and Coulomb’s friction law, which are discretized using special interpolation stencils and enforced via penalty parameters. The resulting nonlinear system depends on displacement unknowns only. To solve it, we use the Newton method but find that it diverges frequently. To understand the divergence pattern, we analyze a simplified 2-cell problem and show that the global convergence of Newton cannot be ensured for any choice of penalty parameters. We thus propose a modified Newton solver, which guarantees convergence for the 2-cell problem and is numerically verified to converge for all the challenging simulations considered herein. While both 1 st - and 2 nd -order variants of IBM, in displacement unknowns, are proposed, the modified Newton solver applies only to the 1 st -order variant.

42 ENGINEERING↗
DOE OSTI2022

Explicit quantification of residence and flushing times in the Salish Sea using a sub-basin scale shoreline resolving model

The Salish Sea, located in the Pacific Northwest region of North America has complex currents and circulation features distributed over numerous interconnected deep basins with islands. Increased risk of exposure to oil spills and untreated wastewater from maritime emergencies and treatment plant failures have led to a need for quantifying residence and flushing characteristics at a sub-basin scale using a shoreline resolving hydrodynamic model. An unstructured grid model of the Salish Sea was developed using FVCOM with a ≈75-100m shoreline resolution. In addition to 7 tides and 23 salinity and temperature monitoring stations, an extensive currents data set from 135 stations collected over a span of three years was used for skill assessment and validation. Explicit forward computations were then conducted to define and quantify residence and flushing times in various sub-basins of interest using (a) Lagrangian particles and (b) Numerical/virtual dye experiments. The results in most basins show expected seasonal variability with longer flushing time associated with summer lower tides and lower freshwater inflows. However, contrary to expectation, flushing time is significantly longer in wintertime in large fjord-like basins such as Hood Canal (≈138 days), likely due to increased stratification and reduced mixing. The flushing time for the Puget Sound region of the Salish Sea is ≈115 days, while Georgia Basin is 240 days when analyzed as stand-alone basins with zero background concentrations. When examined as part of the flushing of the entire Puget Sound filled with virtual dye, the compounded flushing times for embedded sub-basins can be significantly longer in order of magnitudes and largely dictated by the flushing time of Puget Sound. The computed residence and flushing time scales tabulated over 36 sub-basins provide an improved understanding of water renewal in the system, informing pollution management actions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Modeling atmospheric data and identifying dynamics Temporal data-driven modeling of air pollutants

Atmospheric modelling has recently experienced a surge with the advent of deep learning. Most of these models, however, predict concentrations of pollutants following a data-driven approach in which the physical laws that govern their behaviors and relationships remain hidden. Furthermore, with the aid of real-world air quality data collected hourly in different stations throughout Madrid, we present a case study using a series of data-driven techniques with the following goals: (1) Find systems of ordinary differential equations that model the concentration of pollutants and their changes over time; (2) assess the performance and limitations of our model using stability analysis; (3) reconstruct the time series of chemical pollutants not measured in certain stations using delay coordinate embedding results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Controlling interactions between high-frequency phonons and single quantum systems using phononic crystals

The ability to control phonons in solids is key in many fields of quantum science, ranging from quantum information processing to sensing. Phonons often act as a source of noise and decoherence when solid-state quantum systems interact with the phonon bath of their host matrix. In this study, we demonstrate the ability to control the phononic local density of states of the host matrix using phononic crystals and measure its positive impact on single quantum systems. We design and fabricate diamond phononic crystals with features down to around 20 nm, resulting in a high-frequency complete phononic bandgap from 50 to 70 GHz. The engineered local density of states is probed using single silicon-vacancy colour centres embedded in the phononic crystals. We observe an 18-fold reduction in the phonon-induced orbital relaxation rate of the emitters compared to bulk, thereby demonstrating that the phononic crystal suppresses spontaneous single-phonon processes. Furthermore, we show that our approach can efficiently suppress single-phonon-emitter interactions up to 20 K, allowing the investigation of multi-phonon processes in the emitters. Our results represent an important step towards the realization of efficient phonon-emitter interfaces that can be used for quantum acoustodynamics and quantum phononic networks.

quantum information↗
DOE OSTI2024

Multi-resolution partial differential equations preserved learning framework for spatiotemporal dynamics

Traditional data-driven deep learning models often struggle with high training costs, error accumulation, and poor generalizability in complex physical processes. Physics-informed deep learning (PiDL) addresses these challenges by incorporating physical principles into the model. Most PiDL approaches regularize training by embedding governing equations into the loss function, yet this depends heavily on extensive hyperparameter tuning to weigh each loss term. To this end, we propose to leverage physics prior knowledge by “baking” the discretized governing equations into the neural network architecture via the connection between the partial differential equations (PDE) operators and network structures, resulting in a PDE-preserved neural network (PPNN). This method, embedding discretized PDEs through convolutional residual networks in a multi-resolution setting, largely improves the generalizability and long-term prediction accuracy, outperforming conventional black-box models. The effectiveness and merit of the proposed methods have been demonstrated across various spatiotemporal dynamical systems governed by spatiotemporal PDEs, including reaction-diffusion, Burgers’, and Navier-Stokes equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Wet chemistry-synthesized Fe/mixed ferrite soft magnetic composites for high-frequency power conversion

State-of-the-art soft magnetic (SM) alloy systems such as electrical steels and permalloys exhibit large eddy current losses at high-frequency (kHz range and above), limiting their application in fast switching devices. Eddy current losses can be reduced for metallic alloys through a decrease in characteristic length scale or embedding them in an insulation layer. This work proposes the development of fine-scale core materials consisting of metallic SM nano/microparticles with magnetic, inorganic insulation layers, and their synthesis using a wet chemistry-based, scalable method for high frequency and high-power applications. More specifically, a magnetic ferrite coating (Ni 0.5-x Mn x Zn 0.5 Fe 2 O 4 , x = 0.0–0.5) was applied via deposition of ferrite powder produced using a “wet chemistry-based” co-precipitation method; the ferrite was distributed through the Fe micropowder via either ball milling or ultrasonic mixing. Powder cores were prepared by compaction of the synthesized composites. Relative magnetic permeability and core loss were measured at excitation frequencies to 200 kHz. A core loss of 127 kW/m 3 was measured at excitation frequency of 100 kHz and magnetization of 0.02 T. This value improves on the 199.3 kW/m 3 reported in literature at identical excitation conditions for a compact formed from a composite comprising Fe microparticles coated with low-permeability Ni 0.5 Zn 0.5 Fe 2 O 4 . Ultrasonic mixing resulted in slightly lower core loss than ball milling, possibly because ball milling causes loss-increasing deformation of the Fe micropowder. XRD and SEM were used to observe composite composition and core cross-section microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Detecting Masquerade Attacks in Controller Area Networks Using Graph Machine Learning

Modern vehicles rely on a myriad of electronic control units (ECUs) interconnected via controller area networks (CANs) for critical operations. Despite their ubiquitous use and reliability, CANs are susceptible to sophisticated cyberattacks, particularly masquerade attacks, which inject false data that mimic legitimate messages at the expected frequency. These attacks pose severe risks such as unintended acceleration, brake deactivation, and rogue steering. Traditional intrusion detection systems (IDS) often struggle to detect these subtle intrusions due to their seamless integration into normal traffic. This paper introduces a novel framework for detecting masquerade attacks in the CAN bus using graph machine learning (ML). We hypothesize that the integration of shallow graph embeddings with time series features derived from CAN frames enhances the detection of masquerade attacks. We show that by representing CAN bus frames as message sequence graphs (MSGs) and enriching each node with contextual statistical attributes from time series, we can enhance detection capabilities across various attack patterns compared to using graph-based features only. Our method ensures a comprehensive and dynamic analysis of CAN frame interactions, improving robustness and efficiency. Extensive experiments on the ROAD dataset validate the effectiveness of our approach, demonstrating statistically significant improvements in the detection rates of masquerade attacks compared to a baseline that uses graph-based features only as confirmed by Mann-Whitney U and Kolmogorov-Smirnov tests (p < 0.05) .

Marfo, William [Univ. of Texas, El Paso, TX (Unite↗
DOE OSTI2023

Development of Optical Fiber-Based Sensors for Nuclear Microreactor Structural Health Monitoring

This report provides an experimental assessment of two different optical fiber–based acoustic sensors that are being investigated for application in nuclear microreactors to enhance structural health monitoring capabilities. Optical fibers are resilient in high-temperature and high-radiation environments, have a small sensor footprint, are immune to electromagnetic interference, and are capable of spatially distributed sensing. The two sensors investigated here are (1) Fabry–Pérot Cavities (FPCs) between two copper-coated fibers, embedded in nickel capillary tubes and (2) type-I fiber Bragg grating (FBG) arrays contained within metal capillary tubes. These sensors can be interrogated using low-coherence interferometry or swept wavelength interferometry, respectively, to measure the resonant frequencies of the components or systems to which these sensors are bonded. The FPC developed herein has been subjected to temperatures up to nearly 800°C while tack-welded to a tubular test specimen. Even at the highest temperatures, the measured resonant frequencies compared well with those obtained using an accelerometer that was bonded to an unheated portion of the specimen. The FBG array was tested in multiple bonding configurations to a tubular test specimen, all at room temperature, with the understanding that high temperature (i.e., type II) FBGs could be used to obtain similar data at high temperatures; the FBG array data were validated with noncontact laser vibrometry, which is being used at Los Alamos National Laboratory to relate acoustic signatures to component stresses and/or structural defects. The goal of this work is to identify the most promising techniques that are also compatible with operation in a microreactor environment.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Detection of Synchrophasor False Data Injection Attack using Feature Interactive Network

The synchrophasor data recorded by Phasor Measurement Units (PMUs) plays an increasingly critical role in the regulation and situational awareness of power systems. However, the widely installed PMUs are vulnerable to multiple malicious attacks from cyber hackers during data transmission and storage. To address this problem, a Modified Ensemble Empirical Mode Decomposition (MEEMD) is proposed first to extract the intrinsic mode functions of each Synchrophasor Data Attacks (SDA). The frequency-based adaptive screening criterion embedded in MEEMD is used to eliminate the false intrinsic mode functions. Next, a Multivariate Convolutional Neural Network (MCNN) is proposed to identify multiple SDA by utilizing the extracted intrinsic mode functions and original SDA as input vectors. A fusion block as the main structure of MCNN is also leveraged to increase the diversity of features and compress the model parameters. Integrating MEEMD and MCNN, a framework with automatic feature extraction and multi-source information fusion capability, referred to as Feature Interactive Network (FIN), is proposed to detect multiple SDA. Based on the proposed FIN framework, six types of SDA are explored for the first time using actual synchrophasor data in FNET/Grideye that was collected from different locations in the U.S. Eastern Interconnection. Finally, a large quantity of experiments with different attack strengths are used to evaluate the adaptability and classification performance of the proposed FIN.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Retrieving microphysical properties of concurrent pristine ice and snow using polarimetric radar observations

Abstract. Ice and mixed-phase clouds play a key role in our climate system because of their strong controls on global precipitation and radiation budget. Their microphysical properties have been characterized commonly by polarimetric radar measurements. However, there remains a lack of robust estimates of microphysical properties of concurrent pristine ice and aggregates because larger snow aggregates often dominate the radar signal and mask contributions of smaller pristine ice crystals. This paper presents a new method that separates the scattering signals of pristine ice embedded in snow aggregates in scanning polarimetric radar observations and retrieves their respective abundances and sizes for the first time. This method, dubbed ENCORE-ice, is built on an iterative stochastic ensemble retrieval framework. It provides the number concentration, ice water content, and effective mean diameter of pristine ice and snow aggregates with uncertainty estimates. Evaluations against synthetic observations show that the overall retrieval biases in the combined total microphysical properties are within 5 % and that the errors with respect to the truth are well within the retrieval uncertainty. The partitioning between pristine ice and snow aggregates also agrees well with the truth. Additional evaluations against in situ cloud probe measurements from a recent campaign for a stratiform cloud system are promising. Our median retrievals have a bias of 98 % in the total ice number concentration and 44 % in the total ice water content. This performance is generally better than the retrieval from empirical relationships. The ability to separate signals of different ice species and to provide their quantitative microphysical properties will open up many research opportunities, such as secondary ice production studies and model evaluations for ice microphysical processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A Moving Embedded Boundary Approach for the Compressible Navier-Stokes Equations in a Block-Structured Adaptive Refinement Framework

A computational technique has been developed to perform compressible flow simulations involving moving boundaries using an embedded boundary approach within the block-structured adaptive mesh refinement (SAMR) framework of AMReX [1], [91], [92]. We leverage the SAMR capability to obtain quantitatively accurate results whilst using robust, second-order finite volume schemes. A conservative, unsplit, cut-cell approach is utilized and a ghost-cell approach is developed for computing the flux on the moving, embedded boundary faces. A third-order least-squares formulation has been developed to compute the wall velocity gradients, and was found to significantly improve the performance of the solver in terms of the quantitative comparison of surface quantities such as the skin friction coefficient. Various test cases are performed to validate the method, and compared with analytical, experimental, and other numerical results in literature. Inviscid and viscous test cases are performed that span a wide regime of flow speeds - acoustic (harmonically pulsating sphere), smooth flows (expansion fan created by a receding piston) and flows with shocks (shock-cylinder interaction, shock-wedge interaction, pitching NACA 0012 airfoil and shock-cone interaction). A closed system with moving boundaries - an oscillating piston in a cylinder, showed that the percentage error in mass within the system decreases with refinement, demonstrating that the numerical scheme is conservative with grid refinement, but is not discretely conservative. Viscous test cases involve that of a horizontally moving cylinder at Re = 40, an inline oscillating cylinder at Re = 100, and a transversely oscillating cylinder at Re = 185. The judicious use of adaptive mesh refinement with appropriate refinement criteria to capture the regions of interest leads to well-resolved flow features, and good quantitative comparison is observed with the results available in literature.

adaptive refinement↗
DOE OSTI2022

On the broadband vibration isolation performance of nonlocal total-internal-reflection metasurfaces

The concept of a nonlocal elastic metasurface has been recently proposed and experimentally demonstrated in Zhu et al. (2020). When implemented in the form of a total-internal-reflection (TIR) interface, the metasurface can act as an elastic wave barrier that is impenetrable to deep subwavelength waves over an exceptionally wide frequency band. The underlying physical mechanism capable of delivering this broadband subwavelength performance relies on an intentionally nonlocal design that leverages long-range connections between the units forming the fundamental supercell. This paper explores the design and application of a nonlocal TIR metasurface to achieve broadband passive vibration isolation in a structural assembly made of multiple dissimilar elastic waveguides. The specific structural system comprises shell, plate, and beam waveguides, and can be seen as a prototypical structure emulating mechanical assemblies of practical interest for many engineering applications. The study also reports the results of an experimental investigation that confirms the significant vibration isolation capabilities afforded by the embedded nonlocal TIR metasurface. We report these results are particularly remarkable because they show that the performance of the nonlocal metasurface is preserved when applied to a complex structural assembly and under non-ideal incidence conditions of the incoming wave, hence significantly extending the validity of the results presented in Zhu et al. (2020). Results also confirm that, under proper conditions, the original concept of a planar metasurface can be morphed into a curved interface while still preserving full wave control capabilities.

42 ENGINEERING↗
DOE OSTI2020

First-principles studies of strongly correlated states in defect spin qubits in diamond

Using a recently developed quantum embedding theory, we present first-principles calculations of strongly correlated states of spin defects in diamond. Using this theory, effective Hamiltonians are constructed, which can be solved by classical and quantum computers; the latter promise a much more favorable scaling as a function of system size than the former. In particular, we report a study on the neutral group-IV vacancy complexes in diamond, and we discuss their strongly correlated spin-singlet and spin-triplet excited states. Finally, our results provide valuable predictions for experiments aimed at optical manipulation of these defects for quantum information technology applications.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Defect Detection Model Development for Large Scale Thermoplastic Printing

Large-format additive manufacturing (LFAM) offers several advantages, including high throughput, cost-effective pellet-fed extrusion, and the capability to produce large-scale structures. The main pain points of LFAM include start and stops during the printing process, warpage, long layer times that lead to bead freezing, and bead separation due to shrinkage. These issues can lead to overfill, underfill and buildup of material in different sections of a print. This can lead to hidden defects embedded within the printed layers, or even ultimate failure of the printed structure. This ensures these defects can only be identified through nondestructive testing (NDT) inspection methods after printing, which can be timely and costly. Aligned Vision work specializes in 2D projectors with visual inspection systems and machine learning. Traditionally system is used for composite layup and layup inspections. In this work we used the LFAM system at Oak Ridge National Laboratory to create defect rich samples. The Aligned Vision inspection system then performed in-situ monitoring of the print process after each part was printed. This in-situ vision inspection system was used to develop a layer-by-layer inspection model that looks for overfill, underfill, and the buildup of defects using only a camera-based vision system. This leads to the assurance of high-quality production components.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Low power and privacy preserving sensor platform for occupancy detection

A low-cost, low-power, stand-alone sensor platform having a visible-range camera sensor, a thermopile array, a microphone, a motion sensor, and a microprocessor that is configured to perform occupancy detection and counting while preserving the privacy of occupants. The platform is programmed to extract shape/texture from images in spatial domain; motion from video in time domain; and audio features in frequency domain. Embedded binarized neural networks are used for efficient object of interest detection. The platform is also programmed with advanced fusion algorithms for multiple sensor modalities addressing dependent sensor observations. The platform may be deployed for (i) residential use in detecting occupants for autonomously controlling building systems, such as HVAC and lighting systems, to provide energy savings, (ii) security and surveillance, such as to detect loitering and surveil places of interest, (iii) analyzing customer behavior and flows, (iv) identifying high performing stores by retailers.

Velipasalar, Senem↗
DOE OSTI2026

Diabatization with Electrostatic Embedding for Studying Photophysics in Organic Molecular Crystals

Highly emissive organic molecular crystals find applications in several areas, such as organic electronics, solar cells, and sensors. Understanding the excited-state mechanisms underlying these applications is essential for optimizing and controlling them effectively. Exciton models coupled with nonadiabatic dynamics, particularly quantum dynamics, provide crucial insights into photochemical and photophysical processes in molecular crystals. Nevertheless, there remains a lack of general tools and automated workflows to facilitate such simulations. In this paper, we present a computational strategy to investigate the photoactivated dynamics of organic molecular crystals, bridging methodologies traditionally used for molecular systems and materials science, with a particular focus on the interplay between local excitations and charge transfer (CT) processes. We have implemented an interface between the fromage and Overdia programs, enabling the construction of vibronic Hamiltonians for molecular crystals within an excited-state ONIOM(QM:QM′) framework, incorporating long-range electrostatics through a RESP-based Ewald summation. Fragment-based diabatization provides a route to quantum dynamics simulations in weak-to-intermediate coupling regimes. The method was applied to the photophysics of dibenzo[g,p]chrysene (DBC) crystals using time-dependent DFT. The fromage/ Overdia interface was employed to compute the couplings of local excitations and CT states for 18 unique DBC dimers in the crystal and to quantify the influence of electrostatic embedding, which was found to be modest (10−20%). Simulations on π-stacked dimers reproduced the small red shift observed experimentally from solution to crystal, attributed to electronic interactions among fixed monomers rather than crystal electrostatics. Quantum dynamics simulations revealed ultrafast population transfer from bright local excitations to CT states. This approach establishes a robust framework linking molecular and solid-state excited-state dynamics, with potential applications for studying excitations, defects, and impurities in molecular crystals.

Crystals↗
DOE OSTI2026

Semi-Analytical Hierarchical Bayesian Inference of Nonlinear Model Structure in Stochastic Dynamics: Applied to Compartmental Models of Infectious Diseases

A Bayesian computational framework for parsimonious inference in stochastic nonlinear dynamical systems is presented. This framework enables the concurrent estimation of system states, time-varying parameters, time-invariant parameters, and the optimal sparsity structure of the model parameters. Because differential equation-based models are often simplified mechanistic or phenomenological representations, robust inference from noisy measurement data requires explicit treatment of model error and uncertainty. Model error and time-varying parameters can be represented as random processes, enabling inference while making minimal assumptions about the underlying sources of discrepancy and variability. Adopting stochastic differential equation representations affords the model significant flexibility, but can also render it susceptible to overfitting during statistical inversion, where the inferred model may track noise rather than the underlying signal. To alleviate the effects of overfitting and to enable the discovery of the optimal sparse representation of the time-invariant parameters, a Bayesian sparse learning algorithm is embedded within the framework. This sparse learning framework adopts an approximate hierarchical Bayesian setting defined by a series of semi-analytical expressions. The model structure inference framework is validated using a stochastic compartmental model for tracking and forecasting active cases of an infectious disease. Compartmental models describe population-level infectious disease dynamics through interactions among population fractions grouped by disease state. Mathematically, such models consist of a system of coupled ordinary differential equations. This example adopts an expressive compartmental model that includes multiple possible interactions between disease states, motivated by early uncertainty surrounding COVID-19 reinfection dynamics and their implications for long-term epidemic forecasting. The sparse learning exercise permits the inference of a priori unknown epidemiological dynamics from simulated public health data, discovering the nested compartmental model that optimizes the trade-off between average data-fit and model complexity. It is shown that inducing sparsity among the model parameters eliminates redundant interactions between compartments, equivalently revealing the optimal coupling structure between differential equations.

97 MATHEMATICS AND COMPUTING↗