Search NASA⌕ Search

SEARCH · Search NASA

Results for “convention equation”

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 289 records · Page 16

Predicting Ecologically Important Vegetation Variables from Remotely Sensed Optical/Radar Data Using Neural Networks

A number of satellite sensor systems will collect large data sets of the Earth's surface during NASA's Earth Observing System (EOS) era. Efforts are being made to develop efficient algorithms that can incorporate a wide variety of spectral data and ancillary data in order to extract vegetation variables required for global and regional studies of ecosystem processes, biosphere-atmosphere interactions, and carbon dynamics. These variables are, for the most part, continuous (e.g. biomass, leaf area index, fraction of vegetation cover, vegetation height, vegetation age, spectral albedo, absorbed photosynthetic active radiation, photosynthetic efficiency, etc.) and estimates may be made using remotely sensed data (e.g. nadir and directional optical wavelengths, multifrequency radar backscatter) and any other readily available ancillary data (e.g., topography, sun angle, ground data, etc.). Using these types of data, neural networks can: 1) provide accurate initial models for extracting vegetation variables when an adequate amount of data is available; 2) provide a performance standard for evaluating existing physically-based models; 3) invert multivariate, physically based models; 4) in a variable selection process, identify those independent variables which best infer the vegetation variable(s) of interest; and 5) incorporate new data sources that would be difficult or impossible to use with conventional techniques. In addition, neural networks employ a more powerful and adaptive nonlinear equation form as compared to traditional linear, index transformations, and simple nonlinear analyses. These neural networks attributes are discussed in the context of the authors' investigations of extracting vegetation variables of ecological interest.

Kimes, Daniel S.↗

Software Simulates Sight: Flat Panel Mura Detection

In the increasingly sophisticated world of high-definition flat screen monitors and television screens, image clarity and the elimination of distortion are paramount concerns. As the devices that reproduce images become more and more sophisticated, so do the technologies that verify their accuracy. By simulating the manner in which a human eye perceives and interprets a visual stimulus, NASA scientists have found ways to automatically and accurately test new monitors and displays. The Spatial Standard Observer (SSO) software metric, developed by Dr. Andrew B. Watson at Ames Research Center, measures visibility and defects in screens, displays, and interfaces. In the design of such a software tool, a central challenge is determining which aspects of visual function to include while accuracy and generality are important, relative simplicity of the software module is also a key virtue. Based on data collected in ModelFest, a large cooperative multi-lab project hosted by the Optical Society of America, the SSO simulates a simplified model of human spatial vision, operating on a pair of images that are viewed at a specific viewing distance with pixels having a known relation to luminance. The SSO measures the visibility of foveal spatial patterns, or the discriminability of two patterns, by incorporating only a few essential components of vision. These components include local contrast transformation, a contrast sensitivity function, local masking, and local pooling. By this construction, the SSO provides output in units of "just noticeable differences" (JND) a unit of measure based on the assumed smallest difference of sensory input detectable by a human being. Herein is the truly amazing ability of the SSO, while conventional methods can manipulate images, the SSO models human perception. This set of equations actually defines a mathematical way of working with an image that accurately reflects the way in which the human eye and mind behold a stimulus. The SSO is intended for a wide variety of applications, such as evaluating vision from unmanned aerial vehicles, measuring visibility of damage to aircraft and to the space shuttles, predicting outcomes of corrective laser eye surgery, inspecting displays during the manufacturing process, estimating the quality of compressed digital video, evaluating legibility of text, and predicting discriminability of icons or symbols in a graphical user interface.

Source record↗

Revenue prediction for integrated renewable energy and energy storage system using machine learning techniques

