Modeling of the unsteady force for shock–particle interaction

The interaction between a particle and a shock wave leads to unsteady forces that can be an order of magnitude larger than the quasi-steady force in the flow field behind the shock wave. Simple models for the unsteady force have so far not been proposed because of the complicated flow field during the interaction. Here, a simple model is presented based on the work of Parmar et al. (Phil Trans R Soc A 366:2161–2175, 2008). Comparisons with experimental and computational data for both stationary spheres and spheres set in motion by shock waves show good agreement in terms of the magnitude of the peak and the duration of the unsteady force.      

A symbolic algorithm for the automatic computation of multitone-input harmonic balance equations for nonlinear systems

A symbolic algorithm is developed for the automatic generation of harmonic balance equations for multitone input for a class of nonlinear differential systems with polynomial nonlinearities. Generalized expressions are derived for the construction of balance equations for a defined multitone signal form. Procedures are described for determining combinations for a given output frequency from the desired set obtained from box truncated spectra and their permutations to automate symbolic algorithm. An application of method is demonstrated using the well-known Duffing–Van der Pol equation. Then the obtained analytical results are compared with numerical simulations to show the accuracy of the approach. The computation times for both the numerical solutions of equations versus the number of frequency components and the symbolic generation of the equations versus the power of nonlinearity are also investigated.     

On the steady solitary-wave solution of the Green–Naghdi equations of different levels

The steady-state solitary wave solution of high-level Green–Naghdi (GN) equations is obtained by use of the Newton–Raphson method. Four aspects of solitary waves are studied: the wave speed, wave profile, velocity field and particle trajectory. A convergence study is performed for each individual case. Results of the converged model are compared with the existing laboratory experiments and other theoretical solutions for an inviscid and incompressible fluid, including the solutions of the Euler equations. Particle trajectories, predicted by the GN model, show close agreement with the laboratory measurements and provide a new approach to understanding the movement of the particles under a solitary wave. It is further shown that high-level GN equations can predict the solitary wave of the highest height.     

Numerical study of the exact solution of the Navier–Stokes equations describing free-boundary fluid flow

This paper presents a numerical study of the partially invariant solution of the Navier–Stokes equations for the plane case which describes unsteady flow in a layer bounded by a straight solid wall and a free boundary parallel to it. It is found that for different initial flow velocities, a steady state can be established with a decrease or an increase in the initial layer thickness or the layer thickness can be increased infinitely due to fluid inflow from infinity.     

A spectral solver for the Navier–Stokes equations in the velocity–vorticity formulation

A numerical algorithm intended for the study of flows in a cylindrical container under laminar flow conditions is proposed. High resolution of the flow field, governed by the Navier–Stokes equations in velocity–vorticity formulation relative to a cylindrical frame of reference, is achieved through spatial discretisation by means of the spectral method. This method is based on a Fourier expansion in the azimuthal direction and an expansion in Chebyshev polynomials in the (nonperiodic) radial and axial directions. Several regularity constraints are used to take care of the coordinate singularity. These constraints are implemented, together with the boundary conditions at the top, bottom and mantle of the cylinder, via the tau method. The a priori unknown boundary values of the vorticity are evaluated by means of the influence-matrix technique. The compatibility between the mathematical and numerical formulation of the Navier–Stokes equations is established through a tau-correction procedure. The resolved flow field exhibits high-precision satisfaction of the incompressibility constraints for velocity and vorticity and the definition of the vorticity. The performance of the solver is illustrated by resolution of several configurations representative of generic three-dimensional laminar flows.     

A new flux-limiting approach–based kinetic scheme for the Euler equations of gas dynamics

This paper proposes a new kinetic-theory-based high-resolution scheme for the Euler equations of gas dynamics. The scheme uses the well-known connection that the Euler equations are suitable moments of the collisionless Boltzmann equation of kinetic theory. The collisionless Boltzmann equation is discretized using Sweby's flux-limited method and the moment of this Boltzmann level formulation gives a Euler level scheme. It is demonstrated how conventional limiters and an extremum-preserving limiter can be adapted for use in the scheme to achieve a desired effect. A simple total variation diminishing criteria relaxing parameter results in improving the resolution of the discontinuities in a significant way. A 1D scheme is formulated first and an extension to 2D on Cartesian meshes is carried out next. Accuracy analysis suggests that the scheme achieves between first- and second-order accuracy as is expected for any second-order flux-limited method. The simplicity and the explicit form of the conservative numerical fluxes add to the efficiency of the scheme. Several standard 1D and 2D test problems are solved to demonstrate the robustness and accuracy.     

A postprocessing mixed finite element method for the Navier–Stokes equations

A fully discrete postprocessing mixed finite element scheme is considered for solving the time-dependent Navier–Stokes equations. In the PP method, we only consider a non-linear equation in the coarse-level subspace and a linear problem in the fine-level subspace. The analysis shows that the PP scheme can reach the same accuracy as the standard Galerkin method with a very fine mesh size h by an appropriate choice of H. Numerical examples are provided that confirm both the theoretical analysis and the corresponding improvement in computational efficiency.     

A predictor–corrector time integration algorithm for dynamic analysis of nonlinear systems

Nonlinear Dynamics - This paper presents a step-by-step time integration algorithm for efficiently solving second-order nonlinear dynamic problems. The method employs the rewriting of motion as two...     

