Search results

The purpose of this paper is to analyze numerically the blowup in finite time of the solutions to a one-dimensional, bidirectional, nonlinear wave model equation for the propagation of small-amplitude waves in shallow water, as a function of the relaxation time, linear and nonlinear drift, power of the nonlinear advection flux, viscosity coefficient, viscous attenuation, and amplitude, smoothness and width of three types of initial conditions.

An implicit, first-order accurate in time, finite difference method valid for semipositive relaxation times has been used to solve the equation in a truncated domain for three different initial conditions, a first-order time derivative initially equal to zero and several constant wave speeds.

The numerical experiments show a very rapid transient from the initial conditions to the formation of a leading propagating wave, whose duration depends strongly on the shape, amplitude and width of the initial data as well as on the coefficients of the bidirectional equation. The blowup times for the triangular conditions have been found to be larger than those for the Gaussian ones, and the latter are larger than those for rectangular conditions, thus indicating that the blowup time decreases as the smoothness of the initial conditions decreases. The blowup time has also been found to decrease as the relaxation time, degree of nonlinearity, linear drift coefficient and amplitude of the initial conditions are increased, and as the width of the initial condition is decreased, but it increases as the viscosity coefficient is increased. No blowup has been observed for relaxation times smaller than one-hundredth, viscosity coefficients larger than ten-thousandths, quadratic and cubic nonlinearities, and initial Gaussian, triangular and rectangular conditions of unity amplitude.

The blowup of a one-dimensional, bidirectional equation that is a model for the propagation of waves in shallow water, longitudinal displacement in homogeneous viscoelastic bars, nerve conduction, nonlinear acoustics and heat transfer in very small devices and/or at very high transfer rates has been determined numerically as a function of the linear and nonlinear drift coefficients, power of the nonlinear drift, viscosity coefficient, viscous attenuation, and amplitude, smoothness and width of the initial conditions for nonzero relaxation times.

The purpose of this paper is to investigate the two-phase flow problems involving gas–liquid mixture.

The governed equations consist of a range of conservation laws modeling a classification of two-phase flow phenomena subjected to a velocity nonequilibrium for the gas–liquid mixture. Effects of the relative velocity are accounted for in the present model by a kinetic constitutive relation coupled to a collection of specific equations governing mass and volume fractions for the gas phase. Unlike many two-phase models, the considered system is fully hyperbolic and fully conservative. The suggested relaxation approach switches a nonlinear hyperbolic system into a semilinear model that includes a source relaxation term and characteristic linear properties. Notably, this model can be solved numerically without the use of Riemann solvers or linear iterations. For accurate time integration, a high-resolution spatial reconstruction and a Runge–Kutta scheme with decreasing total variation are used to discretize the relaxation system.

The method is used in addressing various nonequilibrium two-phase flow problems, accompanied by a comparative study of different reconstructions. The numerical results demonstrate the suggested relaxation method’s high-resolution capabilities, affirming its proficiency in delivering accurate simulations for flow regimes characterized by strong shocks.

While relaxation methods exhibit notable performance and competitive features, as far as we are aware, there has been no endeavor to address nonequilibrium two-phase flow problems using these methods.

Nonlinear mixed-convective entropy optimized the flow of hyperbolic-tangent nanofluid (HTN) with magnetohydrodynamics (MHD) process is considered over a vertical slendering surface. The impression of activation energy is incorporated in the modeling with the significance of nonlinear radiation, dissipative-function, heat generation/consumption connection and Joule heating. Research in this area has practical applications in the design of efficient heat exchangers, thermal management systems or nanomaterial-based devices.

Suitable set of variables is introduced to transform the PDEs (Partial differential equations) system into required ODEs (Ordinary differential equations) system. The transformed ODEs system is then solved numerically via finite difference method. Graphical artworks are made to predict the control of applicable transport parameters on surface entropy, Bejan number, Sherwood number, skin-friction, Nusselt number, temperature, velocity and concentration fields.

