HomPINNs: homotopy physics-informed neural networks for solving the inverse problems of nonlinear differential equations with multiple solutions.

Due to the complex behavior arising from non-uniqueness, symmetry, and bifurcations in the solution space, solving inverse problems of nonlinear differential equations (DEs) with multiple solutions is a challenging task. To address this, we propose homotopy physics-informed neural networks (HomPINNs), a novel framework that leverages homotopy continuation and neural networks (NNs) to solve inverse problems. The proposed framework begins with the use of NNs to simultaneously approximate unlabeled observations across diverse solutions while adhering to DE constraints. Through homotopy continuation, the proposed method solves the inverse problem by tracing the observations and identifying multiple solutions. The experiments involve testing the performance of the proposed method on one-dimensional DEs and applying it to solve a two-dimensional Gray-Scott simulation. Our findings demonstrate that the proposed method is scalable and adaptable, providing an effective solution for solving DEs with multiple solutions and unknown parameters. Moreover, it has significant potential for various applications in scientific computing, such as modeling complex systems and solving inverse problems in physics, chemistry, biology, etc.

Eigenvibrations of Kirchhoff Rectangular Random Plates on Time-Fractional Viscoelastic Supports via the Stochastic Finite Element Method.

The present work's main objective is to investigate the natural vibrations of the thin (Kirchhoff-Love) plate resting on time-fractional viscoelastic supports in terms of the Stochastic Finite Element Method (SFEM). The behavior of the supports is described by the fractional order derivatives of the Riemann-Liouville type. The subspace iteration method, in conjunction with the continuation method, is used as a tool to solve the non-linear eigenproblem. A deterministic core for solving structural eigenvibrations is the Finite Element Method. The probabilistic analysis includes the Monte-Carlo simulation and the semi-analytical approach, as well as the iterative generalized stochastic perturbation method. Probabilistic structural response in the form of up to the second-order characteristics is investigated numerically in addition to the input uncertainty level. Finally, the probabilistic relative entropy and the safety measure are estimated. The presented investigations can be applied to the dynamics of foundation plates resting on viscoelastic soil.

Selected Problems of Random Free Vibrations of Rectangular Thin Plates with Viscoelastic Dampers.

The main motivation of this work was to present a semi-analytical extension of the correspondence principle in stochastic dynamics. It is demonstrated for the stochastic structural free vibrations of Kirchhoff-Love elastic, isotropic and rectangular plates supported by viscoelastic generalized Maxwell dampers. The ambient temperature of the plate affects the dampers only and is included in a mathematical model using the frequency-temperature correspondence principle. The free vibration problem of the plate-viscoelastic damper system is solved using the continuation method and also the Finite Element Method (FEM). The stochastic approach begins with an initial deterministic sensitivity analysis to detect the most influential parameters and numerical FEM recovery of the polynomial representation for lower eigenfrequencies versus these parameters. A final symbolic integration leads to the first four basic probabilistic characteristics, all delivered as functions of the input uncertainties.

Hyperelastic Microcantilever AFM: Efficient Detection Mechanism Based on Principal Parametric Resonance.

The impetus of writing this paper is to propose an efficient detection mechanism to scan the surface profile of a micro-sample using cantilever-based atomic force microscopy (AFM), operating in non-contact mode. In order to implement this scheme, the principal parametric resonance characteristics of the resonator are employed, benefiting from the bifurcation-based sensing mechanism. It is assumed that the microcantilever is made from a hyperelastic material, providing large deformation under small excitation amplitude. A nonlinear strain energy function is proposed to capture the elastic energy stored in the flexible component of the device. The tip-sample interaction is modeled based on the van der Waals non-contact force. The nonlinear equation governing the AFM's dynamics is established using the extended Hamilton's principle, obeying the Euler-Bernoulli beam theory. As a result, the vibration behavior of the system is introduced by a nonlinear equation having a time-dependent boundary condition. To capture the steady-state numerical response of the system, a developed Galerkin method is utilized to discretize the partial differential equation to a set of nonlinear ordinary differential equations (ODE) that are solved by the combination of shooting and arc-length continuation method. The output reveals that while the resonator is set to be operating near twice the fundamental natural frequency, the response amplitude undergoes a significant drop to the trivial stable branch as the sample's profile experiences depression in the order of the picometer. According to the performed sensitivity analysis, the proposed working principle based on principal parametric resonance is recommended to design AFMs with ultra-high detection resolution for surface profile scanning.

Numerical Analysis of an Unsteady, Electroviscous, Ternary Hybrid Nanofluid Flow with Chemical Reaction and Activation Energy across Parallel Plates.

