SEARCH · Search NASA

Results for “Collocation scheme”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2018

A Müntz wavelets collocation method for solving fractional differential equations

In this paper, we introduce a new set of wavelets, namely Müntz wavelets. These wavelets are defined based on Müntz Legendre polynomials which have the property of orthogonality. We stabilize the wavelet coefficients by utilizing Jacobi polynomials to construct the modified Müntz wavelets. Next, we present a collocation scheme based on Müntz wavelets for solving fractional differential equations. The fractional derivative is described in the Caputo sense. Some error estimates are given and numerical examples are included to demonstrate the applicability and accuracy of the proposed method. This includes the elastic-dissipative mathematical model of linear mechanical devices known as Boltzman operator.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1984

Stresses around holes in pin-loaded orthotropic plates

The effects of pin elasticity, pin/hole friction, and pin/hole clearance on the stresses around a hole in a pin-loaded orthotropic plate were studied. The problem was modeled as a planar problem, the pin and plate being two elastic bodies interacting through contact. Complex variable elasticity, a collocation scheme to enforce the pin/hole boundary conditions, and an iteration procedure to find the contact, slip, and no-slip regions were used to obtain numerical results. The results indicate that pin elasticity is not as important a variable as friction and clearance are. A rigid pin assumption, as opposed to the popular cosinusoidal radial traction, can be used without much error but clearance and friction should not be ignored.

Hyer, M. W.↗

NASA NTRS1984

A spectral element method for fluid dynamics - Laminar flow in a channel expansion

A spectral element method that combines the generality of the finite element method with the accuracy of spectral techniques is proposed for the numerical solution of the incompressible Navier-Stokes equations. In the spectral element discretization, the computational domain is broken into a series of elements, and the velocity in each element is represented as a high-order Lagrangian interpolant through Chebyshev collocation points. The hyperbolic piece of the governing equations is then treated with an explicit collocation scheme, while the pressure and viscous contributions are treated implicitly with a projection operator derived from a variational principle. The implementation of the technique is demonstrated on a one-dimensional inflow-outflow advection-diffusion equation, and the method is then applied to laminar two-dimensional (separated) flow in a channel expansion. Comparisons are made with experiment and previous numerical work.

Patera, A. T.↗

NASA NTRS1985

Higher order methods for convection-diffusion problems

This paper applies C1 cubic Hermite polynomials embedded in an orthogonal collocation scheme to the spatial discretization of the unsteady nonlinear Burgers equation as a model of the equations of fluid mechanics. The temporal discretization is carried out by means of either a noniterative finite difference or an iterative finite difference procedure. Results of this method are compared with those of a second-order finite difference scheme and a splined-cubic Taylor's series scheme. Stability limits are derived and the matrix structure of the several schemes are compared.

Murphy, J. D.↗

NASA NTRS1987

The effects of pin elasticity, clearance, and friction on the stresses in a pin-loaded orthotropic plate

The effects of pin elasticity, clearance, and friction on the stresses in a pin loaded orthotropic plate are studied. The effects are studied by posing the problem as a planar contact elasticity problem, the pin and the plate being two elastic bodies which interact through contact. Coulomb friction is assumed, the pin loads the plate in one of its principal material directions, and the plate is infinite in extent. A collocation scheme and interaction, in conjunction with a complex variable series solution, are used to obtain numerical results. The contact region between the plate and pin is unknown and must be solved for as part of the solution. The same is true of the region of friction induced no slip. Two pin stiffnesses, two clearance levels, two friction levels and two laminates, a (0/+ or - 45/90)s and a (02/+ or - 45)s, are studied. The effects of pin elasticity, clearance, and friction on the load capacity of the plate are assessed by comparing the load capacity of the plate with the capacity when the pin is rigid, perfectly fitting, and frictionless.

Hyer, M. W.↗

NASA NTRS1989

A multidomain spectral collocation method for the Stokes problem

A multidomain spectral collocation scheme is proposed for the approximation of the two-dimensional Stokes problem. It is shown that the discrete velocity vector field is exactly divergence-free and we prove error estimates both for the velocity and the pressure.

Landriani, G. Sacchi↗

NASA NTRS1990