Revenue estimation for integrated renewable energy and energy storage systems is important to support plant owners or operators’ decisions in battery sizing selection that leads to maximized financial performances. A common approach to optimizing revenues of a hybrid hydro and energy storage system is using mixed-integer linear programming (MILP). Although MILP models can provide accurate production cost estimations, they are typically very computationally expensive. To provide a fast yet accurate first-step information to hydropower plant owners or operators who consider integrating energy storage systems, we propose an innovative approach to predicting optimal revenues of an integrated energy generation and storage system. In this study, we examined the performance of two prediction techniques: Generalized Additive Models (GAMs) and machine learning (ML) models developed based on artificial neural networks (ANN). Predictive equations and models are generated based on optimized solutions from a market participation optimization model, the Conventional Hydropower Energy and Environmental Resource System (CHEERS) model. The two predicting techniques reduce the computational time to evaluate annual revenue for one set of battery configurations from 3 h to 1 to 4 min per run while also being implementable with significantly less data. The model validation prediction errors of developed GAMs and ML models are generally below 5%; for model testing predictions, the ML models consistently outperform the regression equations in terms of root mean square errors. This new approach allows plant owners, operators, or potential investors to quickly access multiple battery configurations under different energy generation and market scenarios. This new revenue prediction method will therefore help reduce the barriers, and thereby promoting the deployment of battery hybridization with existing renewable energy sources.

13 HYDRO ENERGY↗

Applying Physics-Informed Neural Networks to Solve Navier–Stokes Equations for Laminar Flow around a Particle

In recent years, Physics-Informed Neural Networks (PINNs) have drawn great interest among researchers as a tool to solve computational physics problems. Unlike conventional neural networks, which are black-box models that “blindly” establish a correlation between input and output variables using a large quantity of labeled data, PINNs directly embed physical laws (primarily partial differential equations) within the loss function of neural networks. By minimizing the loss function, this approach allows the output variables to automatically satisfy physical equations without the need for labeled data. The Navier–Stokes equation is one of the most classic governing equations in thermal fluid engineering. This study constructs a PINN to solve the Navier–Stokes equations for a 2D incompressible laminar flow problem. Flows passing around a 2D circular particle are chosen as the benchmark case, and an elliptical particle is also examined to enrich the research. The velocity and pressure fields are predicted by the PINNs, and the results are compared with those derived from Computational Fluid Dynamics (CFD). Additionally, the particle drag force coefficient is calculated to quantify the discrepancy in the results of the PINNs as compared to CFD outcomes. The drag coefficient maintained an error within 10% across all test scenarios.

Hu, Beichao (ORCID:0009000163215151)↗

Multigrid applications to three-dimensional elliptic coordinate generation

The multigrid algorithm was applied to solve the coupled set of elliptic quasilinear partial differential equations associated with three-dimensional coordinate generation. The results indicate that the multigrid scheme is more than twice as fast as conventional relaxation schemes on moderate-size grids. Convergence factors of order 0.90 per work unit were achieved on 36,000-point grids. The paper covers the form of transformation, develops the set of generation equations, and gives details on the multigrid approach used. Included are a development of the full-approximation storage scheme, details of the smoothing-rate analysis, and a section devoted to rational programming techniques applicable to the multigrid algorithm.

Thames, F. C.↗

High Order Difference Method for Low Mach Number Aeroacoustics

A high order finite difference method with improved accuracy and stability properties for computational aeroacoustics (CAA) at low Mach numbers is proposed. The Euler equations are split into a conservative and a symmetric non- conservative portion to allow the derivation of a generalized energy estimate. Since the symmetrization is based on entropy variables, that splitting of the flux derivatives is referred to as entropy splitting. Its discretization by high order central differences was found to need less numerical dissipation than conventional conservative schemes. Owing to the large disparity of acoustic and stagnation quantities in low Mach number aeroacoustics, the split Euler equations are formulated in perturbation form. The unknowns are the small changes of the conservative variables with respect to their large stagnation values. All nonlinearities and the conservation form of the conservative portion of the split flux derivatives can be retained, while cancellation errors are avoided with its discretization opposed to the conventional conservative form. The finite difference method is third-order accurate at the boundary and the conventional central sixth-order accurate stencil in the interior. The difference operator satisfies the summation by parts property analogous to the integration by parts in the continuous energy estimate. Thus, strict stability of the difference method follows automatically. Spurious high frequency oscillations are suppressed by a characteristic-based filter similar to but without limiter. The time derivative is approximated by a 4-stage low-storage second-order explicit Runge-Kutta method. The method has been applied to simulate vortex sound at low Mach numbers. We consider the Kirchhoff vortex, which is an elliptical patch of constant vorticity rotating with constant angular frequency in irrotational flow. The acoustic pressure generated by the Kirchhoff vortex is governed by the 2D Helmholtz equation, which can be solved analytically using separation of variables.

