Search results

Free and moving boundary problems require the simultaneous solution of unknown field variables and the boundaries of the domains on which these variables are defined. There are many technologically important processes that lead to moving boundary problems associated with fluid surfaces and solid‐fluid boundaries. These include crystal growth, metal alloy and glass solidification, melting and flame propagation. The directional solidification of semi‐conductor crystals by the Bridgman—Stockbarger method^1,2 is a typical example of such a complex process. A numerical model of this growth method must solve the appropriate heat, mass and momentum transfer equations and determine the location of the melt—solid interface. In this work, a Chebyshev pseudospectral collocation method is adapted to the problem of directional solidification. Implementation involves a solution algorithm that combines domain decomposition, a finite‐difference preconditioned conjugate minimum residual method and a Picard type iterative scheme.

This paper aims to reconstruct the spatially varying orthotropic conductivity based on a two-dimensional inverse heat conduction problem described by a partial differential equation (PDE) model with mixed boundary conditions. The proposed discretization uses a highly accurate technique and allows simple implementations. Also, the authors solve the related inverse problem in such a way that smoothness is enforced on the iterations, showing promising results in synthetic examples and real problems with moving heat source.

The discretization procedure applied to the model for the direct problem uses a pseudospectral collocation strategy in the spatial variables and Crank–Nicolson method for the time-dependent variable. Then, the related inverse problem of recovering the conductivity from temperature measurements is solved by a modified version of Levenberg–Marquardt method (LMM) which uses singular scaling matrices. Problems where data availability is limited are also considered, motivated by a face milling operation problem. Numerical examples are presented to indicate the accuracy and efficiency of the proposed method.

The paper presents a discretization for the PDEs model aiming on simple implementations and numerical performance. The modified version of LMM introduced using singular scaling matrices shows the capabilities on recovering quantities with precision at a low number of iterations. Numerical results showed good fit between exact and approximate solutions for synthetic noisy data and quite acceptable inverse solutions when experimental data are inverted.

The paper is significant because of the pseudospectral approach, known for its high precision and easy implementation, and usage of singular regularization matrices on LMM iterations, unlike classic implementations of the method, impacting positively on the reconstruction process.

This paper aims to present thermal and mechanical stresses in solid and hollow thick-walled cylinders and spheres made of functionally graded materials (FGMs) under the effect of heat generation.

Constant internal temperature and convective external conditions in hollow bodies along with internal heat generation with a combination of outer convective conditions in solid bodies are investigated individually. The effect of the heat convection coefficient on solid bodies is additionally discussed. The variation of the FGM properties in the radial direction is adapted to the Mori–Tanaka homogenization schemes, which produces irregular and two-point linear boundary value problems that are numerically solved by the pseudospectral Chebyshev method.

It has been shown that the selection of the mixtures of FGMs has to be made correctly to keep the thermal and mechanical loads acting on objects at low levels.

In this study, both solid and hollow functionally graded cylinders and spheres for different boundary conditions that are as their engineering applications are examined with the proposed method. The results have demonstrated that the pseudospectral Chebyshev method has high accuracy, low calculation costs and ease of application and can be easily adapted to such engineering problems.

The purpose of this paper is to both determine the effects of the nonlinearity on the wave dynamics and assess the temporal and spatial accuracy of five finite difference methods for the solution of the inviscid generalized regularized long-wave (GRLW) equation subject to initial Gaussian conditions.

Two implicit second- and fourth-order accurate finite difference methods and three Runge-Kutta procedures are introduced. The methods employ a new dependent variable which contains the wave amplitude and its second-order spatial derivative. Numerical experiments are reported for several temporal and spatial step sizes in order to assess their accuracy and the preservation of the first two invariants of the inviscid GRLW equation as functions of the spatial and temporal orders of accuracy, and thus determine the conditions under which grid-independent results are obtained.

It has been found that the steepening of the wave increase as the nonlinearity exponent is increased and that the accuracy of the fourth-order Runge-Kutta method is comparable to that of a second-order implicit procedure for time steps smaller than 100th, and that only the fourth-order compact method is almost grid-independent if the time step is on the order of 1,000th and more than 5,000 grid points are used, because of the initial steepening of the initial profile, wave breakup and solitary wave propagation.

This is the first study where an accuracy assessment of wave breakup of the inviscid GRLW equation subject to initial Gaussian conditions is reported.

The purpose of this paper is to apply rectangular Bézier surface patches directly into the mathematical formula used to solve boundary value problems modeled by Laplace's equation. The mathematical formula, called the parametric integral equation systems (PIES), will be obtained through the analytical modification of the conventional boundary integral equations (BIE), with the boundary mathematically described by rectangular Bézier patches.

The paper presents the methodology of the analytic connection of the rectangular patches with BIE. This methodology is a generalization of the one previously used for 2D problems.

In PIES the paper separates the necessity of performing simultaneous approximation of both boundary shape and the boundary functions, as the boundary geometry has been included in its mathematical formalism. The separation of the boundary geometry from the boundary functions enables to achieve an independent and more effective improvement of the accuracy of both approximations. Boundary functions are approximated by the Chebyshev series, whereas the boundary is approximated by Bézier patches.

The originality of the proposed approach lies in its ability to automatic adapt the PIES formula for modified shape of the boundary modeled by the Bézier patches. This modification does not require any dividing the patch into elements and creates the possibility for effective declaration of boundary geometry in continuous way directly in PIES.

