On the usage of NURBS as interface representation in free‐surface flows

When simulating free‐surface flows using the finite element method, there are many cases where the governing equations require information which must be derived from the available discretized geometry. Examples are curvature or normal vectors. The accurate computation of this information directly from the finite element mesh often requires a high degree of refinement—which is not necessarily required to obtain an accurate flow solution. As a remedy and an option to be able to use coarser meshes, the representation of the free surface using non‐uniform rational B‐splines (NURBS) curves or surfaces is investigated in this work. The advantages of a NURBS parameterization in comparison with the standard approach are discussed. In addition, it is explored how the pressure jump resulting from surface tension effects can be handled using doubled interface nodes. Numerical examples include the computation of surface tension in a two‐phase flow as well as the computation of normal vectors as a basis for mesh deformation methods. For these examples, the improvement of the numerical solution compared with the standard approaches on identical meshes is shown. Copyright © 2011 John Wiley & Sons, Ltd.     

The use of Sinc‐collocation method for solving Falkner–Skan boundary‐layer equation

The MHD Falkner–Skan equation arises in the study of laminar boundary layers exhibiting similarity on the semi‐infinite domain. The proposed approach is equipped by the orthogonal Sinc functions that have perfect properties. This method solves the problem on the semi‐infinite domain without truncating it to a finite domain and transforming domain of the problem to a finite domain. In addition, the governing partial differential equations are transformed into a system of ordinary differential equations using similarity variables, and then they are solved numerically by the Sinc‐collocation method. It is shown that the Sinc‐collocation method converges to the solution at an exponential rate. Copyright © 2010 John Wiley & Sons, Ltd.     

影视节目中酸碱指示剂的应用与评价

将影视节目中酸碱指示剂的应用引入到化学教学中，让学生在娱乐中学习化学知识，创造一个有趣的学习情境，提高学生的学习兴趣。     

An improved exponential transformation for nearly singular boundary element integrals in elasticity problems

This paper presents an improved exponential transformation for nearly singular boundary element integrals in elasticity problems. The new transformation is less sensitive to the position of the projection point compared with the original transformation. In our work, the conventional distance function is modified into a new form in the polar coordinate system. Based on the refined distance function, an improved exponential transformation is proposed in the polar coordinate system. Moreover, to perform integrations on irregular elements, an adaptive integration scheme considering both the element shape and the projection point associated with the improved transformation is proposed. Furthermore, when the projection point is located outside the integration element, another nearest point is introduced to subdivide the integration elements into triangular or quadrilateral patches of fine shapes. Numerical examples are presented to verify the proposed method. Results demonstrate the accuracy and efficiency of our method.     

A fixed‐grid b‐spline finite element technique for fluid–structure interaction

We present a fixed‐grid finite element technique for fluid–structure interaction problems involving incompressible viscous flows and thin structures. The flow equations are discretised with isoparametric b‐spline basis functions defined on a logically Cartesian grid. In addition, the previously proposed subdivision‐stabilisation technique is used to ensure inf–sup stability. The beam equations are discretised with b‐splines and the shell equations with subdivision basis functions, both leading to a rotation‐free formulation. The interface conditions between the fluid and the structure are enforced with the Nitsche technique. The resulting coupled system of equations is solved with a Dirichlet–Robin partitioning scheme, and the fluid equations are solved with a pressure–correction method. Auxiliary techniques employed for improving numerical robustness include the level‐set based implicit representation of the structure interface on the fluid grid, a cut‐cell integration algorithm based on marching tetrahedra and the conservative data transfer between the fluid and structure discretisations. A number of verification and validation examples, primarily motivated by animal locomotion in air or water, demonstrate the robustness and efficiency of our approach. Copyright © 2013 John Wiley & Sons, Ltd.     

99Tc NMR as a promising technique for structural investigation of biomolecules: theoretical studies on the solvent and thermal effects of phenylbenzothiazole complex

The phenylbenzothiazole compounds show antitumor properties and are highly selective. In this paper, the ⁹⁹Tc chemical shifts based on the (^99mTc)(CO)₃(NNO) complex conjugated to the antitumor agent 2‐(4′‐aminophenyl)benzothiazole are reported. Thermal and solvent effects were studied computationally by quantum‐chemical methods, using the density functional theory (DFT) (DFT level BPW91/aug‐cc‐pVTZ for the Tc and BPW91/IGLO‐II for the other atoms) to compute the NMR parameters for the complex. We have calculated the ⁹⁹Tc NMR chemical shifts of the complex in gas phase and solution using different solvation models (polarizable continuum model and explicit solvation). To evaluate the thermal effect, molecular dynamics simulations were carried, using the atom‐centered density matrix propagation method at the DFT level (BP86/LanL2dz). The results highlight that the ⁹⁹Tc NMR spectroscopy can be a promising technique for structural investigation of biomolecules, at the molecular level, in different environments. Copyright © 2014 John Wiley & Sons, Ltd.     