Domain identification in impedance computed tomography by spline collocation method

A method for estimating an unknown domain in elliptic boundary value problems is considered. The problem is formulated as an inverse problem of integral equations of the second kind. A computational method is developed using a splice collocation scheme. The results can be applied to the inverse problem of impedance computed tomography (ICT) for image reconstruction.

Kojima, Fumio↗

NASA NTRS2014

Time-Spectral Rotorcraft Simulations on Overset Grids

The Time-Spectral method is derived as a Fourier collocation scheme and applied to NASA's overset Reynolds-averaged Navier-Stokes (RANS) solver OVERFLOW. The paper outlines the Time-Spectral OVERFLOWimplementation. Successful low-speed laminar plunging NACA 0012 airfoil simulations demonstrate the capability of the Time-Spectral method to resolve the highly-vortical wakes typical of more expensive three-dimensional rotorcraft configurations. Dealiasing, in the form of spectral vanishing viscosity (SVV), facilitates the convergence of Time-Spectral calculations of high-frequency flows. Finally, simulations of the isolated V-22 Osprey tiltrotor for both hover and forward (edgewise) flight validate the three-dimensional Time-Spectral OVERFLOW implementation. The Time-Spectral hover simulation matches the time-accurate calculation using a single harmonic. Significantly more temporal modes and SVV are required to accurately compute the forward flight case because of its more active, high-frequency wake.

Rotocraft↗

DOE OSTI2021

Asymptotically Compatible Reproducing Kernel Collocation and Meshfree Integration for Nonlocal Diffusion

Reproducing kernel (RK) approximations are meshfree methods that construct shape functions from sets of scattered data. We present an asymptotically compatible (AC) RK collocation method for nonlocal diffusion models with Dirichlet boundary condition. The numerical scheme is shown to be convergent to both nonlocal diffusion and its corresponding local limit as nonlocal interaction vanishes. The analysis is carried out on a special family of rectilinear Cartesian grids for a linear RK method with designed kernel support. The key idea for the stability of the RK collocation scheme is to compare the collocation scheme with the standard Galerkin scheme, which is stable. In addition, assembling the stiffness matrix of the nonlocal problem requires costly computational resources because high-order Gaussian quadrature is necessary to evaluate the integral. We thus provide a remedy to the problem by introducing a quasi-discrete nonlocal diffusion operator for which no numerical quadrature is further needed after applying the RK collocation scheme. The quasi-discrete nonlocal diffusion operator combined with RK collocation is shown to be convergent to the correct local diffusion problem by taking the limits of nonlocal interaction and spatial resolution simultaneously. The theoretical results are then validated with numerical experiments. We additionally illustrate a connection between the proposed technique and an existing optimization based approach based on generalized moving least squares.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1996

Error and Complexity Analysis for a Collocation-Grid-Projection Plus Precorrected-FFT Algorithm for Solving Potential Integral Equations with LaPlace or Helmholtz Kernels

In this paper we derive error bounds for a collocation-grid-projection scheme tuned for use in multilevel methods for solving boundary-element discretizations of potential integral equations. The grid-projection scheme is then combined with a precorrected FFT style multilevel method for solving potential integral equations with 1/r and e(sup ikr)/r kernels. A complexity analysis of this combined method is given to show that for homogeneous problems, the method is order n natural log n nearly independent of the kernel. In addition, it is shown analytically and experimentally that for an inhomogeneity generated by a very finely discretized surface, the combined method slows to order n(sup 4/3). Finally, examples are given to show that the collocation-based grid-projection plus precorrected-FFT scheme is competitive with fast-multipole algorithms when considering realistic problems and 1/r kernels, but can be used over a range of spatial frequencies with only a small performance penalty.

Phillips, J. R.↗

NASA NTRS1998

Direct Numerical Simulation of a Coolant Jet in a Periodic Crossflow

