Two dimensional aerodynamic interference effects on oscillating airfoils with flaps in ventilated subsonic wind tunnels

The numerical computation of unsteady airloads acting upon thin airfoils with multiple leading and trailing-edge controls in two-dimensional ventilated subsonic wind tunnels is studied. The foundation of the computational method is strengthened with a new and more powerful mathematical existence and convergence theory for solving Cauchy singular integral equations of the first kind, and the method of convergence acceleration by extrapolation to the limit is introduced to analyze airfoils with flaps. New results are presented for steady and unsteady flow, including the effect of acoustic resonance between ventilated wind-tunnel walls and airfoils with oscillating flaps. The computer program TWODI is available for general use and a complete set of instructions is provided

The determination of gravity anomalies from geoid heights using the inverse Stokes' formula, Fourier transforms, and least squares collocation

A numerical method for the determination of gravity anomalies from geoid heights is described using the inverse Stokes formula. This discrete form of the inverse Stokes formula applies a numerical integration over the azimuth and an integration over a cubic interpolatory spline function which approximates the step function obtained from the numerical integration. The main disadvantage of the procedure is the lack of a reliable error measure. The method was applied on geoid heights derived from GEOS-3 altimeter measurements in the calibration area of the GEOS-3 satellite

Lidstone–Euler interpolation and related high even order boundary value problem

AbstractWe consider the Lidstone–Euler interpolation problem and the associated Lidstone–Euler boundary value problem, in both theoretical and computational aspects. After a theorem of existence and uniqueness of the solution to the Lidstone–Euler boundary value problem, we present a numerical method for solving it. This method uses the extrapolated Bernstein polynomials and produces an approximating convergent polynomial sequence. Particularly, we consider the fourth-order case, arising in various physical models. Finally, we present some numerical examples and we compare the proposed method with a modified decomposition method for a tenth-order problem. The numerical results confirm the theoretical and computational ones

Use of Machine Learning for Automated Convergence of Numerical Iterative Schemes

Convergence of a numerical solution scheme occurs when a sequence of increasingly refined iterative solutions approaches a value consistent with the modeled phenomenon. Approximations using iterative schemes need to satisfy convergence criteria, such as reaching a specific error tolerance or number of iterations. The schemes often bypass the criteria or prematurely converge because of oscillations that may be inherent to the solution. Using a Support Vector Machines (SVM) machine learning approach, an algorithm is designed to use the source data to train a model to predict convergence in the solution process and stop unnecessary iterations. The discretization of the Navier Stokes (NS) equations for a transient local hemodynamics case requires determining a pressure correction term from a Poisson-like equation at every time-step. The pressure correction solution must fully converge to avoid introducing a mass imbalance. Considering time, frequency, and time-frequency domain features of its residual’s behavior, the algorithm trains an SVM model to predict the convergence of the Poisson equation iterative solver so that the time-marching process can move forward efficiently and effectively. The fluid flow model integrates peripheral circulation using a lumped-parameter model (LPM) to capture the field pressures and flows across various circulatory compartments. Machine learning opens the doors to an intelligent approach for iterative solutions by replacing prescribed criteria with an algorithm that uses the data set itself to predict convergence

A Finite Element Splitting Extrapolation for Second Order Hyperbolic Equations

Splitting extrapolation is an efficient technique for solving large scale scientific and engineering problems in parallel. This article discusses a finite element splitting extrapolation for second order hyperbolic equations with time-dependent coefficients. This method possesses a higher degree of parallelism, less computational complexity, and more flexibility than Richardson extrapolation while achieving the same accuracy. By means of domain decomposition and isoparametric mapping, some grid parameters are chosen according to the problem. The multiparameter asymptotic expansion of the d-quadratic finite element error is also established. The splitting extrapolation formulas are developed from this expansion. An approximation with higher accuracy on a globally fine grid can be computed by solving a set of smaller discrete subproblems on different coarser grids in parallel. Some a posteriori error estimates are also provided. Numerical examples show that this method is efficient for solving discontinuous problems and nonlinear hyperbolic equations

Numerical methods for computing Casimir interactions

We review several different approaches for computing Casimir forces and related fluctuation-induced interactions between bodies of arbitrary shapes and materials. The relationships between this problem and well known computational techniques from classical electromagnetism are emphasized. We also review the basic principles of standard computational methods, categorizing them according to three criteria---choice of problem, basis, and solution technique---that can be used to classify proposals for the Casimir problem as well. In this way, mature classical methods can be exploited to model Casimir physics, with a few important modifications.Comment: 46 pages, 142 references, 5 figures. To appear in upcoming Lecture Notes in Physics book on Casimir Physic