Towards a Hybrid Adjoint Approach for Arbitrarily Complex Partial Differential Equations

Development and assessment of adjoint sensitivity analysis method for transient two-phase flow simulations

Adjoint-based techniques can provide the sensitivity of an objective function to any number of parameters of a simulation inexpensively at roughly the cost of a single additional flow calculation. This information can be used to perform sensitivity analyses, aerodynamic shape optimization, and to estimate the error in the objective function due to numerical discretization. Existing approaches to derive the numerically discretized adjoint equations involve the so-called discrete and continuous methods, which differ in the order at which discretization and linearization steps are performed. The effect of these contrasting approaches is that they have both strengths and weaknesses over each other in the form of complexity of the formulation and computational expense of the solution. In this paper, we further develop the hybrid approach of Taylor et al. (2012) that combines elements of the continuous and discrete methods with the intention of capturing both their advantages: reducing the time spent on mathematical derivation of the continuous adjoint equations, lowering the computational requirements of the discrete adjoint equations, and generally improving the quality of the adjoint solution. The specific approach investigated in this paper treats the flow conservation equations in a continuous manner and the turbulence transport equations discretely. The framework is designed such that additional transport equations for any turbulence model can be seamlessly included in a discrete fashion for coupling with the mean flow equations. The methodology is demonstrated in an optimization problem of lift-constrained drag minimization of an airfoil in transonic turbulent flow.

Aiming at the heat conduction equation, the heat conduction model is a very typical example of the partial differential equation. This paper focuses on the development and preliminary knowledge of the partial differential equation, the establishment of the heat conduction model and the solution of the heat conduction model. The partial differential equation with fractional Laplace operator is a kind of typical fractional partial differential equation, which has important applications in the fields of science and engineering. Fractional Laplace operator is a kind of nonlocal quasi-differential operator, which is the infinitesimal generator of L é vy steady-state process. It is essentially different from the classical Laplace operator, which leads to the disappearance of some classical properties, which generally brings difficulties to the study of this kind of problems. It is very difficult to solve the explicit solution of partial differential equation with fractional Laplace operator. In order to improve the accuracy and speed of solving partial differential equation, a high precision method for solving partial differential equation based on deep learning network is proposed. A multilayer radial basis function neural network is constructed and its structural layers are determined. By introducing subnetworks into each layer of the multilayer radial basis function neural network, a composite multilayer radial basis function neural network is constructed to improve the real function approximation performance and operation accuracy of the multilayer radial basis function neural network. The high-precision compound multi-layer radial basis function neural network is used to solve the partial differential equation. By giving a specific example of solving the partial differential equation, the solution effect of the method is tested. The results show that the method has very high solution accuracy under the four-layer network structure, which can improve about 1.5 orders of magnitude compared with the first layer. Under different number of training samples, it has higher solution accuracy and speed, and the comprehensive performance is superior. Therefore, it is a practical but challenging work to study the existence of solutions of partial differential equations with fractional Laplace operators. This paper mainly describes the research progress and trends of the existence of solutions of several classes of partial differential equations with fractional Laplace operators, including some of the work done by the author in this field in recent years.

This paper explores the solutions of certain fractional partial differential equations using two methods; the first method involves separation of variables, which is a common technique for solving partial differential equations. However, since many equations cannot be separated in this way, the tensor product of Banach spaces method is applied to find the atomic solutions. To solve the resulting ordinary differential equations, the reproducing Kernel Hilbert space method is used to find numerical solutions, which are then used to find the numerical solution of the partial differential equation. The residual errors indicate that this method is effective and powerful. In summary, this paper presents a study on the solutions of certain fractional partial differential equations using two methods and demonstrates the effectiveness of these methods in finding numerical solutions.

The fractional partial differential equation is an equation with non-integer partial derivative which is now widely used in solving various problems so that in recent years, many researchers are interested to study fractional partial differential equation. The equation which is used in this paper is the non-homogeneous fractional partial differential equation. The aim of this study is to determine general solution of the partial differential equation through three different order equations (α; 1), (β; 2) and (α; β), where 0 < α ≤ 2 and 0 < β ≤ 1 by using Homotopy Analysis Method (HAM). Furthermore, by taking a convergent sequence, it will be analyzed that the convergence of the sequence of differential equations results in the sequence of solutions of partial fractional non-homogeneous partial differential equations converging to the solution function of the fractional partial differential equation.

