Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational algorithms”

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 469 records · Page 26

Geomagnetic field models incorporating physical constraints on the secular variation

This proposal has been concerned with methods for constructing geomagnetic field models that incorporate physical constraints on the secular variation. The principle goal that has been accomplished is the development of flexible algorithms designed to test whether the frozen flux approximation is adequate to describe the available geomagnetic data and their secular variation throughout this century. These have been applied to geomagnetic data from both the early and middle part of this century and convincingly demonstrate that there is no need to invoke violations of the frozen flux hypothesis in order to satisfy the available geomagnetic data.

Constable, Catherine↗

A modified scheil approach for nucleation-dependent solidification pathways

As-solidified microstructures of near-eutectic alloys often contain multiple primary phases that are not expected from equilibrium phase diagrams. Such microstructures are caused by cooling-rate-dependent solidification pathways, a factor not captured by the Scheil–Gulliver model or variations thereof. Here, we present a model and algorithm that incorporate the critical nucleation undercooling for each solid phase into the Scheil–Gulliver model. We hypothesize that the non-equilibrium microstructure formation is primarily governed by a nucleation-competition mechanism. This mechanism accounts for both stable/metastable phase selection and primary-phase formation within eutectic regions driven by asymmetric nucleation barriers. The model is validated against a hypereutectic Al-Fe alloy, where it successfully reproduces the observed microstructural constituents, revealing the key dependencies of solidification microstructure on nucleation kinetics. Applicability to multicomponent systems is demonstrated through a hypereutectic Al–Fe–Si ternary alloy, where the model successfully predicts divorced eutectic microstructures and the associated oscillatory solidification pathways along univariant lines. As a result, the proposed framework establishes a nucleation-dependent computational approach for interpreting and predicting solidification microstructures.

Alloy design↗

Optimal Energy Storage Schedules for Load Leveling and Ramp Rate Control in Distribution Systems

Continued integration of distributed energy resources (DERs) into the grid, such as solar PVs, at a large-scale, contributes into the famous Duck Curve. New DER management algorithms are therefore deemed necessary to alleviate rapid variations within net load profiles of distribution systems. Herein this paper proposes a process to determine the optimal energy storage schedules for leveling the distribution circuit feederhead net load. A series of sensitivity analyses shows how the proposed method can be used to determine the optimal energy storage schedules with different capacities, state of charge requirements, and net load ramp rate limitations.

25 ENERGY STORAGE↗

rBahadur: efficient simulation of structured high-dimensional genotype data with applications to assortative mating

Existing methods for generating synthetic genotype data are ill-suited for replicating the effects of assortative mating (AM). We propose rb_dplr, a novel and computationally efficient algorithm for generating high-dimensional binary random variates that effectively recapitulates AM-induced genetic architectures using the Bahadur order-2 approximation of the multivariate Bernoulli distribution. The rBahadur R library is available through the Comprehensive R Archive Network at https://CRAN.R-project.org/package=rBahadur.

59 BASIC BIOLOGICAL SCIENCES↗

The Tidal Clock Effects of the Lunisolar Gravitational Field and the Earth’s Tidal Deformation

State-of-the-art atomic and optical clocks have the great potential to precisely test fundamental physical assumptions and enhance our understanding of nature. Their widespread applications require us to rigorously deduce the relativistic frequency shift in the framework of general relativity. One interesting question for clocks is tidal field effect which contains variously periodic variations. By introducing tidally deformed Earth, the general algorithms for clock and frequency comparisons are derived. We investigate the effects of external gravitational bodies and Earth’s tidal deformation on the ground- and space-based clocks. The orbital elements of satellite and locations of laboratory are introduced for corresponding clocks, in which the secular and long-period terms and short-period terms are subsequently presented. These effects called tidal clock effects produce the non-negligible contributions in the modern clock experiments and can be directly evaluated from our parameterized formulas. In addition, we also demonstrate for tidal clock effects the position and distance dependences, as well as periodicity. These tidal effects and dependences can provide valuable information for the clock comparison experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

A description of the University of Illinois meteor radar system and some first results

