In this paper we study finite element method for three-dimensional incompressible resistive magnetohydrodynamic equations, in which the velocity, the current density, and the magnetic induction are divergence-free. It...
详细信息
The rapid and successful strides in quantum chemistry in the past decades can be largely credited to a conspicuous synergy between theoretical and computational advancements. However, the architectural computer archet...
详细信息
We establish a new H2−Korn’s inequality and its discrete analog, which greatly simplify the construction of nonconforming elements for a linear strain gradient elastic model. The Specht triangle [41] and the NZT tetr...
详细信息
The direct numerical simulation of metal additive manufacturing processes such as laser powder bed fusion is challenging due to the vast differences in spatial and temporal scales. Classical approaches based on locall...
详细信息
The estimation of quantum entropies and distance measures, such as von Neumann entropy, Rényi entropy, Tsallis entropy, trace distance, and fidelity-induced distances such as the Bures distance, has been a key ar...
详细信息
The final black hole left behind after a binary black hole merger can attain a recoil velocity, or a “kick,” reaching values up to 5000 km/s. This phenomenon has important implications for gravitational wave astron...
详细信息
The final black hole left behind after a binary black hole merger can attain a recoil velocity, or a “kick,” reaching values up to 5000 km/s. This phenomenon has important implications for gravitational wave astronomy, black hole formation scenarios, testing general relativity, and galaxy evolution. We consider the gravitational wave signal from the binary black hole merger GW200129_065458 (henceforth referred to as GW200129), which has been shown to exhibit strong evidence of orbital precession. Using numerical relativity surrogate models, we constrain the kick velocity of GW200129 to vf∼1542−1098+747 km/s or vf≳698 km/s (one-sided limit), at 90% credibility. This marks the first identification of a large kick velocity for an individual gravitational wave event. Given the kick velocity of GW200129, we estimate that there is a less than 0.48% (7.7%) probability that the remnant black hole after the merger would be retained by globular (nuclear star) clusters. Finally, we show that kick effects are not expected to cause biases in ringdown tests of general relativity for this event, although this may change in the future with improved detectors.
Atomistic materials modeling is a critical task with wide-ranging applications, from drug discovery to materials science, where accurate predictions of the target material property can lead to significant advancements...
详细信息
Two-scale simulations are often employed to analyze the effect of the microstructure on a component's macroscopic properties. Understanding these structure-property relations is essential in the optimal design of ...
详细信息
Two-scale simulations are often employed to analyze the effect of the microstructure on a component's macroscopic properties. Understanding these structure-property relations is essential in the optimal design of materials for specific applications. However, these two-scale simulations are typically computationally expensive and infeasible in multi-query contexts such as optimization and material design. To make such analyses amenable, the microscopic simulations can be replaced by inexpensive-to-evaluate surrogate models. Such surrogate models must be able to handle microstructure parameters in order to be used for material design. A previous work focused on the construction of an accurate surrogate model for microstructures under varying loading and material parameters by combining proper orthogonal decomposition and Gaussian process regression. However, that method works only for a fixed geometry, greatly limiting the design space. This work hence focuses on extending the methodology to treat geometrical parameters. To this end, a method that transforms different geometries onto a parent domain is presented, that then permits existing methodologies to be applied. We propose to solve an auxiliary problem based on linear elasticity to obtain the geometrical transformations. The method has a good reducibility and can therefore be quickly solved for many different geometries. Using these transformations, combined with the nonlinear microscopic problem, we derive a fast-to-evaluate surrogate model with the following key features: (1) the predictions of the effective quantities are independent of the auxiliary problem, (2) the predicted stress fields automatically fulfill the microscopic balance laws and are periodic, (3) the method is non-intrusive, (4) the stress field for all geometries can be recovered, and (5) the sensitivities are available and can be readily used for optimization and material design. The proposed methodology is tested on several composite micros
We optimize a numerical time-stabilization routine for a class of phase-field crystal (PFC) models of solidification. By numerical experiments, we demonstrate that our simple approach can improve the accuracy of under...
We optimize a numerical time-stabilization routine for a class of phase-field crystal (PFC) models of solidification. By numerical experiments, we demonstrate that our simple approach can improve the accuracy of underlying time integration schemes by a few orders of magnitude. We investigate different time integration schemes. Moreover, as a prototypical example for applications, we extend our numerical approach to a PFC model of solidification with an explicit temperature coupling.
In this paper, we present a novel penalty model called ExPen for optimization over the Stiefel manifold. Different from existing penalty functions for orthogonality constraints, ExPen adopts a smooth penalty function ...
详细信息
暂无评论