Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conservative schemes”

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.

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

At least 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Error field predictability and consequences for ITER

Abstract ITER coil tolerances are re-evaluated using the modern understanding of coupling to least-stable plasma modes and an updated center-line-traced model of ITER’s coil windings. This reassessment finds the tolerances to be conservative through a statistical, linear study of n = 1 error fields (EFs) due to tilted, shifted misplacements and nominal windings of central solenoid and poloidal field coils within tolerance. We also show that a model-based correction scheme remains effective even when metrology quality is sub-optimal, and compare this to projected empirical correction schemes. We begin with an analysis of the necessity of error field correction (EFC) for daily operation in ITER using scalign laws for the EF penetration threshold. We then consider the predictability of EF dominant mode overlap across early planned ITER scenarios and, as measuring EFs in high power scenarios can pose risks to the device, the potential for extrapolation to the ITER Baseline Scenario (IBS). We find that carefully designing a scenario matching currents proportionally to those of the IBS is far more important than plasma shape or profiles in accurately measuring an optimal correction current set.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1972

Interplanetary shock normals

A scheme has been developed to improve the estimate of an interplanetary shock normal by using both magnetic field and plasma data from a single spacecraft. Calculation of the basic shock model employs a subset of the eight magnetohydrodynamic conservation relations for a shock in an isotropic medium. This subset consists of six equations that are devoid of pressure and temperature terms. A sigma-weighted least squares loss function technique is used to best fit the overdetermination equations with respect to the eleven parameters of the system. This procedure usually yields a normal much more accurate than the one obtained using magnetic field values alone. An example is given of the use of the technique on the data for the Pioneer 7 shock of 29 August 1966.

Lepping, R. P.↗
NASA NTRS1975

A generalized hyperbolic marching method for chemically reacting 3-D supersonic flow using a splitting technique

A generalized hyperbolic marching method employing a nonorthogonal coordinate system and using a split differencing scheme for calculating steady supersonic flow over aerodynamic shapes is presented. It is a second-order-accurate mixed explicit-implicit procedure that solves the inviscid adiabatic and nondiffusive equations for chemically reacting flow in integral conservation-law form. The relationship between the integral and differential forms of the equations are examined and the relative merits of each discussed. The method admits initial Cauchy data situated on any arbitrary surface and integrates them forward along a general curvilinear coordinate, distorting and deforming the surface as it advances. The chemical kinetics term is split from the convective terms which are themselves dimensionally split, thereby freeing the fluid operators from the restricted step-size imposed by the chemical reactions and increasing the computational efficiency. The accuracy of this splitting technique is analyzed, a sufficient stability criterion is established, and comparisons are made with another method.

Rizzi, A. W.↗
NASA NTRS1984

Entropy condition satisfying approximations for the full potential equations of transonic flow

A class of conservative difference approximations for the steady full potential equation was presented. They are, in general, easier to program than the usual density biasing algorithms, and in fact, differ only slightly from them. Rigorous proof indicated that these new schemes satisfied a new discrete entropy inequality, which ruled out expansion shocks, and that they have sharp, steady, discrete shocks. A key tool in the analysis is the construction of a new entropy inequality for the full potential equation itself. Results of some numerical experiments using the new schemes are presented.

Osher, S.↗
NASA NTRS1985

Entropy condition satisfying approximations for the full potential equation of transonic flow

A class of conservative difference approximations for the steady full potential equation was presented. They are, in general, easier to program than the usual density biasing algorithms, and in fact, differ only slightly from them. Rigorous proof indicated that these new schemes satisfied a new discrete entropy inequality, which ruled out expansion shocks, and that they have sharp, steady, discrete shocks. A key tool in the analysis is the construction of a new entropy inequality for the full potential equation itself. Results of some numerical experiments using the new schemes are presented.

Osher, S.↗
NASA NTRS1992

On an origin of numerical diffusion: Violation of invariance under space-time inversion

The invariant properties of the convection equation du/dt + adu/dx = 0 (where d is the partial differential operator) with respect to spatial reflection, time reversal, and space-time inversion are studied. Generally, a finite-difference analog of this equation may possess some or none of these properties. It is shown that, under certain conditions, the von Neumann amplification factor of an analog satisfies a special relation for each invariant property this analog possesses. Particularly, an analog is neutrally stable and thus free of numerical diffusion if it possesses the invariant property related to space-time inversion. It is also explained why generally (1) an upwind scheme possesses neither the invariant property related to spatial reflection nor that related to space-time inversion, and (2) an explicit scheme possesses neither the invariant property related to time reversal nor that related to space-time inversion. Extension to the viscous case and a remarkable connection between the current work and a new numerical framework for solving conservation laws are also discussed.

Chang, Sin-Chung↗
NASA NTRS1976

Split Space-Marching Finite-Volume Method for Chemically Reacting Supersonic Flow

