Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Lessons Learned (3 Years of H2O2 Propulsion System Testing Efforts at NASA's John C. Stennis Space Center)

John C. Stennis Space Center continues to support the Propulsion community in an effort to validate High-Test Peroxide as an alternative to existing/future oxidizers. This continued volume of peroxide test/handling activity at Stennis Space Center (SSC) provides numerous opportunities for the SSC team to build upon previously documented 'lessons learned'. SSC shall continue to strive to document their experience and findings as H2O2 issues surface. This paper is intended to capture all significant peroxide issues that we have learned over the last three years. This data (lessons learned) has been formulated from practical handling, usage, storage, operations, and initial development/design of our systems/facility viewpoint. The paper is intended to be an information type tool and limited in technical rational; therefore, presenting the peroxide community with some issues to think about as the continued interest in peroxide evolves and more facilities/hardware are built. These lessons learned are intended to assist industry in mitigating problems and identifying potential pitfalls when dealing with the requirements for handling high-test peroxide.

Taylor, Gary O.↗
NASA NTRS2011

Should Pruning be a Pre-Processor of any Linear System?

There are many real-world problems whose mathematical models turn out to be linear systems Ax = b, where A is an m x n matrix. Each equation of the linear system is an information. An information, in a physical problem, such as 4 mangoes, 6 bananas, and 5 oranges cost $10, is mathematically modeled as an equation 4x(sub 1) + 6x(sub 2) + 5x(sub 3) = 10 , where x(sub 1), x(sub 2), x(sub 3) are each cost of one mango, that of one banana, and that of one orange, respectively. All the information put together in a specified context, constitutes the physical problem and need not be all distinct. Some of these could be redundant, which cannot be readily identified by inspection. The resulting mathematical model will thus have equations corresponding to this redundant information and hence are linearly dependent and thus superfluous. Consequently, these equations once identified should be better pruned in the process of solving the system. The benefits are (i) less computation and hence less error and consequently a better quality of solution and (ii) reduced storage requirements. In literature, the pruning concept is not in vogue so far although it is most desirable. It is assumed that at least one information, i.e. one equation is known to be correct and which will be our first equation. In a numerical linear system, the system could be slightly inconsistent or inconsistent of varying degree. If the system is too inconsistent, then we should fall back on to the physical problem (PP), check the correctness of the PP derived from the material universe, modify it, if necessary, and then check the corresponding mathematical model (MM) and correct it. In nature/material universe, inconsistency is completely nonexistent. If the MM becomes inconsistent, it could be due to error introduced by the concerned measuring device and/or due to assumptions made on the PP to obtain an MM which is relatively easily solvable or simply due to human error. No measuring device can usually measure a quantity with an accuracy greater that 0.005% or, equivalently with a relative error less than 0.005%. Hence measurement error is unavoidable in a numerical linear system when the quantities are continuous (or even discrete with extremely large number). Assumptions, though not desirable, are usually made when we find the problem sufficiently difficult to be solved within the available means/tools/resources and hence distort the PP and the corresponding MM. The . error thus introduced in the system could (not always necessarily though) make the system somewhat inconsistent. If the inconsistency (contradiction) is too much then one should definitely not proceed to solve the system in terms of getting a least-squares solution or the minimum-norm least-squares solution. All these solutions will be invariably of no real-world use. If, on the other hand, inconsistency is reasonably low, i.e. the system is near-consistent or, equivalently, has near-linearly-dependent rows, then the foregoing solutions are useful. Pruning in such a near-consistent system should be performed based on the desired accuracy and on the definition of near-linear dependence. In this article, we discuss pruning over various kinds of linear systems and strongly suggest its use as a pre-processor or as a part of an algorithm. Ideally pruning should (i) be a part of the solution process (algorithm) of the system, (ii) reduce both computational error and complexity of the process, and (iii) take into account the numerical zero defined in the context. These are precisely what we achieve through our proposed O(mn2) algorithm presented in Matlab, that uses a subprogram of solving a single linear equation and that has embedded in it the pruning.

Sen, Syamal K.↗
NASA NTRS2011

Should Pruning be a Pre-Processor of any Linear System?

There are many real-world problems whose mathematical models turn out to be linear systems Ax = b , where A is an m by x n matrix. Each equation of the linear system is an information. An information, in a physical problem, such as 4 mangoes, 6 bananas, and 5 oranges cost $10, is mathematically modeled as 4x(sub 1) + 6x(sub 2) + 5x (sub 3) = 10, where x(sub 1), x(sub 2), x(sub 3) are each cost of one mango, that of one banana, and that of one orange, respectively. All the information put together in a specified context, constitutes the physical problem and need not be all distinct. Some of these could be redundant, which cannot be readily identified by inspection. The resulting mathematical model will thus have equations corresponding to this redundant information and hence are linearly dependent and thus superfluous. Consequently, these equations once identified should be better pruned in the process of solving the system. The benefits are (i) less computation and hence less error and consequently a better quality of solution and (ii) reduced storage requirements. In literature, the pruning concept is not in vogue so far although it is most desirable. In a numerical linear system, the system could be slightly inconsistent or inconsistent of varying degree. If the system is too inconsistent, then we should fall back on to the physical problem (PP), check the correctness of the PP derived from the material universe, modify it, if necessary, and then check the corresponding mathematical model (MM) and correct it. In nature/material universe, inconsistency is completely nonexistent. If the MM becomes inconsistent, it could be due to error introduced by the concerned measuring device and/or due to assumptions made on the PP to obtain an MM which is relatively easily solvable or simply due to human error. No measuring device can usually measure a quantity with an accuracy greater that 0.005% or, equivalently with a relative error less than 0.005%. Hence measurement error is unavoidable in a numerical linear system when the quantities are continuous (or even discrete with extremely large number). Assumptions, though not desirable, are usually made when we find the problem sufficiently difficult to be solved within the available means/tools/resources and hence distort the PP and the corresponding MM. The error thus introduced in the system could (not always necessarily though) make the system somewhat inconsistent. If the inconsistency (contradiction) is too much then one should definitely not proceed to solve the system in terms of getting a least-squares solution or a minimum norm solution or the minimum-norm least-squares solution. All these solutions will be invariably of no real-world use. If, on the other hand, inconsistency is reasonably low, i.e. the system is near-consistent or, equivalently, has near-linearly-dependent rows, then the foregoing solutions are useful. Pruning in such a near-consistent system should be performed based on the desired accuracy and on the definition of near-linear dependence. In this article, we discuss pruning over various kinds of linear systems and strongly suggest its use as a pre-processor or as a part of an algorithm. Ideally pruning should (i) be a part of the solution process (algorithm) of the system, (ii) reduce both computational error and complexity of the process, and (iii) take into account the numerical zero defined in the context. These are precisely what we achieve through our proposed O(mn2) algorithm presented in Matlab, that uses a subprogram of solving a single linear equation and that has embedded in it the pruning.

Sen, Syamal K.↗
NASA NTRS2002

Numerical Simulation of the Combustion of Fuel Droplets: Finite Rate Kinetics and Flame Zone Grid Adaptation (CEFD)

The NASA Nebraska Space Grant (NSGC) & EPSCoR programs have continued their effort to support outstanding research endeavors by funding the Numerical Simulation of the Combustion of Fuel Droplets study at the University of Nebraska at Lincoln (UNL). This team of researchers has developed a transient numerical model to study the combustion of suspended and moving droplets. The engines that propel missiles, jets, and many other devices are dependent upon combustion. Therefore, data concerning the combustion of fuel droplets is of immediate relevance to aviation and aeronautical personnel, especially those involved in flight operations. The experiments being conducted by Dr. Gogos and Dr. Nayagam s research teams, allow investigators to gather data for comparison with theoretical predictions of burning rates, flame structures, and extinction conditions. The consequent improved hndamental understanding droplet combustion may contribute to the clean and safe utilization of fossil hels (Williams, Dryer, Haggard & Nayagam, 1997, 72). The present state of knowledge on convective extinction of he1 droplets derives fiom experiments conducted under normal gravity conditions. However, any data obtained with suspended droplets under normal gravity are grossly affected by gravity. The need to obtain experimental data under microgravity conditions is therefore well justified and addresses one of the goals of NASA s Human Exploration and Development of Space (HEDS) microgravity combustion experiment.

Gogos, George↗
NASA NTRS2016

The Space-Time Conservative Schemes for Large-Scale, Time-Accurate Flow Simulations with Tetrahedral Meshes

Despite decades of development of unstructured mesh methods, high-fidelity time-accurate simulations are still predominantly carried out on structured, or unstructured hexahedral meshes by using high-order finite-difference, weighted essentially non-oscillatory (WENO), or hybrid schemes formed by their combinations. In this work, the space-time conservation element solution element (CESE) method is used to simulate several flow problems including supersonic jet/shock interaction and its impact on launch vehicle acoustics, and direct numerical simulations of turbulent flows using tetrahedral meshes. This paper provides a status report for the continuing development of the space-time conservation element solution element (CESE) numerical and software framework under the Revolutionary Computational Aerosciences (RCA) project. Solution accuracy and large-scale parallel performance of the numerical framework is assessed with the goal of providing a viable paradigm for future high-fidelity flow physics simulations.

Venkatachari, Balaji Shankar↗
NASA NTRS1989

Downward continuation of the free-air gravity anomalies to the ellipsoid using the gradient solution and terrain correction: An attempt of global numerical computations

The formulas for the determination of the coefficients of the spherical harmonic expansion of the disturbing potential of the earth are defined for data given on a sphere. In order to determine the spherical harmonic coefficients, the gravity anomalies have to be analytically downward continued from the earth's surface to a sphere-at least to the ellipsoid. The goal is to continue the gravity anomalies from the earth's surface downward to the ellipsoid using recent elevation models. The basic method for the downward continuation is the gradient solution (the g sub 1 term). The terrain correction was also computed because of the role it can play as a correction term when calculating harmonic coefficients from surface gravity data. The fast Fourier transformation was applied to the computations.

Wang, Y. M.↗
DOE OSTI2020

A filtering strategy for the numerical convergence of radiation transport through purely absorbing particle clouds

Radiation transport through particle clouds plays a major role in many engineering applications. Here in this paper, we study this problem by means of solving the radiative transport equation (RTE) on a Eulerian continuous domain. The particle clouds are generated through direct numerical simulation of Navier-Stokes equations coupled with Lagrangian particle transport at different Stokes numbers. We confirm the earlier observation noted in the literature that the solution to the RTE on Eulerian mesh diverges when the Eulerian mesh size is of the order of the particle diameter. This observation is often called the homogenization error stemming from relegation of number density onto the Eulerian domain. In order to circumvent this divergence problem, we propose a filtering strategy that spreads the information of the particles to the neighboring cells in a way that the representation of the particles remains the same even when the mesh size is smaller than that of the particle size. We show through our simulations that our filtering strategy solves the issue of homogenization error for the cases where the Beer-Bouger law is valid.

42 ENGINEERING↗
DOE OSTI2023

Control of Fractional Diffusion Problems via Dynamic Programming Equations

In this study, we explore the approximation of feedback control of integro-differential equations containing a fractional Laplacian term. To obtain feedback control for the state variable of this nonlocal equation, we use the Hamilton–Jacobi–Bellman equation. It is well known that this approach suffers from the curse of dimensionality, and to mitigate this problem we couple semi-Lagrangian schemes for the discretization of the dynamic programming principle with the use of Shepard approximation. This coupling enables approximation of high-dimensional problems. Numerical convergence toward the solution of the continuous problem is provided together with linear and nonlinear examples. The robustness of the method with respect to disturbances of the system is illustrated by comparisons with an open-loop control approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Narrow operator models of stellarator equilibria in Fourier Zernike basis

Numerical computation of the ideal magnetohydrodynamic (MHD) equilibrium magnetic field is at the base of stellarator optimisation and provides the starting point for solving more sophisticated partial differential equations like transport or turbulence models. Conventional approaches solve for a single stationary point of the ideal MHD equations, which is fully defined by three invariants and the numerical scheme employed by the solver. We present the first numerical approach that can solve for a continuous distribution of equilibria with fixed boundary and rotational transform, varying only the pressure invariant. This approach minimises the force residual by optimising parameters of multilayer perceptrons that map from a scalar pressure multiplier to the Fourier Zernike basis as implemented in the modern stellarator equilibrium solver DESC.

fusion plasma↗
DOE OSTI2025

Non-Hermitian Quantum Mechanics Approach for Extracting and Emulating Continuum Physics Based on Bound-State-like Calculations

Here, this Letter introduces a unified emulation framework for studying continuum physics in finite quantum systems. Using a reduced basis method, we construct powerful emulators for the inhomogeneous Schrödinger equation that operate in a combined parameter space of complex energy (𝐸) and other inputs (𝜽). Within the space, the emulators simultaneously perform analytical continuation in 𝐸—extracting continuum physics from numerically simpler bound-state-like calculations—and interpolate this entire process across 𝜽. This yields a small, non-Hermitian system whose properties (e.g., resonances and scattering observables) can be rapidly predicted for any 𝜽. Crucially, the complex-𝐸 emulation provides a pathway to compute continuum observables for complex systems where advanced bound-state methods exist but direct continuum calculations are yet to be developed, while the 𝜽 emulation enables rapid parameter-space exploration and can be adapted to accelerate other existing continuum calculations. Demonstrations with two- and three-body systems highlight the method’s effectiveness and suggest its connection to (near-)optimal rational approximation. This Letter presents the key results, with further details reserved for a companion paper.

ab initio calculations↗
DOE OSTI2025

EIC Physics from Lattice QCD: investigations beyond leading twist

This project proposes theoretical studies of Quantum Chromodynamics (QCD), the theory describing the strong nuclear force among the building blocks (quarks and gluons) of the visible matter. These appear only confined within hadrons, that make up more than 99% of the mass of the matter. Understanding QCD will significantly advance many aspects of science, from the sub-nuclear interactions to astrophysics, and a quantitative theoretical description is imperative. However, this is a challenging task because QCD is a highly nonlinear theory. We propose hadron structure calculations within lattice QCD (LQCD), an ideal ab initio approach based on space-time discretization, which allows the study of the properties of fundamental particles numerically. This is done by defining the continuous equations on a discrete four-dimensional lattice, which results in equations with hundreds of billions of degrees of freedom, and must be simulated in powerful computers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Transboundary determinants of avian zoonotic infectious diseases: challenges for strengthening research capacity and connecting surveillance networks

As the climate changes, global systems have become increasingly unstable and unpredictable. This is particularly true for many disease systems, including subtypes of highly pathogenic avian influenzas (HPAIs) that are circulating the world. Ecological patterns once thought stable are changing, bringing new populations and organisms into contact with one another. Wild birds continue to be hosts and reservoirs for numerous zoonotic pathogens, and strains of HPAI and other pathogens have been introduced into new regions via migrating birds and transboundary trade of wild birds. With these expanding environmental changes, it is even more crucial that regions or counties that previously did not have surveillance programs develop the appropriate skills to sample wild birds and add to the understanding of pathogens in migratory and breeding birds through research. For example, little is known about wild bird infectious diseases and migration along the Mediterranean and Black Sea Flyway (MBSF), which connects Europe, Asia, and Africa. Focusing on avian influenza and the microbiome in migratory wild birds along the MBSF, this project seeks to understand the determinants of transboundary disease propagation and coinfection in regions that are connected by this flyway. Through the creation of a threat reduction network for avian diseases (Avian Zoonotic Disease Network, AZDN) in three countries along the MBSF (Georgia, Ukraine, and Jordan), this project is strengthening capacities for disease diagnostics; microbiomes; ecoimmunology; field biosafety; proper wildlife capture and handling; experimental design; statistical analysis; and vector sampling and biology. Here, we cover what is required to build a wild bird infectious disease research and surveillance program, which includes learning skills in proper bird capture and handling; biosafety and biosecurity; permits; next generation sequencing; leading-edge bioinformatics and statistical analyses; and vector and environmental sampling. Creating connected networks for avian influenzas and other pathogen surveillance will increase coordination and strengthen biosurveillance globally in wild birds.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1973

Prediction of overall and blade-element performance for axial-flow pump configurations

A method and a digital computer program for prediction of the distributions of fluid velocity and properties in axial flow pump configurations are described and evaluated. The method uses the blade-element flow model and an iterative numerical solution of the radial equilbrium and continuity conditions. Correlated experimental results are used to generate alternative methods for estimating blade-element turning and loss characteristics. Detailed descriptions of the computer program are included, with example input and typical computed results.

Serovy, G. K.↗
NASA NTRS1974

OSO-7 observations of solar X-rays in the energy range 10-100 keV

Data on 123 hard X-ray bursts observed by the satellite OSO-7 between Oct. 10, 1971 and June 6, 1972 are described and evaluated. Typical duration of a burst is 100 sec. Average spectral indices lie between 3.5 and 5.5 for two-thirds of the 123 bursts, with a median of 4.6. In some events a soft-hard-soft pattern is observed, but there are numerous examples in which the spectrum softens continuously throughout the burst. The mean shape of the hard X-ray time profile as measured by the full width at half maximum does not depend on burst amplitude; nor does the spectral hardness correlate with the flux. The distribution of burst peak fluxes and the observation of large soft X-ray bursts without accompanying hard X rays suggest the existence of a distinct class of solar flares which emit only soft X rays. No center-to-limb variation was found in the frequency of occurrence of bursts or in the fraction with a nonthermal component. Estimates of the energy in the form of nonthermal electrons and in the flare plasma derived from these data indicate that the total amounts in each are comparable.

Datlowe, D. W.↗
NASA NTRS1975

Theoretical studies of massive stars. I - Evolution of a 15-solar-mass star from the zero-age main sequence to neon ignition

The evolution of a star with mass 15 times that of the sun from the zero-age main sequence to neon ignition has been computed by the Henyey method. The hydrogen-rich envelope and all shell sources were explicitly included in the models. An algorithm has been developed for approximating the results of carbon burning, including the branching ratio for the C-12 + C-12 reaction and taking some secondary reactions into account. Penetration of the convective envelope into the core is found to be unimportant during the stages covered by the models. Energy transfer from the carbon-burning shell to the core by degenerate electron conduction becomes important after the core carbon-burning stage. Neon ignition will occur in a semidegenerate core and will lead to a mild 'flash.' Detailed numerical results are given in an appendix. Continuation of the calculations into later stages and variations with the total mass of the star will be discussed in later papers.

Endal, A. S.↗
NASA NTRS1976

A study of efficiency in low resistivity silicon solar cells

A general device-analysis program has been utilized to study the efficiency of silicon solar cells. The analysis is applied to specific geometries of both n(+)-p and n(+)-p-p(+) solar cells, and involves a numerical solution of the basic transport and continuity equations. This approach allows solutions that are free of typical limiting assumptions involved in solving the device equations as well as solutions relating to lifetime, mobility variations, and diffused-region profiles. The analysis includes available empirical information on diffusion length, mobility, and lifetime as a function of doping, as well as a Gaussian profile for the diffused region. Results are presented which illustrate the limitations of efficiency as a function of doping. It is found that the maximum efficiencies for both types of cell converge at lower resistivities to around 16% with air-mass-zero radiation and a single-layer absorbing-SiO antireflecting film. It is also found that the minority-carrier lifetime, both in the n(+) surface and p-type bulk regions, presents serious limitations to conversion efficiency, particularly in the low-resistivity cells.

Dunbar, P. M.↗
NASA NTRS1979

Numerical analysis of turbulent wakes of axial flow turbomachinery rotor blades

The wakes of rotor blades are turbulent, three-dimensional, and are subjected to curvature and rotation effects. The objective of this study is to predict the development of such wakes and compare the predictions with the existing data. A finite difference procedure is employed in the numerical analysis of the wake utilizing the continuity, momentum, and turbulence closure equations in the rotating curvilinear and non-orthogonal coordinate system. The turbulence closure is affected by the two equations model and the simplified Reynolds stress model. The predictions from both the models agree well with the mean velocity and the turbulence data.

Hah, C.↗
NASA NTRS1983

Research on the control of large space structures

The research effort on the control of large space structures at the University of Houston has concentrated on the mathematical theory of finite-element models; identification of the mass, damping, and stiffness matrix; assignment of damping to structures; and decoupling of structure dynamics. The objective of the work has been and will continue to be the development of efficient numerical algorithms for analysis, control, and identification of large space structures. The major consideration in the development of the algorithms has been the large number of equations that must be handled by the algorithm as well as sensitivity of the algorithms to numerical errors.

Denman, E. D.↗