SEARCH · Engineering Papers

Results for “precondition”

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 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS1999

Efficient Multi-Stage Time Marching for Viscous Flows via Local Preconditioning

A new method has been developed to accelerate the convergence of explicit time-marching, laminar, Navier-Stokes codes through the combination of local preconditioning and multi-stage time marching optimization. Local preconditioning is a technique to modify the time-dependent equations so that all information moves or decays at nearly the same rate, thus relieving the stiffness for a system of equations. Multi-stage time marching can be optimized by modifying its coefficients to account for the presence of viscous terms, allowing larger time steps. We show it is possible to optimize the time marching scheme for a wide range of cell Reynolds numbers for the scalar advection-diffusion equation, and local preconditioning allows this optimization to be applied to the Navier-Stokes equations. Convergence acceleration of the new method is demonstrated through numerical experiments with circular advection and laminar boundary-layer flow over a flat plate.

Kleb, William L.↗

NASA NTRSDate not supplied

Implicit Preconditioning for Explicit Multigrid Solvers on Cut-Cell Cartesian Meshes

This work assesses the effectiveness of linearized implicit Euler preconditioning for multigrid solvers using an unpreconditioned, Jacobian-free Newton Krylov method to converge the linear system of equations. Multigrid convergence rates improve to approximately 0.75 across the cases tested including a Mach 2 supersonic wedge, transonic NACA 0012 airfoil, and ONERA M6 wing. While larger Krylov subspaces increase the convergence rate, they also increase the computational cost, such that 4-8 Krylov vectors often offers the fastest turnaround. Further reductions in computational cost are achieved with a sequential hybrid preconditioner that begins with the explicit multigrid solver before transitioning to the preconditioned algorithm later on. In addition, a novel implementation of dual time stepping is extended to include both common BDF methods as well as high-order implicit Runge-Kutta schemes. This particular formulation, which uses A −1 preconditioning, is amenable to matrix-free solvers, and the L-stable methods are especially suited for meshes with arbitrarily small cut-cells. Asymptotic order of convergence is demonstrated for BDF1, BDF2, SDIRK2, and 3rd-order Radau IIA time integration with unsteady 2D vortex simulations.

ARMD↗

DOE OSTI2020

Pyrolysis Vapor and Bio-Oil Preconditioning via Ex Situ Hydrodeoxygenation and Alkylation Using a Heteropolyacid Catalyst

Ex situ catalytic preconditioning of biomass pyrolysis vapors (partial deoxygenation) has the potential for both reducing the oxygen content and beneficially influencing the oxygen speciation of condense bio-oils for enhancement of downstream condensed phase upgrading processes. Here, condensed phase upgrading processes may refer to both coprocessing and hydroprocessing with vacuum gas-oil (VGO). Reduction of oxygen content to produce bio-oils with greater hydrophobicity can enhance VGO-solubility while tailoring of oxygen speciation by removal of reactive carbonyls (e.g., acetic acid) can lessen catalytic deactivation via coking reactions. 1,2 Furthermore, enhanced VGO-solubility and reduced coking leads to greater biogenic carbon incorporation in end-products from condensed phase upgrading processes. By leveraging the hydrodeoxygenation and alkylation activity of a titania-supported molybdenum-heteropolyacid (Mo-HPA/TiO2) catalyst, pyrolysis vapors and resulting condensed bio-oils were effectively preconditioned for downstream upgrading. 3 Data regarding the use of Mo-HPA/TiO2 for pyrolysis vapor preconditioning will be presented on, including method of employment, catalyst characterization, and impacts on bio-oil in relation to downstream condensed phase upgrading processes. Our approach is to: Compare Mo-HPA/TiO2 product distribution to HZSM-5; Investigate Mo-HPA/TiO2 regeneration and characterize changes in acid site character; Assess product differences in terms of enhanced VGO-solubility.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2022

Space-Time Block Preconditioning for Incompressible Flow

