Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiple time stepping”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Temporal averaging of phase measurements in the presence of spurious phase drift - Application to phase-stepped real-time holographic interferometry

A technique that compensates for low spatial frequency spurious phase changes during an interference experiment is developed; it permits temporal averaging of multiple-phase measurements, made before and after object displacement. The method is tested with phase-stepped real-time holographic interferometry applied to cantilever bending of a piezoelectric bimorph ceramic. Results indicate that temporal averaging of the corrected data significantly reduces the white noise in a phase measurement without incurring systematic errors or sacrificing spatial resolution. White noise is reduced from 3 deg to less than 1 deg using these methods.

Ovryn, B.↗
DOE OSTI2022

Robust bidding strategy for aggregation of distributed prosumers in flexiramp market

Distributed prosumers (DPs) are the grid customers that own energy production/storage assets. Due to the flexibility and fast response of their assets, they can procure ancillary service products (ASP) in the wholesale market. An appealing ASP offered by California ISO in the real-time market (RTM) is flexiramp for which market participants do not submit direct offers, and the compensation is based on their energy opportunity costs. Here in this report, we propose a bidding strategy model for DP aggregator participation in the RTM considering energy and flexiramp. First, we develop a risk-averse optimization to determine the optimal energy and reserve product to trade in day-ahead market while considering proper amounts of flexiramp to trade in the RTM. In the RTM, to obtain optimal amounts of energy and flexiramp, the aggregator must submit hourly multi-level price-quantity energy bids for multiple RTM intervals with 15 min time-steps. On this basis, we propose a robust hourly economic bidding strategy model that determines the optimal energy bids in the RTM. We develop an adjustable robust counterpart of the model to address the RTM energy and flexiramp price uncertainties. The simulation results justify the efficacy of our proposed framework in gaining profits from the wholesale market.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Space-Time Block Preconditioning for Incompressible Flow

Parallel-in-time methods have become increasingly popular in the simulation of time-dependent numerical PDEs, allowing for the efficient use of additional message passing interface processes when spatial parallelism saturates. Most methods treat the solution and parallelism in space and time separately. In contrast, all-at-once methods solve the full space-time system directly, largely treating time as simply another spatial dimension. All-at-once methods offer a number of benefits over separate treatment of space and time, most notably significantly increased parallelism and faster time to solution (when applicable). However, the development of fast, scalable all-at-once methods has largely been limited to time-dependent (advection-)diffusion problems. This paper introduces the concept of space-time block preconditioning for the all-at-once solution of incompressible flow. By extending well-known concepts of spatial block preconditioning to the space-time setting, we develop a block preconditioner whose application requires the solution of a space-time (advection-)diffusion equation in the velocity block, coupled with a pressure Schur complement approximation consisting of independent spatial solves at each time-step, and a space-time matrix-vector multiplication. The new method is tested on four classical models in incompressible flow. Finally, the results indicate perfect scalability in refinement of spatial and temporal mesh spacing, perfect scalability in nonlinear Picard iteration count when applied to a nonlinear Navier--Stokes problem, and minimal overhead in terms of number of preconditioner applications compared with sequential time-stepping.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Step-patterned survivorship curves: Mortality and loss of equilibrium responses to high temperature and food restriction in juvenile rainbow trout ( Oncorhynchus mykiss )

While survivorship curves typically exhibit smooth declines over time, step-patterned curves can occur with multiple stressors within a life stage. To explore this process, we examined the effects of heat (24°C) and food restriction on juvenile rainbow trout (Oncorhynchus mykiss Walbaum) in challenge experiments. We observed step-patterned survivorship curves determined by mortality and loss of equilibrium (LOE) endpoints. To examine the cause of heterogeneity in the stress responses from early to late mortality and LOE, we measured indices of energetic reserves. The step transition in the survivorship curves, the peak mortality rates, and start of when individuals reached a critical energetic threshold (14% dry mass; 4.0 kJ∙g -1 energy) all occurred at around days 10–15 of the challenge. The coherence in these temporal patterns suggest heterogeneity in the cohort stress responses, in which an early subgroup died from heat stress and a late subgroup died from starvation. Thus, their endpoint sensitivities resulted in step-patterned survivorship curves. We discuss the implications of the study for understanding effects of multiple stressors on population heterogeneity and note the possible significance of stress response selection under climate change in which heat stress and food limitations occur in concert.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2014

