Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameterization”

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,207 records · Page 67

SU-G-TeP4-03: A Multileaf Collimator Calibration and Quality Assurance Technique Using An Electronic Portal Imaging Device

Purpose: To develop an accurate and quick multileaf collimator (MLC) calibration and quality assurance technique using an electronic portal imaging device (EPID) Methods: The MLC models used include the MLCi and Agility (Elekta Ltd). This technique consists of two 22(L)x10(W) cm{sup 2} fields with 0{sup 0} and 180{sup 0} collimator angles centered to an offset EPID. The MLC opening is estimated by calculating the profile at the image’s center in the image’s horizontal direction. Scans in the image’s vertical direction were calculated every 20 pixels in the inner 70% of estimated MLC opening. The profiles’ edges were fitted with linear equations to determine the image’s rotation angle. Then, crossline profiles were scanned at the center of each leaf taking into account the leaf’s width at isocenter and the rotation angle. The profiles’ edges determine the location of the leaves’ edges and these were subtracted from the reference leaf’s position in order to determine the relative leaf offsets. The edge location of all profiles was determined by using the parameterized gradient of the penumbra region. The technique was tested against an established diode array-based method, and for different MLC systems, patterns, gantry angles, days, energies, beam modalities and MLC openings. Results: The differences between the proposed and established methods were 0.26±0.19mm. The leaf offsets’ deviation was <0.3mm (5 months period). For pattern fields, the differences between predetermined and calculated offsets were 0.18±0.18mm. The leaf offset deviation of measurements with different energies and MLC openings were <0.1mm and <0.3mm, respectively. The differences between offsets of FF and FFF beams were 0.01±0.02mm (<0.07mm). The differences between the offsets at different gantry angles were 0.08±0.15mm. Conclusion: The proposed method proved to be accurate and efficient in calculating the relative leaf offsets. Parameterized field edge is essential to obtain accurate result by eliminating the noise from EPID.

60 APPLIED LIFE SCIENCES↗

WE-H-BRA-02: Radiobiological Modeling of Tumor Control Probability (TCP) and Radiation-Induced Pneumonitis (RP) for Lung Cancer Patients Treated with Monte Carlo-Based Lung SBRT

Purpose: To present radiobiological modeling of TCP using tumor size-adjusted BED(s-BED)and PTV(D99) to lung SBRT patients treated with X-ray Voxel Monte Carlo(XVMC) algorithm, apply parameterized Lyman-NTCP model to predict grade-2 RP and subsequently, compare with clinical outcomes/observations. Methods: Dosimetric parameters and clinical follow-up for XVMC-based lung-SBRT patients were retrospectively evaluated. Patients were treated at Novalis-TX with hybrid(2 non-coplanar partial-arcs plus 3–6 static-beams)plan using HD-MLC/6MV-SRS-beam.For TCP,s-BED modelling was utilized: TCP=EXP[sBED-TCD50]/k/(1.0+EXP[sBED-TCD50]/k), where k=31Gy corresponding to TCD50=0Gy and s-BED was defined as BED10 minus 10 times the tumor diameter(in centimeters)by Ohri et al.(IJROBP,2012). For 2-yr local-control, we used more-realistic MC-computed PTVD99 as a predictive parameter, s-BED(D99).Due to relatively shorter median follow-up interval(12-months),Kaplan-Meier curves were generated to estimate 2-yr observed local-control and compared to predicted-rate by TCP modeling. For NTCP, we employed parameterized Lyman-NTCP model utilizing normal-lung DVH and α/β=3Gy fitted to predict grade-2 RP after lung-SBRT. Results: Total 108 patients (137 tumors) treated for 35–70Gy in 3–5 fractions, either primary-lung(n=74)or metastatic-lung(n=53)tumors were included.F or the given prescription dose with MC-computed MUs, 2-yr local-control rates with s-BED(D99) was 87±8%. Kaplan-Meier generated observed local-control rate at 2-yr was 87.5%,suggesting that PTV(D99) could be a potential predictor (p-value=0.38). Observed vs predicted TCP for primary-lung tumors and metastatic tumors were 97% vs 88±7% and 94% vs 86±9%.NTCP model predicted well for symptomatic-RP with predicted vs observed (3±5% vs 2%). Radiographic and clinically significant RP was observed in 13% and 2% of patients. Higher rates of radiographic change were observed in patients who received >50Gy compared to ≤50Gy(24% vs 10%). Conclusion: Utilizing MC-computed PTVD99, our TCP results were well correlated with clinical outcome. The predicted grade-2 RP rate was comparable to clinical observations. Clinical application of these radiobiological models may potentially allow for target dose escalation and/or lung-toxicity reduction. Further validation of these radiobiological models with longer follow up interval for large cohorts of lung-SBRT patients is anticipated.

