On the numerical solution of ill‐conditioned linear systems by regularization and iteration

We propose to reduce the (spectral) condition number of a given linear system by adding a suitable diagonal matrix to the system matrix, in particular by shifting its spectrum. Iterative procedures are then adopted to recover the solution of the original system. The case of real symmetric positive definite matrices is considered in particular, and several numerical examples are given. This approach has some close relations with Riley's method and with Tikhonov regularization. Moreover, we identify approximately the aforementioned procedure with a true action of preconditioning.     

Development of an improved divergence‐free‐condition compensated coupled framework to solve flow problems with time‐varying geometries

In this article a coupled version of the improved divergence‐free‐condition compensated method will be proposed to simulate time‐varying geometries by direct forcing immersed boundary method. The proposed method can be seen as a quasi‐multi‐moment framework due to the fact that the momentum equations are discretized by both cell‐centered and cell‐face velocity. For simulating time‐varying geometries, a semi‐implicit iterative method is proposed for calculating the direct forcing terms. Treatments for suppressing spurious force oscillations, calculating drag/lift forces, and evaluating velocity and pressure for freshly cells will also be addressed. In order to show the applicability and accuracy, analytical as well as benchmark problems will be investigated by the present framework and compared with other numerical and experimental results.     

An efficient edge based data structure for the compressible Reynolds-averaged Navier–Stokes equations on hybrid unstructured meshes

An efficient edge based data structure has been developed in order to implement an unstructured vertex based finite volume algorithm for the Reynolds-averaged Navier–Stokes equations on hybrid meshes. In the present approach, the data structure is tailored to meet the requirements of the vertex based algorithm by considering data access patterns and cache efficiency. The required data are packed and allocated in a way that they are close to each other in the physical memory. Therefore, the proposed data structure increases cache performance and improves computation time. As a result, the explicit flow solver indicates a significant speed up compared to other open-source solvers in terms of CPU time. A fully implicit version has also been implemented based on the PETSc library in order to improve the robustness of the algorithm. The resulting algebraic equations due to the compressible Navier–Stokes and the one equation Spalart–Allmaras turbulence equations are solved in a monolithic manner using the restricted additive Schwarz preconditioner combined with the FGMRES Krylov subspace algorithm. In order to further improve the computational accuracy, the multiscale metric based anisotropic mesh refinement library PyAMG is used for mesh adaptation. The numerical algorithm is validated for the classical benchmark problems such as the transonic turbulent flow around a supercritical RAE2822 airfoil and DLR-F6 wing-body-nacelle-pylon configuration. The efficiency of the data structure is demonstrated by achieving up to an order of magnitude speed up in CPU times.     

On pointwise decay rates of time‐periodic solutions to the Navier–Stokes equation

We study the existence of a time‐periodic solution with pointwise decay properties to the Navier–Stokes equation in the whole space. We show that if the time‐periodic external force is sufficiently small in an appropriate sense, then there exists a time‐periodic solution $\{u,p\}$ of the Navier–Stokes equation such that $|{?}^j{u(t,x)|=O(|x|}^{1?n?j})$ and $|{?}^j{p(t,x)|=O(|x|}^{?n?j})(j=0,1,\dots )$ uniformly in $t\in \mathbb{R}$ as $|x|\to \infty$. Our solution decays faster than the time‐periodic Stokes fundamental solution and the faster decay of its spatial derivatives of higher order is also described.     

Linear independence of a finite set of time-frequency shifts

This paper introduces an open conjecture in time-frequency analysis on the linear independence of a finite set of time-frequency shifts of a given L² function. Firstly, background and motivation for the conjecture are provided. Secondly, the main progress of this linear independence in the past twenty five years is reviewed. Finally, the partial results of the high dimensional case and other cases for the conjecture are briefly presented.     

Assembling hydrodynamic cloaks to conceal complex objects from drag

Transformation hydrodynamics and the corresponding metamaterials have been proposed as a means to exclude the drag force acting on an object. Here, we report a strategy to deploy the hydrodynamic cloaks in a more practical manner by assembling different-shaped cloaking parts. Our strategy is to first model a square-shaped cloak and a carpet cloak and then combine them to conceal a more complex-shaped space in the three-dimensional hydrodynamic flow. With the derivation of transformation hydrodynamics, the coordinate transformations for each hydrodynamic cloaking are demonstrated with the calculated viscosity tensors. The pressure and velocity fields of the square, triangular (carpet), and exemplary three-dimensional house-shaped cloaks are numerically simulated, thus showing a cloaking effect and reduced drag. This study suggests an efficient way of cloaking complex architectures from fluid-dynamic forces.     

13C chemical shift tensors in MOF α-Mg3(HCOO)6: Which component is more sensitive to host-guest interaction?

Metal–organic frameworks (MOFs) are a class of important porous materials with many current and potential applications. Their applications almost always involve the interaction between host framework and guest species. Therefore, understanding of host–guest interaction in MOF systems is fundamentally important. Solid-state NMR spectroscopy is an excellent technique for investigating host–guest interaction as it provides information complementary to that obtained from X-ray diffraction. In this work, using MOF α-Mg₃(HCOO)₆ as an example, we demonstrated that ¹³C chemical shift tensor of organic linker can be utilized to probe the host–guest interaction in MOFs. Obtaining ¹³C chemical shift tensor components (δ₁₁, δ₂₂, and δ₃₃, where δ₁₁ ≥ δ₂₂ ≥ δ₃₃) in this MOF is particularly challenging as there are six coordinatively equivalent but crystallographically non-equivalent carbons in the unit cell with very similar local coordination environment. Two-dimensional magic-angle-turning experiments were employed to measure the ¹³C chemical shift tensors of each individual crystallographically non-equivalent carbon in three microporous α-Mg₃(HCOO)₆ samples with different guest species. The results indicate that the δ₂₂ component (with its direction approximately being co-planar with the formate anion and perpendicular to the C−H bond) is more sensitive to the adsorbate molecules inside the MOF channel due to the weak C−H···O hydrogen bonding or the ring current effect of benzene. The ¹³C isotropic chemical shift, on the other hand, seems much less sensitive to the subtle changes in the local environment around formate linker induced by adsorption. The approach described in this study may be used in future studies on host–guest interaction within MOFs.     

Calculating nuclear magnetic resonance chemical shifts in solvated systems

The nuclear magnetic resonance (NMR) chemical shift is extremely sensitive to molecular geometry, hydrogen bonding, solvent, temperature, pH, and concentration. Calculated magnetic shielding constants, converted to chemical shifts, can be valuable aids in NMR peak assignment and can also give detailed information about molecular geometry and intermolecular effects. Calculating chemical shifts in solution is complicated by the need to include solvent effects and conformational averaging. Here, we review the current state of NMR chemical shift calculations in solution, beginning with an introduction to the theory of calculating magnetic shielding in general, then covering methods for inclusion of solvent effects and conformational averaging, and finally discussing examples of applications using calculated chemical shifts to gain detailed structural information.     

A discontinuous Galerkin method for Stokes equation by divergence-free patch reconstruction

A discontinuous Galerkin method by patch reconstruction is proposed for Stokes flows. A locally divergence-free reconstruction space is employed as the approximation space, and the interior penalty method is adopted which imposes the normal component penalty terms to cancel out the pressure term. Consequently, the Stokes equation can be solved as an elliptic system instead of a saddle-point problem due to such weak form. The number of degree of freedoms of our method is the same as the number of elements in the mesh for different order of accuracy. The error estimations of the proposed method are given in a classical style, which are then verified by some numerical examples.