Method and Apparatus for Automated Isolation of Nucleic Acids from Small Cell Samples

RNA isolation is a ubiquitous need, driven by current emphasis on microarrays and miniaturization. With commercial systems requiring 100,000 to 1,000,000 cells for successful isolation, there is a growing need for a small-footprint, easy-to-use device that can harvest nucleic acids from much smaller cell samples (1,000 to 10,000 cells). The process of extraction of RNA from cell cultures is a complex, multi-step one, and requires timed, asynchronous operations with multiple reagents/buffers. An added complexity is the fragility of RNA (subject to degradation) and its reactivity to surface. A novel, microfluidics-based, integrated cartridge has been developed that can fully automate the complex process of RNA isolation (lyse, capture, and elute RNA) from small cell culture samples. On-cartridge cell lysis is achieved using either reagents or high-strength electric fields made possible by the miniaturized format. Traditionally, silica-based, porous-membrane formats have been used for RNA capture, requiring slow perfusion for effective capture. In this design, high efficiency capture/elution are achieved using a microsphere-based "microfluidized" format. Electrokinetic phenomena are harnessed to actively mix microspheres with the cell lysate and capture/elution buffer, providing important advantages in extraction efficiency, processing time, and operational flexibility. Successful RNA isolation was demonstrated using both suspension (HL-60) and adherent (BHK-21) cells. Novel features associated with this development are twofold. First, novel designs that execute needed processes with improved speed and efficiency were developed. These primarily encompass electric-field-driven lysis of cells. The configurations include electrode-containing constructs, or an "electrode-less" chip design, which is easy to fabricate and mitigates fouling at the electrode surface; and the "fluidized" extraction format based on electrokinetically assisted mixing and contacting of microbeads in a shape-optimized chamber. A secondary proprietary feature is in the particular layout integrating these components to perform the desired operation of RNA isolation. Apart from a novel functional capability, advantages of the innovation include reduced or eliminated use of toxic reagents, and operator-independent extraction of RNA.

Sundaram, Shivshankar↗
NASA NTRSDate not supplied

Preliminary Considerations for Microwave Consolidation/Sintering of Lunar Regolith Simulant

As NASA prepares to establish permanent habitats on the Moon, a significant first step is to be able to land multiple times in the same area. As was consistently shown during the Apollo program, the very fine granular structure of the lunar regolith (the Moon’s “soil”)poses significant physical and health challenges [1]. One of the most concerning is the hyper-velocity lunar surface ejecta that results from the engine exhaust that exits the rocket as it lands and takes off [2]. It has been determined that to mitigate this, the regolith must be consolidated.

microwave↗
NASA NTRS2016

Quantifying Surface Water Dynamics at 30 Meter Spatial Resolution in the North American High Northern Latitudes 1991-2011

The availability of a dense time series of satellite observations at moderate (30 m) spatial resolution is enabling unprecedented opportunities for understanding ecosystems around the world. A time series of data from Landsat was used to generate a series of three maps at decadal time step to show how surface water has changed from 1991 to 2011 in the high northern latitudes of North America. Previous attempts to characterize the change in surface water in this region have been limited in either spatial or temporal resolution, or both. This series of maps was generated for the NASA Arctic and Boreal Vulnerability Experiment (ABoVE), which began in fall 2015. These maps show a nominal extent of surface water by using multiple observations to make a single map for each time step. This increases the confidence that any detected changes are related to climate or ecosystem changes not simply caused by short duration weather events such as flood or drought. The methods and comparison to other contemporary maps of the region are presented here. Initial verification results indicate 96% producer accuracy and 54% user accuracy when compared to 2-m resolution World View-2 data. All water bodies that were omitted were one Landsat pixel or smaller, hence below detection limits of the instrument.

Carroll, Mark↗
DOE OSTI2023

TransformerG2G: Adaptive time-stepping for learning temporal graph embeddings using transformers