60 APPLIED LIFE SCIENCES↗

Bioaerosols are the dominant source of warm-temperature immersion-mode INPs and drive uncertainties in INP predictability

Ice-nucleating particles (INPs) are rare atmospheric aerosols that initiate primary ice formation, but accurately simulating their concentrations and variability in large-scale climate models remains a challenge. Doing so requires both simulating major particle sources and parameterizing their ice nucleation (IN) efficiency. Validating and improving model predictions of INP concentrations requires measuring their concentrations delineated by particle type. We present a method to speciate INP concentrations into contributions from dust, sea spray aerosol (SSA), and bioaerosol. Field campaign data from Bodega Bay, California, showed that bioaerosols were the primary source of INPs between –12° and –20°C, while dust was a minor source and SSA had little impact. We found that recent parameterizations for dust and SSA accurately predicted ambient INP concentrations. However, the model did not skillfully simulate bioaerosol INPs, suggesting a need for further research to identify major factors controlling their emissions and INP efficiency for improved representation in models.

54 ENVIRONMENTAL SCIENCES↗

The nonsinglet structure function evolution by Laplace method

We derive a general scheme for the evolution of the nonsinglet structure function at the leadingorder (LO) and next-to-leading-order (NLO) by using the Laplace-transform technique. Results for the nonsinglet structure function are compared with MSTW2008, GRV, and CKMT parameterizations and also EMC experimental data in the LO and NLO analysis. The results are in good agreement with the experimental data and other parameterizations in the low- and large-x regions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Neutron star radii, universal relations, and the role of prior distributions

We investigate constraints on neutron star structure arising from the assumptions that neutron stars have crusts, that recent calculations of pure neutron matter limit the equation of state of neutron star matter near the nuclear saturation density, that the high-density equation of state is limited by causality and the largest high-accuracy neutron star mass measurement, and that general relativity is the correct theory of gravity. We explore the role of prior assumptions by considering two classes of equation of state models. In a first, the intermediate- and high-density behavior of the equation of state is parameterized by piecewise polytropes. In the second class, the high-density behavior of the equation of state is parameterized by piecewise continuous line segments. The smallest density at which high-density matter appears is varied in order to allow for strong phase transitions above the nuclear saturation density. We critically examine correlations among the pressure of matter, radii, maximum masses, the binding energy, the moment of inertia, and the tidal deformability, paying special attention to the sensitivity of these correlations to prior assumptions about the equation of state. It is possible to constrain the radii of 1.4 solar mass neutron stars to be larger than 10 km, even without consideration of additional astrophysical observations, for example, those from photospheric radius expansion bursts or quiescent low-mass X-ray binaries. We are able to improve the accuracy of known correlations between the moment of inertia and compactness as well as the binding energy and compactness. Furthermore, we also demonstrate the existence of a correlation between the neutron star binding energy and the moment of inertia.

79 ASTRONOMY AND ASTROPHYSICS↗

TADSim: Discrete Event-based Performance Prediction for Temperature Accelerated Dynamics

