Regularized numerical solution of the nonlinear,two-dimensional,inverse heat-conduction problem

Tomsk. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, No. 1, pp. 106–112, January–February, 1995.     

A general numerical method for solution of gravity wave problems. Part 2: Steady nonlinear gravity waves

A type of numerical scheme for 2D and 3D steady non-linear water wave problems is described. It is based on the finite process method and is insensitive to initial solutions. The relationship between the finite process method and iterative techniques is discussed. As a numerical example the flow past a submerged vortex is solved and the results are compared with those of other authors.     

A simple and efficient fourth-order approximation solution for nonlinear dynamic system

In the present paper, based on the classical Magnus expansion, a simple and efficient fourth-order integrator is given for an arbitrary nonlinear dynamic system, which can preserve the qualitative properties of the exact solution. The proposed method can be considered an averaging technique, and only requires evaluations of exponentials of simple unidimensional integrals. Finally, the numerical examples are given to demonstrate the validity and effectiveness of the method of this paper.     

A coupling difference scheme for the numerical solution of a singular perturbation problem

In this paper, we consider a singularly perturbed problem without turning points. On a special diseretization mesh, a coupling difference scheme, resulting from central difference scheme and Abrahamsson-Keller-Kreiss box scheme, is proposed and the second order convergence, uniform in the small parameter, is proved. Finally, numerical resulls are provided.     

The Riemann problem for nonlinear degenerate wave equations

This paper studies the Riemann problem for a system of nonlinear degenerate wave equations in elasticity. Since the stress function is neither convex nor concave, the shock condition is degenerate. By introducing a degenerate shock under the generalized shock condition, the global solutions are constructively obtained case by case.     

A numerical study on the solution of the Cauchy problem in elasticity

This work deals with the Cauchy problem in two-dimensional linear elasticity. The equations of the problem are discretized through a standard FEM approach and the resulting ill-conditioned discrete problem is solved within the frame of the Tikhonov approach, the choice of the required regularization parameter is accomplished through the Generalized Cross Validation criterion. On this basis a numerical experimentation has been performed and the calculated solutions have been used to highlight the sensitivity to the amount of known data, the noise always present in the data, the regularity of boundary conditions and the choice of the regularization parameter. The aim of the numerical study is to implicitly device some guidelines to be used in the solution of this kind of problems.     

On the numerical solution of the Fokker-Planck equation for nonlinear stochastic systems

The finite element method is applied to the solution of the transient Fokker-Planck equation for several often cited nonlinear stochastic systems accurately giving, for the first time, the joint probability density function of the response for a given initial distribution. The method accommodates nonlinearity in both stiffness and damping as well as both additive and multiplicative excitation, although only the former is considered herein. In contrast to the usual approach of directly solving the backward Kolmogorov equation, when appropriate boundary conditions are prescribed, the probability density function associated with the first passage problem can be directly obtained. Standard numerical methods are employed, and results are shown to be highly accurate. Several systems are examined, including linear, Duffing, and Van der Pol oscillators.     

Development of a numerical wave tank for analysis of nonlinear and irregular wave field

A numerical wave-absorption filter has been developed for an open boundary condition in the analysis of nonlinear and irregular wave evolution. The filter is composed of a simulated sponge layer and Sommerfeld's radiation condition at the outer edge of the layer. The wave-absorption characteristics of the filter have been investigated by applying the linear potential theory and a two-dimensional nonlinear boundary element model. In both cases, the filter is found to he applicable for a wide range of wave parameters. In order to realize an idealized “numerical wave tank”, the present model also incorporates a nonreflective wave generator in the computational domain composed of a series of vertically aligned point sources. Numerous numerical experiments demonstrate that the present approach is effective in generating an arbitrary wave profile without reflection not only at the open boundaries but also at the wave generator.     

A meshless numerical wave tank for simulation of nonlinear irregular waves in shallow water

Time domain simulation of the interaction between offshore structures and irregular waves in shallow water becomes a focus due to significant increase of liquefied natural gas (LNG) terminals. To obtain the time series of irregular waves in shallow water, a numerical wave tank is developed by using the meshless method for simulation of 2D nonlinear irregular waves propagating from deep water to shallow water. Using the fundamental solution of Laplace equation as the radial basis function (RBF) and locating the source points outside the computational domain, the problem of water wave propagation is solved by collocation of boundary points. In order to improve the computation stability, both the incident wave elevation and velocity potential are applied to the wave generation. A sponge damping layer combined with the Sommerfeld radiation condition is used on the radiation boundary. The present model is applied to simulate the propagation of regular and irregular waves. The numerical results are validated by analytical solutions and experimental data and good agreements are observed. Copyright © 2008 John Wiley & Sons, Ltd.     

A solution to the hydraulic fracture problem in the traveling wave form

Solutions to the system of equations describing the propagation of hydraulic fracture cracks in a porous medium are obtained in the traveling wave form. The only sought solution is the separatrix of integral curves on the “penetration depth-crack width” plane. Some necessary dependencies that should be given at the crack inlet are found for the fluid flow rate and the fluid pressure. The crack width and the fluid penetration depth are related by power laws in the limiting cases when the crack propagation processes or the fluid penetration processes are dominant.     

Numerical solution of the dynamic problem of elasticity for bodies with concentrated masses

On the basis of the virtual-displacement and the DAlembert principles, a new initial-boundary-value problem is formulated for the system of dynamic equilibrium equations of two bodies contacting with friction over a constant area. Their masses are assumed to be distributed with certain densities over the contact surfaces. An equivalent generalized problem is formulated. Classes of discontinuous functions of the finite-element method are used to construct computational schemes of high order of accuracy.__________Translated from Prikladnaya Mekhanika, Vol. 40, No. 12, pp. 65–75, December 2004.     

An h-adaptive solution of the spherical blast wave problem

Shock waves and contact discontinuities usually appear in compressible flows, requiring a fine mesh in order to achieve an acceptable accuracy of the numerical solution. The usage of a mesh adaptation strategy is convenient as uniform refinement of the whole mesh becomes prohibitive in three-dimensional (3D) problems. An unsteady h-adaptive strategy for unstructured finite element meshes is introduced. Non-conformity of the refined mesh and a bounded decrease in the geometrical quality of the elements are some features of the refinement algorithm. A 3D extension of the well-known refinement constraint for 2D meshes is used to enforce a smooth size transition among neighbour elements with different levels of refinement. A density-based gradient indicator is used to track discontinuities. The solution procedure is partially parallelised, i.e. the inviscid flow equations are solved in parallel with a finite element SUPG formulation with shock capturing terms while the adaptation of the mesh is sequentially performed. Results are presented for a spherical blast wave driven by a point-like explosion with an initial pressure jump of 10⁵ atmospheres. The adapted solution is compared to that computed on a fixed mesh. Also, the results provided by the theory of self-similar solutions are considered for the analysis. In this particular problem, adapting the mesh to the solution accounts for approximately 4% of the total simulation time and the refinement algorithm scales almost linearly with the size of the problem.