Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical Methods”

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 73 records · Page 4

Numerical method for determination of resonance electron velocities in periodic electric fields

A Monte Carlo method is applied for modeling of a plasma source with a periodic structure of a radiofrequency electric field. The method is modified in order to detect resonance increase of the electron energy at given velocities of the electrons as expected from theoretical considerations. The numerical model is validated by comparison with known analytical results. The applicability of the method for detection of resonance velocities in plasma sources with electric field variation parallel or perpendicular to the electric field vector is proven.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Mathematical Models and Numerical Methods for High-Fidelity Simulation of Ignition of Reactive Mixtures by Nanosecond Plasma Discharges in Realistic Configurations

We present a newly developed framework for the numerical simulation of ignition of reactive mixtures using single or repeated nanosecond discharge pulses. The framework builds upon the AMReX library, using the existing compressible solver PeleC and low-Mach solver PeleLMeX and allowing for adaptive mesh refinement, complex geometries, and execution on next-generation high-performance computing (HPC) systems. High-fidelity elementary models are adopted for weakly-ionised plasma discharges with significant energy deposition, consistent with nanosecond discharge pulses, and then implemented in the solver. The treatment of non-thermal electrons and charged species, thermodynamics of non-equilbrium species, plasma kinetics, limiting time scales, and boundary conditions for charged species are discussed and addressed for computational efficiency. The framework is demonstrated for three relevant applications: single and multi-pulse discharges in air, single pulse ignition of an ethylene/air mixture, and a three-dimensional plasma discharge in air with temperature stratification. The successful application of the framework demonstrates the feasibility of high-fidelity simulation of ignition of air/hydrocarbon mixtures in three-dimensions with multiple discharge pulses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Determining the Minimum Phase‐to‐Phase Gap Distance for Unconventional Transmission Lines Using Numerical Method

With the increasing demand for power transmission, compact, high surge impedance loading (HSIL) high-voltage transmission lines have emerged as a viable solution due to their reduced land acquisition costs and higher power delivery capability. The compactness of a transmission line depends on effective insulation coordination, particularly in determining the phase-to-phase clearance, which is governed by the critical flashover voltage under switching and lightning overvoltage conditions. Traditional methods for phase-to-phase clearance rely on empirical formulas derived from experimental data, which are convenient for conventional high-voltage lines. However, unconventional HSIL lines require a faster and more adaptable evaluation method, as they involve optimized conductor positioning to reduce right-of-way requirements while enhancing natural power loadability. This study presents a simplified numerical approach to determine the minimum phase-to-phase gap, utilizing arc propagation viability curves, and offers an efficient alternative to conventional empirical methods. The proposed method was successfully applied to a 500 kV conventional line as well as 500 and 735 kV unconventional line designs, demonstrating its capability in accurately assessing insulation requirements. Results reveal that the method can support reduced gap clearances while still maintaining reliability, thereby validating its usefulness in optimizing compact transmission line configurations.