Next-generation high-performance computing will require more scalable and flexible performance prediction tools to evaluate software--hardware co-design choices relevant to scientific applications and hardware architectures. Here, we present a new class of tools called application simulators—parameterized fast-running proxies of large-scale scientific applications using parallel discrete event simulation. Parameterized choices for the algorithmic method and hardware options provide a rich space for design exploration and allow us to quickly find well-performing software--hardware combinations. We demonstrate our approach with a TADSim simulator that models the temperature-accelerated dynamics (TAD) method, an algorithmically complex and parameter-rich member of the accelerated molecular dynamics (AMD) family of molecular dynamics methods. The essence of the TAD application is captured without the computational expense and resource usage of the full code. We accomplish this by identifying the time-intensive elements, quantifying algorithm steps in terms of those elements, abstracting them out, and replacing them by the passage of time. We use TADSim to quickly characterize the runtime performance and algorithmic behavior for the otherwise long-running simulation code. We extend TADSim to model algorithm extensions, such as speculative spawning of the compute-bound stages, and predict performance improvements without having to implement such a method. Validation against the actual TAD code shows close agreement for the evolution of an example physical system, a silver surface. Finally, focused parameter scans have allowed us to study algorithm parameter choices over far more scenarios than would be possible with the actual simulation. This has led to interesting performance-related insights and suggested extensions.

97 MATHEMATICS AND COMPUTING↗

Differentiable Quantum Programming with Unbounded Loops

The emergence of variational quantum applications has led to the development of automatic differentiation techniques in quantum computing. Existing work has formulated differentiable quantum programming with bounded loops, providing a framework for scalable gradient calculation by quantum means for training quantum variational applications. However, promising parameterized quantum applications, e.g., quantum walk and unitary implementation, cannot be trained in the existing framework due to the natural involvement of unbounded loops. To fill in the gap, we provide the first differentiable quantum programming framework with unbounded loops, including a newly designed differentiation rule, code transformation, and their correctness proof. Technically, we introduce a randomized estimator for derivatives to deal with the infinite sum in the differentiation of unbounded loops, whose applicability in classical and probabilistic programming is also discussed. We implement our framework with Python and Q# and demonstrate a reasonable sample efficiency. Through extensive case studies, we showcase an exciting application of our framework in automatically identifying close-to-optimal parameters for several parameterized quantum applications.

Computer Science↗

Bond Synergy Model for Bond Energies in Alloy Oxides

In this work we introduce a metal-oxide bond-energy model for alloy oxides based on pure-phase bond energies and bond synergy factors that describe the effect of alloying on the bond energy between cations and oxygen, an important quantity to understand the formation of alloy oxides and their composition. This model is parameterized for binary cation-alloy oxides using density-functional theory energies and is shown to be directly transferable to multi-component alloy oxides. We parameterized the model for alloy oxide energies with metal cations that form the basis of corrosion resistant alloys, including Fe, Ni, Cr, Mo, Mn, W, Co, and Ru. We find that isoelectronic solutes allow quantification of pure-phase bond energies in oxides and the calculated bond energy values give sensible results compared to common experience, including the role of Cr as the passive-layer former in Fe–Ni–Cr alloys for corrosion applications. Additionally, the bond synergy factors give insights into the mutual strengthening and weakening effects of alloying on cation-oxygen bonds and can be related to enthalpy of mixing and charge neutrality constraints. We demonstrate how charge neutrality can be identified and achieved by the oxidation states that the different cations assume depending on alloy composition and the presence of defects.

Electrochemistry↗

Constructing Hardware in a Scale Embedded Language

Chisel is a new open-source hardware construction language developed at UC Berkeley that supports advanced hardware design using highly parameterized generators and layered domain-specific hardware languages. Chisel is embedded in the Scala programming language, which raises the level of hardware design abstraction by providing concepts including object orientation, functional programming, parameterized types, and type inference. From the same source, Chisel can generate a high-speed C++-based cycle-accurate software simulator, or low-level Verilog designed to pass on to standard ASIC or FPGA tools for synthesis and place and route.

Bachan, John↗

ExaSat: An Exascale Static Analysis Tool (ExaSat) v1.0

The ExaSAT framework automates the extraction of parameterized performance models directly from source code using compiler analysis. Its parameterized analytic model enables quantitative evaluation of a broad range of hardware design trade-offs and software optimizations on a variety of different performance metrics, with a primary focus on data movement as a metric. ExaSAT analysis provides insights in the hardware and software tradeoffs and lays the groundwork for exploring a more targeted set of design points using cycle-accurate architectural simulators and emulators.

Chan, Cy↗

BrzostekEcologyLab/CORPSE-soil-jars

