Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear scaling”

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 1,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Interpolated Sounding and Gridded Sounding Value-Added Products

Standard Atmospheric Radiation Measurement (ARM) Climate Research Facility sounding files provide atmospheric state data in one dimension of increasing time and height per sonde launch. Many applications require a quick estimate of the atmospheric state at higher time resolution. The INTERPOLATEDSONDE (i.e., Interpolated Sounding) Value-Added Product (VAP) transforms sounding data into continuous daily files on a fixed time-height grid, at 1-minute time resolution, on 332 levels, from the surface up to a limit of approximately 40 km. The grid extends that high so the full height of soundings can be captured; however, most soundings terminate at an altitude between 25 and 30 km, above which no data is provided. Between soundings, the VAP linearly interpolates atmospheric state variables in time for each height level. In addition, INTERPOLATEDSONDE provides relative humidity scaled to microwave radiometer (MWR) observations.The INTERPOLATEDSONDE VAP, a continuous time-height grid of relative humidity-corrected sounding data, is intended to provide input to higher-order products, such as the Merged Soundings (MERGESONDE; Troyan 2012) VAP, which extends INTERPOLATEDSONDE by incorporating model data. The INTERPOLATEDSONDE VAP also is used to correct gaseous attenuation of radar reflectivity in products such as the KAZRCOR VAP.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Multiresolution GPC-Structured Control of a Single-Loop Cold-Flow Chemical Looping Testbed

Chemical looping is a near-zero emission process for generating power from coal. It is based on a multi-phase gas-solid flow and has extremely challenging nonlinear, multi-scale dynamics with jumps, producing large dynamic model uncertainty, which renders traditional robust control techniques, such as linear parameter varying H ∞ design, largely inapplicable. This process complexity is addressed in the present work through the temporal and the spatiotemporal multiresolution modeling along with the corresponding model-based control laws. Namely, the nonlinear autoregressive with exogenous input model structure, nonlinear in the wavelet basis, but linear in parameters, is used to identify the dominant temporal chemical looping process dynamics. The control inputs and the wavelet model parameters are calculated by optimizing a quadratic cost function using a gradient descent method. The respective identification and tracking error convergence of the proposed self-tuning identification and control schemes, the latter using the unconstrained generalized predictive control structure, is separately ascertained through the Lyapunov stability theorem. The rate constraint on the control signal in the temporal control law is then imposed and the control topology is augmented by an additional control loop with self-tuning deadbeat controller which uses the spatiotemporal wavelet riser dynamics representation. The novelty of this work is three-fold: (1) developing the self-tuning controller design methodology that consists in embedding the real-time tunable temporal highly nonlinear, but linearly parametrizable, multiresolution system representations into the classical rate-constrained generalized predictive quadratic optimal control structure, (2) augmenting the temporal multiresolution loop by a more complex spatiotemporal multiresolution self-tuning deadbeat control loop, and (3) demonstrating the effectiveness of the proposed methodology in producing fast recursive real-time algorithms for controlling highly uncertain nonlinear multiscale processes. The latter is shown through the data from the implemented temporal and augmented spatiotemporal solutions of a difficult chemical looping cold flow tracking control problem.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

FEOTS v0.0.0: a new offline code for the fast equilibration of tracers in the ocean

Abstract. In this paper we introduce a new software framework for the offline calculation of tracer transport in the ocean. The Fast Equilibration of Ocean Tracers Software (FEOTS) is an end-to-end set of tools to efficiently calculate tracer distributions on a global or regional sub-domain using transport operators diagnosed from a comprehensive ocean model. To the best of our knowledge, this is the first application of a transport matrix model to an eddying ocean state. While a Newton–Krylov-based equilibration capability is still under development and not presented here, we demonstrate in this paper the transient modeling capabilities of FEOTS in an application focused on the Argentine Basin, where intense eddy activity and the Zapiola Anticyclone lead to strong mixing of water masses. The demonstration illustrates progress in developing offline passive tracer simulation capabilities, while highlighting the challenges of the impulse response functions approach in capturing tracer transports by a non-linear advection scheme. Our future work will focus on improving the computational efficiency of the code to reduce time-to-solution, using different basis functions to better represent non-linear advection operators, applying FEOTS to a parent model with unstructured grids (Ocean Model for Prediction Across Scales, MPAS-Ocean), and fully implementing a Newton–Krylov steady-state solver.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

[Re] Drivers of evapotranspiration from boreal wildfires