A Direct Numerical Simulation of a coolant jet injected normally into a periodic crossflow is presented. The physical situation simulated represents a periodic module in a coolant hole array with a heated crossflow. A collocated finite difference scheme is used which is fifth-order accurate spatially and second-order accurate temporally. The scheme is based on a fractional step approach and requires the solution of a pressure-Poisson equation. The simulations are obtained for a blowing ratio of 0.25 and a channel Reynolds number of 5600. The simulations reveal the dynamics of several large scale structures including the Counter-rotating Vortex Pair (CVP), the horse-shoe vortex, the shear layer vortex, the wall vortex and the wake vortex. The origins and the interactions of these vortical structures are identified and explored. Also presented are the turbulence statistics and how they relate to the flow structures.

Sharma, Chirdeep↗

NASA NTRS2000

Discontinuous dual-primal mixed finite elements for elliptic problems

We propose a novel discontinuous mixed finite element formulation for the solution of second-order elliptic problems. Fully discontinuous piecewise polynomial finite element spaces are used for the trial and test functions. The discontinuous nature of the test functions at the element interfaces allows to introduce new boundary unknowns that, on the one hand enforce the weak continuity of the trial functions, and on the other avoid the need to define a priori algorithmic fluxes as in standard discontinuous Galerkin methods. Static condensation is performed at the element level, leading to a solution procedure based on the sole interface unknowns. The resulting family of discontinuous dual-primal mixed finite element methods is presented in the one and two-dimensional cases. In the one-dimensional case, we show the equivalence of the method with implicit Runge-Kutta schemes of the collocation type exhibiting optimal behavior. Numerical experiments in one and two dimensions demonstrate the order accuracy of the new method, confirming the results of the analysis.

Bottasso, Carlo L.↗

NASA NTRS2001

Comparison of Implicit Collocation Methods for the Heat Equation

We combine a high-order compact finite difference scheme to approximate spatial derivatives arid collocation techniques for the time component to numerically solve the two dimensional heat equation. We use two approaches to implement the collocation methods. The first one is based on an explicit computation of the coefficients of polynomials and the second one relies on differential quadrature. We compare them by studying their merits and analyzing their numerical performance. All our computations, based on parallel algorithms, are carried out on the CRAY SV1.

Kouatchou, Jules↗

DOE OSTI2022

A hybrid meshfree discretization to improve the numerical performance of peridynamic models

Efficient and accurate calculation of spatial integrals is of major interest in the numerical implementation of peridynamics (PD). The standard way to perform this calculation is a particle-based approach that discretizes the strong form of the PD governing equation. This approach has rapidly been adopted by the PD community since it offers some advantages. Additionally, it is computationally cheaper than other available schemes, can conveniently handle material separation, and effectively deals with nonlinear PD models. Nevertheless, PD models are still computationally very expensive compared with those based on the classical continuum mechanics theory, particularly for large-scale problems in three dimensions. This results from the nonlocal nature of the PD theory which leads to interactions of each node of a discretized body with multiple surrounding nodes. Here, we propose a new approach to significantly boost the numerical efficiency of PD models. We propose a discretization scheme that employs a simple collocation procedure and is truly meshfree; i.e., it does not depend on any background integration cells. In contrast to the standard scheme, the proposed scheme requires a much smaller set of neighboring nodes (keeping the same physical length scale) to achieve a specific accuracy and is thus computationally more efficient. Our new scheme is applicable to the case of linear PD models and within neighborhoods where the solution can be approximated by smooth basis functions. Therefore, to fully exploit the advantages of both the standard and the proposed schemes, a hybrid discretization is presented that combines both approaches within an adaptive framework. The high performance of the developed framework is illustrated by several numerical examples, including brittle fracture and corrosion problems in two and three dimensions.

42 ENGINEERING↗

NASA NTRS1986

Preconditioning for first-order spectral discretization

Efficient solution of the equations from spectral discretizations is essential if the high-order accuracy of these methods is to be realized. Direct solution of these equations is rarely feasible, thus iterative techniques are required. A preconditioning scheme for first-order Chebyshev collocation operators is proposed herein, in which the central finite difference mesh is finer than the collocation mesh. Details of the proper techniques for transferring information between the meshes are given here, and the scheme is analyzed by examination of the eigenvalue spectra of the preconditioned operators. The effect of artificial viscosity required in the inversion of the finite difference operator is examined. A second preconditioning scheme, involving a high-order upwind finite difference operator of the van Leer type is also analyzed to provide a comparison with the present scheme. Finally, the performance of the present scheme is verified by application to several test problems.

