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

Deep learning of accurate force field of Ferroelectric HfO$_2$

学校读者我要写书评

暂无评论

arXiv 2020年

作者： Wu, Jing Zhang, Yuzhi Zhang, Linfeng Liu, Shi Fudan University Shanghai200433 China School of Science Westlake University HangzhouZhejiang310024 China Institute of Natural Sciences Westlake Institute for Advanced Study HangzhouZhejiang310024 China Yuanpei College Peking University Beijing100871 China Program in Applied and Computational Mathematics Princeton University PrincetonNJ08544 United States

The discovery of ferroelectricity in HfO2-based thin films opens up new opportunities for using this silicon-compatible ferroelectric to realize low-power logic circuits and high-density non-volatile memories. The functional performances of ferroelectrics are intimately related to their dynamic responses to external stimuli such as electric fields at finite temperatures. Molecular dynamics is an ideal technique for investigating dynamical processes on large length and time scales, though its applications to new materials is often hindered by the limited availability and accuracy of classical force fields. Here we present a deep neural network-based interatomic force field of HfO2 learned from ab initio data using a concurrent learning procedure. The model potential is able to predict structural properties such as elastic constants, equation of states, phonon dispersion relationships, and phase transition barriers of various hafnia polymorphs with accuracy comparable with density functional theory calculations. The validity of this model potential is further confirmed by the reproduction of experimental sequences of temperature-driven ferroelectric-paraelectric phase transitions of HfO2 with isobaric-isothermal ensemble molecular dynamics simulations. We suggest a general approach to extend the model potential of HfO2 to related material systems including dopants and defects. Copyright © 2020, The Authors. All rights reserved.

关键词： Ferroelectricity

A Nodal Immersed Finite Element-Finite Difference Method

学校读者我要写书评

暂无评论

arXiv 2021年

作者： Wells, David R. Vadala-Roth, Ben Lee, Jae H. Griffith, Boyce E. Department of Mathematics University of North Carolina Chapel HillNC United States Westborough MA United States Department of Mechanical Engineering Institute for Computational Medicine Johns Hopkins University BaltimoreMD United States Department of Mathematics Applied Physical Sciences and 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 Chapel HillNC United States McAllister Heart Institute University of North Carolina Chapel HillNC United States Center for Drug Evaluation and Research U.S. Food and Drug Administration Silver SpringMD United States

The immersed finite element-finite difference (IFED) method is a computational approach to modeling interactions between a fluid and an immersed structure. The IFED method uses a finite element (FE) method to approximate the stresses, forces, and structural deformations on a structural mesh and a finite difference (FD) method to approximate the momentum and enforce the incompressibility of the entire fluid-structure system on a Cartesian grid. The fundamental approach used by this method follows the immersed boundary framework for modeling fluid-structure interaction (FSI), in which a force spreading operator prolongs structural forces to a Cartesian grid, and a velocity interpolation operator restricts a velocity field defined on that grid back onto the structural mesh. With an FE structural mechanics framework, force spreading first requires that the force itself be projected onto the finite element space. Similarly, velocity interpolation requires projecting velocity data onto the FE basis functions. Consequently, evaluating either coupling operator requires solving a matrix equation at every time step. Mass lumping, in which the projection matrices are replaced by diagonal approximations, has the potential to accelerate this method considerably. This paper provides both numerical and computational analyses of the effects of this replacement for evaluating the force projection and for the IFED coupling operators. Constructing the coupling operators also requires determining the locations on the structure mesh where the forces and velocities are sampled. Here we show that sampling the forces and velocities at the nodes of the structural mesh is equivalent to using lumped mass matrices in the IFED coupling operators. A key theoretical result of our analysis is that if both of these approaches are used together, the IFED method permits the use of lumped mass matrices derived from nodal quadrature rules for any standard interpolatory element. This is different from s

关键词： Fluid structure interaction

Unbiased Estimation using a Class of Diffusion Processes

学校读者我要写书评

暂无评论

arXiv 2022年

作者： Ruzayqat, Hamza Beskos, Alexandros Crisan, Dan Jasra, Ajay Kantas, Nikolas Applied Mathematics and Computational Science Program Computer Electrical and Mathematical Sciences and Engineering Division King Abdullah University of Science and Technology Thuwal23955-6900 KSA Saudi Arabia Department of Statistical Science University College London LondonWC1E 6BT United Kingdom Department of Mathematics Imperial College London LondonSW7 2AZ United Kingdom

