Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical mathematics”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Chapter 14: Machine Learning of Combustion LES Models from Reacting Direct Numerical Simulation

In this chapter we demonstrate how supervised deep learning techniques can be used to construct models for the filtered progress variable source term necessary for large eddy simulation (LES). The source data for the model is a direct numerical simulation (DNS) of a reacting flow in a low swirl burner configuration. Filtered quantities taken from the DNS data are used to train a deep neural network (DNN)-based model. An efficient data sampling strategy was devised to ensure that a uniform representation of all the states observed in the filtered DNS data are equally present in the training dataset. A-priori testing of the DNN-based model highlights the representative power of DNN to accurately reproduce the filtered reaction progress variable source term over a range of scales and various flame regimes as seen in an industrial burner.

combustion LES models↗
NASA NTRS1993

Program Computes Flows Of Fluids And Heat

SINDA'85/FLUINT incorporates lumped-parameter-network and one-dimensional-flow mathematical models. System enables analysis of mutual influences of thermal and flow phenomena. Offers two finite-difference numerical solution techniques: forward-difference explicit approximation and Crank-Nicholson approximation. Enables simulation of nonuniform heating and facilitates mathematical modeling of thin-walled heat exchangers. Ability to model nonequilibrium behavior within two-phase volumes included. Recent changes in program improve modeling of real evaporator pumps and other capillary-assist evaporators. Written in FORTRAN 77.

Cullimore, Brent↗
NASA NTRS2001

Radiometry Measurements of Mars at 1064 nm Using the Mars Orbiter Laser Altimeter

Measurements by the Mars Orbiter Laser Altimeter (MOLA) on board the Mars Global Surveyor (MGS) may be used to provides a radiometric measurement of Mars in addition to the topographic measurement. We will describe the principle of operation, a mathematical model, and the receiver calibration in this presentation. MOLA was designed primarily to measure Mars topography, surface roughness end the bidirectional reflectance to the laser beam. To achieve the highest sensitivity the receiver detection threshold is dynamically adjusted to be as low as possible while keeping a predetermined false alarm rate. The average false alarm rate 29 monitored in real time on board MOLA via a noise counter, whose output is fed to the threshold control loop. The false alarm rate at a given threshold is a function of the detector output noise which is the sum of the photo detector, shot noise due to the background light seen by the detector and the dark noise. A mathematical model has been developed that can be used to numerically solve for the optical background power given the MOLA threshold setting and the average noise count. The radiance of Mars can then be determined by dividing the optical power by the solid angle subtended by the MOLA receiver, the receiver optical band-width, end the Mars surface area within the receiver field of view. The phase angle which is the sun-Mars-MOLA angle is available from the MGS database. MOLA also measures simultaneously the bidirectional reflectance of Mars vie its 106-lum loser beam at nadir with nearly zero phase angle. The optical bandwidth of the MOLA receiver is 2um full width at half maximum (FWHM) and centered at 106-lum. The receiver field of view is 0.95mrad FWHM. The nominated spacecraft altitude is 100km and the ground track speed is about 3km/s. Under normal operation, the noise counter are read and the threshold levels are updated at 1Hz. The receiver sensitivity is limited by the detector dark noise to about 0.1nW, which corresponds to less than 2% the maximum radiance during daytime from the brightest area on Mars. The results from the mathematical model agree well with the prelaunch measurements at several calibrated optical power levels. The radiance of sunlit Mars estimated with this technique correlates well with the measurement from the MGS. Thermal Emission Spectrometer (TES) and the Hubble Space Telescope at similar wavelength.

Sun, Xiao-Li↗
DOE OSTI2024

Quantifying local and global mass balance errors in physics-informed neural networks

