检索结果-内蒙古大学图书馆

Learning fluid physics from highly turbulent data using sparse physics-informed discovery of empirical relations (SPIDER)

学校读者我要写书评

暂无评论

arXiv 2021年

作者： Gurevich, Daniel R. Golden, Matthew R. Reinbold, Patrick A.K. Grigoriev, Roman O. Program in Applied and Computational Mathematics Princeton University PrincetonNJ08544 United States School of Physics Georgia Institute of Technology AtlantaGA30332 United States

We show how a complete mathematical model of a physical process can be learned directly from data via a hybrid approach combining three simple and general ingredients: physical assumptions of smoothness, locality, and symmetry, a weak formulation of differential equations, and sparse regression. To illustrate this, we extract a complete system of governing equations of fluid dynamics—the Navier-Stokes equation, the continuity equation, and the boundary conditions—as well as the pressure-Poisson and energy equations, from numerical data describing a highly turbulent channel flow in three dimensions. Whether they represent known or unknown physics, relations discovered using this approach take the familiar form of partial differential equations, which are easily interpretable and readily provide information about the relative importance of different physical effects. The proposed approach offers insight into the quality of the data, serving as a useful diagnostic tool. It is also remarkably robust, yielding accurate results for very high noise levels, and should thus be well-suited for analysis of experimental data. Copyright © 2021, The Authors. All rights reserved.

关键词： Navier Stokes equations

Generating large disordered stealthy hyperuniform systems with ultrahigh accuracy to determine their physical properties

学校读者我要写书评

暂无评论

Physical Review Research 2023年第3期5卷 033190-033190页

作者： Peter K. Morse Jaeuk Kim Paul J. Steinhardt Salvatore Torquato Department of Chemistry Princeton University Princeton New Jersey 08544 USA Department of Physics Princeton University Princeton New Jersey 08544 USA Princeton Institute of Materials Princeton University Princeton New Jersey 08544 USA Princeton Center for Theoretical Science Princeton University Princeton New Jersey 08544 USA Program in Applied and Computational Mathematics Princeton University Princeton New Jersey 08544 USA

Hyperuniform many-particle systems are characterized by a structure factor S(k) that is precisely zero as |k|→0; and stealthy hyperuniform systems have S(k)=0 for the finite range 0<|k|≤K, called the “exclusion region.” Through a process of collective-coordinate optimization, energy-minimizing disordered stealthy hyperuniform systems of moderate size have been made to high accuracy, and their novel physical properties have shown great promise. However, minimizing S(k) in the exclusion region is computationally intensive as the system size becomes large. In this paper, we present an improved methodology to generate such states using double-double precision calculations on graphical processing units (GPUs) that reduces the deviations from zero within the exclusion region by a factor of approximately 1030 for system sizes more than an order of magnitude larger. We further show that this ultrahigh accuracy is required to draw conclusions about their corresponding characteristics, such as the nature of the associated energy landscape and the presence or absence of Anderson localization, which might be masked, even when deviations are relatively small.

关键词： Exotic phases of matter

Heat transport in liquid water from first-principles and deep neural network simulations

学校读者我要写书评

暂无评论

Physical Review B 2021年第22期104卷 224202-224202页

作者： Davide Tisi Linfeng Zhang Riccardo Bertossa Han Wang Roberto Car Stefano Baroni SISSA–Scuola Internazionale Superiore di Studi Avanzati 34136 Trieste Italy Program in Applied and Computational Mathematics Princeton University Princeton New Jersey 08544 USA Laboratory of Computational Physics Institute of Applied Physics and Computational Mathematics Huayuan Road 6 Beijing 100088 People's Republic of China Department of Chemistry Department of Physics and Princeton Institute for the Science and Technology of Materials Princeton University Princeton New Jersey 08544 USA CNR Istituto Officina dei Materiali SISSA unit 34136 Trieste Italy

We compute the thermal conductivity of water within linear response theory from equilibrium molecular dynamics simulations, by adopting two different approaches. In one, the potential energy surface (PES) is derived on the fly from the electronic ground state of density functional theory (DFT) and the corresponding analytical expression is used for the energy flux. In the other, the PES is represented by a deep neural network (DNN) trained on DFT data, whereby the PES has an explicit local decomposition and the energy flux takes a particularly simple expression. By virtue of a gauge invariance principle, established by Marcolongo, Umari, and Baroni, the two approaches should be equivalent if the PES were reproduced accurately by the DNN model. We test this hypothesis by calculating the thermal conductivity, at the GGA (PBE) level of theory, using the direct formulation and its DNN proxy, finding that both approaches yield the same conductivity, in excess of the experimental value by approximately 60%. Besides being numerically much more efficient than its direct DFT counterpart, the DNN scheme has the advantage of being easily applicable to more sophisticated DFT approximations, such as meta-GGA and hybrid functionals, for which it would be hard to derive analytically the expression of the energy flux. We find in this way that a DNN model, trained on meta-GGA (SCAN) data, reduces the deviation from experiment of the predicted thermal conductivity by about 50%, leaving the question open as to whether the residual error is due to deficiencies of the functional, to a neglect of nuclear quantum effects in the atomic dynamics, or, likely, to a combination of the two.