A recently established meteor radar system for observing motions in the upper atmosphere is described, which features a peak transmitter power of 5 megawatts. This large peak power should allow a much greater number of wind determination during a fixed period of time than other meteor radar stations. The essential features of the range algorithm, two velocity algorithms, and the interferometer phase difference algorithm in use are described. A typical measured variation in the radial southward wind as a function of time of day based on data from all altitudes is shown as an example, and some results of spectral analysis of the data, yielding amplitude and phase variation with altitude for significant waves, is discussed.

Geller, M. A.↗

Atmospheric correction of Nimbus-7 Coastal Zone Color Scanner imagery

The Coastal Zone Color Scanner (CZCS) on Nimbus-7 is a scanning radiometer designed to view the ocean in six spectral bands (centered at 443, 520, 550, 670, 750, and 11,500 nm) for the purpose of estimating sea surface chlorophyll and temperature distributions. In the visible bands, the atmosphere obscures the imagery to the extent that at 443 nm, at most, only 20 percent of the observed radiance originates from beneath the sea surface. Retrieving this subsurface radiance from the imagery is complicated by the highly variable nature of the aerosol's contribution. In this paper, an algorithm for the removal of these atmospheric effects from CZCS imagery is described, a preliminary application of the algorithm to an image with very strong horizontal variations in the aerosol optical thickness is presented, and retrieval of the spatial distribution of the aerosol optical thickness is discussed.

Gordon, H. R.↗

Relationships between a potential field-source surface model of the coronal magnetic field and properties of the solar wind at 1 AU

A comparison is made between the properties at the source surface of a potential field model of the solar corona and solar wind properties observed with ISEE 3 at 1 AU and extrapolated back to the sun. This is done for three consecutive Carrington rotations in the summer of 1979. The already known properties of a velocity minimum at the interplanetary current sheet and a tendency for the average velocity to increase with longitudinal distance from the current sheet are recovered. A better correlation between source surface properties and velocity is demonstrated by using the field strength on the source surface rather than longitudinal distance from the current sheet. Applying a simple algorithm to distinguish between transient and corotating interplanetary variations substantially increases the correlation. However, given a reasonable estimate of the number of degrees of freedom for the sample, only the correlation between source surface field strength and flow speed after the transients have been eliminated is statistically significant at the 1 percent level.

Suess, S. T.↗

Structural optimization using optimality criteria methods

Optimality criteria methods take advantage of some concepts as those of statically determinate or indeterminate structures, and certain variational principles of structural dynamics, to develop efficient algorithms for the sizing of structures that are subjected to stiffness-related constraints. Some of the methods and iterative strategies developed over the last decade for calculations of the Lagrange multipliers in stressand displacement-limited problems, as well as for satisfying the appropriate optimality criterion, are discussed. The application of these methods are illustrated by solving problems with stress and displacement constraints.

Khot, N. S.↗

Transformation of two and three-dimensional regions by elliptic systems

Efforts in transferring computational work from the LRC computer to the IRIS Graphics Workstation at MSU are reported and the computation of a conservative solution of a simple hyperbolic equation on an overlapping grid is discussed. Several conclusions concerning computations on overlapping grids are apparent. Problems only occur when there is a major difference in grid spacing on the individual component grids. In the case of hyperbolic equations, it is necessary that both interpolation and extrapolation be applied at the grid boundaries. When interpolated values are used at outflow boundary points, excessive oscillations in the numerical solution may be the result. The same conclusions would be valid for more complicated systems of hyperbolic equations such as the Euler equations for inviscid flow. Some of the solution values would be extrapolated at the overlap boundary, the exact number depending on the number of characteristics pointing out of the overlap region. It is also possible that similar boundary conditions may be needed for some parabolic equations such as high Reynolds number viscous flow equations. Efforts were also expended on the development of three-dimensional conservative interpolation procedures. Finally, the investigation of grid smoothing procedures were initiated during this reporting period. It was decided that the first grid smoothing algorithms will be based on the concepts of variational grid generation.

Mastin, C. W.↗

Parameter learning for performance adaptation

A parameter learning method is introduced and used to broaden the region of operability of the adaptive control system of a flexible space antenna. The learning system guides the selection of control parameters in a process leading to optimal system performance. A grid search procedure is used to estimate an initial set of parameter values. The optimization search procedure uses a variation of the Hooke and Jeeves multidimensional search algorithm. The method is applicable to any system where performance depends on a number of adjustable parameters. A mathematical model is not necessary, as the learning system can be used whenever the performance can be measured via simulation or experiment. The results of two experiments, the transient regulation and the command following experiment, are presented.

