Search NASA⌕ Search

SEARCH · Search NASA

Results for “Continuity equation”

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 361 records · Page 20

Asynchronous distributed-memory task-parallel algorithm for compressible flows on unstructured 3D Eulerian grids

Here, we discuss the implementation of a finite element method, used to numerically solve the Euler equations of compressible flows, using an asynchronous runtime system (RTS). The algorithm is implemented for distributed-memory machines, using stationary unstructured 3D meshes, combining data-, and task-parallelism on top of the Charm++ RTS. Charm++’s execution model is asynchronous by default, allowing arbitrary overlap of computation and communication. Task-parallelism allows scheduling parts of an algorithm independently of, or dependent on, each other. Built-in automatic load balancing enables continuous redistribution of computational load by migration of work units based on real-time CPU load measurement. The RTS also features automatic checkpointing, fault tolerance, resilience against hardware failure, and supports power-, and energy-aware computation. We demonstrate scalability up to 25 x 10 9 cells at $\mathscr{O}$10 4 compute cores and the benefits of automatic load balancing for irregular workloads. The full source code with documentation is available at https://quinoacomputing.org.

42 ENGINEERING↗

A Logarithmic Bayesian Approach to Quantum Error Detection

We consider the problem of continuous quantum error correction from a Bayesian perspective, proposing a pair of digital filters using logarithmic probabilities that are able to achieve near-optimal performance on a three-qubit bit-flip code, while still being reasonable to implement on low-latency hardware. These practical filters are approximations of an optimal filter that we derive explicitly for finite time steps, in contrast with previous work that has relied on stochastic differential equations such as the Wonham filter. By utilizing logarithmic probabilities, we are able to eliminate the need for explicit normalization and can reduce the Gaussian noise distribution to a simple quadratic expression. The state transitions induced by the bit-flip errors are modeled using a Markov chain, which for log-probabilities must be evaluated using a LogSumExp function. We develop the two versions of our filter by constraining this LogSumExp to have either one or two inputs, which favors either simplicity or accuracy, respectively. Using simulated data, we demonstrate that the single-term and two-term filters are able to significantly outperform both a double threshold scheme and a linearized version of the Wonham filter in tests of error detection under a wide variety of error rates and time steps.

97 MATHEMATICS AND COMPUTING↗

Leveraging machine learning to enhance aerosol classification using Single-Particle Mass Spectrometry

Advancing automated classification of atmospheric aerosols from Single-Particle Mass Spectrometry (SPMS) data remains challenging due to overlapping ion signatures, compositional diversity, and limited labeled data. This study evaluates supervised and semi-supervised learning frameworks to enhance aerosol identification by jointly leveraging labeled and unlabeled spectra. Four models were compared: a supervised Support Vector Machine (SVM), a self-training SVM, a stacked autoencoder classifier, and a stacked autoencoder trained using a temporal-ensembling Mean Teacher approach. All models achieved high and stable accuracies (90.0 %–91.1 %), surpassing previous results on the same dataset (87 %) and matching the performance of state-of-the-art deep learning methods. Despite small global metric differences (≤ 1 %), semi-supervised variants yielded up to 5 %–10 % improvements for compositionally rare particle types – such as soot (0.77 % of spectra, F1-score: 0.93–0.97) and hazelnut pollen (0.98 % of spectra, F1-score: 0.97–1.00) – equating to roughly ∼ 187 additional correctly classified spectra. These gains are scientifically significant, as such rare particles exert disproportionate influence on radiative absorption and ice nucleation processes; their improved detection reduces modeled uncertainties in aerosol absorption optical depth and mixed-phase cloud ice nucleation rates. The models' residual misclassifications (≈ 9 %) largely arise from true spectral overlap among chemically adjacent species (e.g., Na- vs. K-feldspar, coated vs. uncoated feldspars), reflecting physical compositional continuity rather than algorithmic error. Collectively, these findings demonstrate that leveraging unlabeled data to learn robust spectral representations and refine classification enhances both fidelity and interpretability, bridging data-driven analysis with aerosol–climate process understanding.

54 ENVIRONMENTAL SCIENCES↗

A note on higher-order and nonlinear limiting approaches for continuously bounds-preserving discontinuous Galerkin methods