Dynamic graph embedding has emerged as a very effective technique for addressing diverse temporal graph analytic tasks (i.e., link prediction, node classification, recommender systems, anomaly detection, and graph generation) in various applications. Such temporal graphs exhibit heterogeneous transient dynamics, varying time intervals, and highly evolving node features throughout their evolution. Hence, incorporating long-range dependencies from the historical graph context plays a crucial role in accurately learning their temporal dynamics. In this paper, we develop a graph embedding model with uncertainty quantification, TransformerG2G, by exploiting the advanced transformer encoder to first learn intermediate node representations from its current state (t) and previous context (over timestamps [t–1,t–l], l is the length of context). Moreover, we employ two projection layers to generate lower-dimensional multivariate Gaussian distributions as each node's latent embedding at timestamp t. We consider diverse benchmarks with varying levels of "novelty" as measured by the TEA (Temporal Edge Appearance) plots. Here, our experiments demonstrate that the proposed TransformerG2G model outperforms conventional multi-step methods and our prior work (DynG2G) in terms of both link prediction accuracy and computational efficiency, especially for high degree of novelty. Furthermore, the learned time-dependent attention weights across multiple graph snapshots reveal the development of an automatic adaptive time stepping enabled by the transformer. Importantly, by examining the attention weights, we can uncover temporal dependencies, identify influential elements, and gain insights into the complex interactions within the graph structure. For example, we identified a strong correlation between attention weights and node degree at the various stages of the graph topology evolution.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1987

DYCAST: A finite element program for the crash analysis of structures

DYCAST is a nonlinear structural dynamic finite element computer code developed for crash simulation. The element library contains stringers, beams, membrane skin triangles, plate bending triangles and spring elements. Changing stiffnesses in the structure are accounted for by plasticity and very large deflections. Material nonlinearities are accommodated by one of three options: elastic-perfectly plastic, elastic-linear hardening plastic, or elastic-nonlinear hardening plastic of the Ramberg-Osgood type. Geometric nonlinearities are handled in an updated Lagrangian formulation by reforming the structure into its deformed shape after small time increments while accumulating deformations, strains, and forces. The nonlinearities due to combined loadings are maintained, and stiffness variation due to structural failures are computed. Numerical time integrators available are fixed-step central difference, modified Adams, Newmark-beta, and Wilson-theta. The last three have a variable time step capability, which is controlled internally by a solution convergence error measure. Other features include: multiple time-load history tables to subject the structure to time dependent loading; gravity loading; initial pitch, roll, yaw, and translation of the structural model with respect to the global system; a bandwidth optimizer as a pre-processor; and deformed plots and graphics as post-processors.

Pifko, A. B.↗
NASA NTRS2020

Image Moment-Based Extended Object Tracking for Complex Motions

A novel image moment-based model for shape estimation and tracking of an extended target moving with a complex trajectory is presented. The proposed extended object tracking algorithm is based on multiple noisy measurement points sampled from the target at each time step. The shape of the object, approximated by an ellipse, is estimated using a combination of image moments. Dynamic models of image moments for constant velocity and coordinated turn motions are mathematically derived. An unscented Kalman filter - interacting multiple model (UKF-IMM) method is used to track the object and estimate its shape. A likelihood function based on average log-likelihood is derived for the IMM filter. Simulation results of the proposed UKF-IMM algorithm with the image momentbased models are presented that show the estimation of the shape of the object moving in a complex trajectory. The intersection over union (IoU), and the root mean square errors (RMSEs) of the position and velocity of the centroid of the ellipse are used as metrics. The comparison results of the proposed algorithm with a benchmark algorithm from literature based on the IoU and RMSE metrics are presented.

Extended object tracking↗
DOE OSTI2022

Introduction of the Adding and Doubling Method for Solving Bateman Equations for Nuclear Fuel Depletion

