Considering the coupled nonlinear Schr?dinger system with multiply components, we provide a novel framework for constructing energy-preserving algorithms. In detail, based on the high order compact finite difference m...
详细信息
Considering the coupled nonlinear Schr?dinger system with multiply components, we provide a novel framework for constructing energy-preserving algorithms. In detail, based on the high order compact finite difference method, Fourier pseudospectral method and wavelet collocation method for spatial discretizations, a series of high accurate conservative algorithms are presented. The proposed algorithms can preserve the corresponding discrete charge and energy conservation laws exactly, which would guarantee their numerical stabilities during long time computations. Furthermore, several analogous multi-symplectic algorithms are constructed as comparison. Numerical experiments for the unstable plane waves will show the advantages of the proposed algorithms over long time and verify the theoretical analysis.
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 tim...
详细信息
Large collections of coupled, heterogeneous agents can manifest complex dynamical behavior presenting difficulties for simulation and analysis. However, if the collective dynamics lie on a low-dimensional manifold the...
详细信息
Large collections of coupled, heterogeneous agents can manifest complex dynamical behavior presenting difficulties for simulation and analysis. However, if the collective dynamics lie on a low-dimensional manifold then the original agent-based model may be approximated with a simplified surrogate model on and near the low-dimensional space where the dynamics live. Analytically identifying such simplified models can be challenging or impossible, but here we present a data-driven coarse-graining methodology for discovering such reduced models. We consider two types of reduced models: globally-based models which use global information and predict dynamics using information from the whole ensemble, and locally-based models that use local information, that is, information from just a subset of agents close (close in heterogeneity space, not physical space) to an agent, to predict the dynamics of an agent. For both approaches we are able to learn laws governing the behavior of the reduced system on the low-dimensional manifold directly from time series of states from the agent-based system. These laws take the form of either a system of ordinary differential equations (ODEs), for the globally-based approach, or a partial differential equation (PDE) in the locally-based case. For each technique we employ a specialized artificial neural network integrator that has been templated on an Euler time stepper (i.e. a ResNet) to learn the laws of the reduced model. As part of our methodology, we utilize the proper orthogonal decomposition (POD) to identify the low-dimensional space of the dynamics. Our globally-based technique uses the resulting POD basis to define a set of coordinates for the agent states in this space, and then seeks to learn the time evolution of these coordinates as a system of ODEs. For the locally-based technique, we propose a methodology for learning a partial differential equation representation of the agents;the PDE law depends on the state variables and
The immersed boundary (IB) method is a non-body conforming approach to fluid-structure interaction (FSI) that uses an Eulerian description of the momentum, viscosity, and incompressibility of a coupled fluid-structure...
详细信息
The immersed boundary (IB) method is a non-body conforming approach to fluid-structure interaction (FSI) that uses an Eulerian description of the momentum, viscosity, and incompressibility of a coupled fluid-structure system and a Lagrangian description of the deformations, stresses, and resultant forces of the immersed structure. Integral transforms with Dirac delta function kernels couple the Eulerian and Lagrangian variables, and in practice, discretizations of these integral transforms use regularized delta function kernels. Many different kernel functions have been proposed, but prior numerical work investigating the impact of the choice of kernel function on the accuracy of the methodology has often been limited to simplified test cases or Stokes flow conditions that may not reflect the method’s performance in applications, particularly at intermediate-to-high Reynolds numbers, or under different loading conditions. This work systematically studies the effect of the choice of regularized delta function in several fluid-structure interaction benchmark tests using the immersed finite element/difference (IFED) method, which is an extension of the IB method that uses a finite element structural discretizations combined with a Cartesian grid finite difference method for the incompressible Navier-Stokes equations. Whereas the conventional IB method spreads forces from the nodes of the structural mesh and interpolates velocities to those nodes, the IFED formulation evaluates the regularized delta function on a collection of interaction points that can be chosen to be denser than the nodes of the Lagrangian mesh. This opens the possibility of using structural discretizations with wide node spacings that would produce gaps in the Eulerian force in nodally coupled schemes (e.g., if the node spacing is comparable to or broader than the support of the regularized delta function). Earlier work with this methodology suggested that such coarse structural meshes can yield imp
This paper analyzes the robust long-term growth rate of expected utility and expected return from holding a leveraged exchange-traded fund (LETF). When the Markovian model parameters in the reference asset are uncerta...
详细信息
We propose a new reionization probe that uses cosmic microwave background (CMB) observations; the cross-correlation between fluctuations in the CMB optical depth which probes the integrated electron density, δτ, and...
详细信息
We propose a new reionization probe that uses cosmic microwave background (CMB) observations; the cross-correlation between fluctuations in the CMB optical depth which probes the integrated electron density, δτ, and the Compton y-map which probes the integrated electron pressure. This cross-correlation is much less contaminated than the y-map power spectrum by late-time cluster contributions. In addition, this cross-correlation can constrain the temperature of ionized bubbles while the optical-depth fluctuations and kinetic Sunyaev Zel’dvich effect can not. We measure this new observable using a Planck y-map as well as a map of optical-depth fluctuations that we reconstruct from Planck CMB temperature data. We use our measurements to derive a first CMB only upper limit on the temperature inside ionized bubbles, Tb≲7.0×105 K (2σ). We also present future forecasts, assuming a fiducial model with characteristic reionization bubble size Rb=5 Mpc and Tb=5×104 K. The signal-to-noise ratio of the fiducial cross-correlation using a signal dominated y-map from the Probe of Inflation and Cosmic Origins (PICO) becomes ≃7 with CMB-S4 δτ and ≃13 with CMB-HD δτ. For the fiducial model, we predict that the CMB-HD—PICO cross-correlation should achieve an accurate measurement of the reionization parameters; Tb≃49800−5100+4500 K and Rb≃5.09−0.79+0.66 Mpc. Since the power spectrum of the electron density fluctuations is constrained by the δτ autospectrum, the temperature constraints should be only weakly model dependent on the details of the electron distributions and should be statistically representative of the temperature in ionized bubbles during reionization. This cross-correlation could, therefore, become an important observable for future CMB experiments.
Node embeddings aim to associate a vector to every vertex of a graph which can then be used for downstream tasks such as clustering, classification, or link prediction. Many popular node embeddings such as no $d$ e2ve...
Node embeddings aim to associate a vector to every vertex of a graph which can then be used for downstream tasks such as clustering, classification, or link prediction. Many popular node embeddings such as no $d$ e2vec and DeepWalk are based upon counting which nodes frequently co-occur in random walks of the graph. In this paper, we show that the performance of such algorithms can be improved by rewiring the edges of the graph through a variety of network indices before running DeepWalk. These rewirings effectively give the random walker an inductive bias and increase the accuracy of a logistic regression classifier applied to the node embedding on several benchmark data sets.
The hierarchical small-world network is a real-world network. It models well the benefit transmission web of the pyramid selling in China and many other countries. In this paper, by applying the spectral graph theory,...
详细信息
暂无评论