Arafat, Easir [The University of Texas at Dallas, ↗

Efficient numerical methods to solve sparse linear equations with application to PageRank

Over the last two decades, the PageRank problem has received increased interest from the academic community as an efficient tool to estimate web-page importance in information retrieval. Despite numerous developments, the design of efficient optimization algorithms for the PageRank problem is still a challenge. Here, we propose three new algorithms with a linear time complexity for solving the problem over a bounded-degree graph. The idea behind them is to set up the PageRank as a convex minimization problem over a unit simplex, and then solve it using iterative methods with small iteration complexity. Our theoretical results are supported by an extensive empirical justification using real-world and simulated data.

97 MATHEMATICS AND COMPUTING↗

Numerical Method for Fractional Advection-Diffusion Equation with Heredity

We propose a method of construction of difference schemes for fractional partial differential equations with delay in time. For the fractional equation with two-sided diffusion, fractional transfer in time, and a functional aftereffect, we construct an implicit difference scheme. We use the shifted Grünwald–Letnikov formulas for the approximation of fractional derivatives with respect to spatial variables and the L1-algorithm for the approximation of fractional derivatives in time. Also we use piecewise constant interpolation and extrapolation by extending the discrete prehistory of the model in time. The algorithm is a fractional analog of a purely implicit method; on each time step, it is reduced to the solution of linear algebraic systems. We prove the stability of the method and find its order of convergence.

97 MATHEMATICS AND COMPUTING↗

Investigation of Numerical Methods for Performance Improvement of MOOSE-based System Analysis Codes

The main objective of this study is to investigate the feasibility of implementation of staggered-grid finite volume method (SG-FVM) using the MOOSE framework to support the development of the advanced system analysis code SAM. This study successfully demonstrated the integration of staggered-grid finite volume method in the application level using the MOOSE framework, although not in the framework level which should be investigated in the future. Several important properties of the implemented SG-FVM, e.g., high-order spatial accuracy and monotonicity preserving, have been demonstrated with selected numerical test cases. The superior performance in execution time was also evident based on a selected one- dimensional flow problem in a loop configuration.

97 MATHEMATICS AND COMPUTING↗

A numerical method for deriving shape functions of nanoparticles for pair distribution function refinements

We report that in the structural refinement of nanoparticles, discrete atomistic modeling can be used for small nanocrystals (< 15 nm), but becomes computationally unfeasible at larger sizes, where instead unit-cell-based small-box modeling is usually employed. However, the effect of the nanocrystal's shape is often ignored or accounted for with a spherical model regardless of the actual shape due to the complexities of solving and implementing accurate shape effects. Recent advancements have provided a way to determine the shape function directly from a pair distribution function calculated from a discrete atomistic model of any given shape, including both regular polyhedra (e.g. cubes, spheres, octahedra) and anisotropic shapes (e.g. rods, discs, ellipsoids) [Olds et al. (2015). J. Appl. Cryst. 48, 1651–1659], although this approach is still limited to small size regimes due to computational demands. In order to accurately account for the effects of nanoparticle size and shape in small-box refinements, a numerical or analytical description is needed. This article presents a methodology to derive numerical approximations of nanoparticle shape functions by fitting to a training set of known shape functions; the numerical approximations can then be employed on larger sizes yielding a more accurate and physically meaningful refined nanoparticle size. Lastly, the method is demonstrated on a series of simulated and real data sets, and a table of pre-calculated shape function expressions for a selection of common shapes is provided.

36 MATERIALS SCIENCE↗

A Numerical Method for Determining the Elements of Circumbinary Orbits and Its Application to Circumbinary Planets and the Satellites of Pluto-Charon

Planets and satellites orbiting a binary system exist in the solar system and extrasolar planetary systems. Their orbits can be significantly different from Keplerian orbits, if they are close to the binary and the secondary-to-primary mass ratio is high. A proper description of a circumbinary orbit is in terms of the free eccentricity e {sub free} at the epicyclic frequency κ {sub 0}, forced eccentricity e {sub forced} at the mean motion n {sub 0}, and oscillations at higher frequencies forced by the non-axisymmetric components of the binary’s potential. We show that accurate numerical values for the amplitudes and frequencies of these terms can be extracted from numerical orbit integrations by applying fast Fourier transformation (FFT) to the cylindrical distance between the circumbinary object and the center of mass of the binary as a function of time. We apply this method to three Kepler circumbinary planets and the satellites of Pluto-Charon. For the satellite Styx of Pluto-Charon, the FFT results for κ {sub 0} and e {sub free} differ significantly from the first-order analytic value and the value reported by Showalter & Hamilton, respectively. We show that the deviation in κ {sub 0} is likely due to the effect of the 3:1 mean-motion resonance and discuss the implications of the lower value for e {sub free}.

79 ASTRONOMY AND ASTROPHYSICS↗

Flame Characterization of a Hydrogen Micromix Burner using Experimental and Numerical Methods - Extended Abstract

Micromix combustion presents a promising solution for integrating hydrogen into gas turbine combustors, offering rapid and compact air-fuel mixing that effectively shortens the combustion zone and significantly reduces nitrogen oxides (NOX) emissions. A major challenge in hydrogen combustion systems is flashback, which is mitigated through short mixing lengths and high flow velocities. This study investigates three fuel injector designs – zero swirl, low swirl, and high swirl – for a 6-mm inner diameter nozzle, where hydrogen is injected through six 0.3 mm diameter holes in a jet-in-crossflow configuration relative to the airstream. Numerical simulations indicate stable flame formation, moderate peak flame temperatures, and subsequent low NOx potential for all designs. These injectors have been fabricated using additive manufacturing, and upcoming experiments will explore flame structure and NOx behavior across varying fuel injection angles, equivalence ratios and pre-mixer velocities. These efforts aim to advance the development of low-emission hydrogen micromix combustion.

hydrogen↗

A parallel adaptive numerical method with generalized curvilinear coordinate transformation for compressible Navier–Stokes equations

A fourth-order finite-volume method for solving the Navier–Stokes equations on a mapped grid with adaptive mesh refinement is proposed, implemented, and demonstrated for the prediction of unsteady compressible viscous flows. The method employs fourth-order quadrature rules for evaluating face-averaged fluxes. Our approach is freestream preserving, guaranteed by the way of computing the averages of the metric terms on the faces of cells. The standard Runge–Kutta marching method is used for time discretization. Solutions of a smooth flow are obtained in order to verify that the method is formally fourth-order accurate when applying the nonlinear viscous operators on mapped grids. Solutions of a shock tube problem are obtained to demonstrate the effectiveness of adaptive mesh refinement in resolving discontinuities. A Mach reflection problem is solved to demonstrate the mapped algorithm on a non-rectangular physical domain. The simulation is compared against experimental results. Future work will consider mapped multiblock grids for practical engineering geometries.

97 MATHEMATICS AND COMPUTING↗

A numerical method for simulating variable density flows in membrane desalination systems

Here, we present a novel method for simulating unsteady, variable density, fluid flows in membrane desalination systems. By assuming the density varies only with concentration and temperature, the scheme decouples the solution of the governing equations into two sequential blocks. The first solves the governing equations for the temperature and concentration fields, which are used to compute all thermophysical properties. The second block solves the conservation of mass and momentum equations for the velocity and pressure. We show that this is computationally more efficient than schemes that iterate over the full coupled equations in one block. We verify that the method achieves second-order spatial-temporal accuracy, and we use the method to investigate buoyancy-driven convection in a desalination process called vacuum membrane distillation. Specifically, we show that with gravity properly oriented, variations in temperature and concentration can trigger a double-diffusive instability that enhances mixing and improves water recovery. We also show that the instability can be strengthened by providing external heating.

97 MATHEMATICS AND COMPUTING↗

THERMAL-GRAVITATIONAL WIND EQUATION FOR THE WIND-INDUCED GRAVITATIONAL SIGNATURE OF GIANT GASEOUS PLANETS: MATHEMATICAL DERIVATION, NUMERICAL METHOD, AND ILLUSTRATIVE SOLUTIONS

The standard thermal wind equation (TWE) relating the vertical shear of a flow to the horizontal density gradient in an atmosphere has been used to calculate the external gravitational signature produced by zonal winds in the interiors of giant gaseous planets. We show, however, that in this application the TWE needs to be generalized to account for an associated gravitational perturbation. We refer to the generalized equation as the thermal-gravitational wind equation (TGWE). The generalized equation represents a two-dimensional kernel integral equation with the Green’s function in its integrand and is hence much more difficult to solve than the standard TWE. We develop an extended spectral method for solving the TGWE in spherical geometry. We then apply the method to a generic gaseous Jupiter-like object with idealized zonal winds. We demonstrate that solutions of the TGWE are substantially different from those of the standard TWE. We conclude that the TGWE must be used to estimate the gravitational signature of zonal winds in giant gaseous planets.

79 ASTRONOMY AND ASTROPHYSICS↗