In Dzanic (2024), a limiting approach for high-order discontinuous Galerkin schemes was introduced which allowed for imposing constraints on the solution continuously (i.e., everywhere within the element). While exact for linear constraint functionals, this approach only imposed a sufficient (but not the minimum necessary) amount of limiting for nonlinear constraint functionals. This short note shows how this limiting approach can be extended to allow exactness for general nonlinear quasiconcave constraint functionals through a nonlinear limiting procedure, reducing unnecessary numerical dissipation. Finally, some examples are shown for nonlinear pressure and entropy constraints in the compressible gas dynamics equations, where both analytic and iterative approaches are used.

97 MATHEMATICS AND COMPUTING↗

Methods for Evaluating DC Arc Incident Energy in PV Systems: Preprint

Renewable energy systems continue to be one of the fastest growing segments of the energy industry. This paper focuses on the understanding of how photovoltaic (PV) technology behaves under dc arc conditions. Emphasis is placed on the electrical safety aspect of DC arc flash incident energy evaluation. Because of the fast proliferation of PV systems and the lack of formal equivalent calculation guidelines such as IEEE 1584 for AC systems, it has been necessary to rely on different equations and models presented by various researchers over the last few years. This paper discusses the behavior of PV systems under arc conditions and presents the results of available methods to estimate the dc arc flash incident energy. This paper provides a comparative analysis of a proposed arc-flash incident energy calculation method against different laboratory tests including those performed by NREL. Detailed explanations are provided regarding the effect of PV module I-V and P-V curves under arcing conditions. Examples of the application of the proposed calculation method to the test measurements are included.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Laser-driven ramp-compression experiments on the national ignition facility

This report details the analyses and related uncertainties in measuring longitudinal-stress–density paths in indirect laser-driven ramp equation-of-state (EOS) experiments. Experiments were conducted at the National Ignition Facility (NIF) located at the Lawrence Livermore National Laboratory. The NIF can deliver up to 2 MJ of laser energy over 30 ns and provide the necessary laser power and control to ramp compress materials to TPa pressures (1 TPa = 10 × 10 6 atmospheres). These data provide low-temperature solid-state EOS data relevant to the extreme conditions found in the deep interiors of giant planets. In these experiments, multi-stepped samples with thicknesses in the range of 40–120 µm experience an initial shock compression followed by a time-dependent ramp compression to peak pressure. Interface velocity measurements from each thickness combine to place a constraint on the Lagrangian sound speed as a function of particle velocity, which in turn allows for the determination of a continuous stress–density path to high levels of compressibility. In this report, we present a detailed description of the experimental techniques and measurement uncertainties and describe how these uncertainties combine to place a final uncertainty in both stress and density. Here, we address the effects of time-dependent deformation and the sensitivity of ramp EOS techniques to the onset of phase transformations.

47 OTHER INSTRUMENTATION↗

DC Arc Incident Energy in Photovoltaic Systems: Methods for Evaluation

Renewable energy systems continue to be one of the fastest growing segments of the energy industry. This article focuses on the understanding of how photovoltaic (PV) technology behaves under dc arc conditions. Emphasis is placed on the electrical safety aspect of dc arc-flash incident energy (IE) evaluation. Because of the fast proliferation of PV systems and lack of formal equivalent calculation guidelines, such as IEEE 1584 for ac systems, it has been necessary to rely on different equations and models presented by various researchers over the last few years. This article discusses the behavior of PV systems under arc conditions and presents the results of the available methods to estimate the dc arc-flash IE. It provides a comparative analysis of a proposed arc-flash IE calculation method against different laboratory tests, including those performed for this article at the National Renewable Energy Laboratory (NREL). Detailed explanations are provided regarding the effect of the PV module current-voltage (I-V) and power-voltage (P-V) curves under arcing conditions. Examples of the application of the proposed calculation method to the test measurements are included.

arc discharges↗

CFD modeling of natural circulation in LiCl-KCl molten salt closed loop