关键词： Thermal conductivity Thermal properties Water Density functional theory Linear response theory Machine learning Molecular dynamics

Machine-learning-based non-Newtonian fluid model with molecular fidelity

学校读者我要写书评

暂无评论

Physical Review E 2020年第4期102卷 043309-043309页

作者： Huan Lei Lei Wu Weinan E Department of Computational Mathematics Science & Engineering and Department of Statistics & Probability Michigan State University East Lansing Michigan 48824 USA Department of Mathematics and Program in Applied and Computational Mathematics Princeton University Princeton New Jersey 08544 USA

We introduce a machine-learning-based framework for constructing continuum a non-Newtonian fluid dynamics model directly from a microscale description. Dumbbell polymer solutions are used as examples to demonstrate the essential ideas. To faithfully retain molecular fidelity, we establish a micro-macro correspondence via a set of encoders for the microscale polymer configurations and their macroscale counterparts, a set of nonlinear conformation tensors. The dynamics of these conformation tensors can be derived from the microscale model, and the relevant terms can be parametrized using machine learning. The final model, named the deep non-Newtonian model (DeePN2), takes the form of conventional non-Newtonian fluid dynamics models, with a generalized form of the objective tensor derivative that retains the microscale interpretations. Both the formulation of the dynamic equation and the neural network representation rigorously preserve the rotational invariance, which ensures the admissibility of the constructed model. Numerical results demonstrate the accuracy of DeePN2 where models based on empirical closures show limitations.

关键词： Learning Neuronal networks Complex fluids Non-Newtonian fluids Rheology Viscoelasticity Multiscale modeling

An analysis of reconstruction noise from undersampled 4D flow MRI

学校读者我要写书评

暂无评论

arXiv 2022年

作者： Partin, Lauren Schiavazzi, Daniele E. Long, Carlos A. Sing Department Of Applied And Computational Mathematics And Statistics University Of Notre Dame Notre DameIN United States Institute For Mathematical And Computational Engineering Institute For Biological And Medical Engineering Pontificia Universidad Católica De Chile Santiago Chile ANID - Millennium Science Initiative Program Millennium Nucleus Center For The Discovery Of Structures In Complex Data Chile ANID - Millennium Science Initiative Program Millennium Nucleus Center For Cardiovascular Magnetic Resonance Chile

Novel Magnetic Resonance (MR) imaging modalities can quantify hemodynamics but require long acquisition times, precluding its widespread use for early diagnosis of cardiovascular disease. To reduce the acquisition times, reconstruction methods from undersampled measurements are routinely used, that leverage representations designed to increase image compressibility. Reconstructed anatomical and hemodynamic images may present visual artifacts. Although some of these artifact are essentially reconstruction errors, and thus a consequence of undersampling, others may be due to measurement noise or the random choice of the sampled frequencies. Said otherwise, a reconstructed image becomes a random variable, and both its bias and its covariance can lead to visual artifacts;the latter leads to spatial correlations that may be misconstrued for visual information. Although the nature of the former has been studied in the literature, the latter has not received as much attention. In this study, we investigate the theoretical properties of the random perturbations arising from the reconstruction process, and perform a number of numerical experiments on simulated and MR aortic flow. Our results show that the correlation length remains limited to two to three pixels when a Gaussian undersampling pattern is combined with recovery algorithms based on 1-norm minimization. However, the correlation length may increase significantly for other undersampling patterns, higher undersampling factors (i.e., 8× or 16× compression), and different reconstruction methods. © 2022, CC BY-NC-ND.

关键词： Magnetic resonance

Earthmover-Based Manifold Learning for Analyzing Molecular Conformation Spaces

学校读者我要写书评

暂无评论

Earthmover-Based Manifold Learning for Analyzing Molecular C...

IEEE International Symposium on Biomedical Imaging

作者： Nathan Zelesko Amit Moscovich Joe Kileel Amit Singer Brown University Program in Applied and Computational Mathematics Princeton University Princeton University

ISBN: (数字)9781538693308

ISBN: (纸本)9781538693315

In this paper, we propose a novel approach for manifold learning that combines the Earthmover's distance (EMD) with the diffusion maps method for dimensionality reduction. We demonstrate the potential benefits of this approach for learning shape spaces of proteins and other flexible macromolecules using a simulated dataset of 3-D density maps that mimic the non-uniform rotary motion of ATP synthase. Our results show that EMD-based diffusion maps require far fewer samples to recover the intrinsic geometry than the standard diffusion maps algorithm that is based on the Euclidean distance. To reduce the computational burden of calculating the EMD for all volume pairs, we employ a wavelet-based approximation to the EMD which reduces the computation of the pairwise EMDs to a computation of pairwise weighted -ℓ 1 distances between wavelet coefficient vectors.