This paper introduces and evaluates the Adding and Doubling Method (ADM) for solving the Bateman equations for depletion systems with varying numbers of nuclides and compares it to the Chebyshev Rational Approximation Method (CRAM), both implemented in the reactor physics analysis application Griffin. ADM, when applied to the Crank-Nicolson Finite Difference method, can produce results comparable in accuracy and precision to CRAM with comparable run times for systems with 35 or 297 nuclides. For systems with more than 300 nuclides, the matrix-matrix operations required by ADM are significantly more costly than the matrix-vector operations required by CRAM, making CRAM the more efficient method for systems with large numbers of nuclides. ADM is an accurate method that maintains other advantages over CRAM in that it does not depend on pre-generated coefficients or require complex number operations. ADM also manages to outperform CRAM by a factor of more than 250 in terms of run time for depletion systems that require multiple Bateman solves while the depletion matrix and time step size remain constant over all depletion intervals.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1989

A time accurate finite volume method for propulsion chamber flows

An implicit three-dimensional time-accurate method for propulsion chamber flows is proposed which uses line Gauss-Seidel relaxation and multiple axial sweeps for the convergence of each time step. The general time-integration algorithm employed includes such schemes as the Euler implicit method. The results of spatial and temporal accuracy tests reveal that Roe's (1981) flux difference splitting provides excellent tracking of acoustic wave speeds. In comparison with other methods, no low mean flow Mach number convergence limitation or Courant number stabilization restriction is observed.

Beddini, R. A.↗
NASA NTRS1992

A comparison of artificial compressibility and fractional step methods for incompressible flow computations

We have applied and compared the efficiency and accuracy of two commonly used numerical methods for the solution of Navier-Stokes equations. The artificial compressibility method augments the continuity equation with a transient pressure term and allows one to solve the modified equations as a coupled system. Due to its implicit nature, one can have the luxury of taking a large temporal integration step at the expense of higher memory requirement and larger operation counts per step. Meanwhile, the fractional step method splits the Navier-Stokes equations into a sequence of differential operators and integrates them in multiple steps. The memory requirement and operation count per time step are low, however, the restriction on the size of time marching step is more severe. To explore the strengths and weaknesses of these two methods, we used them for the computation of a two-dimensional driven cavity flow with Reynolds number of 100 and 1000, respectively. Three grid sizes, 41 x 41, 81 x 81, and 161 x 161 were used. The computations were considered after the L2-norm of the change of the dependent variables in two consecutive time steps has fallen below 10(exp -5).

Chan, Daniel C.↗
DOE OSTI2020

An Overview of the Circular Economy Lifecycle Assessment and VIsualization (CELAVI) Framework

A circular economy emphasizes the efficient use of all resources (e.g., materials, land, water). Despite anticipated overall benefits to society, the transition to a circular economy is likely to create regional differences in impacts. Current tools are unable to fully evaluate these potential externalities, which will be important for informing research prioritization and regional decision making. This presentation will review and contrast nine current impact assessment methods and describe their granularity, scope, data requirements, and capabilities within the context of identifying regional and sectoral transfers of impacts that could result from transitioning to a circular economy for energy systems. The review of current approaches highlights the need for a new hybrid circularity assessment framework that leverages multiple methods. We describe one such approach: the Circular Economy Lifecycle Assessment and VIsualization (CELAVI) framework. CELAVI uses system dynamics to model material flows for multiple circular economy pathways, network theory to track the spatial and sectoral flow of functional units across a graph, and discrete event simulation to to step through time and evaluate lifecycle assessment data at each time step. The framework is designed to be flexible and scalable enough to accommodate multiple energy materials and multiple energy technologies. CELAVI's hybridization of multiple existing methods yields additional capabilities, which may help answer questions about how flows between the technosphere and ecosphere may evolve if the circularity of energy systems is modified.

circular economy↗
DOE OSTI2024

Plant diversity and functional identity drive grassland rhizobacterial community responses after 15 years of CO 2 and nitrogen enrichment

