3D Equations Research Articles

Review of Some Modified Generalized Korteweg–de Vries–Kuramoto–Sivashinsky Equations (Part II)

In part I of this work to appear in Foudations-MDPI 2024, some existence and uniqueness results for the solutions of some equations were reviewed, such as the Korteweg–de Vries equation (KdV), the Kuramoto–Sivashinsky equation (KS), the generalized Korteweg–de Vries–Kuramoto–Sivashinsky equation (gKdV-KS), and the nonhomogeneous boundary value problem for the KdV-KS equation in quarter plane. The main objective of this paper is to review some results of the existence of global attractors for the evolution equations with nonlinearity of the form N(ux), where ux denotes the derivative of u with respect to x, focusing in particular on the Kuramoto–Sivashinsky equation in one and two dimensions. In order to illustrate the general abstract results, we have chosen to discuss in detail the existence of global attractors for the Kuramoto–Sivashinsky (KS) equation in 1D and 2D. Once a global attractor is obtained, the question arises whether it has special regularity properties. Then we give an integrated version of the homogeneous steady state Kuramoto–Sivashinsky equation in Rn. This work ends with a change from rectangular to polar coordinates in the three-dimensional KS equation to give an energy estimate in this case.

Stationary scattering for the nonlinear Schrödinger equation with point-like obstacles: exact solutions

AbstractWe solve the Nonlinear Schrödinger Equation (NLSE) in 1D in presence of one, two and several Dirac delta potentials. With the help of an equivalent central force problem we obtain the analytical solutions in terms of a biparametric family containing the Jacobi functions. Elliptic Jacobi functions are already reported in the literature but they have not been used in the context of a scattering problem under causal boundary conditions. In the simplest examples of one or two Dirac deltas we analyze how the nonlinear term of the equation affects the modulus and phase profiles of the wave function. We also study the transmission curves under the nonlinear modification of the tunneling behavior for the first time. For a Fabri-Perot configuration made of two deltas, we obtain the effect of nonlinear coupling in the positions of the local maxima (resonances). We lay the foundations for nonlinear Anderson localization of 1D BECs in a speckle field. Upon redefinition of parameters these novel results describe the dynamics of a stationary Higgs field in 1D. Finally, we discuss the conditions for soliton formation under the influence of a Dirac comb potential, giving rise to fully correlated locations and intensities of the defects.

Invariant Measures for a Class of Stochastic Third-Grade Fluid Equations in 2D and 3D Bounded Domains

Invariant Measures for a Class of Stochastic Third-Grade Fluid Equations in 2D and 3D Bounded Domains

Calculating traveltimes in 2D general tilted transversely isotropic media using fast sweeping method

Traveltime calculations play an important role in the field of exploration seismology, such as traveltime tomography and seismic imaging and so on. Seismic anisotropy poses a challenge for traveltime calculation, because anisotropic eikonal solvers are more complex than the isotropic counter part. To solve the eikonal equations in 2D tilted transversely isotropic (TTI) media, we have developed a fast algorithm combine with fast sweeping method to compute the first arrival traveltimes of quasi-P (qP)-, quasi-SV (qSV)-, and quasi-SH(qSH)-waves. For the qP- and qSV-waves, we analyzed the quartic coupled slowness surface equation derived from the Christoffel equation. Then, we constructed a local solver to relate traveltime and slowness. We found that in the local solver, one component of the slowness vector is known and the corresponding slowness equation is monotonic. This provides a strong basis for the fast iterative algorithm we proposed, where we use the Newton method to solve the qP- and qSV-wave slowness equation to determine the related traveltimes. For the qSH wave, the slowness equation is quadratic and simple to solve. Numerical experiments demonstrate that the proposed method can obtain accurate traveltimes for simple and complicated 2D TTI models.

Sinc-Galerkin method and a higher-order method for a 1D and 2D time-fractional diffusion equations

In this article, a new numerical algorithm for solving a 1-dimensional (1D) and 2-dimensional (2D) time-fractional diffusion equation is proposed. The Sinc-Galerkin scheme is considered for spatial discretization, and a higher-order finite difference formula is considered for temporal discretization. The convergence behavior of the methods is analyzed, and the error bounds are provided. The main objective of this paper is to propose the error bounds for 2D problems by using the Sinc-Galerkin method. The proposed method in terms of convergence is studied by using the characteristics of the Sinc function in detail with optimal rates of exponential convergence for 2D problems. Some numerical experiments validate the theoretical results and present the efficiency of the proposed schemes.

Modulation instability and rogue waves for two and three dimensional nonlinear Klein-Gordon equation.