The present work is devoted to the numerical study of the stability of shallow jet. The effects of important parameters on the stability behavior for large scale shallow jets are considered and investigated. Connections between the stability theory and observed features reported in the literature are emphasized. The paper aims to discuss these issues.

A linear stability analysis of shallow jet incorporating the effects of bottom topography, bed friction and viscosity has been carried out by using the shallow water stability equation derived from the depth averaged shallow water equations in conjunction with both Chézy and Manning resistance formulae. Effects of the following main factors on the stability of shallow water jets are examined: Rossby number, bottom friction number, Reynolds number, topographic parameters, base velocity profile and resistance model. Special attention has been paid to the Coriolis effects on the jet stability by limiting the rotation number in the range of Ro∈[0, 1.0].

It is found that the Rossby number may either amplify or attenuate the growth of the flow instability depending on the values of the topographic parameters. There is a regime where the near cancellation of Coriolis effects due to other relevant parameters influences is responsible for enhancement of stability. The instability can be suppressed by the bottom friction when the bottom friction number is large enough. The amplification rate may become sensitive to the relatively small Reynolds number. The stability region using the Manning formula is larger than that using the Chézy formula. The combination of these effects may stabilize or destabilize the shallow jet flow. These results of the stability analysis are compared with those from the literature.

Results of linear stability analysis on shallow jets along roughness bottom bed are presented. Different from the previous studies, this paper includes the effects of bottom topography, Rossby number, Reynolds number, resistance formula and bed friction. It is found that the influence of Reynolds number on the stability of the jet is notable for relative small value. Therefore, it is important to experimental investigators that the viscosity should be considered with comparison to the results from inviscid assumption. In contrast with the classical analysis, the use of multi-parameters of the base velocity and topographic profile gives an extension to the jet stability analysis. To characterize the large scale motion, besides the bottom friction as proposed in the related literature, the Reynolds number Re, Rossby number Ro, the topographic parameters and parameters controlling base velocity profile may also be important to the stability analysis of shallow jet flows.

The mechanism of instability in fluid mechanics is one of the topics, which has not been fully understood yet, since there are numerous external and internal agents for a flow field to lose its laminar behavior and to find itself in transition to turbulence. In this study, only the effects of pressure gradients are considered. Commonly used methods by many researchers, in this field, try to connect experimental evidence to numerically obtained results. One of the ways to do this matching is linearized (N‐S) equations under small disturbances which finally yield the so‐called Orr‐Sommerfeld equation. In this paper, Orr‐Sommerfeld equation is solved to eight‐digit accuracy and compared with literature, by using Chebyshev polynomials. After having this accuracy, the flow field has been examined for various pressure gradients. Increasing the pressure gradients either positive or negative, increases the instability region.

The purpose of this paper is to present some numerical methods based on different time stepping and space discretization methods for the Allen‐Cahn equation with non‐periodic boundary conditions.

In space the equation is discretized by the Chebyshev spectral method, while in time the exponential time differencing fourth‐order Runge‐Kutta (ETDRK4) and implicit‐explicit scheme are used. Also, for comparison the finite difference scheme in both space and time is used.

It is found that the use of implicit‐explicit scheme allows use of a large time‐step, since an explicit method has less order of accuracy as compared to implicit‐explicit method. In time‐stepping the proposed ETDRK4 does not behave well for this special kind of partial differential equation.

The paper presents some numerical methods based on different time stepping and space discretization methods for the Allen‐Cahn equation with non‐periodic boundary conditions.

The research is directed at development of an efficient and accurate technique for modelling incompressible free surface flows in which viscous effects may not be neglected. The paper describes the methodology, and gives illustrative results for simple geometries.

The pseudospectral matrix element method of discretisation is selected as the basis for the CFD technique adopted, because of its high spectral accuracy. It is implemented as a means of solving the Navier‐Stokes equations coupled with the modified compressibility method.

The viscous solver has been validated for the benchmark cases of uniform flow past a cylinder at a Reynolds number of 40, and 2D cavity flows. Results for sloshing of a viscous fluid in a tank have been successfully compared with those from a linearised analytical solution. Application of the method is illustrated by the results for the interaction of an impulsive wave with a surface piercing circular cylinder in a cylindrical tank.

The paper demonstrates the viability of the approach adopted. The limitation of small amplitude waves should be tackled in future work.

The results will have particular significance in the context of validating computations from more complex schemes applicable to arbitrary geometries.

The new methodology and results are of interest to the community of those developing numerical models of flow past marine structures.

The purpose of this paper is to study a finite‐thrust orbital transfer optimization problem via a new optimal control method‐Gauss pseudospectral method.

Based on the dynamic equations with the pseudo‐equinoctial elements as state variable, optimality condition is derived and the optimization problem is converted into nonlinear programming problem. Gauss pseudospectral method is used to avoid the two‐point boundary value problem. The dynamic equations are converted into static parameter optimization problem. The state variables and control variables are selected as optimal parameters at all collocation nodes. Two numerical examples of orbital transfer with coplanar and different planes are analyzed, respectively.

The simulation results demonstrate that Gauss pseudospectral method is not sensitive to the initial conditions of orbital transfer. They also show good robustness and control facility.

The precision and efficiency of this trajectory optimization method are demonstrated by applying it to space vehicle orbital transfer with finite thrust optimization problem.