Efforts to manage soils for carbon (C) sequestration remain limited by our understanding of how differences in plant traits and microbial traits mechanistically drive soil organic C (SOC) storage. Addressing this uncertainty is particularly critical in bioenergy agriculture, due to its potential to enhance soil C and provide a C neutral fuel. As such, we examined differences between two contrasting feedstocks, Zea mays (corn) and Miscanthus x giganteus (miscanthus), in the ability of their litter to form new chemically resistant particulate SOC vs. physically protected mineral associated SOC and used this data to improve the parameterization of a microbial SOC model. We tested a hypothesized conceptual model whereby easy to decompose corn litters drive greater microbial carbon use efficiency (CUE) and the formation of more mineral associated SOC over particulate SOC than more complex miscanthus litters. To do this, we performed a soil microcosm experiment where we added 13C enriched aboveground and belowground litters to soils and traced the fate of the 13C into microbial respiration and SOC pools. We found that corn litters promoted higher microbial CUE (0.37) than miscanthus litters (0.24). In turn, corn litter formed approximately 50% more mineral associated SOC than miscanthus litters. Similarly, structurally complex root litters promoted a lower CUE and formed less mineral associated SOC than leaf and shoot litters for both crops. When we used our data to parameterize the SOC model, we found that modeling microbial trait differences uniquely allowed the model to capture the fate of litter C in SOC. Collectively, we found a robust link between litter quality, microbial efficiency, and the formation of SOC. This link bridges the empirical uncertainty in how different crops can form new soil C and provides an empirical basis for modeling SOC transformations.

Ridgeway, Joanna↗

Energy Exascale Earth System Model v2.0

First release of version 2 of the Energy Exascale Earth System Model. The atmosphere component remains EAM. Major changes since version 1 include: all column-physics parameterizations are computed on a separate grid that has approximately half the number of points of the dynamics grid, a new nonhydrostatic dynamical core (running in hydrostatic mode) with semi-Lagrangian tracer transport, CLUBB updated from v1 to v2, a new convective trigger (dCAPE/ULL) based on the dynamic Convective Available Potential Energy (CAPE) (dCAPE) and the Unrestricted Launch Level (ULL) concepts is used in ZM. minimum cloud droplet number changed, gravity wave drag energy conservation fixed and new tunings used, dust emission size distribution changed to emit more coarse dust particles The land component is still ELM. Major changes since version 1 include: using SNICAR-AD for radiation in snow to match the sea-ice model and fixing bugs in snow compaction and water state calculation. The ocean component remains MPAS-ocean. Major change since version 1 include: Redi isopycnal mixing has been updated, tested, and tuned in combination with the Gent-McWilliams parameterization, a sign error was fixed in the 3rd-order flux routines, the mesh used in low-resolution coupled cases was modified and the time steps adjusted, new regionally refined meshes were created, one focused on North America and another on the Southern Ocean, including ice shelf cavities. The sea-ice component remains MPAS-seaice. Major changes since version 1 include: A new heat- and freshwater-conserving coupling of frazil ice, turning off of SSH filtering, addition of SNICAR-AD and snow grain aging. The land-ice component remains MPAS-Albany-landIce (MALI) and is a static ice sheet. There is more out-of-the-box support for cryosphere configurations including the new regionally refined configuration around Antarctica with ice shelf cavities. The river model is MOSART. The half-degree river mesh was redone so it no longer treats Black and Caspian seas as ocean. The coupler remains cpl7/MCT. Major changes since version 1 include handling of ice shelf melt fluxes (heat / freshwater exchange with the ocean), and data icebergs. All components allow regional refinement of their meshes and two separate example of refinement, one in and around North America and one around Antarctica and the Southern Ocean, are provided.

E3SM Project, DOE↗

Energy Exascale Earth System Model v2.1.0

