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

Composite B-Spline Regularized Delta Functions for the Immersed Boundary Method: Divergence-Free Interpolation and Gradient-Preserving Force Spreading

学校读者我要写书评

暂无评论

arXiv 2024年

作者： Gruninger, Cole Griffith, Boyce E. Department of Mathematics University North Carolina Chapel HillNC United States Department of Applied Physical Sciences University of North Carolina Chapel HillNC United States Department of Biomedical Engineering University of North Carolina Chapel HillNC United States Carolina Center for Interdisciplinary Applied Mathematics University of North Carolina Chapel HillNC United States Computational Medicine Program University of North Carolina School of Medicine Chapel HillNC United States McAllister Heart Institute University of North Carolina School of Medicine Chapel HillNC United States

This paper presents an approach to enhance volume conservation in the immersed boundary (IB) method by employing regularized delta functions derived from composite B-splines. These delta functions are constructed using tensor product kernels, similar to the conventional IB method. However, the kernels are B-splines whose polynomial degree varies according to the normal and tangential directions of each velocity component. The conventional IB method, while effective for fluid-structure interaction applications, has long been challenged by poor volume conservation, particularly evident in simulations of pressurized, closed membranes. We demonstrate that composite B-spline regularized delta functions significantly enhance volume conservation through two complementary properties: they provide continuously divergence-free velocity interpolants and maintain the gradient character of forces corresponding to mean pressure jumps across interfaces. By correctly representing these forces as discrete gradients, they eliminate a key source of spurious flows that typically plague immersed boundary computations. Our approach maintains the local nature of the classical IB method, avoiding the computational overhead associated with the non-local Divergence-Free Immersed Boundary (DFIB) method’s construction of an explicit velocity potential which requires additional Poisson solves for interpolation and force spreading operations. Through a series of numerical experiments, we show that sufficiently regular composite B-spline kernels can maintain initial volumes to within machine precision. We provide a detailed analysis of the relationship between kernel regularity and the accuracy of force spreading and velocity interpolation operations. Our findings indicate that composite B-splines of at least C1 regularity produce results comparable to the DFIB method in dynamic simulations, with errors in volume conservation primarily dominated by truncation error of the employed time-stepping s

关键词： Delta functions

Local Divergence-Free Immersed Finite Element-Difference Method Using Composite B-Splines

学校读者我要写书评

暂无评论

arXiv 2024年

作者： Li, Lianxia Gruninger, Cole Lee, Jae H. Griffith, Boyce E. Department of Mathematics University of North Carolina Chapel HillNC United States Department of Biomedical Engineering University of North Carolina Chapel HillNC United States Carolina Center for Interdisciplinary Applied Mathematics University of North Carolina Chapel HillNC United States Computational Medicine Program University of North Carolina School of Medicine Chapel HillNC United States McAllister Heart Institute University of North Carolina School of Medicine Chapel HillNC United States

In the class of immersed boundary (IB) methods, the choice of the regularized delta function plays a crucial role in transferring information between fluid and solid domains through interpolation and spreading operators. Most prior work using the IB method has used isotropic kernels that do not preserve the divergence-free condition of the velocity field, leading to loss of incompressibility of the solid when interpolating the Eulerian velocity to Lagrangian markers. To address this issue, in simulations involving large deformations of incompressible hyperelastic structures immersed in fluid, researchers often use a volumetric stabilization approach such as adding a volumetric energy term and using modified invariants in the constitutive model of the immersed structure. Composite B-spline (CBS) kernels offer an alternative approach by inherently maintaining the discrete divergence-free property. This work evaluates the performance of CBS kernels in terms of their volume conservation and accuracy, comparing them with several traditional isotropic kernel functions using a construction introduced by Peskin (referred to as IB kernels) and B-spline (BS) kernels. Benchmark tests include pressure-loaded and shear-dominated flows, such as an elastic band under differential pressure loads, a pressurized membrane, a compressed block, Cook’s membrane, a slanted channel flow, and a modified Turek-Hron problem. Additionally, we validate our methodology using a complex fluid-structure interaction model of bioprosthetic heart valve dynamics in a pulse duplicator. Results demonstrate that CBS kernels achieve superior volume conservation compared to conventional isotropic kernels, eliminating the need for additional volumetric stabilization techniques typically required to address instabilities arising from volume conservation errors. Further, CBS kernels achieve convergence on coarser fluid grids, while IB and BS kernels need finer grids for comparable accuracy. Unlike IB and BS ke