Mueller, B.↗

gLaSDI: Parametric physics-informed greedy latent space dynamics identification

A parametric adaptive physics-informed greedy Latent Space Dynamics Identification (gLaSDI) method is proposed for accurate, efficient, and robust data-driven reduced-order modeling of high-dimensional nonlinear dynamical systems. In the proposed gLaSDI framework, an autoencoder discovers intrinsic nonlinear latent representations of high-dimensional data, while dynamics identification (DI) models capture local latent-space dynamics. Here, an interactive training algorithm is adopted for the autoencoder and local DI models, which enables identification of simple latent-space dynamics and enhances accuracy and efficiency of data-driven reduced-order modeling. To maximize and accelerate the exploration of the parameter space for the optimal model performance, an adaptive greedy sampling algorithm integrated with a physics-informed residual-based error indicator and random-subset evaluation is introduced to search for the optimal training samples on the fly. Further, to exploit local latent-space dynamics captured by the local DI models for an improved modeling accuracy with a minimum number of local DI models in the parameter space, a -nearest neighbor convex interpolation scheme is employed. The effectiveness of the proposed framework is demonstrated by modeling various nonlinear dynamical problems, including Burgers equations, nonlinear heat conduction, and radial advection. The proposed adaptive greedy sampling outperforms the conventional predefined uniform sampling in terms of accuracy. Compared with the high-fidelity models, gLaSDI achieves 17 to 2,658× speed-up with 1 to 5% relative errors.

97 MATHEMATICS AND COMPUTING↗

Dyn$\mathrm{AMO}$: Multi-agent reinforcement learning for dynamic anticipatory mesh optimization with applications to hyperbolic conservation laws

Here we introduce DynAMO, a reinforcement learning paradigm for Dynamic Anticipatory Mesh Optimization. Adaptive mesh refinement is an effective tool for optimizing computational cost and solution accuracy in numerical methods for partial differential equations. However, traditional adaptive mesh refinement approaches for time-dependent problems typically rely only on instantaneous error indicators to guide adaptivity. As a result, standard strategies often require frequent remeshing to maintain accuracy. In the DynAMO approach, multi-agent reinforcement learning is used to discover new local refinement policies that can anticipate and respond to future solution states by producing meshes that deliver more accurate solutions for longer time intervals. By applying DynAMO to discontinuous Galerkin methods for the linear advection and compressible Euler equations in two dimensions, we demonstrate that this new mesh refinement paradigm can outperform conventional threshold-based strategies while also generalizing to different mesh sizes, remeshing and simulation times, and initial conditions.

97 MATHEMATICS AND COMPUTING↗

Verification of a fully implicit particle-in-cell method for the <!--${MathJax: TeX-AMS-MML_HTMLorMML}--> v &#x2225; -formalism of electromagnetic gyrokinetics in the XGC code

A fully implicit particle-in-cell method for handling the v ∥ -formalism of electromagnetic gyrokinetics has been implemented in XGC. By choosing the v ∥ -formalism, here we avoid introducing the nonphysical skin terms in Ampère's law, which are responsible for the well-known “cancellation problem” in the p ∥ -formalism. The v ∥ -formalism, however, is known to suffer from a numerical instability when explicit time integration schemes are used due to the appearance of a time derivative in the particle equations of motion from the inductive component of the electric field. Here, using the conventional δf scheme, we demonstrate that our implicitly discretized algorithm can provide numerically stable simulation results with accurate dispersive properties. We verify the algorithm using a test case for shear Alfvén wave propagation in addition to a case demonstrating the ion temperature gradient-kinetic ballooning mode (ITG-KBM) transition. The ITG-KBM transition case is compared to results obtained from other δf gyrokinetic codes/schemes, whose verification has already been archived in the literature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electron temperature gradient driven transport model for tokamak plasmas