It is noticed from present numerical examination that Bejan number aggravates for improved estimations of concentration-difference parameter a_2, Eckert number E_c, thermal ratio parameter ?_w and radiation parameter R_d, whereas surface entropy condenses for flow performance index n, temperature-difference parameter a_1, thermodiffusion parameter N_t and mixed convection parameter ?. Sherwood number is enriched with the amplification of pedesis-motion parameter N_b, while opposite development is perceived for thermodiffusion parameter. Lastly, outcomes are matched with formerly published data to authenticate the present numerical investigation.

To the best of the authors' knowledge, no investigation has been reported yet that explains the entropic behavior with activation energy in the flowing of hyperbolic-tangent mixed-convective nanomaterial due to a vertical slendering surface.

This paper aims to propose a new (3+1)-dimensional integrable Hirota bilinear equation characterized by five linear partial derivatives and three nonlinear partial derivatives.

The authors formally use the simplified Hirota's method and lump schemes for determining multiple soliton solutions and lump solutions, which are rationally localized in all directions in space.

The Painlevé analysis shows that the compatibility condition for integrability does not die away at the highest resonance level, but integrability characteristics is justified through the Lax sense.

Multiple-soliton solutions are explored using the Hirota's bilinear method. The authors also furnish a class of lump solutions using distinct values of the parameters via the positive quadratic function method.

The authors also retrieve a bunch of other solutions of distinct structures such as solitonic, periodic solutions and ratio of trigonometric functions solutions.

This work formally furnishes algorithms for extending integrable equations and for the determination of lump solutions.

To the best of the authors’ knowledge, this paper introduces an original work with newly developed Lax-integrable equation and shows new useful findings.

This article offered two meshfree algorithms, namely the local radial basis functions-finite difference (LRBF-FD) approximation and local radial basis functions-differential quadrature method (LRBF-DQM) to simulate the multidimensional hyperbolic wave models and work is an extension of Jiwari (2015).

In the evolvement of the first algorithm, the time derivative is discretized by the forward FD scheme and the Crank-Nicolson scheme is used for the rest of the terms. After that, the LRBF-FD approximation is used for spatial discretization and quasi-linearization process for linearization of the problem. Finally, the obtained linear system is solved by the LU decomposition method. In the development of the second algorithm, semi-discretization in space is done via LRBF-DQM and then an explicit RK4 is used for fully discretization in time.

For simulation purposes, some 1D and 2D wave models are pondered to instigate the chastity and competence of the developed algorithms.

The developed algorithms are novel for the multidimensional hyperbolic wave models. Also, the stability analysis of the second algorithm is a new work for these types of model.

This paper aims to construct a sixth-order weighted essentially nonoscillatory scheme for simulating the nonlinear degenerate parabolic equations in a finite difference framework.

To design this scheme, we approximate the second derivative in these equations in a different way, which of course is still in a conservative form. In this way, unlike the common practice of reconstruction, the approximation of the derivatives of odd order is needed to develop the numerical flux.

The results obtained by the new scheme produce less error compared to the results of other schemes in the literature that are recently developed for the nonlinear degenerate parabolic equations while requiring less computational times.

This research develops a new weighted essentially nonoscillatory scheme for solving the nonlinear degenerate parabolic equations in multidimensional space. Besides, any selection of the constants (sum equals one is the only requirement for them), named the linear weights, will obtain the desired accuracy.

This study aims to examine the impacts of higher memory dependencies on a novel semiconductor material that exhibits generalized photo-piezo-thermo-elastic properties. Specifically, the research focuses on analyzing the behavior of the semiconductor under three distinct temperature models.

The study assumes a homogeneous and orthotropic piezo-semiconductor medium during photo-thermal excitation. The field equations have been devised to encompass higher order parameters, temporal delays and a specifically tailored kernel function to address the problem. The eigenmode technique is used to solve these equations and derive analytical expressions.