关键词： Delta functions

Structural characterization of many-particle systems on approach to hyperuniform states

学校读者我要写书评

暂无评论

arXiv 2021年

作者： Torquato, Salvatore Department of Chemistry Department of Physics Princeton Institute for the Science and Technology of Materials Program in Applied and Computational Mathematics Princeton University PrincetonNJ08544 United States

The study of hyperuniform states of matter is an emerging multidisciplinary field, impinging on topics in the physical sciences, mathematics and biology. The focus of this work is the exploration of quantitative descriptors that herald when a many-particle system in d-dimensional Euclidean space Rd approaches a hyperuniform state as a function of the relevant control parameter. We establish quantitative criteria to ascertain the extent of hyperuniform and nonhyperuniform distance-scaling regimes as well as the crossover point between them in terms of the "volume" coefficient A and "surface-area" coefficient B associated with the local number variance σ2(R) for a spherical window of radius R. The larger the ratio B/A, the larger the hyperuniform scaling regime, which becomes of infinite extent in the limit B/A → ∞. To complement the known direct-space representation of the coefficient B in terms of the total correlation function h(r), we derive its corresponding Fourier representation in terms of the structure factor S(k), which is especially useful when scattering information is available experimentally or theoretically. We also demonstrate that the free-volume theory of the pressure of equilibrium packings of identical hard spheres that approach a strictly jammed state either along the stable crystal or metastable disordered branch dictates that such end states be exactly hyperuniform. Using the ratio B/A, as well as other diagnostic measures of hyperuniformity, including the hyperuniformity index H and the direct-correlation function length scale ξc, we study three different exactly solvable models as a function of the relevant control parameter, either density or temperature, with end states that are perfectly hyperuniform. Specifically, we analyze equilibrium systems of hard rods and "sticky" hard-sphere systems in arbitrary space dimension d as a function of density. We also examine low-temperature excited states of many-particle systems interacting with "steal

关键词： Temperature

ON THE CONVERGENCE OF GRADIENT DESCENT TRAINING FOR TWO-LAYER RELU-NETWORKS IN THE MEAN FIELD REGIME

学校读者我要写书评

暂无评论

arXiv 2020年

作者： Wojtowytsch, Stephan Program in Applied and Computational Mathematics Princeton University PrincetonNJ08544

We describe a necessary and sufficient condition for the convergence to minimum Bayes risk when training two-layer ReLU-networks by gradient descent in the mean field regime with omni-directional initial parameter distribution. This article extends recent results of Chizat and Bach to ReLU-activated networks and to the situation in which there are no parameters which exactly achieve MBR. The condition does not depend on the initalization of parameters and concerns only the weak convergence of the realization of the neural network, not its parameter *** Codes 35Q68, 68T07, 49Q22, 35F20 Copyright © 2020, The Authors. All rights reserved.

关键词： Gradient methods

Local order metrics for many-particle systems across length scales

学校读者我要写书评

暂无评论

arXiv 2024年

作者： Maher, Charles Emmett Torquato, Salvatore Department of Chemistry Princeton University New Jersey08544 United States Department of Chemistry Princeton University PrincetonNJ08544 United States Department of Physics Princeton University PrincetonNJ08544 United States Princeton Institute for the Science and Technology of Materials Princeton University PrincetonNJ08544 United States Program in Applied and Computational Mathematics Princeton University PrincetonNJ08544 United States