First release of version 2.1 of the Energy Exascale Earth System Model. [ATM] The atmosphere component remains EAM. There are no major changes in the default configuration since 2.0. New features include: A semi-lagrangian tracer transport for theta-l dycore, a new algorithm for finding the tropopause, new RRM mesh configurations. Add and update SSP370 and SSP585 cases. Restore the FIDEAL case. [LAND] The land component is ELM. There are no major changes in the default configuration since 2.0. Several option features have been added including: implementation of topography-based subgrid structure (topounits) and accompanying parameterizations and atmospheric forcing downscaling methods; a new plant hydraulics scheme; two-way land-river hydrological coupling through the infiltration of floodplain water; an implementation of perennial crops; updates to the SNICAR-AD snow radiative transfer model; and implementation of soil erosion and sediment yield in ELM-Erosion. Each of these new changes is modular in design and can be turned on or off as the user specifies; they are currently being tested in different “BGC” configurations. [OCEAN] The ocean component remains MPAS-Ocean. Major change since version 2.0 include the addition of the Fox-Kemper et al. 2011 parameterization for submesocale eddies, a correction for barotropic thickness consistency that reduces divergence noise, and the addition of an ocean carbon conservation analysis member. [SEAICE] The sea-ice component remains MPAS-Seaice. Major changes since version 2.0 include: A correction to how shortwave parameters are interpolated in the snicar-ad 5-band radiation scheme, the addition of a sea ice carbon conservation analysis member, updates to the default sea ice biogeochemistry namelist parameters to be consistent with version 2.0 improvements to nitrogen cycling and a correction in the ice-ocean dissolved organic nitrogen coupling. [LAND ICE] The land-ice component remains MPAS-Albany-landIce (MALI). Major changes since 2.0 include an update to the MALI version and the Greenland mesh.[RIVER] The river model is MOSART. There are no major changes in the default configuration since 2.0. A major new optional feature is two-way river-ocean hydrological coupling between MOSART and MPAS-O. This change can be turned on or off as the user specifies, and is being tested in different configurations. [COUPLER] The coupler remains cpl7/MCT. Major changes since version 2.0 include: Carbon budget calculated when heat/water budgets active. Fix a bug in land-atm fluxes for tri-grid configurations. [OTHER] a small bug in the zenith angle calculation was fixed in the data models.

ECP↗

Learning Only on Boundaries: A Physics-Informed Neural Operator for Solving Parametric Partial Differential Equations in Complex Geometries

Recently, deep learning surrogates and neural operators have shown promise in solving partial differential equations (PDEs). However, they often require a large amount of training data and are limited to bounded domains. In this work, we present a novel physics-informed neural operator method to solve parameterized boundary value problems without labeled data. By reformulating the PDEs into boundary integral equations (BIEs), we can train the operator network solely on the boundary of the domain. This approach reduces the number of required sample points from $O(N^d)$ to $O(N^{d-1}$), where $d$ is the domain’s dimension, leading to a significant acceleration of the training process. Additionally, our method can handle unbounded problems, which are unattainable for existing physics-informed neural networks (PINNs) and neural operators. Finally, our numerical experiments show the effectiveness of parameterized complex geometries and unbounded problems.

97 MATHEMATICS AND COMPUTING↗

COSP: Satellite simulation software for model assessment

Errors in the simulation of clouds in general circulation models (GCMs) remain a long-standing issue in climate projections, as discussed in the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report. This highlights the need for developing new analysis techniques to improve our knowledge of the physical processes at the root of these errors. The Cloud Feedback Model Intercomparison Project (CFMIP) pursues this objective, and under that framework the CFMIP Observation Simulator Package (COSP) has been developed. COSP is a flexible software tool that enables the simulation of several satellite-borne active and passive sensor observations from model variables. The flexibility of COSP and a common interface for all sensors facilitates its use in any type of numerical model, from high-resolution cloud-resolving models to the coarser-resolution GCMs assessed by the IPCC, and the scales in between used in weather forecast and regional models. The diversity of model parameterization techniques makes the comparison between model and observations difficult, as some parameterized variables (e.g., cloud fraction) do not have the same meaning in all models. The approach followed in COSP permits models to be evaluated against observations and compared against each other in a more consistent manner. This thus permits a more detailed diagnosis of the physical processes that govern the behavior of clouds and precipitation in numerical models. The World Climate Research Programme (WCRP) Working Group on Coupled Modelling has recommended the use of COSP in a subset of climate experiments that will be assessed by the next IPCC report. Here we describe COSP, present some results from its application to numerical models, and discuss future work that will expand its capabilities.