Characterizing flow within a molten salt closed-loop system is crucial for assessing system requirements, evaluating performance, and identifying potential flaws. Direct flow measurement using instrumentation is challenging due to extreme environmental conditions and the limitations associated with measuring molten salt flow under natural convection. Here, this study aims to provide comprehensive insights into the thermal-hydraulic behavior of a closed loop, with a particular focus on temperature distribution and velocity prediction. The Computational Fluid Dynamics (CFD) model demonstrated the capability to effectively simulate and predict both temperature distributions and flow velocities within the molten salt loop. The CFD model's predictive capability was validated by its ability to replicate temperature measurements under varying boundary conditions. The analysis revealed that the CFD model tends to underpredict temperatures in the cold leg and overpredict them in the hot leg, highlighting the need for continuous model refinement and acknowledging the limitations of using a steady-state approach. Furthermore, the potential of using external temperature measurements to estimate internal molten salt temperatures and predict flow velocity was explored, revealing that this approach could introduce up to a 5.5% error in flow velocity calculations. Line probes mapping temperature distributions across the tube's cross-section and molten salt provided valuable insights into temperature gradients, emphasizing the need for a thermal conductivity equation for molten salt with lower uncertainty to achieve more accurate temperature predictions of the system.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

A First Principles Approach to Spectral Phonon Transport in Heterostructures

Understanding thermal transport across interfaces which give rise to a thermal resistance (also known as Kapitza resistance) is a critical issue affecting the development of nanotechnologies. Much modern and emergent nanotechnology consist of adjacent materials, and phonon mediated heat transfer governs thermal behavior across internal interfaces in these devices. The physics of thermal transport in solids are governed both by phenomena occurring at the atomic scale and interactions with the material's microstructure. The forecasting of fundamental quantities such as temperature, heat flux and thermal conductivity typically employs the semi-classical Boltzmann transport equation to predict the macroscopic behavior of materials in terms of the microscopic dynamics of its heat carriers. Kapitza resistance was first discovered in liquid helium experiments and has led to a fundamental research thrust in micro and nano-scale heat transport, the behavior of thermal carriers across internal interfaces. Thermal interfacial resistance (TIR) is a widely studied phenomenon, first engaged by Swartz and Pohl through their development of the acoustic and diffuse mismatch methods, then continued through myriad efforts with varying methods and approaches in an attempt to resolve carrier behavior at thermal interfaces. Many of the fundamental approaches to TIR have been at the nanoscale, and research is conducted with molecular dynamics (MD) and density functional theory (DFT) methods. The limitations of these methods is system size; atomistic methods tend to be limited to system sizes of 100,000 atoms or less. Larger length-scale methods have also been pursued, based on the principles of acoustic or diffuse mismatch, but not all include simulation of TIR using a full phonon band spectrum, or temperature dependent methods. Our approach to enabling phonon transport in layered materials draws upon our previous work of demonstrating spectrally coupled phonon transport in homogeneous and heterogeneous materials. We use a semi-analytical approach in which the Bose-Einstein (B-E) statistics set the strength of the phonon radiance in a frequency group, but the B-E statistics are informed with information from the transport system. The B-E statistics in a single frequency group feels the influence of all the groups through the spatial temperature. We also include a new field term which is an indicator of the amount of non-equilibrium behavior of the phonon spectrum---this is added to the phonon source term in all groups to ensure closure and conservation of energy, as the phonon groups in the transport system and the analytical systems are coupled. This work builds upon our previous approach by adding a phonon coupling term at an internal interface, using the principles of the DMM through transmission and reflection coefficients. In this work, the coefficients are determined through computing a common temperature at the interface, influenced by the phonon band structure of both materials, in effect, providing mixing between the two material systems and using the common temperature to set the strength of the phonon radiance at the boundaries on either side of the interface. Our approach uses material properties computed along various crystallographic orientations, and while some isotropy is built into the interface condition, the material properties weight the phonon distributions in the proper crystalline direction. Greater resolution of phonon behavior in proximity to an interface, and more accurate predictions of TIR are obtained. While it is true the assumption of diffuse mismatch can yield inconsistent results compared to experiment especially at low temperatures, this work focuses on room temperature and beyond effects, for future applications in nuclear fuel, or thermoelectric devices; a modified mismatch approach may be feasible if applied properly. Additionally, our methods focus on bridging mesoscale to engineering scale

36 MATERIALS SCIENCE↗

A general framework for substructuring‐based domain decomposition methods for models having nonlocal interactions

Abstract A mathematical framework is provided for a substructuring‐based domain decomposition (DD) approach for nonlocal problems that features interactions between points separated by a finite distance. Here, by substructuring it is meant that a traditional geometric configuration for local partial differential equation (PDE) problems is used in which a computational domain is subdivided into non‐overlapping subdomains. In the nonlocal setting, this approach is substructuring‐based in the sense that those subdomains interact with neighboring domains over interface regions having finite volume, in contrast to the local PDE setting in which interfaces are lower dimensional manifolds separating abutting subdomains. Key results include the equivalence between the global, single‐domain nonlocal problem and its multi‐domain reformulation, both at the continuous and discrete levels. These results provide the rigorous foundation necessary for the development of efficient solution strategies for nonlocal DD methods.