Formulating order metrics that sensitively quantify the degree of order/disorder in many-particle systems in d-dimensional Euclidean space Rd across length scales is an outstanding challenge in physics, chemistry, and materials science. Since an infinite set of n-particle correlation functions is required to fully characterize a system, one must settle for a reduced set of structural information, in practice. We initiate a program to use the local number variance σN2 (R) associated with a spherical sampling window of radius R (which encodes pair correlations) and an integral measure derived from it ΣN(Ri, Rj) that depends on two specified radial distances Ri and Rj. Across the first three space dimensions (d = 1, 2, 3), we find these metrics can sensitively describe and categorize the degree of order/disorder of 41 different models of antihyperuniform, nonhyperuniform, disordered hyperuniform, and ordered hyperuniform many-particle systems at a specified length scale R. Using our local variance metrics, we demonstrate the importance of assessing order/disorder with respect to a specific value of R. These local order metrics could also aid in the inverse design of structures with prescribed length-scale-specific degrees of order/disorder that yield desired physical properties. In future work, it would be fruitful to explore the use of higher-order moments of the number of points within a spherical window of radius R [S. Torquato et al., Phys. Rev. X, 11, 021028 (2021)] to devise even more sensitive order metrics. © 2024, CC BY.

关键词： Inverse problems

Ordered and disordered stealthy hyperuniform point patterns across spatial dimensions

学校读者我要写书评

暂无评论

arXiv 2024年

作者： Morse, Peter K. Steinhardt, Paul J. Torquato, Salvatore Department of Chemistry Princeton University PrincetonNJ08544 United States Department of Physics Princeton University PrincetonNJ08544 United States Princeton Institute of Materials Princeton University PrincetonNJ08544 United States Princeton Center for Theoretical Science Princeton University PrincetonNJ08544 United States Program in Applied and Computational Mathematics Princeton University PrincetonNJ08544 United States

In previous work [Phys. Rev. X 5, 021020 (2015)], it was shown that stealthy hyperuniform systems can be regarded as hard spheres in Fourier-space in the sense that the the structure factor is exactly zero in a spherical region around the origin in analogy with the pair-correlation function of real-space hard spheres. In this work, we exploit this correspondence to confirm that the densest Fourier-space hard-sphere system is that of a Bravais lattice. This is in contrast to real-space hard-spheres, whose densest configuration is conjectured to be disordered. We also extend the virial series previously suggested for disordered stealthy hyperuniform systems to higher dimensions in order to predict spatial decorrelation as function of dimension. This prediction is then borne out by numerical simulations of disordered stealthy hyperuniform ground states in dimensions d = 2-8. © 2024, CC BY.

关键词： Ground state

Hyperuniformity of maximally random jammed packings of hyperspheres across spatial dimensions

学校读者我要写书评

暂无评论

Physical Review E 2023年第6期108卷 064602-064602页

作者： Charles Emmett Maher Yang Jiao Salvatore Torquato Department of Chemistry Princeton University Princeton New Jersey 08544 USA Materials Science and Engineering Arizona State University Tempe Arizona 85287 USA Department of Physics Arizona State University Tempe Arizona 85287 USA Department of Physics Princeton University Princeton New Jersey 08544 USA Princeton Materials Institute Princeton University Princeton New Jersey 08544 USA Program in Applied and Computational Mathematics Princeton University Princeton New Jersey 08544 USA

The maximally random jammed (MRJ) state is the most random (i.e., disordered) configuration of strictly jammed (mechanically rigid) nonoverlapping objects. MRJ packings are hyperuniform, meaning their long-wavelength density fluctuations are anomalously suppressed compared to typical disordered systems, i.e., their structure factors S(k) tend to zero as the wave number |k| tends to zero. Here we show that generating high-quality strictly jammed states for Euclidean space dimensions d=3,4, and 5 is of paramount importance in ensuring hyperuniformity and extracting precise values of the hyperuniformity exponent α>0 for MRJ states, defined by the power-law behavior of S(k)∼|k|α in the limit |k|→0. Moreover, we show that for fixed d it is more difficult to ensure jamming as the particle number N increases, which results in packings that are nonhyperuniform. Free-volume theory arguments suggest that the ideal MRJ state does not contain rattlers, which act as defects in numerically generated packings. As d increases, we find that the fraction of rattlers decreases substantially. Our analysis of the largest truly jammed packings suggests that the ideal MRJ packings for all dimensions d≥3 are hyperuniform with α=d−2, implying the packings become more hyperuniform as d increases. The differences in α between MRJ packings and the recently proposed Manna-class random close packed (RCP) states, which were reported to have α=0.25 in d=3 and be nonhyperuniform (α=0) for d=4 and d=5, demonstrate the vivid distinctions between the large-scale structure of RCP and MRJ states in these dimensions. Our paper clarifies the importance of the link between true jamming and hyperuniformity and motivates the development of an algorithm to produce rattler-free three-dimensional MRJ packings.