Several research papers experimentally validated the derived theoretical solution of the partial differential equation (PDE) of Transient Land Drainage (TLD) using laboratory sand tanks, Hele-Shaw models, resistance networks, and field data. In this research, the main objective is to use an Arduino microcontroller linked with the LabVIEW software to derive experimentally-by using an RC analog circuit connected to the Arduino microcontroller - a solution to Kirchhoff’s partial differential equation describing the unsteady electrical flow through the RC circuit. The RC circuit electrically simulates a non-dimensional virtual aquifer for which an experimental Fourier sine series solution in terms of the analog components is to be derived and used for comparison with the theoretical Fourier sine series solution of the transient partial differential equation of land drainage. The derived Fourier sine series solution of the partial differential equation governing the unsteady electrical flow through the RC analog is written in a non-dimensional form to enable its general application to the case study of this research. Isochrones of the water table drawdown based on those derived from the analog solution are similar in trend and behavior to those of the theoretical solution, indicating that both Fourier solutions are physically and mathematically identical.

Numerical convective schemes based on accurate computation of space derivatives

Gain-margin based discrete-continuous method for the stability analysis of real-time hybrid simulation systems

Interface solutions of partial differential equations with point singularity

A comparison of differential quadrature methods for the solution of partial differential equations

We consider the solution of elliptic partial differential equations with random input data. In particular, we are interested in the mean and the correlation of the solution of these partial differential equations. Once the correlation is available, the variance can be computed from its diagonal. If the dependence of the solution on the random input data is linear and the mean of the input data is available, the mean of the solution can be computed whenever the corresponding partial differential equation can be solved. When the correlation of the input data is available, the correlation of the solution is given as the solution of a higher dimensional partial differential equation in the product domain. A Galerkin discretization yields a matrix equation with a typically densely populated right-hand side. This and the higher dimension of the problem makes the equation prohibitively expensive to solve. Since existing approaches to these correlation equations are known to struggle for roughly correlated input data, we use hierarchical matrices and their corresponding arithmetic to represent and solve the matrix equation. The regularity assumptions required for hierarchical matrices can be justified on sufficiently smooth domains. The feasibility of the approach is illustrated for finite element and boundary element discretizations and several specialities of hierarchical matrices for finite element discretizations are discussed. An almost linear scaling with respect to the dimension of the used finite element spaces is verified. The nonlinear dependence of the solution on the random input data is exemplarily discussed for the case of random domains. Extending previous perturbation approaches, we can linearize the problem and can compute the mean and the correlation up to third order accuracy in almost linear time. We discuss that these spaces can become even fourth order accurate under certain circumstances. A full convergence analysis is presented, which shows that higher order boundary elements are required for the solution of the problem. Therefore, we develop a fast multipole method for higher order boundary element methods on parametric surfaces.

HiFlow3 Version 2.0 continues the path as a multi-purpose finite element software, which provides powerful tools for efficient and accurate solution of a wide range of problems modeled by partial differen- tial equations (PDEs). New features and functionalities, which allow to run numerical simulations with more advanced solution algorithms and discretizations in comparison to previous releases of HiFlow3, have been implemented. These comprise Uncertainty Quantification (UQ), energy-efficient multigrid techniques, Schur complement solvers for saddle-point problems, extended third-party library support, and adaptive local mesh refinement in a parallel computing environment. Furthermore, HiFlow3 has been successfully integrated into advanced and state-of-the-art simulation environments by means of the Medical Simulation Markup Language (MSML), for example.The presented new algorithms and features as well as general under-the-hood improvements have enabled excellent and relevant research activities in the fields of both medical engineering and meteorol- ogy and environmental sciences. The described show cases demonstrate the potential and advantages, which HiFlow3 can offer in performing a numerical simulation by means of finite element methods. Especially, the high performance computing capabilities of HiFlow3 – not only in the mentioned fields of applications, but also in general – have been significantly improved in Version 2.0.

Wasserstein generative adversarial uncertainty quantification in physics-informed neural networks

Similarity solutions of nonlinear partial differential equations invariant to a family of affine groups

Moving boundary and boundary value problems occur in many physical and engineering processes involving heat transfer and phase changes, they range from icing problems, for example the Stefan problem, to the melting of metals and plastics using lasers. This book reviews and compares some well known and some recently developed methods to solve these problems. Recent advances in the finite difference solution of linear and non-linear partial differential equations are presented and algorithmic manipulations that enhance the computational efficiency are incorporated into the overall schemes. In this new approach, the the finite difference substitutions are made into the solutions of the partial differential equation rather than classically into the partial differential equation itself. This approach can be exploited for many types of partial differential equation and they may therefore be of interest to researchers studying a wide range of applications in physical problems such as diffusion and heat transfer equations, optoelectronic paraxial wave equations, the quantum mechanical wave equations, time dependent pulse propagation, and similar WKB-derived problemsM