Search results

The objective of this paper is to establish a solution strategy for obtaining dual solutions, namely trivial (conventional) and nontrivial (unconventional) solutions, of coupled pore-fluid flow and chemical dissolution problems in heterogeneous porous media.

Through applying a perturbation to the pore-fluid velocity, original governing partial differential equations of a coupled pore-fluid flow and chemical dissolution problem in heterogeneous porous media are transformed into perturbed ones, which are then solved by using the semi-analytical finite element method. Through switching off and on the applied perturbation terms in the resulting perturbed governing partial differential equations, both the trivial and nontrivial solutions can be obtained for the original governing partial differential equations of the coupled pore-fluid flow and chemical dissolution problem in fluid-saturated heterogeneous porous media.

When a coupled pore-fluid flow and chemical dissolution system is in a stable state, the trivial and nontrivial solutions of the system are identical. However, if a coupled pore-fluid flow and chemical dissolution system is in an unstable state, then the trivial and nontrivial solutions of the system are totally different. This recognition can be equally used to judge whether a coupled pore-fluid flow and chemical dissolution system involving heterogeneous porous media is in a stable state or in an unstable state. The proposed solution strategy can produce dual solutions for simulating coupled pore-fluid flow and chemical dissolution problems in fluid-saturated heterogeneous porous media.

A solution strategy is proposed to obtain the nontrivial solution, which is often overlooked in the computational simulation of coupled pore-fluid flow and chemical dissolution problems in fluid-saturated heterogeneous porous media. The proposed solution strategy provides a useful way for understanding the underlying dynamic mechanisms of the chemical damage effect associated with the stability of structures that are built on soil foundations.

The purpose of this paper is to investigate the behavior of a non-Newtonian nanofluid caused by peristaltic waves along an asymmetric channel. Additionally considered is the production of thermal radiation and activation energy.

The equations of momentum, mass and temperature of Sutterby nanofluids are obtained for long wavelength. By taking into account the velocity, temperature and concentration, the formulation is further finished.

Analyses of the physical variables influencing flow features are represented graphically. The present investigation shows that an enhancement in the temperature ratio parameter results in an increase in both the temperature and concentration. The investigation also shows that the dimensionless reaction rate significantly raises the kinetic energy of the reactant, which permits more particle collisions and as a result, raises the temperature field.

Due to their importance in the treatment of cancer, activation energy and thermal radiation as a route of heat transfer are crucial and exciting phenomena for researchers. So, the cancer cells are killed, and tumors are reduced in size with heat and making hyperthermia therapy a cutting-edge cancer treatment.

The purpose of this paper is to apply the novel instantaneous power flow balance (IPFB)-based identification strategy to a specific practical situation like nonlinear lap joints having single and double bolts. The paper also investigates the identification performance of the proposed power flow method over conventional acceleration-matching (AM) methods and other methods in the literature for nonlinear identification.

A parametric model of the joint assembly formulated using generic beam element is used for numerically simulating the experimental response under sinusoidal excitations. The proposed method uses the concept of substructure IPFB criteria, whereby the algebraic sum of power flow components within a substructure is equal to zero, for the formulation of an objective function. The joint parameter identification problem was treated as an inverse formulation by minimizing the objective function using the Particle Swarm Optimization (PSO) algorithm, with the unknown parameters as the optimization variables.

The errors associated with identified numerical results through the instantaneous power flow approach have been compared with the conventional AM method using the same model and are found to be more accurate. The outcome of the proposed method is also compared with other nonlinear time-domain structural identification (SI) methods from the literature to show the acceptability of the results.

In this paper, the concept of IPFB-based identification method was extended to a more specific practical application of nonlinear joints which is not reported in the literature. Identification studies were carried out for both single-bolted and double-bolted lap joints with noise-free and noise-contamination cases. In the current study, only the zone of interest (substructure) needs to be modelled, thus reducing computational complexity, and only interface sensors are required in this method. If the force application point is outside the substructure, there is no need to measure the forcing response also.

The purpose of this paper is to present the freely available finite element simulation software Pyrit.

In a first step, the design principles and the objective of the software project are defined. Then, the software’s structure is established: The software is organized in packages for which an overview is given. The structure is based on the typical steps of a simulation workflow, i.e., problem definition, problem-solving and post-processing. State-of-the-art software engineering principles are applied to ensure a high code quality at all times. Finally, the modeling and simulation workflow of Pyrit is demonstrated by three examples.

Pyrit is a field simulation software based on the finite element method written in Python to solve coupled systems of partial differential equations. It is designed as a modular software that is easily modifiable and extendable. The framework can, therefore, be adapted to various activities, i.e., research, education and industry collaboration.

The focus of Pyrit are static and quasistatic electromagnetic problems as well as (coupled) heat conduction problems. It allows for both time domain and frequency domain simulations.

In research, problem-specific modifications and direct access to the source code of simulation tools are essential. With Pyrit, the authors present a computationally efficient and platform-independent simulation software for various electromagnetic and thermal field problems.

The dynamic response of the nuclear power plants (NPPs) with pile foundation reinforcement have not yet been systemically investigated in detail. Thus, there is an urgent need to improve evaluation methods for nonlithological foundation reinforcements, as this issue is bound to become an unavoidable task.

A nonlinear seismic wave input method is adopted to consider both a nonlinear viscoelastic artificial boundary and the nonlinear properties of the overburden layer soil. Subsequently, the effects of certain vital parameters on the structural response are analyzed.