Capodaglio, Giacomo↗

Planck 2018 results. VI. Cosmological parameters

In this paper, we present cosmological parameter results from the final full-mission Planck measurements of the cosmic microwave background (CMB) anisotropies, combining information from the temperature and polarization maps and the lensing reconstruction. Compared to the 2015 results, improved measurements of large-scale polarization allow the reionization optical depth to be measured with higher precision, leading to significant gains in the precision of other correlated parameters. Improved modelling of the small-scale polarization leads to more robust constraints on many parameters, with residual modelling uncertainties estimated to affect them only at the 0.5 σ level. We find good consistency with the standard spatially-flat 6-parameter ΛCDM cosmology having a power-law spectrum of adiabatic scalar perturbations (denoted “base ΛCDM” in this paper), from polarization, temperature, and lensing, separately and in combination. A combined analysis gives dark matter density Ω c h 2 = 0.120 ± 0.001, baryon density Ω b h 2 = 0.0224 ± 0.0001, scalar spectral index n s = 0.965 ± 0.004, and optical depth τ = 0.054 ± 0.007 (in this abstract we quote 68% confidence regions on measured parameters and 95% on upper limits). The angular acoustic scale is measured to 0.03% precision, with 100 θ * = 1.0411 ± 0.0003. These results are only weakly dependent on the cosmological model and remain stable, with somewhat increased errors, in many commonly considered extensions. Assuming the base-ΛCDM cosmology, the inferred (model-dependent) late-Universe parameters are: Hubble constant H 0 = (67.4 ± 0.5) km s –1 Mpc –1 ; matter density parameter Ω m = 0.315 ± 0.007; and matter fluctuation amplitude σ 8 = 0.811 ± 0.006. We find no compelling evidence for extensions to the base-ΛCDM model. Combining with baryon acoustic oscillation (BAO) measurements (and considering single-parameter extensions) we constrain the effective extra relativistic degrees of freedom to be N eff = 2.99 ± 0.17, in agreement with the Standard Model prediction N eff = 3.046, and find that the neutrino mass is tightly constrained to Σ m ν < 0.12 eV. The CMB spectra continue to prefer higher lensing amplitudes than predicted in base ΛCDM at over 2 σ , which pulls some parameters that affect the lensing amplitude away from the ΛCDM model; however, this is not supported by the lensing reconstruction or (in models that also change the background geometry) BAO data. The joint constraint with BAO measurements on spatial curvature is consistent with a flat universe, Ω K = 0.001 ± 0.002. Also combining with Type Ia supernovae (SNe), the dark-energy equation of state parameter is measured to be w 0 = –1.03 ± 0.03, consistent with a cosmological constant. We find no evidence for deviations from a purely power-law primordial spectrum, and combining with data from BAO, BICEP2, and Keck Array data, we place a limit on the tensor-to-scalar ratio r 0.002 < 0.06. Standard big-bang nucleosynthesis predictions for the helium and deuterium abundances for the base-ΛCDM cosmology are in excellent agreement with observations. The Planck base-ΛCDM results are in good agreement with BAO, SNe, and some galaxy lensing observations, but in slight tension with the Dark Energy Survey’s combined-probe results including galaxy clustering (which prefers lower fluctuation amplitudes or matter density parameters), and in significant, 3.6 σ , tension with local measurements of the Hubble constant (which prefer a higher value). Simple model extensions that can partially resolve these tensions are not favoured by the Planck data.

79 ASTRONOMY AND ASTROPHYSICS↗

Superconductivity out of a non-Fermi liquid: Free energy analysis

