Search

The eigenfunctions of the Laplace operator are essential in mathematical physics, engineering, and geometry processing. Typically, these are computed by discretizing the domain and performing eigendecomposition, tying the results to a specific mesh. However, this method is unsuitable for continuously-parameterized shapes. We propose a novel representation for eigenfunctions in continuously-parameterized shape spaces, where eigenfunctions are spatial fields with continuous dependence on shape parameters, defined by minimal Dirichlet energy, unit norm, and mutual orthogonality. We implement this with multilayer perceptrons trained as neural fields, mapping shape parameters and domain positions to eigenfunction values. A unique challenge is enforcing mutual orthogonality with respect to causality, where the causal ordering varies across the shape space. Our training method therefore requires three interwoven concepts: (1) learning $n$ eigenfunctions concurrently by minimizing Dirichlet energy with unit norm constraints; (2) filtering gradients during backpropagation to enforce causal orthogonality, preventing earlier eigenfunctions from being influenced by later ones; (3) dynamically sorting the causal ordering based on eigenvalues to track eigenvalue curve crossovers. We demonstrate our method on problems such as shape family analysis, predicting eigenfunctions for incomplete shapes, interactive shape manipulation, and computing higher-dimensional eigenfunctions, on all of which traditional methods fall short.

We propose a hybrid neural network and physics framework for reduced-order modeling of elastoplasticity and fracture. State-of-the-art scientific computing models like the Material Point Method (MPM) faithfully simulate large-deformation elastoplasticity and fracture mechanics. However, their long runtime and large memory consumption render them unsuitable for applications constrained by computation time and memory usage, e.g., virtual reality. To overcome these barriers, we propose a reduced-order framework. Our key innovation is training a low-dimensional manifold for the Kirchhoff stress field via an implicit neural representation. This low-dimensional neural stress field (NSF) enables efficient evaluations of stress values and, correspondingly, internal forces at arbitrary spatial locations. In addition, we also train neural deformation and affine fields to build low-dimensional manifolds for the deformation and affine momentum fields. These neural stress, deformation, and affine fields share the same low-dimensional latent space, which uniquely embeds the high-dimensional simulation state. After training, we run new simulations by evolving in this single latent space, which drastically reduces the computation time and memory consumption. Our general continuum-mechanics-based reduced-order framework is applicable to any phenomena governed by the elastodynamics equation. To showcase the versatility of our framework, we simulate a wide range of material behaviors, including elastica, sand, metal, non-Newtonian fluids, fracture, contact, and collision. We demonstrate dimension reduction by up to 100,000X and time savings by up to 10X.

Linear reduced-order modeling (ROM) simplifies complex simulations by approximating the behavior of a system using a simplified kinematic representation. Typically, ROM is trained on input simulations created with a specific spatial discretization, and then serves to accelerate simulations with the same discretization. This discretization-dependence is restrictive. Becoming independent of a specific discretization would provide flexibility to mix and match mesh resolutions, connectivity, and type (tetrahedral, hexahedral) in training data; to accelerate simulations with novel discretizations unseen during training; and to accelerate adaptive simulations that temporally or parametrically change the discretization. We present a flexible, discretization-independent approach to reduced-order modeling. Like traditional ROM, we represent the configuration as a linear combination of displacement fields. Unlike traditional ROM, our displacement fields are continuous maps from every point on the reference domain to a corresponding displacement vector; these maps are represented as implicit neural fields. With linear continuous ROM (LiCROM), our training set can include multiple geometries undergoing multiple loading conditions, independent of their discretization. This opens the door to novel applications of reduced order modeling. We can now accelerate simulations that modify the geometry at runtime, for instance via cutting, hole punching, and even swapping the entire mesh. We can also accelerate simulations of geometries unseen during training. We demonstrate one-shot generalization, training on a single geometry and subsequently simulating various unseen geometries.

We outline a partial historical summary of the steps through which the nucleosynthesis phenomena induced by {\it slow} neutron captures (the {\it s-process}) were clarified, a scientific achievement in which Franz K\"appeler played a major role. We start by recalling the early phenomenological approach, which yielded a basic understanding of the subject even before models for the parent stellar evolutionary stages were developed. Through such a tool, rough limits for the neutron density and exposure were set, and the crucial fact was understood that more than one nucleosynthesis component is required to account for solar abundances of $s$-process nuclei up to the Pb-Bi region. We then summarize the gradual understanding of the stellar processes actually involved in the production of nuclei from Sr to Pb (the so-called {\it Main Component}, achieved in the last decade of the past century and occurring in red giants of low and intermediate mass, ($M \lesssim$ 8 $M_{\odot}$), populating, in the {\it HR} diagram, the {\it Asymptotic Giant Branch} or {\it AGB} region. We conclude by giving some details on more recent research concerning mixing mechanisms inducing the activation of the main neutron source, $^{13}$C($\alpha$,n)$^{16}$O., Comment: 35 page, 8 figures