The research presents graphical representations of the physical field distribution across different temperatures, higher order plasma heat conduction models and time. The results reveal that the amplitude of the distribution profile is markedly affected by factors such as the memory effect, time, conductive temperature and spatial coordinates. These factors cannot be overlooked in the analysis and design of the semiconductor.

Specific cases are also discussed in detail, offering the potential to advance the creation of precise models and facilitate future simulations.

The research offers valuable information on the physical field distribution across various temperatures, allowing engineers and designers to optimize the design of semiconductor devices. Understanding the impact of memory effect, time, conductive temperature and spatial coordinates enables device performance and efficiency improvement.

This manuscript is the result of the joint efforts of the authors, who independently initiated and contributed equally to this study.

A higher-order implicit shock-capturing scheme is presented for the Euler equations based on time linearization of the implicit flux vector rather than the residual vector.

The flux vector is linearized through a truncated Taylor-series expansion whose leading-order implicit term is an inner product of the flux Jacobian and the vector of differences between the current and previous time step values of conserved variables. The implicit conserved-variable difference vector is evaluated at cell faces by using the reconstructed states at the left and right sides of a cell face and projecting the difference between the left and right states onto the right eigenvectors. Flux linearization also facilitates the construction of implicit schemes with higher-order spatial accuracy (up to third order in the present study). To enhance the diagonal dominance of the coefficient matrix and thereby increase the implicitness of the scheme, wave strengths at cell faces are expressed as the inner product of the inverse of the right eigenvector matrix and the difference in the right and left reconstructed states at a cell face.

The accuracy of the implicit algorithm at Courant–Friedrichs–Lewy (CFL) numbers greater than unity is demonstrated for a number of test cases comprising one-dimensional (1-D) Sod’s shock tube, quasi 1-D steady flow through a converging-diverging nozzle, and two-dimensional (2-D) supersonic flow over a compression corner and an expansion corner.

The algorithm has the advantage that it does not entail spatial derivatives of flux Jacobian so that the implicit flux can be readily evaluated using Roe’s approximate Jacobian. As a result, this approach readily facilitates the construction of implicit schemes with high-order spatial accuracy such as Roe-MUSCL.

A novel finite-volume-based higher-order implicit shock-capturing scheme was developed that uses time linearization of fluxes at cell interfaces.

The purpose of this study is to introduce the top-level design ideas and the overall architecture of earthquake early-warning system for high speed railways in China, which is based on P-wave earthquake early-warning and multiple ways of rapid treatment.

The paper describes the key technologies that are involved in the development of the system, such as P-wave identification and earthquake early-warning, multi-source seismic information fusion and earthquake emergency treatment technologies. The paper also presents the test results of the system, which show that it has complete functions and its major performance indicators meet the design requirements.

The study demonstrates that the high speed railways earthquake early-warning system serves as an important technical tool for high speed railways to cope with the threat of earthquake to the operation safety. The key technical indicators of the system have excellent performance: The first report time of the P-wave is less than three seconds. From the first arrival of P-wave to the beginning of train braking, the total delay of onboard emergency treatment is 3.63 seconds under 95% probability. The average total delay for power failures triggered by substations is 3.3 seconds.

The paper provides a valuable reference for the research and development of earthquake early-warning system for high speed railways in other countries and regions. It also contributes to the earthquake prevention and disaster reduction efforts.

The purpose of this paper is to study large deviations for the solution processes of a stochastic equation incorporated with the effects of nonlocal condition.

A weak convergence approach is adopted to establish the Laplace principle, which is same as the large deviation principle in a Polish space. The sufficient condition for any family of solutions to satisfy the Laplace principle formulated by Budhiraja and Dupuis is used in this work.

Freidlin–Wentzell type large deviation principle holds good for the solution processes of the stochastic functional integral equation with nonlocal condition.

The asymptotic exponential decay rate of the solution processes of the considered equation towards its deterministic counterpart can be estimated using the established results.