In this paper, we present an in-depth analysis of the condensation energy E c for a superconductor in a situation when superconductivity emerges out of a non-Fermi liquid due to pairing mediated by a massless boson. This is the case for electronic-mediated pairing near a quantum-critical point in metal, for pairing in SYK-type models, and for phonon-mediated pairing in the properly defined limit, when the dressed Debye frequency vanishes. We consider a subset of these quantum-critical models, in which the pairing in a channel with a proper spatial symmetry is described by an effective 0 + 1 dimensional model with the effective dynamical interaction V(Ω m ) = g¯ γ /|Ω m | γ , where γ is model-specific (the γ model). In previous papers, we argued that the pairing in the γ model is qualitatively different from that in a Fermi liquid, and the gap equation at T = 0 has an infinite number of topologically distinct solutions, Δ n (ω m ), where an integer n, running between 0 and infinity, is the number of zeros of Δ n (ω m ) on the positive Matsubara axis. This gives rise to the set of extrema of E c at E c,n , of which E c,0 is the global minimum. Here, the spectrum E c,n is discrete for a generic γ < 2 but becomes continuous at γ = 2–0. Here, we discuss in more detail the profile of the condensation energy near each E c,n and the transformation from a discrete to a continuous spectrum at γ → 2. We also discuss the free energy and the specific heat of the γ model in the normal state.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A novel physics-based and data-supported microstructure model for part-scale simulation of laser powder bed fusion of Ti-6Al-4V

The elasto-plastic material behavior, material strength and failure modes of metals fabricated by additive manufacturing technologies are significantly determined by the underlying process-specific microstructure evolution. In this work a novel physics-based and data-supported phenomenological microstructure model for Ti-6Al-4V is proposed that is suitable for the part-scale simulation of laser powder bed fusion processes. The model predicts spatially homogenized phase fractions of the most relevant microstructural species, namely the stable β-phase, the stable α s -phase as well as the metastable Martensite α m -phase, in a physically consistent manner. In particular, the modeled microstructure evolution, in form of diffusion-based and non-diffusional transformations, is a pure consequence of energy and mobility competitions among the different species, without the need for heuristic transformation criteria as often applied in existing models. The mathematically consistent formulation of the evolution equations in rate form renders the model suitable for the practically relevant scenario of temperature- or time-dependent diffusion coefficients, arbitrary temperature profiles, and multiple coexisting phases. Due to its physically motivated foundation, the proposed model requires only a minimal number of free parameters, which are determined in an inverse identification process considering a broad experimental data basis in form of time-temperature transformation diagrams. Subsequently, the predictive ability of the model is demonstrated by means of continuous cooling transformation diagrams, showing that experimentally observed characteristics such as critical cooling rates emerge naturally from the proposed microstructure model, instead of being enforced as heuristic transformation criteria. Eventually, the proposed model is exploited to predict the microstructure evolution for a realistic selective laser melting application scenario and for the cooling/quenching process of a Ti-6Al-4V cube of practically relevant size. Numerical results confirm experimental observations that Martensite is the dominating microstructure species in regimes of high cooling rates, e.g., due to highly localized heat sources or in near-surface domains, while a proper manipulation of the temperature field, e.g., by preheating the base-plate in selective laser melting, can suppress the formation of this metastable phase.

Inverse parameter identification↗

Is the Adoption of Electric Vehicles (EVs) and Solar Photovoltaics (PVs) Interdependent or Independent? An Integrated EVs-PVs Modeling Framework

Transportation and residential energy consumption account for more than one-half of the overall energy consumption in the United States. Adoption of electric vehicles (EVs) can play a key role in decarbonizing the transportation sector, while the adoption of renewable energy sources (e.g., solar photovoltaics or PVs) could bring similar benefits to the residential energy sector and in turn support transport electrification. Although the market share for both EVs and PVs continue to grow, both these emerging technologies are deployed rather disjoint without considering the existence of potential similarities among users who own (or aspire to own) EVs and PVs. This might be due to lack of understanding of the behavioral interdependence in consumer preferences towards these technologies. To fill this gap in knowledge, this study utilizes data from the 2018 WholeTraveler Transportation Behavior Study to develop an integrated model system that explores EV and PV adoption behaviors. A structural equations model (SEM) is employed that incorporates direct effects as well as error correlations among the adoption behaviors for EVs and PVs. Model results indicate that the adoption behavior for both these technologies is indeed interconnected and significantly influenced by attitudes, values and personality traits. Findings from this research suggest that incentives (e.g., subsidies) that drive 'bundled' adoption of EVs and PVs could accelerate the adoption of both these sustainable technologies. The results highlights the need to consider transport and building energy efficient technology adoption behavior in a single integrated structure.

ADVANCED PROPULSION SYSTEMS,ENERGY PLANNING, POLIC↗

The electric vehicles-solar photovoltaics Nexus: Driving cross-sectoral adoption of sustainable technologies