Peek, Mark D.↗

Ice tracking techniques, implementation, performance, and applications

Present techniques of ice tracking make use both of cross-correlation and of edge tracking, the former being more successful in heavy pack ice, the latter being critical for the broken ice of the pack margins. Algorithms must assume some constraints on the spatial variations of displacements to eliminate fliers, but must avoid introducing any errors into the spatial statistics of the measured displacement field. We draw our illustrations from the implementation of an automated tracking system for kinematic analyses of ERS-1 and JERS-1 SAR imagery at the University of Alaska - the Alaska SAR Facility's Geophysical Processor System. Analyses of the ice kinematic data that might have some general interest to analysts of cloud-derived wind fields are the spatial structure of the fields, and the evaluation and variability of average deformation and its invariants: divergence, vorticity and shear. Many problems in sea ice dynamics and mechanics can be addressed with the kinematic data from SAR.

Rothrock, D. A.↗

Unified approach for incompressible flows

An unified approach for solving both compressible and incompressible flows was investigated in this study. The difference in CFD code development between incompressible and compressible flows is due to the mathematical characteristics. However, if one can modify the continuity equation for incompressible flows by introducing pseudocompressibility, the governing equations for incompressible flows would have the same mathematical characters as compressible flows. The application of a compressible flow code to solve incompressible flows becomes feasible. Among numerical algorithms developed for compressible flows, the Centered Total Variation Diminishing (CTVD) schemes possess better mathematical properties to damp out the spurious oscillations while providing high-order accuracy for high speed flows. It leads us to believe that CTVD schemes can equally well solve incompressible flows. In this study, the governing equations for incompressible flows include the continuity equation and momentum equations. The continuity equation is modified by adding a time-derivative of the pressure term containing the artificial compressibility. The modified continuity equation together with the unsteady momentum equations forms a hyperbolic-parabolic type of time-dependent system of equations. The continuity equation is modified by adding a time-derivative of the pressure term containing the artificial compressibility. The modified continuity equation together with the unsteady momentum equations forms a hyperbolic-parabolic type of time-dependent system of equations. Thus, the CTVD schemes can be implemented. In addition, the boundary conditions including physical and numerical boundary conditions must be properly specified to obtain accurate solution. The CFD code for this research is currently in progress. Flow past a circular cylinder will be used for numerical experiments to determine the accuracy and efficiency of the code before applying this code to more specific applications.

Chang, Tyne-Hsien↗

Simulation of Stagnation Region Heating in Hypersonic Flow on Tetrahedral Grids

Hypersonic flow simulations using the node based, unstructured grid code FUN3D are presented. Applications include simple (cylinder) and complex (towed ballute) configurations. Emphasis throughout is on computation of stagnation region heating in hypersonic flow on tetrahedral grids. Hypersonic flow over a cylinder provides a simple test problem for exposing any flaws in a simulation algorithm with regard to its ability to compute accurate heating on such grids. Such flaws predominantly derive from the quality of the captured shock. The importance of pure tetrahedral formulations are discussed. Algorithm adjustments for the baseline Roe / Symmetric, Total-Variation-Diminishing (STVD) formulation to deal with simulation accuracy are presented. Formulations of surface normal gradients to compute heating and diffusion to the surface as needed for a radiative equilibrium wall boundary condition and finite catalytic wall boundary in the node-based unstructured environment are developed. A satisfactory resolution of the heating problem on tetrahedral grids is not realized here; however, a definition of a test problem, and discussion of observed algorithm behaviors to date are presented in order to promote further research on this important problem.

Gnoffo, Peter A.↗

Data Fusion for Earth Science Remote Sensing

