Search NASA⌕ Search

SEARCH · Search NASA

Results for “Projection scheme”

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,063 records · Page 59

Group Projected subspace pursuit for IDENTification of variable coefficient differential equations (GP-IDENT)

We propose an effective and robust algorithm for identifying partial differential equations (PDEs) with space-time varying coefficients from the noisy observation of a single solution trajectory. Identifying unknown differential equations from noisy data is a difficult task, and it is even more challenging with space and time varying coefficients in the PDE. The proposed algorithm, GP-IDENT, has three ingredients: (i) we use B-spline bases to express the unknown space and time varying coefficients, (ii) we propose Group Projected Subspace Pursuit (GPSP) to find a sequence of candidate PDEs with various levels of complexity, and (iii) we propose a new criterion for model selection using the Reduction in Residual (RR) to choose an optimal one among a pool of candidates. The new GPSP considers group projected subspaces which is more robust than existing methods in distinguishing correlated group features. We test GP-IDENT on a variety of PDEs and PDE systems, and compare it with the state-of-the-art parametric PDE identification algorithms under different settings to illustrate its outstanding performance. Furthermore, our experiments show that GP-IDENT is effective in identifying the correct terms from a large dictionary, and our model selection scheme is robust to noise.

Data-driven method↗

Impacts of Suppressing Excessive Light Rain on Aerosol Radiative Effects and Health Risks

Global climate models (GCMs) have been used widely to study radiative forcing and health risks of aerosols. A recent study using two GCMs found that light rain plays a dominant role in controlling aerosol loading. However, “too much light rain and too little heavy rain” is a longstanding bias in GCMs. It is unclear how much light rain affects aerosol-cloud-radiation interactions and health risks from air pollution. Here, in this work, we show that, with the correction of the rainfall intensity spectrum in the National Center for Atmospheric Research Community Atmosphere Model version 5.3 by introducing a stochastic deep convection scheme, the reduced frequency of light rain (1–20 mm d -1 ) results in changes of aerosol direct radiative effects (DRE) of up to -0.5 ± 0.03 W/m 2 and aerosol cloud radiative effects (CRE) of up to -0.9 ± 0.03 W/m 2 . The total (CRE + DRE) radiative effects of light rain-mediated aerosol changes exceed the present-day anthropogenic forcing of aerosols relative to preindustrial levels from the Coupled Model Intercomparison Project (CMIP5&6) models. However, the correction of the rainfall intensity spectrum has little effect on anthropogenic aerosol forcing (defined as the radiative perturbation due to changes in aerosol concentrations between the industrial era and preindustrial levels). Due to increased exposure to fine particulates (PM 2.5 ), the estimated global total premature mortality is much higher than previously estimated, by 300,000 ± 60,000 deaths per year, and is more severe in populous regions such as India and China. The findings in this study highlight the need to understand uncertainties in radiative effects and health risks of aerosols due to simulation biases of precipitation in GCMs.

54 ENVIRONMENTAL SCIENCES↗

Determination of nanoparticle collision cross section distribution functions in low pressure plasma synthesis reactors via ion mobility spectrometry

Low pressure (< 20 Torr) non-thermal plasma flow tube reactors are often used for the synthesis of high purity nanoparticles (NPs) via nucleation and subsequent growth from vapor phase precursors. In spite of their utility, process monitoring (i.e. output NP size characterization) in such reactors is difficult, as there is a dearth of techniques available for online nanoparticle size distribution function measurement. In this study, we developed and applied an ion mobility spectrometry (IMS) system consisting of a low pressure differential mobility analyzer (LPDMA) and electrical detector to determine the collision cross section distribution functions of Si NPs synthesized in a radiofrequency non-thermal SiH4-Ar plasma operated at ~2 Torr (266 Pa). The collision cross section, roughly proportional to projected area, is a parameter quantifying the size and structure of nanometer scale species in the vapor phase. We introduce the collision cross section distribution function as a metric for potential use in online process monitoring in NP synthesis. Proper inversion of the collision cross section distribution function requires a priori knowledge of the LPDMA transfer function. We utilized a tandem differential mobility analyzer approach, coupled with a Twomey-Markwoski based scheme to determine LPDMA transfer functions. Subsequent application of these transfer functions in collision cross section distribution function determination showed that at the outlet of the non-thermal plasma flow tube reactor, both negatively and positively charged Si NPs persisted with nearly identical collision cross section distribution functions. NPs were found to have mode “mobility equivalent” diameters near 10 nm, and via TEM analysis were found to persist at the reactor outlet as small aggregates composed of ~5 nm diameter primary particles. Size distributions functions inferred from collision cross section measurements were compared to size distribution functions inferred from TEM images; excellent agreement was found between IMS and TEM for both mean mobility equivalent diameter and distribution function width. In total, this study shows that (1) IMS is a viable approach for process monitoring in non-thermal plasma NP flow tube synthesis systems, (2) although NPs are modestly aggregated at low-pressure plasma reactor outlets, the extent of aggregation is considerably less than observed in most atmospheric pressure and equilibrium NP synthesis systems (e.g. flames), and (3) after exiting the plasma reactor, the decharging of NPs from highly negatively charge states to a bipolar charge distribution likely drives aggregation on the plasma boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ice-Sheet / Ocean Interaction Model for Greenland Fjords Using High-Order Discontinuous Galerkin Methods (Final Report)