Despite the recycling challenges in ionic fluids, they have a significant advantage over traditional solvents. Ionic liquids make it easier to separate the end product and recycle old catalysts, particularly when the reaction media is a two-phase system. In the current analysis, the properties of transient, electroviscous, ternary hybrid nanofluid flow through squeezing parallel infinite plates is reported. The ternary hybrid nanofluid is synthesized by dissolving the titanium dioxide (TiO2), aluminum oxide (Al2O3), and silicon dioxide (SiO2) nanoparticles in the carrier fluid glycol/water. The purpose of the current study is to maximize the energy and mass transfer rate for industrial and engineering applications. The phenomena of fluid flow is studied, with the additional effects of the magnetic field, heat absorption/generation, chemical reaction, and activation energy. The ternary hybrid nanofluid flow is modeled in the form of a system of partial differential equations, which are subsequently simplified to a set of ordinary differential equations through resemblance substitution. The obtained nonlinear set of dimensionless ordinary differential equations is further solved, via the parametric continuation method. For validity purposes, the outcomes are statistically compared to an existing study. The results are physically illustrated through figures and tables. It is noticed that the mass transfer rate accelerates with the rising values of Lewis number, activation energy, and chemical reaction. The velocity and energy transfer rate boost the addition of ternary NPs to the base fluid.

Computational Valuation of Darcy Ternary-Hybrid Nanofluid Flow across an Extending Cylinder with Induction Effects.

The flow of an electroconductive incompressible ternary hybrid nanofluid with heat conduction in a boundary layer including metallic nanoparticles (NPs) over an extended cylindrical with magnetic induction effects is reported in this research. The ternary hybrid nanofluid has been synthesized with the dispersion of titanium dioxide, cobalt ferrite, and magnesium oxide NPs in the base fluid water. For a range of economical and biological applications, a computational model is designed to augment the mass and energy conveyance rate and promote the performance and efficiency of thermal energy propagation. The model has been written as a system of partial differential equations. Which are simplified to the system of ordinary differential equations through similarity replacements. The computing approach parametric continuation method is used to further process the resultant first order differential equations. The results are validated with the bvp4c package for accuracy and validity. The outcomes are displayed and analyzed through Figures and Tables. It has been observed that the inverse Prandtl magnetic number and a larger magnetic constant reduce the fluid flow and elevate the energy profile. The variation of ternary hybrid NPs significantly boosts the thermophysical features of the base fluid.

Numerical Simulation of a Time-Dependent Electroviscous and Hybrid Nanofluid with Darcy-Forchheimer Effect between Squeezing Plates.

In this article, the behavior of transient electroviscous fluid flow is investigated through squeezing plates containing hybrid nanoparticles. A hybrid nanofluid MoS2+Au/C2H6O2-H2O was formulated by dissolving the components of an inorganic substance such as molybdenum disulfide (MoS2) and gold (Au) in a base fluid of ethylene glycol/water. This hybrid non-liquid flow was modeled by various nonlinear mathematical fluid flow models and subsequently solved by numerical as well as analytical methods. For the numerical solution of nonlinear ODEs, a built-in function BVP4C was used in MATLAB, and the same problem was solved in MATHEMATICA by HAM. The result of the present problem related to the results obtained from the existing literature under certain conditions. The outcomes revealed that the concentration profiles were more sensitive to homogeneity diversity parameters. The simulation of the various physical parameters of the model indicated that the heat transfer through a mixture of hybrid nanofluids was greater than a simple nanofluid. In addition, the phenomenon of mixed convection was considered to improve the velocity of simple nanofluids and hybrid nanofluids, when both cases have low permeability. A rise in the volume fraction of the nanomaterials, Φ, was associated with an increase in the heat transfer rate. It was observed that the heat transfer rate of the hybrid nanofluids MoS2+Au/C2H6O2-H2O was higher than that of the single nanofluids MoS2/C2H6O2-H2O.

Steady Squeezing Flow of Magnetohydrodynamics Hybrid Nanofluid Flow Comprising Carbon Nanotube-Ferrous Oxide/Water with Suction/Injection Effect.

The main purpose of the current article is to scrutinize the flow of hybrid nanoliquid (ferrous oxide water and carbon nanotubes) (CNTs + Fe3O4/H2O) in two parallel plates under variable magnetic fields with wall suction/injection. The flow is assumed to be laminar and steady. Under a changeable magnetic field, the flow of a hybrid nanofluid containing nanoparticles Fe3O4 and carbon nanotubes are investigated for mass and heat transmission enhancements. The governing equations of the proposed hybrid nanoliquid model are formulated through highly nonlinear partial differential equations (PDEs) including momentum equation, energy equation, and the magnetic field equation. The proposed model was further reduced to nonlinear ordinary differential equations (ODEs) through similarity transformation. A rigorous numerical scheme in MATLAB known as the parametric continuation method (PCM) has been used for the solution of the reduced form of the proposed method. The numerical outcomes obtained from the solution of the model such as velocity profile, temperature profile, and variable magnetic field are displayed quantitatively by various graphs and tables. In addition, the impact of various emerging parameters of the hybrid nanofluid flow is analyzed regarding flow properties such as variable magnetic field, velocity profile, temperature profile, and nanomaterials volume fraction. The influence of skin friction and Nusselt number are also observed for the flow properties. These types of hybrid nanofluids (CNTs + Fe3O4/H2O) are frequently used in various medical applications. For the validity of the numerical scheme, the proposed model has been solved by another numerical scheme (BVP4C) in MATLAB.

