Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical integration”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Implementation of a High-Speed Multichannel Data Acquisition System for Magnetic Diagnostics and Plasma Centroid Position Control in ISTTOK

In tokamak and other fusion devices, magnetic control is the main tool that allows to regulate the plasma current, position and shape; it is in charge of actuating the desired plasma current waveform, steering the plasma position to a given set point and maintain the plasma shape close to a prescribed plasma equilibrium. This work describes the application of several physics concepts and computational tools in order to obtain a novel optimal controller for the plasma centroid position, which has been implemented and tested in the real-time plasma control system at the ISTTOK tokamak. A key point for the development of the new control system was the installation of a recently upgraded hardware, that numerically integrates in real-time the magnetic probes signals.

poloidal field coils↗
DOE OSTI2025

Gaussian FLOWERS: Wind-rose-based analytical integration of Gaussian wake model for extremely fast AEP estimation

A major cost in the study of wind farm layout optimization is the repeated evaluation of the annual energy production (AEP). The current approach to estimating AEP requires a large set of flow simulations to be performed that cover each discrete wind speed and direction combination contained within the wind rose, followed by a probability-weighted sum of the power production resulting from each simulation. Even with inexpensive engineering wake models, this numerical integration scheme can lead to high computational costs. In this paper, we derive an analytical formulation for estimating farm AEP across every wind direction, based on a Gaussian wake velocity model, which reduces the number of wind farm simulations to a single function evaluation. As a result, we find that the Gaussian-FLOWERS approach reduces the time for AEP calculations by more than two orders of magnitude with a small trade-off in accuracy when compared to a conventional approach. This massive reduction in computation cost is useful to reduce overall costs in wind farm layout optimization studies.

17 WIND ENERGY↗
NASA NTRS1984

A numerical transport scheme which avoids negative mixing ratios

A new scheme for numerically integrating the transport equation has been developed. This square root method avoids the problem of negative mixing ratios by using the square root of the concentration instead of the concentration itself as an advective time step. Conservation of total mass at every time step is achieved by extending the concept of quadratic conservation to the numerical time integration method. A few time schemes that fulfill the requirement of step-by-step quadratic conservation are described, the simplest of which is the modified Lax-Wendroff method.

Schneider, H.-R.↗
DOE OSTI2023

A conservative Galerkin solver for the quasilinear diffusion model in magnetized plasmas

We propose a conservative Galerkin scheme for the quasilinear model in three-dimensional momentum space and three-dimensional spectral space, with cylindrical symmetry. We construct an unconditionally conservative weak form and use a discretization that preserves conservation properties independent of the wave emission probability. The discrete operators, combined with a consistent quadrature rule, preserve all the conservation laws rigorously. The proposed scheme is quite general: it works for both relativistic and non-relativistic systems, for both magnetized and unmagnetized plasmas, and even for problems with time-dependent dispersion relations. We represent the particle distribution by continuous basis functions and use discontinuous basis functions for the wave spectral energy density, which enables the application of a positivity-preserving technique. We adopt the marching simplex algorithm, designed initially for computer graphics, for numerical integration on the resonance manifold. Furthermore, the numerical examples with a bump-on-tail initial configuration show how the unstable waves produce strong momentum space diffusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1977

Unsteady supersonic aerodynamic theory for interfering surfaces by the method of potential gradient

A generalized solution of the hyperbolic wave equation was further developed to relate the velocity components at a field point to the potential gradient distribution in the dependence domain. Singular integrals were evaluated in closed form, with numerical integration methods for more complex but analytic functions. Idealization of the lifting surfaces by trapezoidal elements with two sides parallel to the streamlines is computationally efficient. Streamwise integrals were performed analytically, and spanwise integrals were neccessary only on element leading and trailing sides. All integrands vanish on the Mach cone. Pressure distribution on a double delta wing and generalized aerodynamic coefficients for three AGARD planforms were calculated and compared with available results.

Jones, W. P.↗
NASA NTRS1991

Effects of anisotropic conduction and heat pipe interaction on minimum mass space radiators