This report summarizes the development and validation of the Non-hydrostatic Unified Model of the Ocean (NUMO). The overarching science goal of the project was to advance Earth System science by proposing a more realistic representation of physical processes in a Greenland fjord. To achieve this objective, we proposed to build a new, separate, high-resolution ocean model to realistically represent the complex three-dimensional physics of the fjord, its bathymetry, coastlines using continuous/discontinuous Galerkin (CG/DG) methods. The report describes the technical details of model components, i.e., choice of equations set, time integration scheme, boundary conditions, spatial discretization methods, mesh generation. We show a suite of validation results for test cases ranging from the Kovasznay flow, through the lock-exchange, simplified two-dimensional fjord, to a three-dimensional density current on a slope. Results are compared with the literature. The report is concluded with a discussion of efforts taken to perform a simulation of Sermilik Fjord with realistic bathymetry data.

58 GEOSCIENCES↗

Compiler-Assisted Detection of Transient Memory Errors

The probability of bit flips in hardware memory systems is projected to increase significantly as memory systems continue to scale in size and complexity. Effective hardware-based error detection and correction requires that the complete data path, involving all parts of the memory system, be protected with sufficient redundancy. First, this may be costly to employ on commodity computing platforms and second, even on high-end systems, protection against multi-bit errors may be lacking. Therefore, augmenting hardware error detection schemes with software techniques is of consider- able interest. In this paper, we consider software-level mechanisms to comprehensively detect transient memory faults. We develop novel compile-time algorithms to instrument application programs with checksum computation codes so as to detect memory errors. Unlike prior approaches that employ checksums on computational and architectural state, our scheme verifies every data access and works by tracking variables as they are produced and consumed. Experimental evaluation demonstrates that the proposed comprehensive error detection solution is viable as a completely software-only scheme. We also demonstrate that with limited hardware support, overheads of error detection can be further reduced.

def-use analysis↗

Global Energy and Water Cycle Experiment (GEWEX) and the Continental-scale International Project (GCIP)

A discussion of the objectives of the Global Energy and Water Cycle Experiment (GEWEX) and the Continental-scale International Project (GCIP) is presented in vugraph form. The objectives of GEWEX are as follows: determine the hydrological cycle by global measurements; model the global hydrological cycle; improve observations and data assimilation; and predict response to environmental change. The objectives of GCIP are as follows: determine the time/space variability of the hydrological cycle over a continental-scale region; develop macro-scale hydrologic models that are coupled to atmospheric models; develop information retrieval schemes; and support regional climate change impact assessment.

Vane, Deborah↗

HSR Model Deformation Measurements from Subsonic to Supersonic Speeds

This paper describes the video model deformation technique (VMD) used at five NASA facilities and the projection moire interferometry (PMI) technique used at two NASA facilities. Comparisons between the two techniques for model deformation measurements are provided. Facilities at NASA-Ames and NASA-Langley where deformation measurements have been made are presented. Examples of HSR model deformation measurements from the Langley Unitary Wind Tunnel, Langley 16-foot Transonic Wind Tunnel, and the Ames 12-foot Pressure Tunnel are presented. A study to improve and develop new targeting schemes at the National Transonic Facility is also described. The consideration of milled targets for future HSR models is recommended when deformation measurements are expected to be required. Finally, future development work for VMD and PMI is addressed.

Burner, A. W.↗

Computational Astrophysics Consortium (2006 - 2013)