Parallel-in-time methods have become increasingly popular in the simulation of time-dependent numerical PDEs, allowing for the efficient use of additional message passing interface processes when spatial parallelism saturates. Most methods treat the solution and parallelism in space and time separately. In contrast, all-at-once methods solve the full space-time system directly, largely treating time as simply another spatial dimension. All-at-once methods offer a number of benefits over separate treatment of space and time, most notably significantly increased parallelism and faster time to solution (when applicable). However, the development of fast, scalable all-at-once methods has largely been limited to time-dependent (advection-)diffusion problems. This paper introduces the concept of space-time block preconditioning for the all-at-once solution of incompressible flow. By extending well-known concepts of spatial block preconditioning to the space-time setting, we develop a block preconditioner whose application requires the solution of a space-time (advection-)diffusion equation in the velocity block, coupled with a pressure Schur complement approximation consisting of independent spatial solves at each time-step, and a space-time matrix-vector multiplication. The new method is tested on four classical models in incompressible flow. Finally, the results indicate perfect scalability in refinement of spatial and temporal mesh spacing, perfect scalability in nonlinear Picard iteration count when applied to a nonlinear Navier--Stokes problem, and minimal overhead in terms of number of preconditioner applications compared with sequential time-stepping.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Scalable preconditioning for the stabilized contact mechanics problem

We present a family of preconditioning strategies for the contact problem in fractured and faulted porous media. We combine low-order continuous finite elements to simulate the bulk deformation with piecewise constant Lagrange multipliers to impose the frictional contact constraints. This formulation is not uniformly inf-sup stable and requires stabilization. We improve previous work by Franceschini et al. (2020) by introducing a novel jump stabilization technique that requires only local geometrical and mechanical properties. We then design scalable preconditioning strategies that take advantage of the block structure of the Jacobian matrix using a physics-based partitioning of the unknowns by field type, namely displacement and Lagrange multipliers. The key to the success of the proposed preconditioners is a pseudo-Schur complement obtained by eliminating the Lagrange multiplier degrees of freedom, which can then be efficiently solved using an optimal multigrid method. Numerical results, including complex real-world problems, are presented to illustrate theoretical properties, scalability and robustness of the preconditioner. A comparison with other approaches available in the literature is also provided.

58 GEOSCIENCES↗

DOE OSTI2022

Accelerating astronomical and cosmological inference with preconditioned Monte Carlo

ABSTRACT We introduce preconditioned Monte Carlo (PMC), a novel Monte Carlo method for Bayesian inference that facilitates efficient sampling of probability distributions with non-trivial geometry. PMC utilizes a Normalizing Flow (NF) in order to decorrelate the parameters of the distribution and then proceeds by sampling from the preconditioned target distribution using an adaptive Sequential Monte Carlo (SMC) scheme. The results produced by PMC include samples from the posterior distribution and an estimate of the model evidence that can be used for parameter inference and model comparison, respectively. The aforementioned framework has been thoroughly tested in a variety of challenging target distributions achieving state-of-the-art sampling performance. In the cases of primordial feature analysis and gravitational wave inference, PMC is approximately 50 and 25 times faster, respectively, than nested sampling (NS). We found that in higher dimensional applications, the acceleration is even greater. Finally, PMC is directly parallelisable, manifesting linear scaling up to thousands of CPUs.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

End-to-end GPU acceleration of low-order-refined preconditioning for high-order finite element discretizations