Physics-informed neural networks (PINN) have recently become attractive for solving partial differential equations (PDEs) that describe physics laws. By including PDE-based loss functions, physics laws such as mass balance are enforced softly in PINN. This paper investigates how mass balance constraints are satisfied when PINN is used to solve the resulting PDEs. We investigate PINN’s ability to solve the 1D saturated groundwater flow equations (diffusion equations) for homogeneous and heterogeneous media and evaluate the local and global mass balance errors. We compare the obtained PINN’s solution and associated mass balance errors against a two-point finite volume numerical method and the corresponding analytical solution. We also evaluate the accuracy of PINN in solving the 1D saturated groundwater flow equation with and without incorporating hydraulic heads as training data. We demonstrate that PINN’s local and global mass balance errors are significant compared to the finite volume approach. Tuning the PINN’s hyperparameters, such as the number of collocation points, training data, hidden layers, nodes, epochs, and learning rate, did not improve the solution accuracy or the mass balance errors compared to the finite volume solution. Mass balance errors could considerably challenge the utility of PINN in applications where ensuring compliance with physical and mathematical properties is crucial.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1987

A generalized procedure for constructing an upwind based TVD scheme

A generalized formulation for constructing second- and higher-order accurate TVD (total variation diminishing) schemes is presented. A given scheme is made TVD by limiting antidiffusive flux differences with some linear functions, so-called limiters. The general idea of the formulation and its mathematical proof of Harten's TVD conditions is shown by applying the Lax-Wendroff method to scalar nonlinear equations and a constant-coefficient system of conservation laws. For the system of equations, several definitions are derived for the argument used in the limiter function and present their performance in numerical experiments. The formulation is extended to the nonlinear system. It is demonstrated that the present procedure can easily convert existing central or upwind, and second- or higher-order differencing schemes to preserve monotonicity and yield physically admissible solutions. The formulation is simple mathematically as well as numerically; both matrix-vector multiplication and Riemann solver are avoided. Although the notion of TVD is based on the initial value problem, application to the steady Euler equations of the formulation is also made.

Liou, Meng-Sing↗
NASA NTRS1982

A mathematical model of post-instability in fluid mechanics

Postinstability of fluids is eliminated in numerical models by introducing multivalued velocity fields after discarding the principle of impenetrability. Smooth functions are shown to be incapable of keeping the derivatives from going towards infinity when iterating solutions for the governing equations such as those defined by Navier-Stokes. Enlarging the class of functions is shown to be necessary to eliminate the appearance of imaginary characteristic roots in the systems of arbitrary partial differential equations, a condition which leads to physically impossible motions. The enlarging is demonstrated to be achievable by allowing several individual particles with different velocities to appear at the same point of space, and the subsequent multivaluedness of the solutions is purely a mathematical concern, rather than one of actual physical existence. Applications are provided for an inviscid fluid and for turbulence.

Zak, M. A.↗
DOE OSTI2024

Can mesoscale models capture the effect from cluster wakes offshore?

Long wakes from offshore wind turbine clusters can extend tens of kilometers downstream, affecting the wind resource of a large area. Given the ability of mesoscale numerical weather prediction models to capture important atmospheric phenomena and mechanisms relevant to wake evolution, they are often used to simulate wakes behind large wind turbine clusters and their impact over a wider region. Yet, uncertainty persists regarding the accuracy of representing cluster wakes via mesoscale models and their wind turbine parameterizations. Here, we evaluate the accuracy of the Fitch wind farm parameterization in the Weather Research and Forecasting model in capturing cluster-wake effects using two different options to represent turbulent mixing in the planetary boundary layer. To this end, we compare operational data from an offshore wind farm in the North Sea that is fully or partially waked by an upstream array against high-resolution mesoscale simulations. In general, we find that mesoscale models accurately represent the effect of cluster wakes on front-row turbines of a downstream wind farm. However, the same models may not accurately capture cluster-wake effects on an entire downstream wind farm, due to misrepresenting internal-wake effects.

17 WIND ENERGY↗
DOE OSTI2022

WORM (Write One, Read Many)