The Computational Astrophysics Consortium was a SciDAC Science Application Project that operated during the period 2006 - 2013. The collaboration was headed by Stan Woosley at UC Santa Cruz, and (initially) included partners at UC Berkeley (UCB), LLNL, LANL, LBNL, Stonybrook, Arizona, Stanford, and Johns Hopkins. This report summarizes activities at UCB over the project period. More detailed information about the project accomplishments can be found in the collaboration final report submitted by PI Stan Woosley, and attached to this document. The CAC effort at UCB was originally led by PI Chris Mckee, and included participating faculty members Jon Arons, Peter Nugent, and Eliot Quataert. In 2010, Daniel Kasen moved from UCSC to a faculty position at UCB, and took over as UCB institution PI. The project supported several graduate students (Michael Childress, Siva Dharba, Andrew Myers, Janos Botyanszki, Jeffery Silverman) and postdocs (Dovi Poznanski, Antonino Cucchiarra, Jeffrey Oishi, Bradley Cenko, Eric Hsiao and Jason Dexter). Several codes were further developed by the UCB CAC group. The Orion code is a 3-dimensional radiation hydrodynamics code, used by McKee’s group and often applied to problems in star formation. As part of the CAC, an adaptive ray tracing ionization scheme was implemented, and a chemistry solver was incorporated for coupling the ionization state of the gas to its thermodynamics. The Sedona code is 3-dimensional monte carlo radiation transport code for calculation the light curves and spectra. As part of the CAC, the physics of the code was improved and domain decomposition techniques were implemented for scaling on massively parallel systems.

79 ASTRONOMY AND ASTROPHYSICS↗

An extraction scheme for future CEBAF FFA based energy upgrade

Jefferson lab is considering an energy increase from current 12 GeV to 22 GeV for its CEBAF accelerator. This will be accomplished by recirculating 5-6 additional turns through two parallel CEBAF LINACs using an FFA arc at each end of the racetrack. The total recirculation turns would be 10 times, the first four turns use present conventional arcs to make the 180-degree bends from one LINAC to the other. However, the last 5-6 turns will all share a single beam line inside two FFA arcs. This reduces the footprint and the cost of the project significantly. On the other hand, having the trajectories of last 5-6 recirculating beams close to each other makes it challenging to extract beams from different passes with different energies. In this paper we will explain our present extraction system for 12 GeV, our challenges and limitations, and a possible extraction solution for the 22 GeV upgrade with the goal of extracting beam at different turns/energies to different experimental halls.

Accelerator Physics↗

Coupling Mechanistically the Convective Motions and Cloud Macrophysics in a Climate Model (Final Technical Report for CM4 Project)

Global climate models fail to adequately represent shallow cumulus conditions. The activities at the University of Texas at Austin focus on deploying advanced statistical approaches to identify linkages between numerical models and observations using the E3SM-Single Column Model (SCM) and different ARM campaigns focused on shallow cumulus conditions at SGP. The current version of the SCM uses CLUBB to represent shallow cumulus. The parameterization represents a significant evolution from older generation parameterizations insofar as it contains a number of prognostic equations that keep track of the turbulent nature of the environment that are carried forward in time. The parameters associated with this scheme address the rate at which different prognostic quantities that determine vertical turbulent fluxes of heat and moisture decay with time. The large number of parameters gives the CLUBB scheme a lot of flexibility. The expectation of its developers is that it should represent any shallowly convective environment. The concern among developers and other theorists was not whether it could represent any one environment, but whether a single judicious choice of parameters would allow it to predict the turbulent characteristics of a range of environments. Through several exhaustive studies that made use of Bayesian statistical calibration, tens of thousands of SCM experiments, and ensembles of observations of shallow cumulus at SGP that were also modeled using Large Eddy Simulations (LES), we have only achieved partial success in using CLUBB within SCM to capture important aspects of what has been observed and modeled with LES at SGP. Even as we worked to improve observational tests and make a more scientifically focused attempt at modifying subsets of CLUBB parameters, model performance continues to be lacking in some important respects. In one calibration, we were able to achieve a near perfect match to observed properties of liquid water path, cloud radiative effects, and cloud fraction but with clouds that were trapped within a shallow layer and without the highly turbulent fluxes and vertical extent observed with ARM data and modeled with LES (Figure 1). The conclusion we are reaching is that shallow cumulus conditions are much more challenging to represent than was previously thought. This result is also prompting us to consider some aspects of the experiment design that may be limiting the performance of SCM in this case. The work has garnered the attention of Vincent Larson (CLUBB developer at Univ. Wisconsin – Milwaukee) and model developers Steven Klein and Yunyan Zhang (LLNL) with whom this work and will continue.