Streett, C. L.↗

NASA NTRSDate not supplied

Explicit Discontinuous Galerkin Methods for Conservation Laws

The two explicit DG methods in this study are based on a ‘predictor-corrector’ formulation, the first introduced by Lörcher, Gassner, and Munz (2007, 2008) called space–time expansion discontinuous Galerkin or STE-DG scheme, and the second, introduced independently by the author (Huynh 2006, 2013) called the upwind moment scheme. The predictor step of the two methods is essentially identical using a Cauchy-Kovalevsky (CK) procedure, which involves no interaction of the data among neighboring cells. The corrector step also shares the same space-time integration formulation and is where interaction of the data among neighboring cells takes place; the difference, however, is in how the resulting space-time volume integral is estimated. As a consequence of the different estimates, for the case of advection in one spatial dimension (1D), the moment scheme has a CFL (Courant-Friedrichs-Lewy) condition of 1 for all p and is accurate to order 2p+1, i.e., it possesses the super accuracy property, whereas the STE-DG method has a more restrictive CFL condition and is accurate to the expected order of p+1. For 1D advection, compared with the CFL conditions of 1/(2p+1) of standard RK-DG (Runge-Kutta) scheme where space and time discretization are of the same order, the moment scheme allows a significantly larger time step size. It also turns out that the scheme yields a result identical to Van Leer’s scheme III (1977), which amounts to shifting the data a distance of advection corresponding to the time step and projecting the result onto the space of polynomial solutions. Contrary to Van Leer’s approach, however, the space-time ‘predictor-corrector’ formulation facilitates extensions to the case of systems of equations. Concerning 2D extensions, in the case of advection, when the flow is along the diagonal direction, the CFL conditions for the moment schemes become restrictive as will be shown by Fourier (Von Neumann) stability and accuracy analyses. Since the moment scheme employs the right Radau points as collocation points in time, the method is closely related to the implicit Radau IIA scheme, which is stable for any time step size. The role of Radau IIA in relieving stability restriction for these explicit DG schemes remains to be explored

Discontinuous Galerkin↗

NASA NTRS2000

Assimilation of (A) TOVS data at the NASA Goddard Data Assimilation Office

At the NASA Goddard Data Assimilation Office (DAO), a 1D variational radiance assimilation system has been developed. This system, called DAOTOVS (Tiros operational vertical sounder (TOVS)), uses (A)TOVS level 1b radiances. It has been implemented within the DAO's semi-operational system as well as within the next generation data assimilation system that uses a finite-volume dynamical core. We will show results from (A)TOVS assimilation, including stratospheric analyses and validation. We will also describe our systematic error correction scheme which is based on collocated radiosondes.

Joiner, Joanna↗

NASA NTRS2013

High-Order Space-Time Methods for Conservation Laws

Current high-order methods such as discontinuous Galerkin and/or flux reconstruction can provide effective discretization for the spatial derivatives. Together with a time discretization, such methods result in either too small a time step size in the case of an explicit scheme or a very large system in the case of an implicit one. To tackle these problems, two new high-order space-time schemes for conservation laws are introduced: the first is explicit and the second, implicit. The explicit method here, also called the moment scheme, achieves a Courant-Friedrichs-Lewy (CFL) condition of 1 for the case of one-spatial dimension regardless of the degree of the polynomial approximation. (For standard explicit methods, if the spatial approximation is of degree p, then the time step sizes are typically proportional to 1/p(exp 2)). Fourier analyses for the one and two-dimensional cases are carried out. The property of super accuracy (or super convergence) is discussed. The implicit method is a simplified but optimal version of the discontinuous Galerkin scheme applied to time. It reduces to a collocation implicit Runge-Kutta (RK) method for ordinary differential equations (ODE) called Radau IIA. The explicit and implicit schemes are closely related since they employ the same intermediate time levels, and the former can serve as a key building block in an iterative procedure for the latter. A limiting technique for the piecewise linear scheme is also discussed. The technique can suppress oscillations near a discontinuity while preserving accuracy near extrema. Preliminary numerical results are shown

Huynh, H. T.↗