Equations are formulated for the two dimensional, anisotropic conduction of heat in space radiator fins. The transverse temperature field was obtained by the integral method, and the axial field by numerical integration. A shape factor, defined for the axial boundary condition, simplifies the analysis and renders the results applicable to general heat pipe/conduction fin interface designs. The thermal results are summarized in terms of the fin efficiency, a radiation/axial conductance number, and a transverse conductance surface Biot number. These relations, together with those for mass distribution between fins and heat pipes, were used in predicting the minimum radiator mass for fixed thermal properties and fin efficiency. This mass is found to decrease monotonically with increasing fin conductivity. Sensitivities of the minimum mass designs to the problem parameters are determined.

Baker, Karl W.↗
NASA NTRS1994

Algorithms For Integrating Nonlinear Differential Equations

Improved algorithms developed for use in numerical integration of systems of nonhomogenous, nonlinear, first-order, ordinary differential equations. In comparison with integration algorithms, these algorithms offer greater stability and accuracy. Several asymptotically correct, thereby enabling retention of stability and accuracy when large increments of independent variable used. Accuracies attainable demonstrated by applying them to systems of nonlinear, first-order, differential equations that arise in study of viscoplastic behavior, spread of acquired immune-deficiency syndrome (AIDS) virus and predator/prey populations.

Freed, A. D.↗
NASA NTRS1982

A numerical method for integrating the kinetic equations of droplet spectra evolution by condensation/evaporation and by coalescence/breakup processes

An extension of the method of moments is developed for the numerical integration of the kinetic equations of droplet spectra evolution by condensation/evaporation and by coalescence/breakup processes. The number density function n sub k (x,t) in each separate droplet packet between droplet mass grid points (x sub k, x sub k+1) is represented by an expansion in orthogonal polynomials with a given weighting function. In this way droplet number concentrations, liquid water contents and other moments in each droplet packet are conserved and the problem of solving the kinetic equations is replaced by one of solving a set of coupled differential equations for the number density function moments. The method is tested against analytic solutions of the corresponding kinetic equations. Numerical results are obtained for different coalescence/breakup and condensation/evaporation kernels and for different initial droplet spectra. Also droplet mass grid intervals, weighting functions, and time steps are varied.

Emukashvily, I. M.↗
NASA NTRS1966

Numerical Studies of the Non-linear Vlasov Equation

The subject of this dissertation is the numerical integration of the initial-value problem for the non-linear Vlasov equation. The Vlasov equation is used to describe the dynamics of a "collisionless", one-dimensional, classical electron gas confined between two perfectly reflecting boundaries. Only the long-range Coulomb interactions of the electrons are considered; effects associated with the discrete structure are neglected. The numerical results obtained for non-linear Landau damping compare well with similar results obtained by Knorr. A general statement of the results on stable initial conditions is: As the degree of non-linearity of the initial conditions is increased, the deviation from linear Landau damping appears sooner and is more severe. In some cases damping was observed to cease. Curves showing the time dependence of the damping decrement are derived and compared with predictions of non-linear theories. New results obtained in this study include the observation that for strongly non-linear cases, the damping of the electric field causes an initially Maxwellian fo (v, O) to develop a peak in the neighborhood of the phase velocity; strong growth of the second harmonic is seen after fo (v, t) develops such a peak. Also new in this study is the interpretation of the development of a certain class of strongly unstable initial conditions as approaching an inhomogeneous equilibrium.

NUMERICAL INTEGRATION↗
DOE OSTI2025

Single-shot femtosecond laser-induced breakdown spectroscopy in low-density gases

Single-shot fs laser-induced breakdown spectroscopy (LIBS) has been used to interrogate low-density gases of H 2 ⁠, CH 4 ⁠, CO 2 ⁠, and H 2 O ranging from 0.5 × 10 15 to 9.1 × 10 15 cm −3 (15–275 mTorr). Strong spectral signatures were observed above the noise baseline for each test gas at all densities considered. Unique molecular signatures were observed for both CH 4 and CO 2 along with atomic features for each gas. The LIBS intensity response with respect to density was determined for each species using Voigt lineshape fitting or numerical integration. Finally, the collection and analysis of LIBS spectra of composite mixtures of CH 4 –CO 2 and CH 4 –H 2 demonstrates that the extraction of individual densities from gas mixtures is feasible at low densities.