Computational reproducibility is a difficult challenge across science. I attempted to use R 3.6.1 to reproduce linear model fits, done originally using v2.6.0 for a 2009 paper on the drivers of large-scale forest evapotranspiration after wildfire. Model outputs were largely identical, aside from minor formatting changes, except for one–out of 12 total– regression in which the median residual value changed very slightly (in the sixth decimal place). I suggest that this essentially successful reproducibility is due to the relative simplicity of the script, its use of only base R functions, and R’s historically conservative approach to breaking changes.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1978

Dynamic model verification of large structural systems

The objective of the present methodology is two-fold: (1) to process test data obtained from either modal survey tests, or slow sine-sweep tests, to extract a set of orthogonal modes best matching the test data while being commensurate with the dynamic model, and (2) to modify submatrices of the dynamic model mass and stiffness matrices to adjust the model to best fit the test data. The method has been implemented using a linear statistical sequential estimator for computation on a CDC computer. Demonstration problems involving Space Shuttle quarter-scale vibration test data and dynamic models have been run. This paper will discuss the general methodology and experience to date.

Lee, L. T.↗
NASA NTRS1984

Numerical experiments with a general circulation model concerning the distribution of ozone in the stratosphere

The distribution of ozone below 60 km altitude has been simulated in two experiments employing a nine-layer quasi-geostrophic spectral model and linear parameterization of ozone photochemistry, the first of which included thermal and orographic forcing of the planetary scale waves, while the second omitted it. The first experiment exhibited a high latitude winter ozone buildup which was due to a Brewer-Dodson circulation forced by large amplitude (planetary scale) waves in the winter lower stratosphere. Photochemistry was also found to be important down to lower altitudes (20 km) in the summer stratosphere than had previously been supposed.

Kurzeja, R. J.↗
NASA NTRS1985

CO2 impact ionization-driven plasma instability observed by Pioneer Venus Orbiter at Periapsis

Observations of enhanced ac electric field noise about Pioneer Venus periapsis are shown to be related to spacecraft-generated impact ionization of the ambient CO2. The frequency of the electric field noise is found to peak in the vicinity of the CO2(+) ion plasma frequency and to closely follow the form of the neutral CO2 density profile. When the electric field noise in all channels is normalized by the square root of the CO2 number density, the ratio is constant. Since the impact electron density measured by the Pioneer Venus Langmuir probe, is observed to scale directly with the neutral CO2, the growth of the electric field amplitude is found to be linear in time with a growth rate proportional to the CO2(+) ion plasma frequency. On the basis of these results the impact ionization-driven instability is shown to be the ion acoustic instability. Implications for the lack of observations by Pioneer Venus of reflected-O(+)-driven instabilities, as have been proposed for the space shuttle, are discussed.

Curtis, S. A.↗
NASA NTRS1986

Toward the application of the Kalman filter to regional open ocean modeling

A partial differential equation model is defined for ocean meteorological prediction and synoptic analysis. The Kalman filter used for data assimilation is described and applied to the one-dimensional linear barotropic quasi-geostrophic model with periodic and open boundary conditions. The model accounts for eddy scale dynamics in the ocean. The assumptions made in the forecast model are discussed, along with comparisons of the error variances expected with the filter and from an objective analysis method. The effectiveness of the Kalman filter is demonstrated and subsequent efforts to extend the filter to two dimensions are indicated.

Miller, R. N.↗
NASA NTRS1989

Accretion by magnetic neutron stars

The capture of material from the wind or from the atmosphere or envelope of a binary companion star is described and the resulting types of accretion flows are discussed. An explanation is given for the formation of a magnetosphere around the neutron star. Conditions for stable flow and for angular and linear momentum conservation are discussed within the framework of neutron star accretion and rough estimates of the scale of the magnetosphere are obtained. The beat-frequency magnetospheric model of quasi-periodic intensity oscillations is presented.

Lamb, F. K.↗
NASA NTRS1994

Subpixel resolution from multiple images

Multiple images taken from similar locations and under similar lighting conditions contain similar, but not identical, information. Slight differences in instrument orientation and position produces mismatches between the projected pixel grids. These mismatches ensure that any point on the ground is sampled differently in each image. If all the images can be registered with respect to each other to a small fraction of a pixel accuracy, then the information from the multiple images can be combined to increase linear resolution by roughly the square root of the number of images. In addition, the gray-scale resolution of the composite image is also improved. We describe methods for multiple image registration and combination, and discuss some of the problems encountered in developing and extending them. We display test results with 8:1 resolution enhancement, and Viking Orbiter imagery with 2:1 and 4:1 enhancements.

Cheeseman, Peter↗
NASA NTRS1994

On the effects of tidal interaction on thin accretion disks: An analytic study

We calculate tidal effects on two-dimensional thin accretion disks in binary systems. We apply a perturbation expansion to obtain an analytic solution of the tidally induced waves. We obtain spiral waves that are stronger at the inner parts of the disks, in addition to a local disturbance which scales like the strength of the local tidal force. Our results agree with recent calculations of the linear response of the disk to tidal interaction.