A new model for electron temperature gradient (ETG) modes is developed as a component of the multi-mode anomalous transport module [Rafiq et al., Phys Plasmas 20, 032506 (2013)] to predict a time-dependent electron temperature profile in conventional and low aspect ratio tokamaks. This model is based on two-fluid equations that govern the dynamics of low-frequency short- and long-wavelength electromagnetic toroidal ETG driven drift modes. A low collisionality NSTX discharge is used to scan the plasma parameter dependence on the ETG real frequency, growth rate, and electron thermal diffusivity. Electron thermal transport is discovered in the deep core region where modes are more electromagnetic in nature. Several previously reported gyrokinetic trends are reproduced, including the dependencies of density gradients, magnetic shear, β and gradient of β (β'), collisionality, safety factor, and toroidicity, where β is the ratio of the plasma pressure to the magnetic pressure. Furthermore, the electron heat diffusivity associated with the ETG mode is discovered to be on a scale consistent with the experimental diffusivity determined by power balance analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electron Temperature Gradient Driven Transport Model for Tokamak Plasmas

A new model for electron temperature gradient (ETG) modes is developed as a component of the Multi-Mode anomalous transport module [T. Rafiq \textit{et al.,} Phys Plasmas \textbf{20}, 032506 (2013)] to predict a time dependent electron temperature profile in conventional and low aspect ratio tokamaks. This model is based on two-fluid equations that govern the dynamics of low-frequency short- and long-wavelength electromagnetic toroidal ETG driven drift modes. A low collisionality NSTX discharge is used to scan the plasma parameter dependence on the ETG real frequency, growth rate, and electron thermal diffusivity. Electron thermal transport is discovered in the deep core region where modes are more electromagnetic in nature. Several previously reported gyrokinetic trends are reproduced, including the dependencies of density gradients, magnetic shear, $\beta$ and gradient of $\beta$ $(\betap)$, collisionality, safety factor, and toroidicity, where $\beta$ is the ratio of plasma pressure to the magnetic pressure. The electron heat diffusivity associated with the ETG mode is discovered to be on a scale consistent with the experimental diffusivity determined by power balance analysis.

ETG↗

Analysis of a theoretically optimized transonic airfoil

Numerical optimization was used in conjunction with an inviscid, full potential equation, transonic flow analysis computer code to design an upper surface contour for a conventional airfoil to improve its supercritical performance. The modified airfoil was tested in a compressible flow wind tunnel. The modified airfoil's performance was evaluated by comparison with test data for the baseline airfoil and for an airfoil developed by optimization of leading edge of the baseline airfoil. While the leading edge modification performed as expected, the upper surface re-design did not produce all of the expected performance improvements. Theoretical solutions computed using a full potential, transonic airfoil code corrected for viscosity were compared to experimental data for the baseline airfoil and the upper surface modification. These correlations showed that the theory predicted the aerodynamics of the baseline airfoil fairly well, but failed to accurately compute drag characteristics for the upper surface modification.

Lores, M. E.↗

Airfoil shape and thickness effects on transonic airloads and flutter

A transient pulse technique is used to obtain harmonic forces from a time-marching solution of the complete unsteady transonic small perturbation potential equation. The unsteady pressures and forces acting on a model of the NACA 64A010 conventional airfoil and the MBB A-3 supercritical airfoil over a range of Mach numbers are examined in detail. Flutter calculations at constant angle of attack show a similar flutter behavior for both airfoils, except for a boundary shift in Mach number associated with corresponding Mach number shift in the unsteady aerodynamic forces. Differences in the static aeroelastic twist behavior for the two airfoils are significant.

Bland, S. R.↗

SYMBOD - A computer program for the automatic generation of symbolic equations of motion for systems of hinge-connected rigid bodies

A computer program is described that can automatically generate symbolic equations of motion for systems of hinge-connected rigid bodies with tree topologies. The dynamical formulation underlying the program is outlined, and examples are given to show how a symbolic language is used to code the formulation. The program is applied to generate the equations of motion for a four-body model of the Galileo spacecraft. The resulting equations are shown to be a factor of three faster in execution time than conventional numerical subroutines.

Macala, G. A.↗

Reassessment of the theory of stimulated Raman scattering

A modification of the standard theory of stimulated Raman scattering (SRS) first proposed by Sparks (1974, 1975) is analyzed and shown to incorporate a possibly important physical effect; however, its original formulation is incorrect. The analysis is based on an exact numerical integration of the coupled equations of the modified theory, the results of which are compared with both the conventional theory of SRS and with one set of experimental data. A reformulation of the modified theory is suggested that leads to a gain which is in somewhat better agreement with the data than is the conventional theory.

Fralick, G. C.↗

The measurement of precipitation with synthetic aperture radar

The radar equation for the measurement of precipitation by SAR is identical to that for a conventional radar. The achievable synthetic beamwidth, beta(s), is proportional to sigma(v)/U, the ratio of the spread of the precipitation Doppler spectrum to the platform velocity. Thus, a small beta(s) can be achieved only with small sigma(v) or from a fast-moving vehicle such as a spacecraft. Also, the along-track resolution is variable with sigma(v) and is not known. Nevertheless, the reflectivity is measured correctly. A possible approach to the measurement of sigma(v) is noted. The C-band SAR proposed for the SIR-C mission is capable of detecting a rain rate as small as 0.5 mm/h at nadir when the beam is filled. Because the cross-track beam dimension is about 20 km wide, the use of a high-resolution microwave radiometer to correct for the unfilled beam and the variation of gain across it is suggested. Alternatively, the cross-track dimension should be decreased to no more than about 5 km by increasing the antenna width and/or decreasing the wavelength.

Atlas, David↗

Equipotential transformation of multipole systems to dipole systems

The vector sum of fields produced by a set of magnetic dipoles with different magnetic moments that are both tilted from the planetary spin axis and offset from the planetary center by different amounts can be used to completely represent a planetary magnetic field in the conventional form of spherical harmonic expansion. The scalar equipotential transformation analytically results in 24 equations that may subsequently be solved for the 24 adjustable parameters in dipole systems with the predetermined main dipole. Attention is given to this method's application to the Jovian magnetic field.

Wang, W. X.↗

A new non-eddy viscosity subgrid-scale model and its application to channel flow

To date, most large-eddy simulations (LES) have been carried out with eddy viscosity subgrid scale (SGS) models, with only a few exceptions that used the mixed model. Even though the assumptions behind Smagorinsky's model are rather stringent, it has been applied successfully to a variety of turbulent flows. This success is attributed to the ability of eddy viscosity models to drain energy from large scales, thus simulating the dissipative nature of turbulence. Most SGS models are absolutely dissipative, i.e. they remove energy from the large scales at every instant. However, SGS stresses may transfer energy back to the large scales intermittently; this reverse transfer or backscatter is especially important in geophysical flows and in transition. In a fully developed channel flow, there is reverse flow of energy from small to large scales near the walls, but eddy viscosity models are unable to account for this important feature. The dynamic localization eddy viscosity model of Ghosal et al. (1995) allows backscatter by co-evolving an auxiliary equation for the SGS energy; however, the computational cost is considerably larger than for conventional SGS models (Cabot 1994). In this report, a new non-eddy viscosity model based on local approximation of total quantities in terms of filtered ones is introduced; the scale similarity model of Bardina (1983) is a special case of this model. This procedure does not require the assumption of homogeneity, permits backscatter of energy from small to large scales, and is readily implemented in finite difference codes. The results of applying the proposed model to second order finite volume simulation of plane channel flow at high Reynolds numbers (Re(sub b) = 38,000) is described in this report. Greater emphasis is placed on the high Reynolds number flow since it provides a more rigorous test of the SGS model and its potential application. The results are compared to ones produced by the conventional and dynamic Smagorinsky models and the spectral LES of Piomelli (1993).

Shah, K. B.↗