58 GEOSCIENCES↗

Extraction scheme for future CEBAF FFA based energy upgrade

Jefferson lab is considering an energy increase from current 12 GeV to 22 GeV for its CEBAF accelerator. This will be accomplished by recirculating 5-6 additional turns through two parallel CEBAF LINACs using an FFA arc at each end of the racetrack. The total recirculation turns would be 10 times, the first four turns use present conventional arcs to make the 180-degree bends from one LINAC to the other. However, the last 5-6 turns will all share a single beam line inside two FFA arcs. This reduces the footprint and the cost of the project significantly. On the other hand, having the trajectories of last 5-6 recirculating beams close to each other makes it challenging to extract beams from different passes with different energies. In this paper we will explain our present extraction system for 12 GeV, our challenges and limitations, and a possible extraction solution for the 22 GeV upgrade with the goal of extracting beam at different turns/energies to different experimental halls.

Kazimi, Reza↗

Scalable FBP decomposition for cone-beam CT reconstruction

Filtered Back-Projection (FBP) is a fundamental compute intense algorithm used in tomographic image reconstruction. Cone-Beam Computed Tomography (CBCT) devices use a cone-shaped X-ray beam, in comparison to the parallel beam used in older CT generations. Distributed image reconstruction of cone-beam datasets typically relies on dividing batches of images into different nodes. This simple input decomposition, however, introduces limits on input/output sizes and scalability.We propose a novel decomposition scheme and reconstruction algorithm for distributed FPB. This scheme enables arbitrarily large input/output sizes, eliminates the redundancy arising in the end-to-end pipeline and improves the scalability by replacing two communication collectives with only one segmented reduction. Finally, we implement the proposed decomposition scheme in a framework that is useful for all current-generation CT devices (7th gen). In our experiments using up to 1024 GPUs, our framework can construct 40963 volumes, for real-world datasets, in under 16 seconds (including I/O).

Chen, Peng↗

FACE-MDS Phase 2: Model Output

Model output from Phase 2 of the Free Air CO2 Enrichment Model Data Synthesis (FACE-MDS) project. Results from 12 models simulating six elevated CO2 experiments using the protocols, met data, and other data supplied in the accompanying dataset: Phase 2 Free Air CO2 Enrichment Model-Data Synthesis: Meteorological Data and Protocols. File naming convention can be found in the file 'facemds-file-table-modeldata.xml' and information about variable naming convention can be found in the file 'facemds-standard-name-table-modvars.xml'. We followed the ALMA v3 convention for meteorological variables and developed a hierarchical naming scheme for ecosystem variables. Additional info can be found in the protocols in the accompanying meteorological dataset (doi:10.15485/1480328). We replaced the original naming convention mentioned in the protocol with the one described in 'facemds-standard-name-table-modvars.xml'.

54 ENVIRONMENTAL SCIENCES↗

Unsupervised Learning Based Interaction Force Model for Nonspherical Particles in Incompressible Flows

This project provides a neural network-based interaction force model for gas-solid flows from low to intermediate Reynolds numbers and concentration, which can be linked to MFiX-DEM. We have constructed a database of the interaction force between the irregular-shaped particles using a spherical harmonic method and the fluid phase based on the particle-resolved direct numerical simulation (PR-DNS) with immersed boundary-based gas kinetic scheme. Unsupervised learning method, i.e., variational auto-encoder (VAE) has been applied to extract the primitive shape factors determining the drag force, lifting forces, and torque. The interaction force model has been trained and validated with a simple but effective multi-layer feed-forward neural network: multi-layer perceptron (MLP), which will be concatenated after the encoder of the previously trained VAE for geometry feature extraction for single, irregular particles. We have trained transpose convolutional neural networks with the PR-DNS data to predict the velocity and pressure gradient of the single particle systems and utilized them to calculate drag force of multi-particle systems. This model can provide high computational efficiency because it does not require collecting multiparticle system data from PR-DNS.

99 GENERAL AND MISCELLANEOUS↗

Development of a standard connector for orbital replacement units for serviceable spacecraft