A space-marching finite-volume method employing a nonorthogonal coordinate system and using a split differencing scheme for calculating steady supersonic flow over aerodynamic shapes is presented. It is a second-order-accurate mixed explicit-implicit procedure that solves the inviscid adiabatic and nondiffusive equations for chemically reacting flow in integral conservation-law form. The relationship between the finite-volume and differential forms of the equations is examined and the relative merits of each discussed. The method admits initial Cauchy data situated on any arbitrary surface and integrates them forward along a general curvilinear coordinate, distorting and deforming the surface as it advances. The chemical kinetics term is split from the convective terms which are themselves dimensionally split, thereby freeing the fluid operators from the restricted step size imposed by the chemical reactions and increasing the computational efficiency. The accuracy of this splitting technique is analyzed, a sufficient stability criterion is established, a representative flow computation is discussed, and some comparisons are made with another method.

Rizzi, Arthur W.↗
DOE OSTI2020

A fully implicit, scalable, conservative nonlinear relativistic Fokker–Planck 0D-2P solver for runaway electrons

Upon application of a sufficiently strong electric field, electrons break away from thermal equilibrium and approach relativistic speeds. These highly energetic ‘runaway’ electrons (~ MeV) play a significant role in tokamak disruption physics, and therefore their accurate understanding is essential to develop reliable mitigation strategies. As such, we have developed a fully implicit solver for the 0D-2P (i.e., including two momenta coordinates) relativistic nonlinear Fokker–Planck equation (rFP). As in earlier implicit rFP studies (NORSE, CQL3D), electron–ion interactions are modeled using the Lorentz operator, and synchrotron damping using the Abraham–Lorentz–Dirac reaction term. However, our implementation improves on these earlier studies by (1) ensuring exact conservation properties for electron collisions, (2) strictly preserving positivity, and (3) being scalable algorithmically and in parallel. Key to our proposed approach is an efficient multigrid preconditioner for the linearized rFP equation, a multigrid elliptic solver for the Braams–Karney potentials, and a novel adaptive technique to determine the associated boundary values. We verify the accuracy and efficiency of the proposed scheme with numerical results ranging from small electric-field electrical conductivity measurements to the accurate reproduction of runaway tail dynamics when strong electric fields are applied.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1988

A solution method for the unsteady incompressible Navier-Stokes equations in generalized coordinate systems

A solution method based on a fractional step approach is developed for obtaining time-dependent solutions of the three-dimensional, incompressible Navier-Stokes equations in generalized coordinate systems. The governing equations are discretized conservatively by finite volumes using a staggered mesh system. The primitive variable formulation uses the volume fluxes across the faces of each computational cell as dependent variables. This procedure, combined with accurate and consistent approximations of geometric parameters, is done to satisfy the discretized mass conservation equation to machine accuracy as well as to gain favorable convergence properties of the Poisson solver. The discretized equations are second-order-accurate in time and space and no smoothing terms are added. An approximate-factorization scheme is implemented in solving the momentum equations. A novel ZEBRA scheme with four-color ordering is devised for the efficient solution of the Poisson equation. Several two and three-dimensional solutions are compared with other numerical and experimental results to validate the present method.

Rosenfeld, Moshe↗
NASA NTRS1975

Three-dimensional finite element analysis for high velocity impact

A finite element algorithm for solving unsteady, three-dimensional high velocity impact problems is presented. A computer program was developed based on the Eulerian hydroelasto-viscoplastic formulation and the utilization of the theorem of weak solutions. The equations solved consist of conservation of mass, momentum, and energy, equation of state, and appropriate constitutive equations. The solution technique is a time-dependent finite element analysis utilizing three-dimensional isoparametric elements, in conjunction with a generalized two-step time integration scheme. The developed code was demonstrated by solving one-dimensional as well as three-dimensional impact problems for both the inviscid hydrodynamic model and the hydroelasto-viscoplastic model.

Chan, S. T. K.↗
NASA NTRS1982

An efficient full potential implicit method based on characteristics for analysis of supersonic flows

A nonlinear aerodynamic prediction technique based on the full potential equation in conservation form has been developed for the treatment of supersonic flows. The method uses the theory of characteristic signal propagation to accurately simulate the flow structure, which includes shock waves and mixed elliptic-hyperbolic crossflow. An implicit approximate factorization scheme is employed to solve the finite-differenced equation. The necessary body-fitted grid system in every marching plane is generated numerically, using an elliptic grid solver. Results are shown for conical and nonconical wing-body combinations and compared with experimental data and Euler calculations. The method demonstrates an enormous savings in execution time and memory requirements over Euler methods.

Shankar, V.↗
NASA NTRS1985

Finite element algorithms for compressible flow computation on a supercomputer

Methods of applying computational fluid dynamics techniques to engineering problems are discussed. The linking factors between theoretical and applied research are the shape and weight functions which produce the Galerkin finite element schemes. The flow is unsteady, viscous, heat conducting and compressible, and steady-state flows are the asymptotic limit of unsteady flows. All flows are described by the time-averaged Navier-Stokes equations (NSE) with mass and energy conservation. Local curvilinear intrinsic coordinates are applied to discretize the NSE in arbitrary geometric domains. A modified weighted residuals approach defines the discrete analogs of the physical systems modeled and a general interpolants method is used to derive families of numerical models, both implicit and explicit, finite difference and finite elemnt, from a single point of departure. Time-dependent element approximations are achieved with a progressive assembly of generalized elements method. Hyperbolic steady-state Euler algorithms and quasi-parabolic are spatial marching algorithms for solving the discretized equations. Application of the techniques is illustrated through calculation of the compressible flow around a hypersonic flight vehicle at 100,000 ft at zero angle of attack at speeds of Mach 4-10.