In this article, we present algorithms and implementations for the end-to-end GPU acceleration of matrix-free low-order-refined preconditioning of high-order finite element problems. The methods described here allow for the construction of effective preconditioners for high-order problems with optimal memory usage and computational complexity. The preconditioners are based on the construction of a spectrally equivalent low-order discretization on a refined mesh, which is then amenable to, for example, algebraic multigrid preconditioning. The constants of equivalence are independent of mesh size and polynomial degree. For vector finite element problems in H(curl) and H(div) (e.g., for electromagnetic or radiation diffusion problems), a specially constructed interpolation–histopolation basis is used to ensure fast convergence. Detailed performance studies are carried out to analyze the efficiency of the GPU algorithms. The kernel throughput of each of the main algorithmic components is measured, and the strong and weak parallel scalability of the methods is demonstrated. The different relative weighting and significance of the algorithmic components on GPUs and CPUs is discussed. Results on problems involving adaptively refined nonconforming meshes are shown, and the use of the preconditioners on a large-scale magnetic diffusion problem using all spaces of the finite element de Rham complex is illustrated.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1983

A preconditioned formulation of the Cauchy-Riemann equations

A preconditioning of the Cauchy-Riemann equations which results in a second-order system is described. This system is shown to have a unique solution if the boundary conditions are chosen carefully. This choice of boundary condition enables the solution of the first-order system to be retrieved. A numerical solution of the preconditioned equations is obtained by the multigrid method.

Phillips, T. N.↗

NASA NTRS1986

Preconditioned methods for solving the incompressible and low speed compressible equations

Acceleration methods are presented for solving the steady state incompressible equations. These systems are preconditioned by introducing artificial time derivatives which allow for a faster convergence to the steady state. The compressible equations in conservation form with slow flow are also considered. Two arbitrary functions, alpha and beta, are introduced in the general preconditioning. An analysis of this system is presented and an optimal value for beta is determined given a constant, alpha. It is further shown that the resultant incompressible equations form a symmetric hyperbolic system and so are well posed. Several generalizations to the compressible equations are presented which generalize previous results.

Turkel, E.↗

NASA NTRS1987

Preconditioned methods for solving the incompressible and low speed compressible equations

The introduction of artificial time derivatives which permit faster convergence to the steady state is used to precondition steady state incompressible equation-solving acceleration methods. An optimal value is obtained for beta, given constant alpha, where alpha and beta are arbitrary functions in the general preconditioning stage. A symmetric hyperbolic system is formed by the resulting incompressible equations. Compressible equations in conservation form with slow flow are also considered.

Turkel, Eli↗

NASA NTRS1990

Fourier analysis of finite element preconditioned collocation schemes

The spectrum of the iteration operator of some finite element preconditioned Fourier collocation schemes is investigated. The first part of the paper analyses one-dimensional elliptic and hyperbolic model problems and the advection-diffusion equation. Analytical expressions of the eigenvalues are obtained with use of symbolic computation. The second part of the paper considers the set of one-dimensional differential equations resulting from Fourier analysis (in the tranverse direction) of the 2-D Stokes problem. All results agree with previous conclusions on the numerical efficiency of finite element preconditioning schemes.

Deville, Michel O.↗

NASA NTRS1989

Preconditioning matrices for the pseudospectral approximation of first-order operators

The behavior of the eigenvalues of preconditioning matrices for the pseudospectral approximation to the derivative operator has been analyzed in one and two dimensions. The one-dimensional analysis resulted in real and positive eigenvalues for the selected tridiagonal matrices. In the two-dimensional analysis, the eigenvalues of the selected block-diagonal matrices behaved well, but the preconditioner is full and therefore not suitable for applications. The Richardson scheme has been applied in the unpreconditioned as well as the preconditioned version to find the solution of the model problem.

Funaro, D.↗

NASA NTRS1990

Preconditioned upwind methods to solve 3-D incompressible Navier-Stokes equations for viscous flows

A computational method for calculating low-speed viscous flowfields is developed. The method uses the implicit upwind-relaxation finite-difference algorithm with a nonsingular eigensystem to solve the preconditioned, three-dimensional, incompressible Navier-Stokes equations in curvilinear coordinates. The technique of local time stepping is incorporated to accelerate the rate of convergence to a steady-state solution. An extensive study of optimizing the preconditioned system is carried out for two viscous flow problems. Computed results are compared with analytical solutions and experimental data.