The current trend for spacecraft to be serviceable and repairable in orbit has led to a modular approach to satellite subsystem design. Spacecraft equipment, such as sensors, tape recorders, computers, transponders, batteries, etc., housed in discrete modular units, (Orbital Replacement Units(ORUs)) can be attached and detached to the spacecraft as needed. The interface between the ORU and the spacecraft is crucial. The projected use of robotics and the need for a common mechanism capable of performing several functions puts many constraints on the design of the interface. Such an interface mechanisms was designed and developed called the Standard Interface Connector (SIC) that mates the ORU to the spacecraft. The SIC also provides for the flow of fluids, data, and power between the module and spacecraft. The baseline design presented can be configured to provide various attachment schemes. Tests on SIC models have demonstrated the functionality of the design and its compatibility with current robotics.

Heath, Ellen F.↗

Theoretical Analysis of Photoelectron Spectra of Pure and Mixed Metal Clusters: Disentangling Size, Structure, and Composition Effects

A theoretical/computational description and analysis of the spectra of electron binding energies of Al 12 - , Al 13 - and Al 12 Ni- clusters, which differ in size and/or composition by a single atom yet possess strikingly different measured photoelectron spectra, is presented. It is shown that the measured spectra can not only be reproduced computationally with quantitative fidelity – this is achieved through a combination of state-of-the-art density functional theory with a highly accurate scheme for conversion of the Kohn-Sham eigenenergies into electron binding energies – but also explained in terms of the effects of size, structure/symmetry and composition. Furthermore, a new methodology is developed and applied that provides for disentanglement and differential assignment of the separate roles played by size, structure/symmetry and composition in defining the observed differences in the measured spectra. The methodology is general and applicable to any finite system, homogeneous or heterogeneous. Finally, we project that in combination with advances in synthesis techniques this methodology will become an indispensable computation-based aid in the design of controlled synthesis protocols for manufacture of nanosystems and nanodevices with precisely desired electronic and other characteristics.

36 MATERIALS SCIENCE↗

High-Order Preconditioners for Modern Computing Architectures: Line-Based Discontinuous Galerkin Method on GPUs

The goals of our project are to be able to overcome issues that plague computing in modern architecture such as limited memory bandwidth, badly conditioned Jacobians, and restrictive time steps. We begin by describing the line-based discontinuous Galerkin (line DG) method for a scalar two-dimensional problem. Line DG has been shown to provide a significant increase in sparsity for Jacobian matrix, especially in the diagonal block elements. In addition to this, the implicit time integration scheme helps alleviate the time step restriction but requires solving a (potentially badly conditioned) linear system. We can precondition the linear system to facilitate the solution. Since we want to solve the system on the GPU, the increase in sparsity from the line DG method is helpful to alleviate the memory issue. We successfully implemented a computation of the Jacobian on the CPU and are currently working on the GPU implementation. Finally, I discuss how this opportunity to intern at Lawrence Livermore National Laboratory through the National Science Foundation Mathematical Sciences Graduate Internship has affected my personal and professional growth.

97 MATHEMATICS AND COMPUTING↗

Inclusion of a cold hardening scheme to represent frost tolerance is essential to model realistic plant hydraulics in the Arctic–boreal zone in CLM5.0-FATES-Hydro

As temperatures decrease in autumn, vegetation of temperate and boreal ecosystems increases its tolerance to freezing. This process, known as hardening, results in a set of physiological changes at the molecular level that initiate modifications of cell membrane composition and the synthesis of anti-freeze proteins. Together with the freezing of extracellular water, anti-freeze proteins reduce plant water potentials and xylem conductivity. To represent the responses of vegetation to climate change, land surface schemes increasingly employ “hydrodynamic” models that represent the explicit fluxes of water from soil and through plants. The functioning of such schemes under frozen soil conditions, however, is poorly understood. Nonetheless, hydraulic processes are of major importance in the dynamics of these systems, which can suffer from, e.g., winter “frost drought” events. In this study, we implement a scheme that represents hardening into CLM5.0-FATES-Hydro. FATES-Hydro is a plant hydrodynamics module in FATES, a cohort model of vegetation physiology, growth, and dynamics hosted in CLM5.0. We find that, in frozen systems, it is necessary to introduce reductions in plant water loss associated with hardening to prevent winter desiccation. This work makes it possible to use CLM5.0-FATES-Hydro to model realistic impacts from frost droughts on vegetation growth and photosynthesis, leading to more reliable projections of how northern ecosystems respond to climate change.

58 GEOSCIENCES↗