A suitable range for the size of the overburden foundation is suggested. Then, when piles are used to reinforce the overburden foundation, the peak frequencies in the floor response spectra (FRS) in the horizontal direction becomes higher (38%). Finally, the Poisson ratio of the foundation soil has a significant influence on the FRS peak frequency in the vertical direction (reduce 35%–48%).

The quantifiable results are performed to demonstrate the seismic responses with respect to key design parameters, including foundational dimensions, the Poisson Ratio of the soil and the depth of the foundation. The results can help guide the development of seismic safety requirements for NPPs.

The purpose of this paper is to design a simple and accurate a-posteriori Lagrangian-based error estimator is developed for the class of backward differentiation formula (BDF) algorithms with variable time step size, and the adaptive time-stepping in BDF algorithms is demonstrated for efficient time-dependent simulations in fluid flow and heat transfer.

The Lagrange interpolation polynomial is used to predict the time derivative, and then the accurate primary result is obtained by the Gauss integral, which is applied to evaluate the local error. Not only the generalized formula of the proposed error estimator is presented but also the specific expression for the widely applied BDF1/2/3 is illustrated. Two essential executable MATLAB functions to implement the proposed error estimator are appended for practical applications. Then, the adaptive time-stepping is demonstrated based on the newly proposed error estimator for BDF algorithms.

The validation tests show that the newly proposed error estimator is accurate such that the effectivity index is always close to unity for both linear and nonlinear problems, and it avoids under/overestimation of the exact local error. The applications for fluid dynamics and coupled fluid flow and heat transfer problems depict the advantage of adaptive time-stepping based on the proposed error estimator for time-dependent simulations.

In contrast to existing error estimators for BDF algorithms, the present work is more accurate for the local error estimation, and it can be readily extended to practical applications in engineering with a few changes to existing codes, contributing to efficient time-dependent simulations in fluid flow and heat transfer.

The purpose of this study is to further advance the multiple space/time subdomain framework with model reduction. Existing linear multistep (LMS) methods that are second-order time accurate, and useful for practical applications, have a significant limitation. They do not account for separable controllable numerical dissipation of the primary variables. Furthermore, they have little or no significant choices of altogether different algorithms that can be integrated in a single analysis to mitigate numerical oscillations that may occur. In lieu of such limitations, under the generalized single-step single-solve (GS4) umbrella, several of the deficiencies are circumvented.

The GS4 framework encompasses a wide variety of LMS schemes that are all second-order time accurate and offers controllable numerical dissipation. Unlike existing state-of-art, the present framework permits implicit–implicit and implicit–explicit coupling of algorithms via differential algebraic equations (DAE). As further advancement, this study embeds proper orthogonal decomposition (POD) to further reduce model sizes. This study also uses an iterative convergence check in acquiring sufficient snapshot data to adequately capture the physics to prescribed accuracy requirements. Simple linear/nonlinear transient numerical examples are presented to provide proof of concept.

The present DAE-GS4-POD framework has the flexibility of using different spatial methods and different time integration algorithms in altogether different subdomains in conjunction with the POD to advance and improve the computational efficiency.

The novelty of this paper is the addition of reduced order modeling features, how it applies to the previous DAE-GS4 framework and the improvement of the computational efficiency. The proposed framework/tool kit provides all the needed flexibility, robustness and adaptability for engineering computations.

The numerical methods are of great importance for approximating the solutions of a system of nonlinear singular ordinary differential equations. In this paper, the authors present the biorthogonal flatlet multiwavelet collocation method (BFMCM) as a numerical scheme for a class of system of Lane–Emden equations with initial or boundary or four-point boundary conditions.

The approach is involved in combining the biorthogonal flatlet multiwavelet (BFM) with the collocation method. The authors investigate the properties and procedure of the BFMCM for first time on this class of equations. By using the BFM and the collocation points, the method is constructed and it transforms the nonlinear differential equations problem into a system of nonlinear algebraic equations. The unknown coefficients of the assuming solution are determined by solving the obtained system. Additionally, convergence analysis and numerical stability of the suggested method are provided.

According to the attained results, the proposed BFMCM has more accurate results in comparison with results of other methods. The maximum absolute errors are calculated by using the BFMCM for comparison purposes provided.

The key desirable properties of BFMCM are its efficiency, simple applicability and minimizes errors. Therefore, the proposed method can be used to solve nonlinear problems or problems with singular points.

This study aims to achieve the post-stall pitching maneuver (PSPM) and decrease the deflection frequency of aircraft actuators controlled by the robust backstepping method based on event-triggered mechanism (ETM), nonlinear disturbance observer (NDO) and dynamic surface control (DSC) techniques.

To estimate unsteady aerodynamic disturbances (UADs) to suppress their adverse effects, the NDO is designed. To avoid taking the derivative of the virtual control law directly and eliminate the coupling term of the system states and dynamic surface errors in the stability analysis, an improved DSC is developed. Combined with the NDO and DSC techniques, a robust backstepping method is proposed to achieve the PSPM. Furthermore, to decrease the deflection frequency of the aircraft actuators, a state-dependent ETM is introduced.

An ETM-and-NDO-based backstepping method with an improved DSC technique is developed to achieve the PSPM and decrease the deflection frequency of aircraft actuators. And simulation results are presented to verify the effectiveness of the proposed paper.

Few studies have been conducted on the control of the PSPM in which the lateral and longitudinal attitude dynamics are coupled with each other considering the UADs. Moreover, the mechanism that can decrease the deflection frequency of aircraft actuators is rarely developed in existing research. This study proposes an ETM-and-NDO-based backstepping scheme to address these problems with satisfactory performance of the PSPM.

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.