Hsu, C.-H.↗

NASA NTRS1992

Preconditioning for stratospheric sudden warmings - Sensitivity studies with a numerical model

The phenomenon of preconditioning for sudden stratospheric warming events is investigated using a 3D primitive equation model adapted from the model used by Rose and Brasseur (1989). In order to identify the conditions that make the flow preconditioned, the model was initialized alternately with observed and with artificially constructed 'hybrid' stratospheric flow patterns. Two observed warmings were investigated, February 1979 and February 1989. The artificially constructed hybrid flow fields used upper-level winds (above 10 mb) from one period and lower-level winds (below 10 mb) from another. The results of model simulations indicate that the flow in the lower and middle stratosphere is capable of serving as a predictor of the potential for a major sudden warming, with the occurrence of a warming depending only on the lower-level winds.

Smith, Anne K.↗

NASA NTRS1992

Propulsion-related flowfields using the preconditioned Navier-Stokes equations

A previous time-derivative preconditioning procedure for solving the Navier-Stokes is extended to the chemical species equations. The scheme is implemented using both the implicit ADI and the explicit Runge-Kutta algorithms. A new definition for time-step is proposed to enable grid-independent convergence. Several examples of both reacting and non-reacting propulsion-related flowfields are considered. In all cases, convergence that is superior to conventional methods is demonstrated. Accuracy is verified using the example of a backward facing step. These results demonstrate that preconditioning can enhance the capability of density-based methods over a wide range of Mach and Reynolds numbers.

Venkateswaran, S.↗

NASA NTRS1992

Shape reanalysis and sensitivities utilizing preconditioned iterative boundary solvers

The computational advantages associated with the utilization of preconditined iterative equation solvers are quantified for the reanalysis of perturbed shapes using continuum structural boundary element analysis (BEA). Both single- and multi-zone three-dimensional problems are examined. Significant reductions in computer time are obtained by making use of previously computed solution vectors and preconditioners in subsequent analyses. The effectiveness of this technique is demonstrated for the computation of shape response sensitivities required in shape optimization. Computer times and accuracies achieved using the preconditioned iterative solvers are compared with those obtained via direct solvers and implicit differentiation of the boundary integral equations. It is concluded that this approach employing preconditioned iterative equation solvers in reanalysis and sensitivity analysis can be competitive with if not superior to those involving direct solvers.

Guru Prasad, K.↗

NASA NTRS1993

Preconditioned domain decomposition scheme for three-dimensional aerodynamic sensitivity analysis

A preconditioned domain decomposition scheme is introduced for the solution of the 3D aerodynamic sensitivity equation. This scheme uses the iterative GMRES procedure to solve the effective sensitivity equation of the boundary-interface cells in the sensitivity analysis domain-decomposition scheme. Excluding the dense matrices and the effect of cross terms between boundary-interfaces is found to produce an efficient preconditioning matrix.

Eleshaky, Mohammed E.↗

NASA NTRS1993

The multigrid preconditioned conjugate gradient method

A multigrid preconditioned conjugate gradient method (MGCG method), which uses the multigrid method as a preconditioner of the PCG method, is proposed. The multigrid method has inherent high parallelism and improves convergence of long wavelength components, which is important in iterative methods. By using this method as a preconditioner of the PCG method, an efficient method with high parallelism and fast convergence is obtained. First, it is considered a necessary condition of the multigrid preconditioner in order to satisfy requirements of a preconditioner of the PCG method. Next numerical experiments show a behavior of the MGCG method and that the MGCG method is superior to both the ICCG method and the multigrid method in point of fast convergence and high parallelism. This fast convergence is understood in terms of the eigenvalue analysis of the preconditioned matrix. From this observation of the multigrid preconditioner, it is realized that the MGCG method converges in very few iterations and the multigrid preconditioner is a desirable preconditioner of the conjugate gradient method.

Tatebe, Osamu↗