skip to main content


Title: Statistical Metamodeling of Dynamic Network Loading
Dynamic traffic assignment models rely on a network performance module known as dynamic network loading (DNL), which expresses flow propagation, flow conservation, and travel delay at a network level. The DNL defines the so-called network delay operator , which maps a set of path departure rates to a set of path travel times (or costs). It is widely known that the delay operator is not available in closed form, and has undesirable properties that severely complicate DTA analysis and computation, such as discontinuity, nondifferentiability, nonmonotonicity, and computational inefficiency. This paper proposes a fresh take on this important and difficult issue, by providing a class of surrogate DNL models based on a statistical learning method known as Kriging . We present a metamodeling framework that systematically approximates DNL models and is flexible in the sense of allowing the modeler to make trade-offs among model granularity, complexity, and accuracy. It is shown that such surrogate DNL models yield highly accurate approximations (with errors below 8%) and superior computational efficiency (9 to 455 times faster than conventional DNL procedures such as those based on the link transmission model). Moreover, these approximate DNL models admit closed-form and analytical delay operators, which are Lipschitz continuous and infinitely differentiable, with closed-form Jacobians. We provide in-depth discussions on the implications of these properties to DTA research and model applications.  more » « less
Award ID(s):
1662968
NSF-PAR ID:
10122278
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Transportation research. Part B, Methodological
Volume:
117
ISSN:
1879-2367
Page Range / eLocation ID:
740-756
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Dynamic user equilibrium (DUE) is the most widely studied form of dynamic traffic assignment (DTA), in which road travelers engage in a non-cooperative Nash-like game with departure time and route choices. DUE models describe and predict the time-varying traffic flows on a network consistent with traffic flow theory and travel behavior. This paper documents theoretical and numerical advances in synthesizing traffic flow theory and DUE modeling, by presenting a holistic computational theory of DUE, which is numerically implemented in a MATLAB package. In particular, the dynamic network loading (DNL) sub-problem is formulated as a system of differential algebraic equations based on the Lighthill-Whitham-Richards fluid dynamic model, which captures the formation, propagation and dissipation of physical queues as well as vehicle spillback on networks. Then, the fixed-point algorithm is employed to solve the DUE problems with simultaneous route and departure time choices on several large-scale networks. We make openly available the MATLAB package, which can be used to solve DUE problems on user-defined networks, aiming to not only facilitate benchmarking a wide range of DUE algorithms and solutions, but also offer researchers a platform to further develop their own models and applications. The MATLAB package and computational examples are available online. 
    more » « less
  2. Additive manufacturing has been recognized as an industrial technological revolution for manufacturing, which allows fabrication of materials with complex three-dimensional (3D) structures directly from computer-aided design models. Using two or more constituent materials with different physical and mechanical properties, it becomes possible to construct interpenetrating phase composites (IPCs) with 3D interconnected structures to provide superior mechanical properties as compared to the conventional reinforced composites with discrete particles or fibers. The mechanical properties of IPCs, especially response to dynamic loading, highly depend on their 3D structures. In general, for each specified structural design, it could take hours or days to perform either finite element analysis (FEA) or experiments to test the mechanical response of IPCs to a given dynamic load. To accelerate the physics-based prediction of mechanical properties of IPCs for various structural designs, we employ a deep neural operator (DNO) to learn the transient response of IPCs under dynamic loading as surrogate of physics-based FEA models. We consider a 3D IPC beam formed by two metals with a ratio of Young’s modulus of 2.7, wherein random blocks of constituent materials are used to demonstrate the generality and robustness of the DNO model. To obtain FEA results of IPC properties, 5000 random time-dependent strain loads generated by a Gaussian process kennel are applied to the 3D IPC beam, and the reaction forces and stress fields inside the IPC beam under various loading are collected. Subsequently, the DNO model is trained using an incremental learning method with sequence-to-sequence training implemented in JAX, leading to a 100X speedup compared to widely used vanilla deep operator network models. After an offline training, the DNO model can act as surrogate of physics-based FEA to predict the transient mechanical response in terms of reaction force and stress distribution of the IPCs to various strain loads in one second at an accuracy of 98%. Also, the learned operator is able to provide extended prediction of the IPC beam subject to longer random strain loads at a reasonably well accuracy. Such superfast and accurate prediction of mechanical properties of IPCs could significantly accelerate the IPC structural design and related composite designs for desired mechanical properties. 
    more » « less
  3. Abstract

    Groundwater discharge zones connect aquifers to surface water, generating baseflow and serving as ecosystem control points across aquatic ecosystems. The influence of groundwater discharge on surface flow connectivity, fate and transport of contaminants and nutrients, and thermal habitat depends strongly on hydrologic characteristics such as the spatial distribution, age, and depth of source groundwater flow paths. Groundwater models have the potential to predict spatial discharge characteristics within river networks, but models are often not evaluated against these critical characteristics and model equifinality with respect to discharge processes is a known challenge. We quantify discharge characteristics across a suite of groundwater models with commonly used frameworks and calibration data. We developed a base model (MODFLOW‐NWT) for a 1,570‐km2watershed in the northeastern United States and varied the calibration data, control of river‐aquifer exchange directionality, and resolution. Most models (n = 11 of 12) fit similarly to calibration metrics, but patterns in discharge location, flow path depth, and subsurface travel time varied substantially. We found (1) a 15% difference in the percent of discharge going to first‐order streams, (2) threefold variations in flow path depth, and (3) sevenfold variations in the subsurface travel times among the models. We recalibrated three models using a synthetic discharge location data set. Calibration with discharge location data reduced differences in simulated discharge characteristics, suggesting an approach to improved equifinality based on widespread field‐based mapping of discharge zones. Our work quantifying variation across common modeling approaches is an important step toward characterizing and improving predictions of groundwater discharge characteristics.

     
    more » « less
  4. This paper presents methods to obtain analytical solutions to a class of continuous traffic equilibrium problems, where continuously distributed customers from a bounded two-dimensional service region seek service from one of several discretely located facilities via the least congested travel path. We show that under certain conditions, the traffic flux at equilibrium, which is governed by a set of partial differential equations, can be decomposed with respect to each facility and solved analytically. This finding paves the foundation for an efficient solution scheme. Closed-form solution to the equilibrium problem can be obtained readily when the service region has a certain regular shape, or through an additional conformal mapping if the service region has an arbitrary simply connected shape. These results shed light on some interesting properties of traffic equilibrium in a continuous space. This paper also discusses how service facility locations can be easily optimized by incorporating analytical formulas for the total generalized cost of spatially distributed customers under congestion. Examples of application contexts include gates or booths for pedestrian traffic, as well as launching sites for air vehicles. Numerical examples are used to show the superiority of the proposed optimization framework, in terms of both solution quality and computation time, as compared with traditional approaches based on discrete mathematical programming and partial differential equation solution methods. An example with the metro station entrances at the Beijing Railway Station is also presented to illustrate the usefulness of the proposed traffic equilibrium and location design models. 
    more » « less
  5. : In order to evaluate urban earthquake resilience, reliable structural modeling is needed. However, detailed modeling of a large number of structures and carrying out time history analyses for sets of ground motions are not practical at an urban scale. Reduced-order surrogate models can expedite numerical simulations while maintaining necessary engineering accuracy. Neural networks have been shown to be a powerful tool for developing surrogate models, which often outperform classical surrogate models in terms of scalability of complex models. Training a reliable deep learning model, however, requires an immense amount of data that contain a rich input-output relationship, which typically cannot be satisfied in practical applications. In this paper, we propose model-informed symbolic neural networks (MiSNN) that can discover the underlying closed-form formulations (differential equations) for a reduced-order surrogate model. The MiSNN will be trained on datasets obtained from dynamic analyses of detailed reinforced concrete special moment frames designed for San Francisco, California, subject to a series of selected ground motions. Training the MiSNN is equivalent to finding the solution to a sparse optimization problem, which is solved by the Adam optimizer. The earthquake ground acceleration and story displacement, velocity, and acceleration time histories will be used to train 1) an integrated SNN, which takes displacement and velocity states and outputs the absolute acceleration response of the structure; and 2) a distributed SNN, which distills the underlying equation of motion for each story. The results show that the MiSNN can reduce computational cost while maintaining high prediction accuracy of building responses. 
    more » « less