The long runtime of high-fidelity partial differential equation (PDE) solvers makes them unsuitable for time-critical applications. We propose to accelerate PDE solvers using reduced-order modeling (ROM). Whereas prior ROM approaches reduce the dimensionality of discretized vector fields, our continuous reduced-order modeling (CROM) approach builds a low-dimensional embedding of the continuous vector fields themselves, not their discretization. We represent this reduced manifold using continuously differentiable neural fields, which may train on any and all available numerical solutions of the continuous system, even when they are obtained using diverse methods or discretizations. We validate our approach on an extensive range of PDEs with training data from voxel grids, meshes, and point clouds. Compared to prior discretization-dependent ROM methods, such as linear subspace proper orthogonal decomposition (POD) and nonlinear manifold neural-network-based autoencoders, CROM features higher accuracy, lower memory consumption, dynamically adaptive resolutions, and applicability to any discretization. For equal latent space dimension, CROM exhibits 79$\times$ and 49$\times$ better accuracy, and 39$\times$ and 132$\times$ smaller memory footprint, than POD and autoencoder methods, respectively. Experiments demonstrate 109$\times$ and 89$\times$ wall-clock speedups over unreduced models on CPUs and GPUs, respectively. Videos and codes are available on the project page: https://crom-pde.github.io

This work proposes a model-reduction approach for the material point method on nonlinear manifolds. Our technique approximates the $\textit{kinematics}$ by approximating the deformation map using an implicit neural representation that restricts deformation trajectories to reside on a low-dimensional manifold. By explicitly approximating the deformation map, its spatiotemporal gradients -- in particular the deformation gradient and the velocity -- can be computed via analytical differentiation. In contrast to typical model-reduction techniques that construct a linear or nonlinear manifold to approximate the (finite number of) degrees of freedom characterizing a given spatial discretization, the use of an implicit neural representation enables the proposed method to approximate the $\textit{continuous}$ deformation map. This allows the kinematic approximation to remain agnostic to the discretization. Consequently, the technique supports dynamic discretizations -- including resolution changes -- during the course of the online reduced-order-model simulation. To generate $\textit{dynamics}$ for the generalized coordinates, we propose a family of projection techniques. At each time step, these techniques: (1) Calculate full-space kinematics at quadrature points, (2) Calculate the full-space dynamics for a subset of `sample' material points, and (3) Calculate the reduced-space dynamics by projecting the updated full-space position and velocity onto the low-dimensional manifold and tangent space, respectively. We achieve significant computational speedup via hyper-reduction that ensures all three steps execute on only a small subset of the problem's spatial domain. Large-scale numerical examples with millions of material points illustrate the method's ability to gain an order of magnitude computational-cost saving -- indeed $\textit{real-time simulations}$ -- with negligible errors.

Universal meshes have recently appeared in the literature as a compu- tationally efficient and robust paradigm for the generation of conforming simpli- cial meshes for domains with evolving boundaries. The main idea behind a univer- sal mesh is to immerse the moving boundary in a background mesh (the universal mesh), and to produce a mesh that conforms to the moving boundary at any given time by adjusting a few of elements of the background mesh. In this manuscript we present the application of universal meshes to the simulation of brittle fracturing. To this extent, we provide a high level description of a crack propagation algorithm and showcase its capabilities. Alongside universal meshes for the simulation of brit- tle fracture, we provide other examples for which universal meshes prove to be a powerful tool, namely fluid flow past moving obstacles. Lastly, we conclude the manuscript with some remarks on the current state of universal meshes and future directions.

The use of the interaction integral to compute stress intensity factors around a crack tip requires selecting an auxiliary field and a material variation field. We formulate a family of these fields accounting for the curvilinear nature of cracks that, in conjunction with a discrete formulation of the interaction integral, yield optimally convergent stress intensity factors. We formulate three pairs of auxiliary and material variation fields chosen to yield a simple expression of the interaction integral for different classes of problems. The formulation accounts for crack face tractions and body forces. Distinct features of the fields are their ease of construction and implementation. The resulting stress intensity factors are observed converging at a rate that doubles the one of the stress field. We provide a sketch of the theoretical justification for the observed convergence rates, and discuss issues such as quadratures and domain approximations needed to attain such convergent behavior. Through two representative examples, a circular arc crack and a loaded power function crack, we illustrate the convergence rates of the computed stress intensity factors. The numerical results also show the independence of the method on the size of the domain of integration.

The relativistic precession can be quickly inferred from the nonlinear polar orbit equation without actually solving it., Comment: Accepted for publication in Astrophysics & Space Science

We study the general relativistic orbital equation and using a straightforward perturbation method and a mathematical device first introduced by d'Alembert, we work out approximate expressions of a bound planetary orbit in the form of trigonometrical polynomials and the first three terms of the power series development of the perihelion advance. The results are applied to a more precise determination of the total mass of the double pulsar J0737-3039., Comment: 8 pages. Accepted for publication in "Astrophysics & Space Science"