Beginning in 2004, NASA has supported the development of an international network of ground-based remote sensing installations for the measurement of greenhouse gas columns. This collaboration has been successful and is currently used in both carbon cycle investigations and in the efforts to validate the GOSAT space-based column observations of CO2 and CH4. With the support of a grant, this research group has established a network of ground-based column observations that provide an essential link between the satellite observations of CO2, CO, and CH4 and the extensive global in situ surface network. The Total Carbon Column Observing Network (TCCON) was established in 2004. At the time of this report seven sites, employing modern instrumentation, were operational or were expected to be shortly. TCCON is expected to expand. In addition to providing the most direct means of tying the in situ and remote sensing data sets together, TCCON provides a means of testing the retrieval algorithms of SCIAMACHY and GOSAT over the broadest variation in atmospheric state. TCCON provides a critically maintained and long timescale record for identification of temporal drift and spatial bias in the calibration of the space-based sensors. Finally, the global observations from TCCON are improving our understanding of how to use column observations to provide robust estimates of surface exchange of C02 and CH4 in advance of the launch of OCO and GOSAT. TCCON data are being used to better understand the impact of both regional fluxes and long-range transport on gradients in the C02 column. Such knowledge is essential for identifying the tools required to best use the space-based observations. The technical approach and methodology of retrieving greenhouse gas columns from near-IR solar spectra, data quality and process control are described. Additionally, the impact of and relevance to NASA of TCCON and satellite validation and carbon science are addressed.

data fusion↗

Observability and Estimation of Distributed Space Systems via Local Information-Exchange Networks

In this work, we develop an approach to formation estimation by explicitly characterizing formation's system-theoretic attributes in terms of the underlying inter-spacecraft information-exchange network. In particular, we approach the formation observer/estimator design by relaxing the accessibility to the global state information by a centralized observer/estimator- and in turn- providing an analysis and synthesis framework for formation observers/estimators that rely on local measurements. The noveltyof our approach hinges upon the explicit examination of the underlying distributed spacecraft network in the realm of guidance, navigation, and control algorithmic analysis and design. The overarching goal of our general research program, some of whose results are reported in this paper, is the development of distributed spacecraft estimation algorithms that are scalable, modular, and robust to variations inthe topology and link characteristics of the formation information exchange network. In this work, we consider the observability of a spacecraft formation from a single observation node and utilize the agreement protocol as a mechanism for observing formation states from local measurements. Specifically, we show how the symmetry structure of the network, characterized in terms of its automorphism group, directly relates to the observability of the corresponding multi-agent system The ramification of this notion of observability over networks is then explored in the context of distributed formation estimation.

observability↗

Lunar calibration and performance assessments of the NOAA-20 VIIRS reflective solar bands

The Moon has played a vital role in the radiometric stability monitoring of NASAs Earth Observing System sensors, such as MODIS on the Terra and Aqua spacecrafts. The lunar calibration methodologies developed for the two MODIS instruments were later extended to apply to the VIIRS instrument on the Suomi-NPP spacecraft. The follow-on VIIRS instrument on the NOAA-20 (N20) spacecraft (previously JPSS-1) has been successfully operating since its launch in November 2017. Fifteen VIIRS spectral bands are in the reflective solar spectral region, covering wavelengths from 0.4 to 2.3 μm. Similar to the previous MODIS and VIIRS instruments, the N20 VIIRS views the Moon via a spacecraft roll maneuver on a near-monthly basis at a confined phase angle range (-51.5° to -50.5°) to facilitate radiometric stability monitoring with reduced uncertainties. In this paper, we briefly present the lunar calibration methodology and also report the derived instrument gain parameters from the regularly scheduled lunar calibrations. The USGS Robotic Lunar Observatory (ROLO) model is used to provide the predicted lunar irradiance for each lunar measurement to account for the variation in the observation geometry. The spatial characterization algorithms using the Moon, previously developed for MODIS and SNPP VIIRS, have also been extended to N20 VIIRS to track its on-orbit band-to-band registration (BBR).

VIIRS↗

Application of an efficient generator-coordinate subspace-selection algorithm to neutrinoless double- β decay

The generator coordinate method begins with the variational construction of a set of nonorthogonal mean-field states that span a subspace of the full many-body Hilbert space. These states are then often projected onto states with good quantum numbers to restore symmetries, leading to a set with members that can be similar to one another, and it is sometimes possible to reduce this set without greatly affecting results. Here, we propose a greedy algorithm that we call the energy-transition-orthogonality procedure (ENTROP) to select subsets of important states. As applied here, the approach selects on the basis of diagonal energy, orthogonality, and contribution to the matrix element that governs neutrinoless double-β decay. We present both shell-model and preliminary ab initio calculations of this matrix element for the decay of 76Ge, with quadrupole deformation parameters and the isoscalar pairing strength as generator coordinates. ENTROP converges quickly, reducing significantly the number of basis states needed for an accurate calculation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