We study the problem of unbiased estimation of expectations with respect to (w.r.t.) π a given, general probability measure on (d, B(d)) that is absolutely continuous with respect to a standard Gaussian measure. We focus on simulation associated to a particular class of diffusion processes, sometimes termed the Schrödinger-Föllmer Sampler, which is a simulation technique that approximates the law of a particular diffusion bridge process {Xt}tϵ[0,1]on d, d ϵ 0. This latter process is constructed such that, starting at X0= 0, one has X1∼ π. Typically, the drift of the diffusion is intractable and, even if it were not, exact sampling of the associated diffusion is not possible. As a result, [10, 16] consider a stochastic Euler-Maruyama scheme that allows the development of biased estimators for expectations w.r.t. π. We show that for this methodology to achieve a mean square error of O(ϵ2), for arbitrary ϵ > 0, the associated cost is O(ϵ-5). We then introduce an alternative approach that provides unbiased estimates of expectations w.r.t. π, that is, it does not suffer from the time discretization bias or the bias related with the approximation of the drift function. We prove that to achieve a mean square error of O(ϵ2), the associated cost (which is random) is, with high probability, O(ϵ-2| log(ϵ)|2+δ), for any δ > 0. We implement our method on several examples including Bayesian inverse *** Codes 60J60, 62D05, 65C40 © 2022, CC BY.

关键词： Diffusion

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction

学校读者我要写书评

暂无评论

arXiv 2022年

作者： Kolahdouz, Ebrahim M. Wells, David R. Rossi, Simone Aycock, Kenneth I. Craven, Brent A. Griffith, Boyce E. Center for Computational Biology Flatiron Institute Simons Foundation New YorkNY United States Department of Mathematics University of North Carolina Chapel HillNC United States Division of Applied Mechanics Office of Science and Engineering Laboratories Center for Devices and Radiological Health United States Food and Drug Administration Silver SpringMD United States Departments of Mathematics and 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 Chapel HillNC United States McAllister Heart Institute University of North Carolina Chapel HillNC United States

This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressibl

关键词： Fluid structure interaction

A Nodal Immersed Finite Element-Finite Difference Method

学校读者我要写书评

暂无评论

SSRN

SSRN 2022年

作者： Wells, David Vadala-Roth, Ben Lee, Jae Ho Griffith, Boyce E. Department of Mathematics University of North Carolina Chapel HillNC United States U.S. Army Corps of Engineers Engineer Research and Development Center Coastal and Hydraulic Laboratory VicksburgMS United States Department of Mechanical Engineering Institute for Computational Medicine Johns Hopkins University BaltimoreMD United States Departments of Mathematics Applied Physical Sciences and 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 Chapel HillNC United States McAllister Heart Institute University of North Carolina Chapel HillNC United States

The immersed finite element-finite difference (IFED) method is a computational approach to modeling interactions between a fluid and an immersed structure. The IFED method uses a finite element (FE) method to approximate the stresses and forces on a structural mesh and a finite difference (FD) method to approximate the momentum of the entire fluid-structure system on a Cartesian grid. The fundamental approach used by this method follows the immersed boundary framework for modeling fluid-structure interaction (FSI), in which a force spreading operator prolongs structural forces to a Cartesian grid, and a velocity interpolation operator restricts a velocity field defined on that grid back onto the structural mesh. With an FE structural mechanics framework, force spreading first requires that the force itself be projected onto the finite element space. Similarly, velocity interpolation requires projecting velocity data onto the FE basis functions. Consequently, evaluating either coupling operator requires solving a matrix equation at every time step. Mass lumping, in which the projection matrices are replaced by diagonal approximations, has the potential to accelerate this method considerably. This paper provides both numerical and computational analyses of the effects of this replacement for evaluating the force projection and for the IFED coupling operators. Constructing the coupling operators also requires determining the locations on the structure mesh where the forces and velocities are sampled. Here we show that sampling the forces and velocities at the nodes of the structural mesh is equivalent to using lumped mass matrices in the IFED coupling operators. A key theoretical result of our analysis is that if both of these approaches are used together, the IFED method permits the use of lumped mass matrices derived from nodal quadrature rules for any standard interpolatory element. This is different from standard FE methods, which require specialized treatments to

关键词： Fluid structure interaction

An immersed peridynamics model of fluid-structure interaction accounting for material damage and failure

学校读者我要写书评

暂无评论

arXiv 2022年

作者： Kim, Keon Ho Bhalla, Amneet P.S. Griffith, Boyce E. Department of Mathematics University of North Carolina Chapel HillNC United States Department of Mechanical Engineering San Diego State University San DiegoCA United States Departments of Mathematics Applied Physical Sciences and 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 Chapel HillNC United States McAllister Heart Institute University of North Carolina Chapel HillNC United States