An immersed boundary method wall model for high‐Reynolds‐number channel flow over complex topography

High‐Reynolds‐number channel flows regularly encounter topographies composed of multiple length scales and that protrude into the boundary layer. Physically, the presence of immersed obstacles leads to increased velocity gradients, turbulence production, and manifestation of wakes. Considerable challenges are associated with numerically describing the presence of obstacles in channel flows. Common approaches include generation of a computational mesh that is uniquely designed for the flow and obstacle, the immersed boundary method, and terrain‐following coordinates. There are challenges and limitations associated with each of these techniques. Specification of boundary conditions representing the perimeter of solid obstacles is a primary challenge of the immersed boundary method. In this document, a simplistic canopy stress‐like wall model is used to impose boundary conditions. The model isolates aerodynamically relevant local frontal areas through evaluation of the gradient of the topographic height field. The gradient of the height field describes both the surface‐normal direction and the frontal area, making it ideal for detecting areas on which the flow impinges. The model is tested in numerical simulations of turbulent half‐channel flow over topographies with different obstacles affixed–right prisms, rectangular prisms, ellipsoidal mounds, and sinusoids. In all cases, the performance is strong relative to datasets presented in the literature. Results are finally presented for numerical simulation of flow over complex synthetic fractal‐like topography and a synthetic city. These results show interesting trends in how the turbulent multiscale flow field responds to multiscale topography. Copyright © 2012 John Wiley & Sons, Ltd.     

A moving particle semi‐implicit method for free surface flow: Improvement in inter‐particle force stabilization and consistency restoring

The moving particle semi‐implicit (MPS) method has been widely applied in free surface flows. However, the implementation of MPS remains limited because of compressive instability occurred when the particles are under compressive stress states. This study proposed an inter‐particle force stabilization and consistency restoring MPS (IFS‐CR‐MPS) method to overcome this numerical instability. For inter‐particle force stabilization, a hyperbolic‐shaped quintic kernel function is developed with a non‐negative and smooth second order derivative to satisfy the stability criterion under compressive stress state. Then, a contrastive study is conducted on the contradiction between the common understanding of the conventional MPS hyperbolic‐shaped kernel function and its performance. The result shows that the conventional MPS hyperbolic‐shaped kernel function can easily cause violent repulsive inter‐particle force and then lead to the compressive instability. Therefore, the first order derivative of the modified hyperbolic‐shaped quintic kernel function is recommended as the form of the contribution of the neighbor particles to achieve a more stable inter‐particle repulsive force. For consistency restoring, the Taylor series expansion and the hyperbolic‐shaped quintic kernel are combined to improve the accuracy of the viscosity and pressure calculation. The IFS‐CR‐MPS algorithm is subsequently verified by the inviscid hydrostatic pressure, jet impacting, and viscous droplet impacting problems. These results can be used for choosing kernel function and the contribution of neighbor particles in particle methods. Copyright © 2016 John Wiley & Sons, Ltd.     

Dark solitons of the Biswas–Milovic equation by the first integral method

This paper studies the Biswas–Milovic equation that is a generalized version of the familiar nonlinear Schrodinger's equation describing the propagation of solitons through optical fibers for trans-continental and trans-oceanic distances with Kerr law nonlinearity by the aid of the first integral method. The dark 1-soliton solution is retrieved by the aid of this method and a couple of other singular periodic solutions are also obtained.     

An anti‐diffusive volume of fluid method for interfacial fluid flows

A volume of fluid (VOF) method is developed combining a first‐order limited downwind scheme with higher order accurate schemes. The method is characterized by retaining a sharp fluid interface and a reduction in numerical diffusion near the interface, but avoids complicated geometrical reconstruction as occurs in most volume tracing algorithms. To demonstrate the accuracy and robustness of the method, a selection of numerical experiments are presented involving a pure advection problem, a water wave impact caused by a dam breaking and liquid sloshing in a partially filled tank. Copyright © 2011 John Wiley & Sons, Ltd.