Residential and transportation energy consumption account for more than one-half of the overall energy consumption in the United States. Adoption of electric vehicles (EVs) can play a key role in decarbonizing the transportation sector, while the adoption of renewable energy sources (e.g., solar photovoltaics [PVs]) could bring similar benefits to the residential energy sector and in turn support transport electrification. Although the market shares for both EVs and PVs continue to grow, both of these emerging technologies are deployed rather disjointly, without considering the existence of potential similarities among users who own (or aspire to own) these technologies. This might be due to lack of understanding of the behavioral interdependence in consumer preferences toward these technologies. To fill this gap in knowledge, this study utilizes data from the 2018 WholeTraveler Transportation Behavior Study to develop an integrated model system that explores interactive EV and PV adoption behaviors. A structural equation model is employed that incorporates direct effects as well as error correlations among the adoption behaviors for EVs and PVs. Model results indicate that the adoption behavior for both these technologies is indeed interconnected and significantly influenced by attitudes, values, and personality traits. Findings from this research suggest that incentives (e.g., subsidies) that drive bundled adoption of EV-PV systems could accelerate the adoption of both of these sustainable technologies. In conclusion, this study highlights the need to consider transport and building energy-efficient technology adoption behavior in a single integrated structure.

14 SOLAR ENERGY↗

Scattering using real-time path integrals

Background: Path integrals are a powerful tool for solving problems in quantum theory that are not amenable to a treatment by perturbation theory. Most path integral computations require an analytic continuation to imaginary time. While imaginary time treatments of scattering are possible, imaginary time is not a natural framework for treating scattering problems. More importantly, quantum algorithms for calculating path integrals require real-time evolution. Purpose: Here, we test a recently introduced method for performing direct calculations of scattering observables using real-time path integrals in order to understand the challenges facing real-time path integral calculations of scattering observables. Method: The computations are based on a new interpretation of the path integral as the expectation value of a potential functional on cylinder sets of continuous paths with respect to a complex probability distribution. The method can in principle be applied to arbitrary short-range potentials. Results: The method is applied to compute matrix elements of Møller wave operators applied to narrow wave packets. These are used to calculate half-shell sharp-momentum transition matrix elements for one-dimensional potential scattering. The calculations for half-shell transition operator matrix elements converge to the numerical solution of the Lippmann-Schwinger equation. Conclusions: This work presents a proof in principle that scattering observables can be computed using real-time Feynman path integrals. While the computational method is not efficient, it can be improved. It provides a laboratory for studying quantum computational algorithms that are applicable to scattering problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Sliding Mode Network Perimeter Control

Urban traffic congestion is a chronic problem faced by many cities in the US and worldwide. It results in inefficient infrastructure use as well as increased vehicle fuel consumption and emission levels. Congestion is intertwined with delay, as road users waste precious hours on the road, which in turn reduces productivity. Researchers have developed, and continue to design, tools and systems to alleviate this problem. Network perimeter control is one such tool that has been studied extensively. It attempts to control the flow of vehicles entering a protected area to ensure that the congested regime predetermined by the Network Fundamental Diagram (NFD) is not reached. In this paper, a method derived from sliding mode control theory is presented. Its main advantages over proportional-integral controllers include (1) minimal tuning, (2) no linearization of the governing equations, (3) no assumptions with regard to the shape of the NFD, and (4) ability to handle various demand profiles without the need to retune the controller. A sliding mode controller was implemented and tested on a congested grid network. The results show that the proposed controller produces network-wide delay savings and disperses congestion effectively.

42 ENGINEERING↗

A general framework for substructuring-based domain decomposition methods for models having nonlocal interactions [including apppendix B]

A rigorous mathematical framework is provided for a substructuring-based domain-decomposition approach for nonlocal problems that feature interactions between points separated by a finite distance. Here, by substructuring it is meant that a traditional geometric configuration for local partial differential equation problems is used in which a computational domain is subdivided into non-overlapping subdomains. In the nonlocal setting, this approach is substructuring-based in the sense that those subdomains interact with neighboring domains over interface regions having finite volume, in contrast to the local PDE setting in which interfaces are lower dimensional manifolds separating abutting subdomains Key results include the equivalence between the global, single-domain nonlocal problem and its multi-domain reformulation, both at the continuous and discrete levels. These results provide the rigorous foundation necessary for the development of efficient solution strategies for nonlocal domain-decomposition methods.

97 MATHEMATICS AND COMPUTING↗