In this paper, the G2 interpolation by Pythagorean-hodograph (PH) quintic curves in ℝd, d ≥ 2, is considered. The obtained results turn out as a useful tool in practical applications. Independently of the dimension d, they supply a G2 quintic PH spline that locally interpolates two points, two tangent directions and two curvature vectors at these points. The interpolation problem considered is reduced to a system of two polynomial equations involving only tangent lengths of the interpolating curve as unknowns. Although several solutions might exist, the way to obtain the most promising one is suggested based on a thorough asymptotic analysis of the smooth data case. The numerical algorithm traces this solution from a particular set of data to the general case by a homotopy continuation method. Numerical examples confirm the efficiency of the proposed method.

Computational modeling and simulation are presented on the motion of red blood cells behind a moving interface in a capillary. The methodology is based on an immersed boundary method and the skeleton structure of the red blood cell (RBC) membrane is modeled as a spring network. As by the nature of the problem, the computational domain is moving with either a designated RBC or an interface in an infinitely long two-dimensional channel with an undisturbed flow field in front of the computational domain. The tanking-treading and the inclination angle of a cell in a simple shear flow are briefly discussed for the validation purpose. We then present and discuss the results of the motion of red blood cells behind a moving interface in a capillary, which show that the RBCs with higher velocity than the interface speed form a concentrated slug behind the moving interface.

Image fusion is an imaging technique to visualize information from multiple imaging sources in one single image, which is widely used in remote sensing, medical imaging etc. In this work, we study two variational approaches to image fusion which are closely related to the standard TV-L2 and TV-L1 image approximation methods. We investigate their convex optimization formulations, under the perspective of primal and dual, and propose their associated new image decomposition models. In addition, we consider the TV-L1 based image fusion approach and study the specified problem of fusing two discrete-constrained images and where and are the sets of linearly-ordered discrete values. We prove that the TV-L1 based image fusion actually gives rise to the exact convex relaxation to the corresponding nonconvex image fusion constrained by the discrete-valued set This extends the results for the global optimization of the discrete-constrained TV-L1 image approximation [8, 36] to the case of image fusion. As a big numerical advantage of the two proposed dual models, we show both of them directly lead to new fast and reliable algorithms, based on modern convex optimization techniques. Experiments with medical images, remote sensing images and multi-focus images visibly show the qualitative differences between the two studied variational models of image fusion. We also apply the new variational approaches to fusing 3D medical images.

In this paper, we study shifted restated full orthogonalization method with deflation for simultaneously solving a number of shifted systems of linear equations. Theoretical analysis shows that with the deflation technique, the new residual of shifted restarted FOM is still collinear with each other. Hence, the new approach can solve the shifted systems simultaneously based on the same Krylov subspace. Numerical experiments show that the deflation technique can significantly improve the convergence performance of shifted restarted FOM.

Resorting to recent results on subperiodic trigonometric quadrature, we provide three product Gaussian quadrature formulas exact on algebraic polynomials of degree n on circular lunes. The first works on any lune, and has cardinality. The other two have restrictions on the lune angular intervals, but their cardinality is

Experimental observations show that a strong magnetic field has a dramatic influence on the sedimentation of RBCs, which motivates us to model the sedimentation of red blood cell (RBC) under strong external magnetic body force. To model the sedimentation of a RBC in a square duct and a circular pipe, a recently developed technique derived from the lattice Boltzmann and the distributed Lagrange multiplier/fictitious domain methods (LBM-DLM/FD) is extended to employ the mesoscopic network model for simulations of the sedimentation of a RBC in flow. The flow is simulated by the LBM with a strong magnetic body force, while the network model is used for modeling RBC deformation. The fluid-RBC interactions are enforced by the Lagrange multiplier. The sedimentation of RBC in a square duct and a circular pipe is simulated, which demonstrates the developed method's capability to model the sedimentation of RBCs in various flows. Numerical results illustrate that the terminal settling velocity increases incrementally with the exerted body force. The deformation of RBC has a significant effect on the terminal settling velocity due to the change in the frontal area. The larger the exerted force, the smaller the frontal area and the larger the RBC deformation become. Additionally, the wall effect on the motion and deformation of RBC is also investigated.

A computational fluid dynamics solver based on homogeneous cavitation model is employed to compute the two-phase cavitating flow. The model treats the two-phase regime as the homogeneous mixture of liquid and vapour which are locally assumed to be under both kinetic and thermodynamic equilibrium. As our focus is on pressure wave formation, propagation and its impact on cavitation bubble, the compressibility effects of liquid water have to be accounted for and hence the flow is considered to be compressible. The cavitating flow disturbed by the introduced pressure wave is simulated to investigate the unsteady features of cavitation due to the external perturbations. It is observed that the cavity becomes unstable, locally experiencing deformation or collapse, which depends on the shock wave intensity and freestream flow speed.

We consider in this paper numerical approximation of the linear Fluid-Structure Interaction (FSI). We construct a new class of pressure-correction schemes for the linear FSI problem with a fixed interface, and prove rigorously that they are unconditionally stable. These schemes are computationally very efficient, as they lead to, at each time step, a coupled linear elliptic system for the velocity and displacement in the whole region and a discrete Poisson equation in the fluid region.

We have developed a numerical method for simulating viscous flow through a compliant closed tube, driven by a pair of fluid source and sink. As is natural for tubular flow simulations, the problem is formulated in axisymmetric cylindrical coordinates, with fluid flow described by the Navier-Stokes equations. Because the tubular walls are assumed to be elastic, when stretched or compressed they exert forces on the fluid. Since these forces are singularly supported along the boundaries, the fluid velocity and pressure fields become unsmooth. To accurately compute the solution, we use the velocity decomposition approach, according to which pressure and velocity are decomposed into a singular part and a remainder part. The singular part satisfies the Stokes equations with singular boundary forces. Because the Stokes solution is unsmooth, it is computed to second-order accuracy using the immersed interface method, which incorporates known jump discontinuities in the solution and derivatives into the finite difference stencils. The remainder part, which satisfies the Navier-Stokes equations with a continuous body force, is regular. The equations describing the remainder part are discretized in time using the semi-Lagrangian approach, and then solved using a pressure-free projection method. Numerical results indicate that the computed overall solution is second-order accurate in space, and the velocity is second-order accurate in time.