A parallel gas-kinetic Bhatnagar–Gross–Krook method for the solution of viscous flows on two-dimensional hybrid grids

In this study, a parallel implementation of gas-kinetic Bhatnagar–Gross–Krook method on two-dimensional hybrid grids is presented. Boundary layer regions in wall bounded viscous flows are discretised with quadrilateral grid cells stretched in the direction normal to the solid surface while the rest of the flow domain is discretised by triangular cells. The parallel solution algorithm on hybrid grids is based on the domain decomposition using METIS, a graph partitioning software. The flow solutions obtained in parallel significantly improve the computation time, a significant deficiency of gas-kinetic methods. Several validation test cases presented show the accuracy and robustness of the method developed.     

Functional and generalized separable solutions to unsteady Navier–Stokes equations

The paper studies unsteady Navier–Stokes equations with two space variables. It shows that the non-linear fourth-order equation for the stream function with three independent variables admits functional separable solutions described by a system of three partial differential equations with two independent variables. The system is found to have a number of exact solutions, which generate new classes of exact solutions to the Navier–Stokes equations. All these solutions involve two or more arbitrary functions of a single argument as well as a few free parameters. Many of the solutions are expressed in terms of elementary functions, provided that the arbitrary functions are also elementary; such solutions, having relatively simple form and presenting significant arbitrariness, can be especially useful for solving certain model problems and testing numerical and approximate analytical hydrodynamic methods. The paper uses the obtained results to describe some model unsteady flows of viscous incompressible fluids, including flows through a strip with permeable walls, flows through a strip with extrusion at the boundaries, flows onto a shrinking plane, and others. Some blow-up modes, which correspond to singular solutions, are discussed.     

A semi-analytical solution for the confined compression of hydrated soft tissue

Confined compression is a common experimental technique aimed at gaining information on the properties of biphasic mixtures comprised of a solid saturated by a fluid, a typical example of which are soft hydrated biological tissues. When the material properties (elastic modulus, permeability) are assumed to be homogeneous, the governing equation in the axial displacement reduces to a Fourier equation which can be solved analytically. For the more realistic case of inhomogeneous material properties, the governing equation does not admit, in general, a solution in closed form. In this work, we propose a semi-analytical alternative to Finite Element analysis for the study of the confined compression of linearly elastic biphasic mixtures. The partial differential equation is discretised in the space variable and kept continuous in the time variable, by use of the Finite Difference Method, and the resulting system of ordinary differential equations is solved by means of the Laplace Transform method.     

On the long-time limit of the Garvin–Alterman–Loewenthal solution for a buried blast load

The Garvin–Alterman–Loewenthal solution refers to the problem of a line blast load suddenly applied in the interior of an elastic half-space. It is expected that the long-time asymptotic limit of this solution should be equal to the solution of a related static problem. This expectation is justified here. First, the solution of the static problem is constructed. Then, the asymptotic limit of the transient problem is found, correcting previously published results.     

A time splitting fictitious domain algorithm for fluid–structure interaction problems (A fictitious domain algorithm for FSI)

A Finite Element Method in mixed Eulerian and Lagrangian formulation is developed to allow direct numerical simulations of dynamical interaction between an incompressible fluid and a hyper-elastic incompressible solid. A Fictitious Domain Method is applied so that the fluid is extended inside the deformable solid volume and the velocity field in the entire computational domain is resolved in an Eulerian framework. Solid motion, which is tracked in a Lagrangian framework, is imposed through the body force acting on the fluid within the solid boundaries. Solid stress smoothing on the Lagrangian mesh is performed with the Zienkiewicz–Zhu patch recovery method. High-order Gaussian integration quadratures over cut elements are used in order to avoid sub-meshing within elements in the Eulerian mesh that are intersected by the Lagrangian grid. The algorithm is implemented and verified in two spatial dimensions by comparing with the well validated simulations of solid deformation in a lid driven cavity and periodic elastic wall deformation driven by a time-dependent flow. It shows good agreement with the numerical results reported in the literature. In 3-D the method is validated against previously reported numerical simulations of 3-D rhythmically contracting alveolated ducts.     

Exact solution of the Muskat–Leibenzon problem for a growing elliptic bubble

The exact solution of the two-phase time-dependent Hele–Shaw problem (in other words, the plane Muskat–Leibenzon problem) in which a fluid occupied an unbounded channel is displaced by another fluid incoming through a slitted cut in the channel. In this case the interface between the phases, namely, fluids of different viscosity, evolves as an ellipse whose area and eccentricity vary continuously.     

Finite volume element method for analysis of unsteady reaction–diffusion problems

A finite volume element method is developed for analyzing unsteady scalar reaction-diffusion problems in two dimensions. The method combines the concepts that are employed in the finite volume and the finite element method together. The finite volume method is used to discretize the unsteady reaction-diffusion equation, while the finite element method is applied to estimate the gradient quantities at cell faces. Robustness and efficiency of the combined method have been evaluated on uniform rectangular grids by using available numerical solutions of the two-dimensional reaction-diffusion problems. The numerical solutions demonstrate that the combined method is stable and can provide accurate solution without spurious oscillation along the high-gradient boundary layers.     

A new type of multiple-lump and interaction solution of the Kadomtsev–Petviashvili I equation

Nonlinear Dynamics - In this paper, the solution in the form of Grammian of the Kadomtsev–Petviashvili I equation is employed to investigate a new type of multiple-lump solution. The bound...