Spradley, L. W.↗
NASA NTRS1994

A mapping method for the gravitational few-body problem with dissipation

Recently a new class of numerical integration methods - 'mixed variable symplectic integrators' - has been introduced for studying long-term evolution in the conservative gravitational few-body problem. These integrators are an order of magnitude faster than conventional ordinary differential equations (ODE) integration methods. Here we present a simple modification of this method to include small non-gravitational forces. The new scheme provides a similar advantage of computational speed for a larger class of problems in Solar System dynamics.

Malhotra, Renu↗
NASA NTRS1994

Nonlinear evolution of the coronal magnetic field under reconnective relaxation

Recently, Vekstein et al. (Vekstein, Priest, & Steele 1993) have developed a model for coronal heating in which the corona responds to photospheric footpoint motions by small-scale reconnection events that bring about a relaxed state while conserving magnetic helicity but not field-line connectivity. Vekstein et al. consider a partially open field configuration in which magnetic helicity is ejected to infinity on open field lines but retained in the closed-field region. Under this scheme, they describe the evolution of an initially potential field, in response to helicity injection, in the linear regime. The present work uses numerical calculations to extend the model of Vekstein et al. into the fully nonlinear regime. The results show a rise and bulging of the field lines of the closed-field region with increasing magnetic helicity, to a point where further solutions are impossible. We interpret these solution-sequence endpoints as indicating a possible loss of equilibrium, in the sense that a relaxed equilibrium state may no longer be available to the corona when sufficient helicity has been injected. The rise and bulging behavior is reminiscent of what is observed in a helmet streamer just before the start of a coronal mass ejection (CME), and so our model suggests that a catastrophic loss of magnetic equilibrium might be the initiation mechanism for CMEs. We also find that some choices of boundary conditions can result in qualitative changes in the magnetic topology, with the appearance of magnetic islands. Whether or not this behavior occurs depends on the relative strengths of the fields in the closed- and open-field regions; in particular, island formation is most likely when the open field (which is potential) is strong and thus acts to confine the force-free closed field. Finally, we show that the energy released through reconnective relaxation can be a substantial fraction of the magnetic energy injected into the corona through footpoint motions and may be sufficient for heating the corona above active regions.

Wolfson, R.↗
DOE OSTI2020

Computational model for predicting the performance of fluid-powered longiudinal piston modules and motors.

A Matlab-based computational model line code has been developed for predicting the performance of fluid-powered longitudinal piston modules and motors. A system of equations has been developed based upon applying the conservation of mass to a control volume defined by a reciprocating piston. The method results in a differential equation that predicts the pressure gradient in the control volume as a function of rotor rotation angle. The model uses an integration scheme to predict the pressure state in the control volume. The orifice equation is used to predict the flowrate in and out of the control volume. From the resulting solution the torque and speed characteristics of the motor configuration can be predicted, optimized, or tailored for a specific requirement.SAND2019-13983 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Taylor, Michael↗
NASA NTRS1981

An implicit marching procedure for the treatment of supersonic flow fields using the conservative full potential equation

A new aerodynamic prediction technique based on the conservative full potential equation is developed for the treatment of supersonic flow fields. This new technique bridges the gap between simplistic linear theory methods and complex Euler solvers. A novel local density linearization concept and a second order accurate retarded density scheme, both producing the correct artificial viscosity, are introduced in developing an implicit marching scheme for solving the scalar phi. The method produces results that compare well with Euler solvers and requires an order of magnitude less computer time and significantly less computer memory over existing nonlinear codes. The scalar phi formulation can be extended to handle subsonic pockets in the marching direction and also is suitable for developing inverse procedures where the shape corresponding to a prescribed loading is sought.

Shankar, V.↗
NASA NTRS1981

A PANSONIC Navier-Stokes solver

A finite-difference formulation of the full Navier-Stokes equations which demonstrates a capability to economically solve two-dimensional problems has been developed. The basic algorithm was derived from the full, Reynolds-averaged, conservative, Navier-Stokes equations expressed in curvilinear coordinates. Eddy viscosity was determined by the Baldwin and Lomax algebraic turbulence model. This non-iterative, second-order accurate, implicit, numerical algorithm is based on the approximate factorization finite-difference scheme of Beam and Warming. Results indicate a facility for solving subsonic, transonic, and supersonic (hence PANSONIC) flows about arbitrary airfoils for a wide range of Reynolds numbers, Mach numbers, and angles of attack. Current computations demonstrate that vectorized implementations of this algorithm can solve steady-state, two-dimensional problems in five to ten minutes of computer time.

Cooper, G. K.↗