This paper develops and benchmarks an immersed peridynamics method to simulate the deformation, damage, and failure of hyperelastic materials within a fluid-structure interaction framework. The immersed peridynamics method describes an incompressible structure immersed in a viscous incompressible fluid. It expresses the momentum equation and incompressibility constraint in Eulerian form, and it describes the structural motion and resultant forces in Lagrangian form. Coupling between Eulerian and Lagrangian variables is achieved by integral transforms with Dirac delta function kernels, as in standard immersed boundary methods. The major difference between our approach and conventional immersed boundary methods is that we use peridynamics, instead of classical continuum mechanics, to determine the structural forces. We focus on non-ordinary state-based peridynamic material descriptions that allow us to use a constitutive correspondence framework that can leverage well-characterized nonlinear constitutive models of soft materials. The convergence and accuracy of our approach are compared to both conventional and immersed finite element methods using widely used benchmark problems of nonlinear incompressible elasticity. We demonstrate that the immersed peridynamics method yields comparable accuracy with similar numbers of structural degrees of freedom for several choices of the size of the peridynamic horizon. We also demonstrate that the method can generate grid-converged simulations of fluid-driven material damage growth, crack formation and propagation, and rupture under large deformations. © 2022, CC BY-SA.

关键词： Continuum mechanics

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

Hard convex lens-shaped particles: Characterization of dense disordered packings

学校读者我要写书评

暂无评论

Physical Review E 2019年第6期100卷 062902-062902页

作者： Giorgio Cinacchi Salvatore Torquato Departamento de Física Teórica de la Materia Condensada Instituto de Física de la Materia Condensada (IFIMAC) Instituto de Ciencias de Materiales “Nicolás Cabrera” Universidad Autónoma de Madrid Ciudad Universitaria de Cantoblanco E-28049 Madrid Spain Department of Chemistry and Department of Physics Institute for the Science and Technology of Materials Program for Applied and Computational Mathematics Princeton University Princeton New Jersey 08544 USA

Among the family of hard convex lens-shaped particles (lenses), the one with aspect ratio equal to 2/3 is “optimal” in the sense that the maximally random jammed (MRJ) packings of such lenses achieve the highest packing fraction ϕMRJ≃0.73 [G. Cinacchi and S. Torquato, Soft Matter 14, 8205 (2018)]. This value is only a few percent lower than ϕDKP=0.76210⋯, the packing fraction of the corresponding densest-known crystalline (degenerate) packings [G. Cinacchi and S. Torquato, J. Chem. Phys. 143, 224506 (2015)]. By exploiting the appreciably reduced propensity that a system of such optimal lenses has to positionally and orientationally order, disordered packings of them are progressively generated by a Monte Carlo method–based procedure from the dilute equilibrium isotropic fluid phase to the dense nonequilibrium MRJ state. This allows us to closely monitor how the (micro)structure of these packings changes in the process of formation of the MRJ state. The gradual changes undergone by the many structural descriptors calculated here can coherently and consistently be traced back to the gradual increase in contacts between the hard particles until the isostatic mean value of ten contact neighbors per lens is reached at the effectively hyperuniform MRJ state. Compared to the MRJ state of hard spheres, the MRJ state of such optimal lenses is denser (less porous), more disordered, and rattler-free. This set of characteristics makes them good glass formers. It is possible that this conclusion may also hold for other hard convex uniaxial particles with a correspondingly similar aspect ratio, be they oblate or prolate, and that, by using suitable biaxial variants of them, that set of characteristics might further improve.

关键词： Granular packing Jamming Colloids Disordered systems Granular materials Lyotropic liquid crystals Monte Carlo methods

An Immersed Peridynamics Model of Fluid-Structure Interaction Accounting for Material Damage and Failure

学校读者我要写书评

暂无评论

SSRN

SSRN 2022年

作者： Kim, Keon Ho Bhalla, Amneet Pal Singh Griffith, Boyce E. Department of Mathematics University of North Carolina Chapel HillNC United States Department of Mechanical Engineering San Diego State University San DiegoCA United States Departments of Mathematics Applied Physical Sciences and 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 Chapel HillNC United States McAllister Heart Institute University of North Carolina Chapel HillNC United States

This paper develops and benchmarks an immersed peridynamics method to simulate the deformation, damage, and failure of hyperelastic materials within a fluid-structure interaction framework. The immersed peridynamics method describes an incompressible structure immersed in a viscous incompressible fluid. It expresses the momentum equation and incompressibility constraint in Eulerian form, and it describes the structural motion and resultant forces in Lagrangian form. Coupling between Eulerian and Lagrangian variables is achieved by integral transforms with Dirac delta function kernels, as in standard immersed boundary (IB) methods. The major difference between our approach and conventional IB methods is that we use peridynamics, instead of classical continuum mechanics, to determine the structural forces. We focus on non-ordinary state-based peridynamic material descriptions that allow us to use a constitutive correspondence framework that can leverage well characterized nonlinear constitutive models of soft materials. The convergence and accuracy of our approach are compared to both conventional and immersed finite element methods using widely used benchmark problems of nonlinear incompressible elasticity. We demonstrate that the immersed peridynamics method yields comparable accuracy with similar numbers of structural degrees of freedom for several choices of the size of the peridynamic horizon. We also demonstrate that the method can generate grid-converged simulations of fluid-driven material damage growth, crack formation and propagation, and rupture under large deformations. © 2022, The Authors. All rights reserved.

关键词： Fluid structure interaction