WORM (Write One, Run Many) is an easy to use, cross platform, embedded and extensible, functional programming language designed to facilitate the creation of input-decks for computer codes that use standard ASCII text files for input. WORM makes it easy to create generic (yet, complex and powerful) reusable models. Additionally its nature allows for complex calculations and routines to be coded once and easily reused, further simplifying the creation of input decks. WORM (Write One, Run Many) is a powerful and versatile tool designed to improve the efficiency of today’s criticality safety analyst by allowing: + input decks for parametric studies to be created quickly and easily, + calculations and variables to be imbedded into any input deck, thus allowing for meaningful parameter specifications, + problems to be specified using any combination of units, and + complex mathematically defined models to be created. A very simple syntax is employed, and therefore the WORM is easy to learn. A WORM model is essentially a standard input deck with some of its numerical values replaced by WORM code. WORM code may include and evaluate the following mathematical operators and functions: addition, subtraction, multiplication, division, exponentiation, modulus, sine, cosine, tangent, arcsine, arccosine, arctangent, the natural logarithm, logarithm base 10, integer truncation, absolute value, and random number. Several common constants, e.g., pi, e, and Avogadros’s Number (both as 6.022e23 and 0.6022), are predefined in WORM. Additionally, many unit conversion factors are also predefined: millimeters, meters, inches, feet, yards, and mils to centimeters; kilograms, pounds, and ounces to grams; liters, milliliters, gallons, and fluid ounces to cubic centimeters; and angular degrees to radians. For parametric studies, WORM supports various shorthand list specifications: the explicit step size, linear interpolation, and logarithmic interpolation. The list notation sequentially assigns multiple values to a name. WORM creates an input deck for each value of the name. If multiple lists are used, WORM steps through each list individually, i.e., WORM creates input decks corresponding to each and every permutation of the list values. Additionally, a library of standard material definitions and Perl subroutines are included. Any one of these files can be incorporated into the subject model with a simple WORM read command. WORM is completely written in Perl, the Practical Extraction and Reporting Language. Perl is one of the most portable programming languages available today. As such, the WORM works on practically any computer platform.

Sartor, Raymond↗
DOE OSTI2026

Assessing the Impact of a Forest Canopy on Near-Surface Wind Statistics

Representing the forest canopy in atmospheric numerical models should improve simulated winds within and above the canopy up to a few hundred meters above the ground. Here, in this study, we implement a forest canopy parameterization into the Weather Research and Forecasting (WRF) Model in a large-eddy simulation (LES) mode by applying drag forces across multiple layers within the canopy height. We use unique observations from the Lidar Experiments for Assessing Flow over Forests (LEAFF) field campaign at the Wind River Experimental Forest (WREF) in the U.S. Pacific Northwest to evaluate model performance. In a 2-day case study, the canopy parameterization improved wind predictions both within and above the canopy, particularly during the daytime and at finer grid resolution. Without it, winds were frequently overpredicted above the canopy. Similarly, derived quantities such as the wind shear index also yielded estimates closer to observations with the canopy parameterization implemented. These findings suggest that representing the canopy using drag forces alone can improve simulated mean winds up to 200 m above the surface. Furthermore, second-order statistical moments of wind were more sensitive to canopy density than first-order moments, especially during the daytime. This increased sensitivity and the improved daytime performance in wind speed—evidenced by the lowest bias from observations (3% compared to 20% over diurnal cycle)—imply that winds above the canopy layer are strongly influenced by how well turbulence above the canopy is modeled. The results of this study can serve as a foundation for parameterizing forest canopy effects in coarser weather forecast models.

Energy - Wind↗
NASA NTRS1978

Application of semianalytical satellite theories to precision orbit determination

Those factors which limit the usefulness of current implementations of the semianalytical approach are discussed. Numerical and analytical enhancements to the semianalytical approach are considered. A simple mathematical model is provided to estimate the computational speed of a semianalytical theory employing the suggested enhancements. The model can factor in current experience with semianalytical theories (integration stepsizes, quadrature orders, speed of recursive formulations, etc.) and the characteristics of the particular output requirement (observation span (or orbit determination interval), observation rate, observation model, etc.). Comparisons with numerical integration are suggested.