We perform the modulation instability analysis of the 2D and 3D nonlinear Klein-Gordon equation. The instability region depends on dispersion and wavenumbers of the plane wave. The N-breathers of the nonlinear Klein-Gordon equation are constructed directly from its 2N-solitons obtained in history. The regularity conditions of breathers are established. The dynamic behaviors of breathers of the 2D nonlinear Klein-Gordon equation are consistent with modulation instability analysis. Furthermore, by means of the bilinear method together with improved long-wave limit technique, we obtain general high order rogue waves of the 2D and 3D nonlinear Klein-Gordon equation. In particular, the first- and second-order rogue waves and lumps of the 2D nonlinear Klein-Gordon equation are investigated by using their explicit expressions. We find that their dynamic behaviors are similar to the nonlinear Schrödinger equation. Finally, the first-order rational solutions are illustrated for the 3D nonlinear Klein-Gordon equation. It is demonstrated that the rogue waves of the 2D and 3D nonlinear Klein-Gordon equation always exist by choosing dispersion and wavenumber of plane waves.

The pressure-wired Stokes element: a mesh-robust version of the Scott–Vogelius element

The Scott–Vogelius finite element pair for the numerical discretization of the stationary Stokes equation in 2D is a popular element which is based on a continuous velocity approximation of polynomial order k and a discontinuous pressure approximation of order k-1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$k-1$$\\end{document}. It employs a “singular distance” (measured by some geometric mesh quantity Θz≥0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\Theta \\left( \ extbf{z}\\right) \\ge 0$$\\end{document} for triangle vertices z\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ extbf{z}$$\\end{document}) and imposes a local side condition on the pressure space associated to vertices z\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ extbf{z}$$\\end{document} with Θz=0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Theta \\left( \ extbf{z}\\right) =0$$\\end{document}. The method is inf-sup stable for any fixed regular triangulation and k≥4\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$k\\ge 4$$\\end{document}. However, the inf-sup constant deteriorates if the triangulation contains nearly singular vertices 0<Θz≪1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$0<\\Theta \\left( \ extbf{z}\\right) \\ll 1$$\\end{document}. In this paper, we introduce a very simple parameter-dependent modification of the Scott–Vogelius element with a mesh-robust inf-sup constant. To this end, we provide sharp two-sided bounds for the inf-sup constant with an optimal dependence on the “singular distance”. We characterise the critical pressures to guarantee that the effect on the divergence-free condition for the discrete velocity is negligibly small, for which we provide numerical evidence.

Quantum Physics-Informed Neural Networks.

In this study, the PennyLane quantum device simulator was used to investigate quantum and hybrid, quantum/classical physics-informed neural networks (PINNs) for solutions to both transient and steady-state, 1D and 2D partial differential equations. The comparative expressibility of the purely quantum, hybrid and classical neural networks is discussed, and hybrid configurations are explored. The results show that (1) for some applications, quantum PINNs can obtain comparable accuracy with less neural network parameters than classical PINNs, and (2) adding quantum nodes in classical PINNs can increase model accuracy with less total network parameters for noiseless models.

Improving the conditioning of the Method of Fundamental Solutions for the Helmholtz equation on domains in polar or elliptic coordinates

A new approach to overcome the ill-conditioning of the Method of Fundamental Solutions (MFS) combining Singular Value Decomposition (SVD) and an adequate change of basis was introduced in [1] as MFS-SVD. The original formulation considered polar coordinates and harmonic polynomials as basis functions and is restricted to the Laplace equation in 2D. In this work, we start by adapting the approach to the Helmholtz equation in 2D and later extending it to elliptic coordinates. As in the Laplace case, the approach in polar coordinates has very good numerical results both in terms of conditioning and accuracy for domains close to a disk but does not perform so well for other domains, such as an eccentric ellipse. We therefore consider the MFS-SVD approach in elliptic coordinates with Mathieu functions as basis functions for the latter. We illustrate the feasibility of the approach by numerical examples in both cases.

Element-based finite volume method applied to autogenous bead-on-plate GTAW process using the enthalpy energy equation

The welding process is the most used technique for metal joining. Understanding the temperature variation along the welded part during the process can prevent the appearance of failures. The experimental investigation process is quite time-consuming and costly. Therefore, numerical simulation processes based on the finite difference method, finite element method, finite volume method, or meshless Element-Free Galerkin (EFG) methods are important tools to optimize the welding process. The main goal of the present study is to show the feasibility of the Element-based Finite-Volume Method (EbFVM) approach for actual engineering applications. To solve the unsteady 2D and 3D thermal energy equation using enthalpy as an independent variable, an in-house Fortran code has been developed based on the EbFVM approach in conjunction with unstructured and structured meshes. The numerical simulations, with four types of different heat sources, were performed for applications of real welding processes with variations in density and enthalpy as a function of temperature. The results are presented in terms of thermal cycles and temperature fields. Furthermore, the developed code was confronted against experimental works from the literature, simulated and lab-controlled experiments with AISI 409 ferritic workpieces, and exact analytical solutions with thermocouples fixed in different positions. In general, the numerical results from the current investigation are in close agreement with the results from the literature and the experimental results performed by the authors. The numerical results also highlighted the differences between the 2D and 3D models for thermal cycles near the bead weld.