关键词： Manifolds Wavelet transforms Euclidean distance Eigenvalues and eigenfunctions Axles Geometry

On the Particle Approximation of Lagged Feynman-Kac Formulae

学校读者我要写书评

暂无评论

arXiv 2024年

作者： Awadelkarim, Elsiddig Caffarel, Michel Moral, Pierre Del Jasra, Ajay Applied Mathematics and Computational Science Program Computer Electrical and Mathematical Sciences and Engineering Division King Abdullah University of Science and Technology Thuwal23955-6900 Saudi Arabia University of Toulouse CNRS-Lab. de Chimie et Physique Quantiques Toulouse31062 France Centre de Recherche Inria Bordeaux Sud-Ouest Talence33405 France School of Data Science The Chinese University Hong Kong Shenzhen China

In this paper we examine the numerical approximation of the limiting invariant measure associated with Feynman-Kac formulae. These are expressed in a discrete time formulation and are associated with a Markov chain and a potential function. The typical application considered here is the computation of eigenvalues associated with non-negative operators as found, for example, in physics or particle simulation of rare-events. We focus on a novel lagged approximation of this invariant measure, based upon the introduction of a ratio of time-averaged Feynman-Kac marginals associated with a positive operator iterated l ∈ N times;a lagged Feynman-Kac formula. This estimator and its approximation using Diffusion Monte Carlo (DMC) have been extensively employed in the physics literature. In short, DMC is an iterative algorithm involving N ∈ N particles or walkers simulated in parallel, that undergo sampling and resampling operations. In this work, it is shown that for the DMC approximation of the lagged Feynman-Kac formula, one has an almost sure characterization of the L1-error as the time parameter (iteration) goes to infinity and this is at most of O(exp{−κl}/N), for_ κ > 0. In addition a non-asymptotic in time, and time uniform L1−bound is proved which is O(l/√N). We also prove a novel central limit theorem to give a characterization of the exact asymptotic in time variance. This analysis demonstrates that the strategy used in physics, namely, to run DMC with N and l small and, for long time enough, is mathematically justified. Our results also suggest how one should choose N and l in practice. We emphasize that these results are not restricted to physical applications;they have broad relevance to the general problem of particle simulation of the Feynman-Kac formula, which is utilized in a great variety of scientific and engineering fields. © 2024, CC BY.

关键词： Eigenvalues and eigenfunctions

THE COMPLEX INTERPLAY BETWEEN RISK TOLERANCE AND THE SPREAD OF INFECTIOUS DISEASES

学校读者我要写书评

暂无评论

arXiv 2024年

作者： Nguyen, Maximilian Freedman, Ari Cheung, Matthew Saad-Roy, Chadi Espinoza, Baltazar Grenfell, Bryan Levin, Simon Lewis-Sigler Institute Princeton University PrincetonNJ08544 United States Department of Ecology and Evolutionary Biology Princeton University PrincetonNJ United States Program in Applied and Computational Mathematics Princeton University PrincetonNJ United States Miller Institute for Basic Research in Science Department of Integrative Biology University of California Berkeley BerkeleyCA United States Biocomplexity Institute University of Virginia CharlottesvilleVA United States

Risk-driven behavior provides a feedback mechanism through which individuals both shape and are collectively affected by an epidemic. We introduce a general and flexible compartmental model to study the effect of heterogeneity in the population with regards to risk tolerance. The interplay between behavior and epidemiology leads to a rich set of possible epidemic dynamics. Depending on the behavioral composition of the population, we find that increasing heterogeneity in risk tolerance can either increase or decrease the epidemic size. We find that multiple waves of infection can arise due to the interplay between transmission and behavior, even without the replenishment of susceptibles. We find that increasing protective mechanisms such as the effectiveness of interventions, the number of risk-averse people in the population, and the duration of intervention usage reduces the epidemic overshoot. When the protection is pushed past a critical threshold, the epidemic dynamics enter an underdamped regime where the epidemic size exactly equals the herd immunity threshold. Lastly, we can find regimes where epidemic size does not monotonically decrease with a population that becomes increasingly risk-averse. Copyright © 2024, The Authors. All rights reserved.

关键词： Epidemiology

On the local well-posedness for the relativistic euler equations for a liquid body

学校读者我要写书评

暂无评论

arXiv 2021年

作者： Ginsberg, Daniel Lindblad, Hans Program in Applied and Computational Mathematics Princeton University PrincetonNJ08544 United States Johns Hopkins University Department of Mathematics 3400 N. Charles St. BaltimoreMD21218 United States

We prove a local existence theorem for the free boundary problem for a relativistic fluid in a fixed spacetime. Our proof involves an a priori estimate which only requires control of derivatives tangential to the boundary, which holds also in the Newtonian compressible case. Copyright © 2021, The Authors. All rights reserved.

关键词：