Artificial neural networks↗
NASA NTRS1975

Integrated potential formulation of unsteady supersonic aerodynamics for interacting wings

A numerical integrated velocity potential method for the determination of unsteady aerodynamic forces on arbitrary interacting wings and tails in supersonic flow has been developed. Constant Mach number has been assumed throughout the flow field. Normalwash and sidewash integrals have been derived. The upwash integral remains to be derived. Singular integrals in the expressions for the velocity components have been evaluated in closed form. Lifting surfaces are represented by triangular elements defined by arbitrarily spaced characteristic lines and the true surface edges. The wake field is represented by rectangular strip elements. Velocity potential distributions and generalized aerodynamic coefficients have been compared.

Appa, K.↗
NASA NTRS1979

Shuttle program. Onorbit navigation integrator results for typical shuttle orbits

Three types of navigation onorbit numerical integrators were evaluated: (1) power integrators with no delta-V incorporation, just coasting (using Taylor series expansion integrators); (2) coasting integrators using the Cowell method of special perturbations; and (3) coasting integrator using the Pines variation of parameter perturbation method. Results show that the super G integrator is a very simple and effective for 2 and 4 second time steps. Since IMU delta-V data can be easily incorporated in the integration scheme, its use as the standard onorbit navigation propagator for the maintenance of the current state was implemented in the onboard navigation software. The Pines formulation method with a Runge-Kutta-Gill fourth-order integrator method produces excellent results up to 300 second time steps. On orbit prediction with this method was implemented in the onboard onorbit navigation scheme. The Runge-Kutta third order, using Cowell's method, is an excellent general purpose determination integrator for time steps up to a 60 second duration.

Olszewski, O. W.↗
DOE OSTI2020

Quantum Monte Carlo benchmarking of large noncovalent complexes in the L7 benchmark set

We have used diffusion Monte Carlo (DMC) to perform calculations on the L7 benchmark set. DMC is a stochastic numerical integration scheme in real-space and part of a larger set of quantum Monte Carlo methods. The L7 set was designed to test the ability of electronic structure methods to include dispersive interactions. While the agreement between DMC and quantum-chemical state-of-the-art methods is excellent for some of the structures, there are significant differences in others. In contrast to wavefunction-based quantum chemical methods, DMC is a first-principle many-body method with the many-body wavefunction evolving in real space. It includes explicitly all electron–electron interactions and is relatively insensitive to the size of the basis set.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1974

Nonlinearities in analyses of unsteady flow around oscillating wings

An analysis of nonlinearities involved in unsteady subsonic viscous flow around wings at high Reynolds numbers is presented. The finite element method was applied for the numerical integration of the vorticity transport equations. Physical and numerical aspects of the nonlinearities are discussed. Magnification of nonlinearities in the convection terms near wing surfaces due to changing geometry at leading and trailing edges of the wing chord is determined. As a bubble is generated near the leading edge, the timewise variation of vorticity and pressure distributions over the surface is shown. The different procedures for the integration of the vorticity transport equations in three dimensions along spanwise, chordwise, and normal directions are discussed.

Bratanow, T.↗
DOE OSTI2025

A three-dimensional laser ray-tracing methodology for radiation-hydrodynamics simulations

We report on a methodology for performing laser ray-tracing in three spatial dimensions for radiation-hydrodynamics simulation codes. Our method, which is an extension of that developed in Haines et al., Comput. Fluids 201, 104478 (2020), utilizes an automatically generated separate mesh for the laser ray-tracing from the radiation-hydrodynamics mesh. This enables the laser mesh to be tailored to minimize ray noise with significantly fewer rays than would be required when the ray-tracing is performed on the radiation-hydrodynamics mesh, primarily by allowing the use of high-aspect-ratio cells that are not suitable for hydrodynamics solvers. For a planar target, we show that our method provides a ≈ 100× reduction in computational expense to achieve a fixed level of ray noise relative to ray-tracing directly on the radiation-hydrodynamics mesh. The relatively low ray requirement also enables efficient computation of cross-beam energy transfer. Each cell in the logically cubic laser mesh is a non-convex dodecahedron with triangular sides, and numerical integration of the ray trajectories and inverse bremsstrahlung is performed by mapping each cell to the unit cube. We will describe our methodology in detail as well as its implementation in the xRAGE radiation-hydrodynamics code, discuss performance, and present the results from applying the methodology to test problems with analytic solutions for laser ray-tracing through a quadratic density gradient with an analytic solution as well as for a laser-driven heat front. In 3D radiation-hydrodynamics simulations of laser-driven experiments performed on the National Ignition Facility, laser ray-tracing with our methodology uses less than 1% of total computational time while introducing acceptably low levels of ray noise.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A User Guide to PARET/ANL