Optimal decay estimate and asymptotic profile for solutions to the generalized Zakharov–Kuznetsov–Burgers equation in 2D

We consider the Cauchy problem for the generalized Zakharov–Kuznetsov–Burgers equation in 2D. This is one of the nonlinear dispersive–dissipative equations, which has a spatial anisotropic dissipative term −μuxx. In this paper, we prove that the solution to this problem decays at the rate of t−34 in the L∞-sense, provided that the initial data u0(x,y) satisfies u0∈L1(R2) and some appropriate regularity assumptions. Moreover, we investigate the more detailed large time behavior and obtain a lower bound of the L∞-norm of the solution. As a result, we prove that the given decay rate t−34 of the solution to be optimal. Furthermore, combining the techniques used for the parabolic equations and for the Schrödinger equation, we derive the explicit asymptotic profile for the solution.

Supercloseness of linear finite element method on Vulanović–Bakhvalov mesh for singularly perturbed convection–diffusion equation in 1D

Supercloseness of linear finite element method is studied on the Vulanović–Bakhvalov mesh for a singularly perturbed convection–diffusion problem. A novel interpolation is proposed to tackle challenges arising from the subinterval adjacent to the transition point. This innovative interpolation method is formulated by leveraging the distinctive features of the boundary layer and incorporating the specific attributes of the Vulanović–Bakhvalov mesh. On the basis of that, supercloseness of 2th order can be obtained.

Energy-conserving neural network for turbulence closure modeling

In turbulence modeling, we are concerned with finding closure models that represent the effect of the subgrid scales on the resolved scales. Recent approaches gravitate towards machine learning techniques to construct such models. However, the stability of machine-learned closure models and their abidance by physical structure (e.g. symmetries, conservation laws) are still open problems. To tackle both issues, we take the ‘discretize first, filter next’ approach. In this approach we apply a spatial averaging filter to existing fine-grid discretizations. The main novelty is that we introduce an additional set of equations which dynamically model the energy of the subgrid scales. Having an estimate of the energy of the subgrid scales, we can use the concept of energy conservation to derive stability. The subgrid energy containing variables is determined via a data-driven technique. The closure model is used to model the interaction between the filtered quantities and the subgrid energy. Therefore the total energy should be conserved. Abiding by this conservation law yields guaranteed stability of the system. In this work, we propose a novel skew-symmetric convolutional neural network architecture that satisfies this law. The result is that stability is guaranteed, independent of the weights and biases of the network. Importantly, as our framework allows for energy exchange between resolved and subgrid scales it can model backscatter. To model dissipative systems (e.g. viscous flows), the framework is extended with a diffusive component. The introduced neural network architecture is constructed such that it also satisfies momentum conservation. We apply the new methodology to both the viscous Burgers' equation and the Korteweg-De Vries equation in 1D. The novel architecture displays superior stability properties when compared to a vanilla convolutional neural network.

Electrostatic conductive disc singularity resolved

The conventional electrostatic solutions for two-dimensional (2D) electrodes possess edge singularities for the surface charge density σ and the normal component of the electric field En. These singularities are generally non-physical because they admit infinite gradients of the concentration of free charge carriers. In particular, they are unacceptable in the studies of the local field sensitive effects, such as the electric breakdown and the ferroelectric domain reversal. We claim that account for diffusion of free charge carriers leads to the disappearance of the edge singularities. This generalization occurs consistently within the same basic concept of conduction. Specifically, we consider the case of U-biased circular disc electrode of radius a. Account for diffusion leads here to a strongly nonlinear integral 2D equation for the electrostatic potential φ(r). Numerical solution of this equation shows that the law σ(a)∝U2/a takes place. Outside a close vicinity of the disc edge, we stay close to the conventional electrostatic solution for φ and σ.

Numerical study of the Serre–Green–Naghdi equations in 2D *

A detailed numerical study of solutions to the Serre–Green–Naghdi (SGN) equations in 2D with vanishing curl of the velocity field is presented. The transverse stability of line solitary waves, 1D solitary waves being exact solutions of the 2D equations independent of the second variable, is established numerically. The study of localized initial data as well as crossing 1D solitary waves does not give an indication of existence of stable structures in SGN solutions localized in two spatial dimensions. For the numerical experiments, an approach based on a Fourier spectral method with a Krylov subspace technique is applied.