54 ENVIRONMENTAL SCIENCES↗

Thermodynamics of the Madden-Julian Oscillation in a Regional Model with Constrained Moisture

In order to identify the main thermodynamic processes that sustain the Madden Julian Oscillation, an eddy available potential energy budget analysis is performed on a WRF simulation with moisture constrained by observations. The model realistically simulates the two MJO episodes observed during the winter of 2007-2008. The analysis shows that instabilities and damping associated with variations in diabatic heating and energy transport work in concert to provide the MJO with its observed characteristics. The results are used to construct a simplified paradigm of MJO thermodynamics. Furthermore, the effect of moisture nudging on the simulation is analyzed to understand the limitations of the model cumulus parameterization. Without moisture nudging, the parameterization fails to provide adequate low-level (upper-level) moistening during the early (late) stage of the MJO active phase. The moistening plays a critical role in providing stratiform heating variability that is an important source of eddy available potential energy for the model MJO.

54 ENVIRONMENTAL SCIENCES↗

Challenges and Prospects for Reducing Coupled Climate Model SST Biases in the Eastern Tropical Atlantic and Pacific Oceans: The U.S. CLIVAR Eastern Tropical Oceans Synthesis Working Group

Well-known problems trouble coupled general circulation models of the eastern Atlantic and Pacific Ocean basins. Model climates are significantly more symmetric about the equator than is observed. Model sea surface temperatures are biased warm south and southeast of the equator, and the atmosphere is too rainy within a band south of the equator. Near-coastal eastern equatorial SSTs are too warm, producing a zonal SST gradient in the Atlantic opposite in sign to that observed. The U.S. Climate Variability and Predictability Program (CLIVAR) Eastern Tropical Ocean Synthesis Working Group (WG) has pursued here an updated assessment of coupled model SST biases, focusing on the surface energy balance components, on regional error sources from clouds, deep convection, winds, and ocean eddies; on the sensitivity to model resolution; and on remote impacts. Motivated by the assessment, the WG makes the following recommendations: 1) encourage identification of the specific parameterizations contributing to the biases in individual models, as these can be model dependent; 2) restrict multimodel intercomparisons to specific processes; 3) encourage development of high-resolution coupled models with a concurrent emphasis on parameterization development of finer-scale ocean and atmosphere features, including low clouds; 4) encourage further availability of all surface flux components from buoys, for longer continuous time periods, in persistently cloudy regions; and 5) focus on the eastern basin coastal oceanic upwelling regions, where further opportunities for observational–modeling synergism exist.

58 GEOSCIENCES↗

The Second Wind Forecast Improvement Project (WFIP 2): General Overview

In 2015 the U.S. Department of Energy initiated a four-year study, the second Wind Forecast Improvement Project (WFIP 2), to improve the representation of boundary-layer physics and related processes in mesoscale models for better wind and wind power forecasts. This goal challenges numerical weather prediction models in complex terrain due to inherent assumptions underlying their boundary-layer parameterizations. The WFIP 2 effort involved the wind industry, universities, NOAA, and DOE's national laboratories in an integrated observational and modeling study. Observations spanned 18-months to assure a full annual cycle of continuously recorded observations from remote-sensing and in situ measurement systems. The study area comprised the Columbia Basin of eastern Washington and Oregon, containing more than 4 GW of installed wind capacity; nests of observational systems captured important atmospheric scales from mesoscale to NWP sub-grid scale. Model improvements targeted NOAA's High-Resolution Rapid Refresh (HRRR) model to facilitate transfer of improvements to NWS operational forecast models. This paper describes the general WFIP 2 scope and objectives, the particular scientific challenges of improving wind forecasts in complex terrain, early successes of the project, and an integrated approach to archiving observations and model output. It provides an introduction for a set of more detailed BAMS papers addressing WFIP 2 observational science, modeling challenges and solutions, practical transfer of boundary layer advances to industry use, and advances in coupling improved mesoscale models to microscale models that can represent interactions between wind plants and the atmosphere.

WFIP2, boundary layer, complex terrain, forecastin↗