General partial differential equation solvers for saturated-unsaturated seepage

Is the neural tangent kernel of PINNs deep learning general partial differential equations always convergent?

To gain an understanding of the electroelastic properties of tactile piezoelectric sensors used in the characterization of soft tissue, the frequency-dependent electric impedance response of thick piezoelectric disks has been calculated using finite element modeling. To fit the calculated to the measured response, a new method was developed using harmonic overtones for tuning of the calculated effective elastic, piezoelectric, and dielectric parameters. To validate the results, the impedance responses of 10 piezoelectric disks with diameter-to-thickness ratios of 20, 6, and 2 have been measured from 10 kHz to 5 MHz. A two-dimensional, general purpose finite element partial differential equation solver with adaptive meshing capability run in the frequency-stepped mode, was used. The equations and boundary conditions used by the solver are presented. Calculated and measured impedance responses are presented, and resonance frequencies have been compared in detail. The comparison shows excellent agreement, with average relative differences in frequency of 0.27%, 0.19%, and 0.54% for the samples with diameter-to-thickness ratios of 20, 6, and 2, respectively. The method of tuning the effective elastic, piezoelectric, and dielectric parameters is an important step toward a finite element model that describes the properties of tactile sensors in detail.

Explicit high-order energy-preserving methods for general Hamiltonian partial differential equations

In this paper we shall introduce new constructions of approximate solutions of general linear partial differential equations with constant coefficients on the whole spaces, and establish fundamental estimates of the solutions depending on the inhomogeneous terms. This will be done by combining general ideas of the Tikhonov regularization and discretization of bounded linear operator equations on reproducing kernel Hilbert spaces. Furthermore, we will provide approximate solutions for the related inverse source problems.

Flash memory-based computing-in-memory (CIM) architectures have gained popularity due to their remarkable performance in various computation tasks of data processing, including machine learning, neuron networks, and scientific calculations. Especially in the partial differential equation (PDE) solver that has been widely utilized in scientific calculations, high accuracy, processing speed, and low power consumption are the key requirements. This work proposes a novel flash memory-based PDE solver to implement PDE with high accuracy, low power consumption, and fast iterative convergence. Moreover, considering the increasing current noise in nanoscale devices, we investigate the robustness against the noise in the proposed PDE solver. The results show that the noise tolerance limit of the solver can reach more than five times that of the conventional Jacobi CIM solver. Overall, the proposed flash memory-based PDE solver offers a promising solution for scientific calculations that require high accuracy, low power consumption, and good noise immunity, which could help to develop flash-based general computing.

Increasingly, ground water hydrologists and engineers must solve saturated/unsaturated flow and transport problems. The unsaturated component of flow is especially important when considering transport of contaminants that originate from surface or near‐surface sources such as waste spill sites, landfills, and hazardous waste storage sites. Numerical solution of such problems has been accomplished with special purpose computer programs developed for solution of unsaturated flow problems. However, general purpose partial differential equation (PDE) solvers also exist, which often provide greater model flexibility. In this paper, the performance of PDE2D, a general PDE solver, and SWMS_2D, a special purpose program, are compared for solution of unsaturated flow and transport problems. Both programs accurately calculated changes in water volume in the flow domain for four examples, but PDE2D did not calculate boundary fluxes as well as SWMS_2D, resulting in high mass balance errors. However, PDE2D yielded a smoother solution of the example solute transport problem, despite a large mass balance error. SWMS_2D offers ease of use and excellent mass balances while PDE2D offers the user greater flexibility in terms of the types of problems solved, and a grid generator that ensures that the element size and distribution of elements are appropriate for a given problem.