Scattering problems for the wave equation in 1D: D’Alembert-type representations and a reconstruction method

We derive the extension of the classical d’Alembert formula for the wave equation, which provides the analytical solution for the direct scattering problem for a medium with constant refractive index. Analogous formulae exist already in the literature, but in the current work this is derived in a natural way for general incident field, by employing results obtained via the Fokas method. This methodology is further extended to a medium with piecewise constant refractive index, providing the apparatus for the solution of the associated inverse scattering problem. Hence, we provide an exact reconstruction method which is also valid for phaseless data.

Estimation of the concentration boundary layer adjacent to a flat surface using computational fluid dynamics

Dissolution-permeation (D/P) experiments are widely used during preclinical development due to producing results with better predictability than traditional monophasic experiments. However, it is difficult to compare absorption across in vitro setups given the propensity to only report apparent permeability. We therefore developed an approach to predict the concentration boundary layer for any D/P device by using computational fluid dynamics (CFD). The Navier-Stokes and continuity equation in 2D were solved numerically in MATLAB and by finite element methods in COMSOL v6.1 to predict the momentum (ηf′) and concentration ηg boundary layer for a flow over a flat plate, i.e. the classical Blasius boundary layer flow. A MATLAB algorithm was developed to calculate the edge of either boundary layer. The methodology to determine the concentration boundary layer based on Blasius’s analysis provided an accurate estimate for both ηf′ and ηg, resulting in, ηf′/ηg, at high Schmidt numbers (Sc ∼ 1000) within 14 % of the Blasius solution and 6.6 % of the accepted Schmidt number correlation (Sc1/3=ηf′/ηg). The methodology based on the Blasius analysis of the concentration boundary layer using velocity and concentration profiles computed using CFD presented herein will enable characterization/analysis of complex D/P apparatuses used in preclinical development, where an analytical solution may not be available.

Airborne Transmission of SARS-CoV-2: The Contrast between Indoors and Outdoors

COVID-19 is an airborne disease, with the vast majority of infections occurring indoors. In comparison, little transmission occurs outdoors. Here, we investigate the airborne transmission pathways that differentiate the indoors from outdoors and conclude that profound differences exist, which help to explain why SARS-CoV-2 transmission is much more prevalent indoors. Near- and far-field transmission pathways are discussed along with factors that affect infection risk, with aerosol concentration, air entrainment, thermal plumes, and occupancy duration all identified as being influential. In particular, we present the fundamental equations that underpin the Wells–Riley model and show the mathematical relationship between inhaled virus particles and quanta of infection. A simple model is also presented for assessing infection risk in spaces with incomplete air mixing. Transmission risk is assessed in terms of aerosol concentration using simple 1D equations, followed by a description of thermal plume–ceiling interactions. With respect to this, we present new experimental results using Schlieren visualisation and computational fluid dynamics (CFD) based on the Eulerian–Lagrangian approach. Pathways of airborne infection are discussed, with the key differences identified between indoors and outdoors. In particular, the contribution of thermal and exhalation plumes is evaluated, and the presence of a near-field/far-field feedback loop is postulated, which is absent outdoors.

On the Interplay Between Vortices and Harmonic Flows: Hodge Decomposition of Euler’s Equations in 2d

On the Interplay Between Vortices and Harmonic Flows: Hodge Decomposition of Euler’s Equations in 2d

A mimetic interpolation-free cell-centered finite volume scheme for the 2D and 3D heterogeneous anisotropic diffusion equations

We propose a mimetic interpolation-free cell-centered finite volume scheme for diffusion equations on arbitrary two-dimensional (2D) and three dimensional (3D) unstructured, possibly non star-shaped or nonplanar meshes. This scheme borrows ideas from the well known mimetic finite difference methods. This interpolation-free means that we do not need to introduce auxiliary unknowns in the process of constructing the scheme, and the final scheme has only cell-centered unknowns. After carefully handling continuity of solutions and normal fluxes, we construct piecewise linear approximations on each cell. This makes our scheme applicable to arbitrary discontinuities and arbitrary mesh topologies. Interpolation-free makes the program easy to implement and the 2D and 3D programs have a unified framework. It is worth mentioning that the traditional star-shaped mesh assumption has been abolished, and the new scheme has more relaxed restrictions on meshes. Some representative and even challenging numerical examples verify that the scheme is robust and almost has the optimal convergence order on some benchmark polygonal and polyhedral meshes.