Cefola, P. J.↗
NASA NTRS1987

A generalized procedure for constructing an upwind-based TVD scheme

A generalized formulation for constructing second- and higher-order accurate TVD (total variation diminishing) schemes is presented. A given scheme is made TVD by limiting antidiffusive flux differences with some nonlinear functions, so-called limiters. The general idea of the formulation and its mathematical proof of Harten's TVD conditions is shown by applying the Lax-Wendroff method to a scalar nonlinear equation and constant-coefficient system of conservation laws. For the system of equations, several definitions are derived for the argument used in the limiter function and present their performance to numerical experiments. Then the formulation is formally extended to the nonlinear system of equations. It is demonstrated that use of the present procedure allows easy conversion of existing central or upwind, and second- or higher-order differencing schemes so as to preserve monotonicity and to yield physically admissible solutions. The formulation is simple mathematically as well as numerically; neither matrix-vector multiplication nor Riemann solver is required. Roughly twice as much computational effort is needed as compared to conventional scheme. Although the notion of TVD is based on the initial value problem, application to the steady Euler equations of the formulation is also made. Numerical examples including various ranges of problems show both time- and spatial-accuracy in comparison with exact solutions.

Liou, Meng-Sing↗
DOE OSTI2021

Qualitative properties of mathematical model for data flow

In this paper, properties of a recently proposed mathematical model for data flow in large-scale asynchronous computer systems are analyzed. In particular, the existence of special weak solutions based on propagating fronts is established. Qualitative properties of these solutions are investigated, both theoretically and numerically.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Numerical fluid dynamics for FRG flow equations: Zero-dimensional QFTs as numerical test cases. II. Entropy production and irreversibility of RG flows

Herein we demonstrate that the reformulation of renormalization group (RG) flow equations as nonlinear heat equations has severe implications on the understanding of RG flows in general. We demonstrate by explicitly constructing an entropy function for a zero-dimensional Z 2 -symmetric model that the dissipative character of generic nonlinear diffusion equations is also hard-coded in the functional RG equation. This renders RG flows manifestly irreversible, revealing the semigroup property of RG transformations on the level of the flow equation itself. Additionally, we argue that the dissipative character of RG flows, its irreversibility and the entropy production during the RG flow may be linked to the existence of a so-called C– / A-function. In total, this introduces an asymmetry in the so-called RG time—in complete analogy to the thermodynamic arrow of time—and allows for an interpretation of infrared actions as equilibrium solutions of dissipative RG flows equations. The impossibility of resolving microphysics from macrophysics is evident in this framework. Furthermore, we directly link the irreversibility and the entropy production in RG flows to an explicit numerical entropy production, which is manifest in diffusive and non-linear partial differential equations (PDEs) and a standard mathematical tool for the analysis of PDEs. Using exactly solvable zero-dimensional Z 2 -symmetric models, we explicitly compute the (numerical) entropy production related to the total variation nonincreasing property of the PDE during RG flows toward the infrared limit. Finally, we discuss generalizations of our findings and relations to the C– / A-theorem as well as how our work may help to construct truncations of RG flow equations in the future, including numerically stable schemes for solving the corresponding PDEs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

The reflection of a blast wave by a very intense explosion

We demonstrate that the geometric similarity of Taylor’s blast wave persists beyond reflection from an ideal surface. Upon impacting the surface, the spherical symmetry of the blast wave is lost but its cylindrical symmetry endures. As the flow acquires dependence on a second spatial dimension, an analytic solution of the Euler equations becomes elusive. However, the preservation of axisymmetry, geometric similarity and planar symmetry in the presence of a mirror-like surface causes all flow solutions to collapse when scaled by the height of burst (HOB) and the shock arrival time at the surface. The scaled blast volume for any yield, HOB and ambient air density follows a single universal trajectory for all scaled time, both before and after reflection.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