Optimal Placement of Viscoelastic Vibration Dampers for Kirchhoff Plates Based on PSO Method.

The main subject of this study is to determine the optimal position of a fixed number of viscoelastic dampers on the surface of a thin (Kirchhoff-Love) plate. It is assumed that the dampers are described according to the generalized Maxwell model. In order to determine the optimal position of the dampers, a metaheuristic optimization method is used, called the particle swarm optimization method. The non-dimensional damping ratio of the first mode of the plate vibrations is assumed as an objective function in the task. The dynamic characteristics of the plate with dampers are determined by solving the non-linear eigenproblem using the continuation method. The finite element method is used to determine the stiffness matrix and the mass matrix occurring in the considered eigenproblem. The results of exemplary numerical calculations are also presented, where the final optimal arrangement of dampers on the surface of sample plates with different boundary conditions is shown graphically.

A continuation method for spatially discretized models with nonlocal interactions conserving size and shape of cells and lattices.

In this paper, we introduce a continuation method for the spatially discretized models, while conserving the size and shape of the cells and lattices. This proposed method is realized using the shift operators and nonlocal operators of convolution types. Through this method and using the shift operator, the nonlinear spatially discretized model on the uniform and nonuniform lattices can be systematically converted into a spatially continuous model; this renders both models point-wisely equivalent. Moreover, by the convolution with suitable kernels, we mollify the shift operator and approximate the spatially discretized models using the nonlocal evolution equations, rendering suitable for the application in both experimental and mathematical analyses. We also demonstrate that this approximation is supported by the singular limit analysis, and that the information of the lattice and cells is expressed in the shift and nonlocal operators. The continuous models designed using our method can successfully replicate the patterns corresponding to those of the original spatially discretized models obtained from the numerical simulations. Furthermore, from the observations of the isotropy of the Delta-Notch signaling system in a developing real fly brain, we propose a radially symmetric kernel for averaging the cell shape using our continuation method. We also apply our method for cell division and proliferation to spatially discretized models of the differentiation wave and describe the discrete models on the sphere surface. Finally, we demonstrate an application of our method in the linear stability analysis of the planar cell polarity model.

Nonlinear Behavior of Electrostatically Actuated Microbeams with Coupled Longitudinal⁻Transversal Vibration.

Microelectromechanical switch has become an essential component in a wide variety of applications, ranging from biomechanics and aerospace engineering to consumer electronics. Electrostatically actuated microbeams and microplates are chief parts of many MEMS instruments. In this study, the nonlinear characteristics of coupled longitudinal-transversal vibration are analyzed, while an electrostatically actuated microbeam is designed considering that the frequency ratio is two to one between the first longitudinal vibration and transversal vibration. The nonlinear governing equations are truncated into a set of coupled ordinary differential equations by the Galerkin method. Then the equations are solved using the multiple-scales method and the nonlinear dynamics of the internal resonance is investigated. The influence of bias voltage, longitudinal excitation and frequency detuning parameters are mainly analyzed. Results show that using the pseudo-arclength continuation method, the nonlinear amplitude-response curves can be plotted continuously. The saturation and jump phenomena are greatly affected by the bias voltage and the detuning frequency. Beyond the critical excitation amplitude, the response energy will transfer from the longitudinal motion to the transversal motion, even the excitation is employed on the longitudinal direction. The large-amplitude jump of the low-order vibration mode can be used to detect the variation of the conditions or parameters, which shows great potential in improving precision of MEMS switches.

Speed-up hyperspheres homotopic path tracking algorithm for PWL circuits simulations.

In the present work, we introduce an improved version of the hyperspheres path tracking method adapted for piecewise linear (PWL) circuits. This enhanced version takes advantage of the PWL characteristics from the homotopic curve, achieving faster path tracking and improving the performance of the homotopy continuation method (HCM). Faster computing time allows the study of complex circuits with higher complexity; the proposed method also decrease, significantly, the probability of having a diverging problem when using the Newton-Raphson method because it is applied just twice per linear region on the homotopic path. Equilibrium equations of the studied circuits are obtained applying the modified nodal analysis; this method allows to propose an algorithm for nonlinear circuit analysis. Besides, a starting point criteria is proposed to obtain better performance of the HCM and a technique for avoiding the reversion phenomenon is also proposed. To prove the efficiency of the path tracking method, several cases study with bipolar (BJT) and CMOS transistors are provided. Simulation results show that the proposed approach can be up to twelve times faster than the original path tracking method and also helps to avoid several reversion cases that appears when original hyperspheres path tracking scheme was employed.