关键词： Granular solids

Structural properties of hyperuniform Voronoi networks

学校读者我要写书评

暂无评论

Physical Review E 2025年第3期111卷 034123-034123页

作者： Eli Newby Wenlong Shi Yang Jiao Reka Albert Salvatore Torquato Department of Physics Pennsylvania State University University Park Pennsylvania 16802 USA Materials Science and Engineering Arizona State University Tempe Arizona 85287 USA Department of Physics Arizona State University Tempe Arizona 85287 USA 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 Program in Applied and Computational Mathematics Princeton University Princeton New Jersey 08544 USA

Disordered hyperuniform many-particle systems are recently discovered exotic states of matter, characterized by the complete suppression of normalized infinite-wavelength density fluctuations, as in perfect crystals, while lacking conventional long-range order, as in liquids and glasses. In this work, we begin a program to quantify the structural properties of nonhyperuniform and hyperuniform networks. In particular, large two-dimensional (2D) Voronoi networks (graphs) containing approximately 10,000 nodes are created from a variety of different point configurations, including the antihyperuniform hyperplane intersection process (HIP), nonhyperuniform Poisson process, nonhyperuniform random sequential addition (RSA) saturated packing, and both non-stealthy and stealthy hyperuniform point processes. We carry out an extensive study of the Voronoi-cell area distribution of each of the networks by determining multiple metrics that characterize the distribution, including their average areas and corresponding variances as well as higher-order cumulants (i.e., skewness γ1 and excess kurtosis γ2). We show that the HIP distribution is far from Gaussian, as evidenced by a high skewness (γ1=3.16) and large positive excess kurtosis (γ2=16.2). The Poisson (with γ1=1.07 and γ2=1.79) and non-stealthy hyperuniform (with γ1=0.257 and γ2=0.0217) distributions are Gaussian-like distributions, since they exhibit a small but positive skewness and excess kurtosis. The RSA (with γ1=0.450 and γ2=−0.0384) and the highest stealthy hyperuniform distributions (with γ1=0.0272 and γ2=−0.0626) are also non-Gaussian because of their low skewness and negative excess kurtosis, which is diametrically opposite of the non-Gaussian behavior of the HIP. The fact that the cell-area distributions of large, finite-sized RSA and stealthy hyperuniform networks (e.g., with N≈10,000 nodes) are narrower, have larger peaks, and smaller tails than a Gaussian distribution implies that in the thermodynamic limit th

关键词： Gaussian distribution

Mutually beneficial confluence of structure-based modeling of protein dynamics and machine learning methods

学校读者我要写书评

暂无评论

Current Opinion in Structural Biology 2023年 78卷 102517页

作者： Banerjee, Anupam Saha, Satyaki Tvedt, Nathan C. Yang, Lee-Wei Bahar, Ivet Computational and Systems Biology University of Pittsburgh School of Medicine Pittsburgh 15261 PA United States Computational and Applied Mathematics and Statistics The College of William and Mary Williamsburg 23185 VA United States Institute of Bioinformatics and Structural Biology and PhD Program in Biomedical Artificial Intelligence National Tsing Hua University Hsinchu 300044 Taiwan Physics Division National Center for Theoretical Sciences Taipei 106319 Taiwan

Proteins sample an ensemble of conformers under physiological conditions, having access to a spectrum of modes of motions, also called intrinsic dynamics. These motions ensure the adaptation to various interactions in the cell, and largely assist in, if not determine, viable mechanisms of biological function. In recent years, machine learning frameworks have proven uniquely useful in structural biology, and recent studies further provide evidence to the utility and/or necessity of considering intrinsic dynamics for increasing their predictive ability. Efficient quantification of dynamics-based attributes by recently developed physics-based theories and models such as elastic network models provides a unique opportunity to generate data on dynamics for training ML models towards inferring mechanisms of protein function, assessing pathogenicity, or estimating binding affinities. © 2022 The Authors

关键词：