Dgani, R.↗
NASA NTRS1995

What determines the direction of minimum variance of the magnetic field fluctuations in the solar wind?

The solar wind is not an isotropic medium; two symmetry axis are provided, first the radial direction (because the mean wind is radial) and second the spiral direction of the mean magnetic field, which depends on heliocentric distance. Observations show very different anisotropy directions, depending on the frequency waveband; while the large-scale velocity fluctuations are essentially radial, the smaller scale magnetic field fluctuations are mostly perpendicular to the mean field direction, which is not the expected linear (WkB) result. We attempt to explain how these properties are related, with the help of numerical simulations.

Grappin, R.↗
NASA NTRS2015

Use of North American and European Air Quality Networks to Evaluate Global Chemistry-Climate Modeling of Surface Ozone

We test the current generation of global chemistry-climate models in their ability to simulate observed, present-day surface ozone. Models are evaluated against hourly surface ozone from 4217 stations in North America and Europe that are averaged over 1 degree by 1 degree grid cells, allowing commensurate model-measurement comparison. Models are generally biased high during all hours of the day and in all regions. Most models simulate the shape of regional summertime diurnal and annual cycles well, correctly matching the timing of hourly (approximately 15:00 local time (LT)) and monthly (mid-June) peak surface ozone abundance. The amplitude of these cycles is less successfully matched. The observed summertime diurnal range (25 ppb) is underestimated in all regions by about 7 parts per billion, and the observed seasonal range (approximately 21 parts per billion) is underestimated by about 5 parts per billion except in the most polluted regions, where it is overestimated by about 5 parts per billion. The models generally match the pattern of the observed summertime ozone enhancement, but they overestimate its magnitude in most regions. Most models capture the observed distribution of extreme episode sizes, correctly showing that about 80 percent of individual extreme events occur in large-scale, multi-day episodes of more than 100 grid cells. The models also match the observed linear relationship between episode size and a measure of episode intensity, which shows increases in ozone abundance by up to 6 parts per billion for larger-sized episodes. We conclude that the skill of the models evaluated here provides confidence in their projections of future surface ozone.

diurnal variations↗
DOE OSTI2026

Hydrogen Production System Scaling Using a High-Fidelity Simulation-Optimization Framework

Proton exchange membrane (PEM) electrolyzers are widely used for hydrogen production, yet few validated, high-fidelity tools can reliably guide scale-up. Using measured performance from a 50-hour hardware-in-the-loop pilot test, a physics-based, plant-level model of a 1.25 MW PEM electrolyzer and its balance-of-plant (BoP) subsystems is developed and validated. The model couples electrochemistry and thermal/flow submodels and is calibrated against pilot test data via a genetic algorithm (GA) workflow. Validation yields a mean absolute percentage error (APE) of 0.43% for cell voltage and stack power. Two scale-out strategies are then benchmarked under a common 7-day wind-and-photovoltaic (PV) profile: (i) linear duplication of 1.25 MW blocks and (ii) shared-BoP architectures. Sharing BoP between stacks reduces BoP energy by 27% at 10 MW and 34% at 100 MW (vs. linear duplication) and improves system specific energy consumption (SEC) to 52.9 and 52.6 kWh/kg, respectively (from 54.0 kWh/kg with linear duplication). Partial-load studies (25-100% set-point) show that cumulative hydrogen production remains nearly constant down to 50% load because all cases use the same weekly renewable-energy input. Below 50%, the power cap limits how much energy can be used within 168 h, which reduces hydrogen output. The model further indicates that the practical operating optimum lies between 50% and 85% load, where efficiency gains begin to appear without significant loss in hydrogen output. Moreover, the efficiency gains at lower loads are offset by reduced production. The validated framework supports scenario-based engineering trade-off studies for large configurations (10-100 MW) and for operating policies under variable renewables.

08 HYDROGEN↗
NASA NTRS1992

Wrinkle ridges on Venusian plains: Indicators of shallow crustal stress orientations at local and regional scales