usual boundary conditions, and constraints u~U(t, z)~Em, where the control space U = U(t, z) is a given subset of Era. Here z =(z 1 . . . . . : ) are state variables and u =(u 1 .. . . , u m) control variables. As stated by P. K. RASHEVSKY [6], very general partial differential equations and systems can be written locally in the form (2). For instance, second order partial differential equations such as z~,+z~s= 1 u) can be immediately written in the form f ( t , s, z 1, zlt, z~,

A very general partial differential equation in space and time satisfied by the gene frequency in a monoecious population distributed continuously over an arbitrary habitat is derived. The treatment is restricted to a single diallelic locus in the absence of mutation and random drift, and it is supposed that time is continuous, births and deaths occur at random, and migration is independent of genotype. With the further assumptions that migration is isotropic and homogeneous, the population density is constant and uniform (as permitted by the population regulation mechanism included in the formulation), and Hardy-Weinberg proportions obtain locally, this partial differential equation reduces to the simplest multidimensional generalization of the classical Fisher-Haldane cline model. The efficacy of migration and selection in maintaining genetic variability at equilibrium in this model is investigated by deducing conditions for the existence of clines under various circumstances. The effects of the degree of dominance, a neutral belt between the regions where a particular allele is advantageous and deleterious, finiteness of the habitat, and habitat dimensionality are evaluated. Provided at least one of the alleles is favored only in a finite region, excluding the special case in which its total effective selective coefficient is zero, if conditions for supporting a cline are too unfavorable because migration is too strong, selection is too weak, or both, a cline cannot exist at all. Thus, unless there is overdominance, the population must be monomorphic. It is possible for a cline which can barely exist under the prevailing ecological circumstances to show a large amount of variation in gene frequency.

The primary aim is to develop a new class of pseudo-differential operators involving Linear Canonical Transform (LCT) and some of its essential results in Schwartz space and tempered distribution space. The secondary aim is to find applications of the LCT while constructing new general partial differential equations (Heat and Wave) and derive their solutions. Further, we have discussed particular cases of these partial differential equations. Furthermore, we have also shown the resolution of the solution of these equations graphically.

In Chapter 1 we presented basic ideas for the reduction of boundary value problems of the Laplacian to various forms of boundary integral equations based on the direct approach. This reduction can be easily extended to more general partial differential equations. Here we will consider, in particular, the Helmholtz equation, the Lame system, the Stokes equations and the biharmonic equation.

A system of two time‐dependent, first order differential equations of electrical membrane resistance and of hydrostatic pressure is derived. The equations are based on a set of flux equations for a 1,–1 valent electrolyte, combined with the equation of continuity, and the condition for conservation of mass. New insight into the functioning of the membrane oscillator is gained by the strict derivation of the relaxation equation of the averaged salt concentration with convection from the general partial differential equation. This relaxation equation is successfully transformed into the useful time‐dependent first order differential equation of the function of the electrical membrane resistance. Thereby the postulated relation for the concomitant changes of the electrical membrane resistance and the averaged salt concentration in the membrane phase is employed. The equations describe both the stationary and the non‐stationary states of Teorell's membrane oscillator. Particularly, the stationary pressure‐dependent, non‐monotonous current‐voltage curves and the current‐dependent, non‐monotonous pressure‐volume flux relations are calculated. These characteristics form the basis for mechano‐electrical transduction. Furthermore, the relaxation oscillations of membrane voltage, volume flux and hydrostatic pressure and the corresponding electroosmotic and hydrodynamic phase planes, representing stable limit cycles, are described. The stability criteria of the time‐dependent states are discussed.

Assessing basic research : Some partial indicators of scientific progress in radio astronomy

PDE blending surfaces with C2 continuity

ABSTRACTFor most systems, the appearance of Turing bifurcation means that it is possible to excite Turing–Turing bifurcation, thus inducing superimposed spatial patterns, multistable spatial patterns co‐existing, and others.It is undoubtedly of interest to qualitatively analyze the structures and stability of these spatially heterogeneous solutions generated by Turing–Turing bifurcation. Therefore, in this paper, with the aid of the center manifold theory and the normal form method, for general partial functional differential equations with nonlinear diffusion, the third‐order normal form is firstly derived. It is locally topologically equivalent to the primitive partial functional differential equations at the Turing–Turing bifurcation point. And then the explicit formulae for the coefficients in the normal form associated with three different spatial modes are presented. As an application, a predator–prey model with predator–taxis is considered. It is theoretically revealed that the system admits the coexistence of a pair of stable steady states with a single characteristic wavelength and the coexistence of four stable steady states with different single characteristic wavelengths. Further, the parameter regions in which these phenomena will arise are quantitatively given. It is illustrated that the self‐diffusion of prey and predator–taxis complement each other to promote the formation for the spatially homogeneous distributions of populations.

The partial differential equation that describes the flow, of non-Newtonian, power-law, slightly compressible fluids in porous media is derived. An approximate solution, in closed form, is developed for the unsteady-state flow behavior and verified by. two different methods. Using the unsteady-state solution, a method for analyzing injection test data is formulated and used to analyze four injection tests. Theoretical results were used to derive steady-state equations of flow, equivalent transient drainage radius, and a method for analyzing isochronal test data. The theoretical fundamentals of the flow, of non-Newtonian power-law fluids in porous media are established. Introduction Non-Newtonian power-law fluids are those that obey the relation = constant. Here, is the viscosity, e is the shear rate at which the viscosity is measured, and n is a constant. Examples of such fluids are polymers. This paper establishes the theoretical foundation of the flow of such fluids in porous media. The partial differential equation describing this flow is derived and solved for unsteady-state flow. In addition, a method for interpreting isochronal tests and an equation for calculating the equivalent transient drainage radius are presented. The unsteady-state flow solution provides a method for interpreting flow tests (such as injection tests).Non-Newtonian power-law fluids are injected into the porous media for mobility control, necessitating a basic porous media for mobility control, necessitating a basic understanding of the flow behavior of such fluids in porous media. Several authors have studied the porous media. Several authors have studied the rheological properties of these fluids using linear flow experiments and standard viscometers. Van Poollen and Jargon presented a theoretical study of these fluids. They described the flow by the partial differential equation used for Newtonian fluids and accounted for the effect of shear rate on viscosity by varying the viscosity as a function of space. They solved the equation numerically using finite difference. The numerical results showed that the pressure behavior vs time differed from that for Newtonian fluids. However, no methods for analyzing flow-test data (such as injection tests) were offered. This probably was because of the lack of analytic solution normally required to understand the relationship among the variables.Recently, injectivity tests were conducted using a polysaccharide polymer (biopolymer). The data showed polysaccharide polymer (biopolymer). The data showed anomalies when analyzed using methods derived for Newtonian fluids. Some of these anomalies appeared to be fractures. However, when the methods of analysis developed here were applied, the anomalies disappeared. Field data for four injectivity tests are reported and used to illustrate our analysis methods. Theoretical Consideration General Consideration The partial differential equation describing the flow of a non-Newtonian, slightly compressible power-law fluid in porous media derived in Appendix A is ..........(1) where the symbols are defined in the nomenclature. JPT P. 155