Abstract Improved understanding of bacterial community responses to multiple environmental filters over long time periods is a fundamental step to develop mechanistic explanations of plant–bacterial interactions as environmental change progresses. This is the first study to examine responses of grassland root‐associated bacterial communities to 15 years of experimental manipulations of plant species richness, functional group and factorial enrichment of atmospheric CO 2 (eCO 2 ) and soil nitrogen (+N). Across the experiment, plant species richness was the strongest predictor of rhizobacterial community composition, followed by +N, with no observed effect of eCO 2 . Monocultures of C 3 and C 4 grasses and legumes all exhibited dissimilar rhizobacterial communities within and among those groups. Functional responses were also dependent on plant functional group, where N 2 ‐fixation genes, NO 3− ‐reducing genes and P‐solubilizing predicted gene abundances increased under resource‐enriched conditions for grasses, but generally declined for legumes. In diverse plots with 16 plant species, the interaction of eCO 2 +N altered rhizobacterial composition, while +N increased the predicted abundance of nitrogenase‐encoding genes, and eCO 2 +N increased the predicted abundance of bacterial P‐solubilizing genes. Synthesis : Our findings suggest that rhizobacterial community structure and function will be affected by important global environmental change factors such as eCO 2 , but these responses are primarily contingent on plant species richness and the selective influence of different plant functional groups.

Revillini, Daniel↗
DOE OSTI2022

Stability analysis of Monte Carlo burnup calculation based on PWR fuel rods

The existing Monte Carlo burnup algorithms may give some unstable results when using multiple burnable cells during burnup calculations. A very small time step is required to maintain the stability of the system, which pose some limits for the realistic applications. In this paper, we firstly repeat the instability issue by using a Monte Carlo code based on a typical PWR fuel rod model. Then, a theoretical stability analysis was carried out by using the perturbation theory in order to get a better understanding about the property of system stability. We determine a perturbation matrix which relates some main influence factors, such as the time step and the height of fuel rod, it can help people estimate the stability of burnup calculation system and the results are quantitatively consistent with the Monte Carlo code. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRSDate not supplied

Permanent Magnetic Synchronous Motor (PMSM) Model Development

The Advanced Air Transport Technology (AATT) Project seeks to enhance the capabilities of fixed-wing subsonic transport through improved energy efficiency and environmental compatibility. One element of this effort is the use of high efficiency PMSM systems. Understanding the behavior and limitations of the se motors allows the National Aeronautics and Space Administration (NASA) to collaborate with and inform its partners worldwide. PMSM modeling is achieved using MATLAB®-Simulink® (The Math Works, Inc., Natick, Massachusetts). Through the analysis of an existing high-fidelity NASA Electrical Aircraft Testbed (NEAT) model, simplified models of varying fidelity can be developed with the goal of building a PMSM model capable of running in real-time for use in a piloted simulation environment. By utilizing system identification methods, it has been demonstrated that internal electrical components like the inverter can be represented by a simple variable gain thus greatly reducing the required simulation step size and subsequently decreasing model run-time by multiple orders of magnitude

motor controls↗
NASA NTRS2017

Interactive Visualization of High-Dimensional Petascale Ocean Data

We describe an application for interactive visualization of 5 petabytes of time-varying multivariate data from a high-resolution global ocean circulation model. The input data are 10311 hourly (ocean time) time steps of various 2D and 3D fields from a 22-billion point 1/48- degree “lat-lon cap” configuration of the MIT General Circulation Model (MITgcm). We map the global horizontal model domain onto our 128-screen (8x16) tiled display wall to produce a canonical tiling with approximately one MITgcm grid point per display pixel, and using this tiling we encode the entire time series for multiple native and computed scalar quantities at a collection of ocean depths. We reduce disk bandwidth requirements by converting the model’s floating point data to 16-bit fixed point values, and compressing those values with a lossless video encoder, which together allow synchronized playback at 24 time steps per second across all 128 displays. The application allows dynamic assignment of any two encoded tiles to any display, and has multiple interfaces for quickly specifying various orderly arrangements of tiles. All subsequent rendering is done on the fly, with run time control of colormaps, transfer functions, histogram equalization, and labeling. The two data streams on each screen can be rendered independently and combined in various ways, including blending, differencing, horizontal/ vertical wipes, and checkerboarding. The two data streams on any screen can optionally be displayed as a scatterplot in their joint attribute space. All scatterplots and map-view plots from the same x/y location and depth are linked so they all show the current brushable selection. Ocean scientists have used the system, and have found previously unidentified features in the data.

Interactive↗