The plains regions of Venus exhibit a complex array of structural features, including deformation belts of various types, wrinkle ridges, grabens, and enigmatic radar-bright linears. Probably the most pervasive of these structures are the wrinkle ridges, which appear to be morphologically identical to their counterparts on the Moon and Mars. Almost all workers agree that wrinkle ridges result from horizontal compressive stresses in the crust; they either are explained as flexural fold structures, or alternatively as scarps or folds related to reverse faults. Wrinkle ridges generally are narrow, have small amplitudes, and commonly are closely spaced as well, characteristics that imply a shallow crustal origin. If wrinkle ridges are due to horizontally directed compressive stresses in the shallow crust, as generally has been inferred, then the trends of these features provide a means to map both local and regional orientations of principal stresses in the uppermost part of the venusian crust: maximum compressive stress is normal to the ridges, minimum compressive stress is normal to the topographic surface, and thus the wrinkle ridge trends trace the orientation of the intermediate principal stress. Because there are few plains areas on Venus totally devoid of wrinkle ridges, it should be possible to establish a number of interesting relationships on a near-global scale by mapping the trends of wrinkle ridges wherever they occur. The present study is addressing three questions: (1) Do the trends of wrinkle ridges define domains that are large relative to the sizes of individual plains regions? If so, can these domains be related to large-scale topographic or geologic features? (2) Are regional trends of wrinkle ridges affected by local features such as coronae? If so, is it possible to determine the relative ages of the far-field and local stresses from detailed study of trend inheritance or superposition relationships? (3) What is the relationship between wrinkle ridges and the larger ridges that make up ridge belts?

Mcgill, George E.↗
DOE OSTI2025

Demonstration of Vertical 2T-nC FeRAM Hybrid Cell and Its Scalability for High-Density 3-D Ferroelectric Capacitor Memory

In this work, we present a comprehensive experimental and modeling study on the scaling of vertical 2T-nC ferroelectric random access memory (FeRAM) hybrid cells, comprising n metal-ferroelectric–metal (MFM) capacitors, to demonstrate a high-performance and high-density 3-D capacitor memory. Our contributions include: 1) successful process integration of vertical 2T-3C FeRAM cells by stacking MFM structures on top of Si CMOS transistors; 2) experimental validation of memory cell functionality, confirming the feasibility of the vertical 2T-nC FeRAM architecture; 3) an analysis of scaling effects on parasitic capacitance in densely integrated 3-D arrays, using 3-D technology computer-aided design (TCAD) simulations; 4) exploration of aggressive stacking of write bitlines (WBLs) to enhance memory density, where ferroelectric linear capacitance ( C FE ) enables self-boosted inhibition under the V W /2 scheme, but renders the V W /3 scheme ineffective due to intolerable write disturbances; and 5) assessment of horizontal scaling, revealing significant increases in read disturbances caused by interplane capacitance between adjacent WBLs ( C Z ). This work represents an early exploration into the potential of 2T-nC FeRAM as a scalable and efficient 3-D memory solution.

42 ENGINEERING↗
DOE OSTI2020

Reducing communication in algebraic multigrid with multi-step node aware communication

Algebraic multigrid (AMG) is often viewed as a scalable [Formula: see text] solver for sparse linear systems. Yet, AMG lacks parallel scalability due to increasingly large costs associated with communication, both in the initial construction of a multigrid hierarchy and in the iterative solve phase. This work introduces a parallel implementation of AMG that reduces the cost of communication, yielding improved parallel scalability. It is common in Message Passing Interface (MPI), particularly in the MPI-everywhere approach, to arrange inter-process communication, so that communication is transported regardless of the location of the send and receive processes. Performance tests show notable differences in the cost of intra- and internode communication, motivating a restructuring of communication. In this case, the communication schedule takes advantage of the less costly intra-node communication, reducing both the number and the size of internode messages. Node-centric communication extends to the range of components in both the setup and solve phase of AMG, yielding an increase in the weak and strong scaling of the entire method.

Computer Science↗
DOE OSTI2021

Numerical solution of large scale Hartree–Fock–Bogoliubov equations

The Hartree–Fock–Bogoliubov (HFB) theory is the starting point for treating superconducting systems. However, the computational cost for solving large scale HFB equations can be much larger than that of the Hartree–Fock equations, particularly when the Hamiltonian matrix is sparse, and the number of electrons N is relatively small compared to the matrix size N b . We first provide a concise and relatively self-contained review of the HFB theory for general finite sized quantum systems, with special focus on the treatment of spin symmetries from a linear algebra perspective. We then demonstrate that the pole expansion and selected inversion (PEXSI) method can be particularly well suited for solving large scale HFB equations. For a Hubbard-type Hamiltonian, the cost of PEXSI is at most $\mathcal{O}$(N b 2 ) for both gapped and gapless systems, which can be significantly faster than the standard cubic scaling diagonalization methods. We show that PEXSI can solve a two-dimensional Hubbard-Hofstadter model with N b up to 2.88 × 10 6 , and the wall clock time is less than 100 s using 17 280 CPU cores. Finally, this enables the simulation of physical systems under experimentally realizable magnetic fields, which cannot be otherwise simulated with smaller systems.

97 MATHEMATICS AND COMPUTING↗