PARET was originally created in 1969 at what is now Idaho National Laboratory (INL), to analyze reactivity insertion events in research and test reactor cores cooled by light or heavy water, with fuel composed of either plates or pins. The use of PARET is also appropriate for fuel assemblies with curved fuel plates when their radii of curvatures are large with respect to the fuel plate thickness. The PARET/ANL version of the code has been developed at Argonne National Laboratory (ANL) under the sponsorship of the U.S. Department of Energy/NNSA since the inception of the Reactor Conversion Program. Since Reduced Enrichment for Research and Test Reactors (RERTR) began in 1978, PARET/ANL has been benchmarked to experimental data including SPERT testing, and used to determine the expected transient behavior of a large number of reactors regardless of enrichment both inside and outside the Reactor Conversion Program. This document provides the pertinent information for the use of PARET/ANL Version 7.6. PARET/ANL models the various fueled regions of a reactor core as channels. Each of these channels consists of a single flat fuel plate/pin (including cladding and, optionally, a gap) with water coolant on each side. In slab geometry the coolant channels for a given fuel plate are of identical dimensions (mirror symmetry), but they can be of different thickness in each channel. There can be many channels, but each channel is independent and coupled only through reactivity feedback effects to the whole core. The time-dependent differential equations that represent the system are replaced by an equivalent set of finite-difference equations in space and time, which are integrated numerically. PARET/ANL uses fundamentally the same numerical scheme as RELAP5 for the time-integration of the point-kinetics equations. The one-dimensional thermal-hydraulic model includes temperature-dependent thermal properties of the solid materials, such as heat capacity and thermal conductivity, as well as the transient heat production and heat transfer from the fuel meat to the coolant. Temperature- and pressure-dependent thermal properties of the coolant such as enthalpy, density, thermal conductivity, and viscosity are also used in determining parameters such as friction factors and heat transfer coefficients. The code first determines the steady-state solution for the initial state. Then the solution of the transient is obtained by integration in time and space. Multiple heat transfer, DNB and flow instability correlations are available. The code was originally developed to model reactors cooled by an open loop, which was adequate for rapid transients in pool-type cores. An external loop model appropriate for Miniature Neutron Source Reactors (MNSR’s) was also added to PARET/ANL to model natural circulation within the vessel, heat transfer from the vessel to pool and heat loss by evaporation from the pool. PARET/ANL also contains models for decay heat after shutdown, control rod reactivity versus time or position, time-dependent pump flow, and loss-of-flow event with flow reversal as well as logic for trips on period, power, and flow. Feedback reactivity effects from coolant density changes and temperature changes are represented by tables. Feedback reactivity from fuel heat-up (Doppler Effect) is represented by a four-term polynomial in powers of fuel temperature. Photo-neutrons produced in beryllium or in heavy water may be included in the point-kinetics equations by using additional delayed neutron groups.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1975

A three-dimensional viscoelastic model of a strike slip fault

An analytic approximation to the Green's function for the displacements due to a strike slip point source in an elastic layer over a viscoelastic half-space is developed. This approximate Green's function is useful because it can be analytically integrated over the fault surface. Comparison with a numerical integration of the exact solution integral indicates that the approximation is quite good. The approximate Green's function is integrated analytically to obtain the displacements due to a finite rectangular strike slip fault in an elastic layer over a viscoelastic half-space. Ground displacements and angle changes from a model survey net are computed to illustrate the viscoelastic relaxation which follows a fracture in the elastic region.

Rundle, J. B.↗