Abstract In this paper we present a reconstruction technique for the reduction of unsteady flow data based on neural representations of time‐varying vector fields. Our approach is motivated by the large amount of data typically generated in numerical simulations, and in turn the types of data that domain scientists can generatein situthat are compact, yet useful, for post hoc analysis. One type of data commonly acquired during simulation are samples of the flow map, where a single sample is the result of integrating the underlying vector field for a specified time duration. In our work, we treat a collection of flow map samples for a single dataset as a meaningful, compact, and yet incomplete, representation of unsteady flow, and our central objective is to find a representation that enables us to best recover arbitrary flow map samples. To this end, we introduce a technique for learning implicit neural representations of time‐varying vector fields that are specifically optimized to reproduce flow map samples sparsely covering the spatiotemporal domain of the data. We show that, despite aggressive data reduction, our optimization problem — learning a function‐space neural network to reproduce flow map samples under a fixed integration scheme — leads to representations that demonstrate strong generalization, both in the field itself, and using the field to approximate the flow map. Through quantitative and qualitative analysis across different datasets we show that our approach is an improvement across a variety of data reduction methods, and across a variety of measures ranging from improved vector fields, flow maps, and features derived from the flow map.
more »
« less
Representing and Computing the B-Derivative of the Piecewise-Differentiable Flow of a Class of Nonsmooth Vector Fields
Abstract This paper concerns first-order approximation of the piecewise-differentiable flow generated by a class of nonsmooth vector fields. Specifically, we represent and compute the Bouligand (or B-)derivative of the piecewise-differentiable flow generated by a vector field with event-selected discontinuities. Our results are remarkably efficient: although there are factorially many “pieces” of the derivative, we provide an algorithm that evaluates its action on a tangent vector using polynomial time and space, and verify the algorithm's correctness by deriving a representation for the B-derivative that requires “only” exponential time and space to construct. We apply our methods in two classes of illustrative examples: piecewise-constant vector fields and mechanical systems subject to unilateral constraints.
more »
« less
- Award ID(s):
- 2045014
- PAR ID:
- 10345345
- Date Published:
- Journal Name:
- Journal of Computational and Nonlinear Dynamics
- Volume:
- 17
- Issue:
- 9
- ISSN:
- 1555-1415
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We introduce SLANG.D, an extension to the Slang shading language that incorporates first-class automatic differentiation support. The new shading language allows us to transform a Direct3D-based path tracer to be fully differentiable with minor modifications to existing code. SLANG.D enables a shared ecosystem between machine learning frameworks and pre-existing graphics hardware API-based rendering systems, promoting the interchange of components and ideas across these two domains. Our contributions include a differentiable type system designed to ensure type safety and semantic clarity in codebases that blend differentiable and non-differentiable code, language primitives that automatically generate both forward and reverse gradient propagation methods, and a compiler architecture that generates efficient derivative propagation shader code for graphics pipelines. Our compiler supports differentiating code that involves arbitrary control-flow, dynamic dispatch, generics and higher-order differentiation, while providing developers flexible control of checkpointing and gradient aggregation strategies for best performance. Our system allows us to differentiate an existing real-time path tracer, Falcor, with minimal change to its shader code. We show that the compiler-generated derivative kernels perform as efficiently as handwritten ones. In several benchmarks, the SLANG.D code achieves significant speedup when compared to prior automatic differentiation systems.more » « less
-
Infinitesimal contraction analysis, wherein global convergence results are obtained from properties of local dynamics, is a powerful analysis tool. In this paper, we generalize infinitesimal contraction analysis to hybrid systems in which state-dependent guards trigger transitions defined by reset maps between modes that may have different norms and need not be of the same dimension. In contrast to existing literature, we do not restrict mode sequence or dwell time. We work in settings where the hybrid system flow is differentiable almost everywhere and its derivative is the solution to a jump-linear-time-varying differential equation whose jumps are defined by a saltation matrix determined from the guard, reset map, and vector field. Our main result shows that if the vector field is infinitesimally contracting, and if the saltation matrix is non-expansive, then the intrinsic distance between any two trajectories decreases exponentially in time. When bounds on dwell time are available, our approach yields a bound on the intrinsic distance between trajectories regardless of whether the dynamics are expansive or contractive. We illustrate our results using wo examples: a constrained mechanical system and an electrical circuit with an ideal diode.more » « less
-
Infinitesimal contraction analysis, wherein global convergence results are obtained from properties of local dynamics, is a powerful analysis tool. In this letter, we generalize infinitesimal contraction analysis to hybrid systems in which state-dependent guards trigger transitions defined by reset maps between modes that may have different norms and need not be of the same dimension. In contrast to existing literature, we do not restrict mode sequence or dwell time. We work in settings where the hybrid system flow is differentiable almost everywhere and its derivative is the solution to a jump-linear-time-varying differential equation whose jumps are defined by a saltation matrix determined from the guard, reset map, and vector field. Our main result shows that if the vector field is infinitesimally contracting, and if the saltation matrix is non-expansive, then the intrinsic distance between any two trajectories decreases exponentially in time. When bounds on dwell time are available, our approach yields a bound on the intrinsic distance between trajectories regardless of whether the dynamics are expansive or contractive. We illustrate our results using two examples: a constrained mechanical system and an electrical circuit with an ideal diode.more » « less
-
Abstract We showed earlier that the level set function of a monotonic advancing front is twice differentiable everywhere with bounded second derivative and satisfies the equation classically. We show here that the second derivative is continuous if and only if the flow has a single singular time where it becomes extinct and the singular set consists of a closedC1manifold with cylindrical singularities. © 2017 Wiley Periodicals, Inc.more » « less
An official website of the United States government

