Search NASA⌕ Search

SEARCH · Search NASA

Results for “CASE METHOD”

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 505 records · Page 28

Ab-Initio Description and Prediction of Properties of Carbon-Based and Other Non-Metallic Materials

We have resolved the long-standing problem consisting of 30%-50% theoretical underestimates of the band gaps of non-metallic materials. We describe the Bagayoko, Zhao, and Williams (BZW) method that rigorously circumvents the basis-set and variational effect presumed to be a cause of these underestimates. We present ab-initio, computational results that are in agreement with experiment for diamond (C), silicon (Si), silicon carbides (3C-SiC and 4H-SiC), and other semiconductors (GaN, BaTiO3, AlN, ZnSe, ZnO). We illustrate the predictive capability of the BZW method in the case of the newly discovered cubic phase of silicon nitride (c-Si3N4) and of selected carbon nanotabes [(10,0), and (8,4)]. Our conclusion underscores the inescapable need for the BZW method in ab-initio calculations that employ a basis set in a variational approach. Current nanoscale trends amplify this need. We estimate that the potential impact of applications of the BZW method in advancing our understanding of nonmetallic materials, in informing experiment, and particularly in guiding device design and fabrication is simply priceless.

Bagayoko, D.↗

Land surface model parameter optimisation using in situ flux data: comparison of gradient-based versus random search algorithms (a case study using ORCHIDEE v1.9.5.2)

Land surface models (LSMs), which form the land component of earth system models, rely on numerous processes for describing carbon, water and energy budgets, often associated with highly uncertain parameters. Data assimilation (DA) is a useful approach for optimising the most critical parameters in order to improve model accuracy and refine future climate predictions. In this study, we compare two different DA methods for optimising the parameters of seven plant functional types (PFTs) of the ORCHIDEE LSM using daily averaged eddy-covariance observations of net ecosystem exchange and latent heat flux at 78 sites across the globe. We perform a technical investigation of two classes of minimisation methods – local gradient-based (the L-BFGS-B algorithm, limited memory Broyden–Fletcher–Goldfarb–Shanno algorithm with bound constraints) and global random search (the genetic algorithm) – by evaluating their relative performance in terms of the model–data fit and the difference in retrieved parameter values. We examine the performance of each method for two cases: when optimising parameters at each site independently (“single-site” approach) and when simultaneously optimising the model at all sites for a given PFT using a common set of parameters (“multi-site” approach). We find that for the single site case the random search algorithm results in lower values of the cost function (i.e. lower model–data root mean square differences) than the gradient-based method; the difference between the two methods is smaller for the multi-site optimisation due to a smoothing of the cost function shape with a greater number of observations. The spread of the cost function, when performing the same tests with 16 random first-guess parameters, is much larger with the gradient-based method, due to the higher likelihood of being trapped in local minima. When using pseudo-observation tests, the genetic algorithm results in a closer approximation of the true posterior parameter value in the L-BFGS-B algorithm. We demonstrate the advantages and challenges of different DA techniques and provide some advice on using it for the LSM parameter optimisation.

58 GEOSCIENCES↗

Addressing delayed case reporting in infectious disease forecast modeling

Infectious disease forecasting is of great interest to the public health community and policymakers, since forecasts can provide insight into disease dynamics in the near future and inform interventions. Due to delays in case reporting, however, forecasting models may often underestimate the current and future disease burden. In this paper, we propose a general framework for addressing reporting delay in disease forecasting efforts with the goal of improving forecasts. We propose strategies for leveraging either historical data on case reporting or external internet-based data to estimate the amount of reporting error. We then describe several approaches for adapting general forecasting pipelines to account for under- or over-reporting of cases. We apply these methods to address reporting delay in data on dengue fever cases in Puerto Rico from 1990 to 2009 and to reports of influenza-like illness (ILI) in the United States between 2010 and 2019. Through a simulation study, we compare method performance and evaluate robustness to assumption violations. Our results show that forecasting accuracy and prediction coverage almost always increase when correction methods are implemented to address reporting delay. Some of these methods required knowledge about the reporting error or high quality external data, which may not always be available. Provided alternatives include excluding recently-reported data and performing sensitivity analysis. This work provides intuition and guidance for handling delay in disease case reporting and may serve as a useful resource to inform practical infectious disease forecasting efforts.

59 BASIC BIOLOGICAL SCIENCES↗

A new database of building-space-specific internal loads and load schedules for performance based code compliance modeling of commercial buildings

Building-level loads and load profiles prescribed by current modeling rules save modelers time and avoid gaming during whole building performance modeling. However, recent studies show that they sometimes insufficiently capture the entire building performance due to the varied loads and load profiles for different space types. As a solution to this issue, this paper develops a database of building-space-specific loads and load profiles used in code compliance modeling. The existing sets of loads and load profiles are reviewed and the challenges behind using them for specific research topics are discussed. Then, the proposed method to develop the building-space-specific loads and load profiles is introduced. After that, the database for these building-space-specific loads and load profiles is presented. In addition, one case is studied to demonstrate the applications of these loads and load profiles. In this case study, three methods are used to develop building energy models: space-specific (using knowledge of the distribution and location of space types and applying the space-specific data in the developed database), building-level (assuming a lack of knowledge of the space types and using the building-level data in the developed database), and calculated-ratio (assuming knowledge of the distribution of space types but not their locations and calculating weighted average values based on the space-specific data in the developed database). Finally, the energy results simulated by using these three methods are compared, which show building-level methods can produce energy results up to 20% different than the space-specific methods. Finally, this paper discusses the application scope and maintenance of this new database.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Hermite analogs of the lowest order Raviart–Thomas mixed method for convection–diffusion equations

The Raviart–Thomas mixed finite-element method of the lowest order (Raviart and Thomas in mixed finite-element methods for second-order elliptic problems, Lecture Notes in mathematics, Springer, New York, 1977) commonly known as the RT{sub 0} method, is a well-established and popular numerical tool to solve diffusion-like problems providing flux continuity across inter-element boundaries. Douglas and Roberts extended the method to the case of more general second-order boundary-value problems including the convection–diffusion equations (cf. this journal Douglas in Comput Appl Math 1:91–103; 1982). The main drawback of these methods, however, is the poor representation of the primal variable by piecewise constant functions. The Hermite analog of the RT{sub 0} method for treating pure diffusion phenomena proposed in Ruas (J Comput Appl Math, 246:234–242; 2013) proved to be a valid alternative to attain higher order approximation of the primal variable while keeping intact the matrix structure and the quality of the discrete flux variable of the original RT{sub 0} method. Non-trivial extensions of this method are studied here that can be viewed as Hermite analogs of the two Douglas and Roberts’ versions of the RT{sub 0} method, to solve convection–diffusion equations. A detailed convergence study is carried out for one of the Hermite methods, and numerical results illustrate the performance of both of them, as compared to each other and to the corresponding mixed methods.

97 MATHEMATICS AND COMPUTING↗

A performance analysis method for distributed real-time robotic systems: A case study of remote teleoperation

Robot coordination and control systems for remote teleoperation applications are by necessity implemented on distributed computers. Modeling and performance analysis of these distributed robotic systems is difficult, but important for economic system design. Performance analysis methods originally developed for conventional distributed computer systems are often unsatisfactory for evaluating real-time systems. The paper introduces a formal model of distributed robotic control systems; and a performance analysis method, based on scheduling theory, which can handle concurrent hard-real-time response specifications. Use of the method is illustrated by a case of remote teleoperation which assesses the effect of communication delays and the allocation of robot control functions on control system hardware requirements.

Lefebvre, D. R.↗

Integrated approaches for reducing sample size for measurements of trace elemental impurities in plutonium by ICP-OES and ICP-MS

This study has demonstrated that bulk plutonium chemical analysis can be performed at small scales (\50 mg material) through three case studies. Analytical methods were developed for ICP-OES and ICP-MS instruments to measure trace impurities and gallium content in plutonium metals with comparable or improved detection limits, measurement accuracy and precision. In two case studies, the sample size has been reduced by 109, and in the third case study, by as much as 50009, so that the plutonium chemical analysis can be performed in a facility rated for lower-hazard and lower-security operations.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Comparison of approximate solutions to the phonon Boltzmann transport equation with the relaxation time approximation: Spherical harmonics expansions and the discrete ordinates method

In this study, the phonon Boltzmann transport equation is used to analyze model problems in one and two spatial dimensions, under transient and steady-state conditions. New, explicit solutions are obtained by using the P 1 and P 3 approximations, based on expansions in spherical harmonics, and are compared with solutions from the discrete ordinates method. For steady-state energy transfer, it is shown that analytic expressions derived using the P 1 and P 3 approximations agree quantitatively with the discrete ordinates method, in some cases for large Knudsen numbers, and always for Knudsen numbers less than unity. However, for time-dependent energy transfer, the P N solutions differ qualitatively from converged solutions obtained by the discrete ordinates method. Although they correctly capture the wave-like behavior of energy transfer at short times, the P 1 and P 3 approximations rely on one or two wave velocities, respectively, yielding abrupt, step-changes in temperature profiles that are absent when the angular dependence of the phonon velocities is captured more completely. It is shown that, with the gray approximation, the P 1 approximation is formally equivalent to the so-called “hyperbolic heat equation.” In conclusion, these results support the use of the P N approximation to find solutions to the phonon Boltzmann transport equation for steady-state conditions. Such solutions can be useful in the design and analysis of devices that involve heat transfer at nanometer length scales, where continuum-scale approaches become inaccurate.

36 MATERIALS SCIENCE↗

Comparative Architecture Analysis: Using Laminar Structure to Unify Multiple Grid Architectures

A great many architectures for various aspects of electric grids have been proposed over recent years and new ones continue to be suggested. Generally speaking, there is no means except ad hoc analysis by which one can reduce these architectures to some common basis, or extract common elements. This leads to confusion in choosing an architecture for implementation, as well as making it difficult to identify core product or platform opportunities with sufficient potential markets to encourage vendors to provide or develop them. The discipline of Grid Architecture as defined for the US Department of Energy (DOE) in 2014 has multiple purposes. Among them are two that are relevant here: Providing rigorous bases for architectural structures and Identifying common aspects of apparently divergent architectural structures. For the first of these, mathematical methods are applied where possible to provide the rigorous basis for structure. For the second, mapping of architectural structures to a canonical structure puts apparently disparate architectures onto a common basis. In the best case, these two methods coincide. In this paper, Laminar Coordination structure serves as the rigorously defined common basis for analysis of several real and proposed grid architectures.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Interrelation Between Safety Factors and Reliability

An evaluation was performed to establish relationships between safety factors and reliability relationships. Results obtained show that the use of the safety factor is not contradictory to the employment of the probabilistic methods. In many cases the safety factors can be directly expressed by the required reliability levels. However, there is a major difference that must be emphasized: whereas the safety factors are allocated in an ad hoc manner, the probabilistic approach offers a unified mathematical framework. The establishment of the interrelation between the concepts opens an avenue to specify safety factors based on reliability. In cases where there are several forms of failure, then the allocation of safety factors should he based on having the same reliability associated with each failure mode. This immediately suggests that by the probabilistic methods the existing over-design or under-design can be eliminated. The report includes three parts: Part 1-Random Actual Stress and Deterministic Yield Stress; Part 2-Deterministic Actual Stress and Random Yield Stress; Part 3-Both Actual Stress and Yield Stress Are Random.

Elishakoff, Isaac↗

Impact of Consolidation Radiation Therapy in Stage III-IV Diffuse Large B-cell Lymphoma With Negative Post-Chemotherapy Radiologic Imaging

Purpose: While consolidation radiation therapy (i.e., RT administered after chemotherapy) is routine treatment for patients with early-stage diffuse large B-cell lymphoma (DLBCL), the role of consolidation RT in stage III-IV DLBCL is controversial. Methods and Materials: Cases of patients with stage III-IV DLBCL treated from 1991 to 2009 at Duke University, who achieved a complete response to chemotherapy were reviewed. Clinical outcomes were calculated using the Kaplan-Meier method and were compared between patients who did and did not receive RT, using the log-rank test. A multivariate analysis was performed using Cox proportional hazards model. Results: Seventy-nine patients were identified. Chemotherapy (median, 6 cycles) consisted of anti-CD20 antibody rituximab combined with cyclophosphamide, doxorubicin, vincristine, and prednisone (R-CHOP; 65%); cyclophosphamide, doxorubicin, vincristine, and prednisone (CHOP; 22%); or other (13%). Post-chemotherapy imaging consisted of positron emission tomography (PET)/computed tomography (CT) (73%); gallium with CT (14%); or CT only (13%). Consolidation RT (median, 25 Gy) was given to involved sites of disease in 38 (48%) patients. Receipt of consolidation RT was associated with improved in-field control (92% vs. 69%, respectively, p = 0.028) and event-free survival (85% vs. 65%, respectively, p = 0.014) but no difference in overall survival (85% vs. 78%, respectively, p = 0.15) when compared to patients who did not receive consolidation RT. On multivariate analysis, no RT was predictive of increased risk of in-field failure (hazard ratio [HR], 8.01, p = 0.014) and worse event-free survival (HR, 4.3, p = 0.014). Conclusions: Patients with stage III-IV DLBCL who achieve negative post-chemotherapy imaging have improved in-field control and event-free survival with low-dose consolidation RT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Approaching the complete basis set limit of CCSD(T) for large systems by the third-order incremental dual-basis set zero-buffer F12 method

The third-order incremental dual-basis set zero-buffer approach was combined with CCSD(T)-F12x (x = a, b) theory to develop a new approach, i.e., the inc3-db-B0-CCSD(T)-F12 method, which can be applied as a black-box procedure to efficiently obtain the near complete basis set (CBS) limit of the CCSD(T) energies also for large systems. We tested this method for several cases of different chemical nature: four complexes taken from the standard benchmark sets S66 and X40, the energy difference between isomers of water hexamer and the rotation barrier of biphenyl. The results show that our method has an error relative to the best estimation of CBS energy of only 0.2 kcal/mol or less. By parallelization, our method can accomplish the CCSD(T)-F12 calculations of about 60 correlated electrons and 800 basis functions in only several days, which by standard implementation are impossible for ordinary hardware. We conclude that the inc3-db-B0-CCSD(T)-F12a/AVTZ method, which is of CCSD(T)/AV5Z quality, is close to the limit of accuracy that one can achieve for large systems currently.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A fourth-order box method for solving the boundary layer equations

A fourth order box method for calculating high accuracy numerical solutions to parabolic, partial differential equations in two variables or ordinary differential equations is presented. The method is the natural extension of the second order Keller Box scheme to fourth order and is demonstrated with application to the incompressible, laminar and turbulent boundary layer equations. Numerical results for high accuracy test cases show the method to be significantly faster than other higher order and second order methods.

Wornom, S. F.↗

L2 Milestone: Neutron Capture Cross Sections from Surrogate (p, d) Measurements: Determination of the Unknown 87 Y(n, g) Cross Section and Assessment of the Method Via the 90 Zr(n, g) Benchmark Case: Theory Report

Cross sections for compound-nuclear reactions involving unstable targets are important for many applications, but can often not be measured directly. Here we describe a method for extracting cross sections for neutron-capture on unstable isotopes from indirect (surrogate) measurements. The surrogate reaction, which produces the compound nucleus of interest, has to be described and the decay of the nucleus has to be modeled. We outline the approach for one-neutron pickup and report on the determination of the 90 Zr(n, γ ) reaction from surrogate 92 Zr(p,d) data, which is compared to the directly-measured capture cross section and thus provides a benchmark for the method. We then apply the method to determine the 87 Y(n, γ ) cross section, which has not been measured directly. The work was carried out in the context of an LLNL L2 Milestone. This report addresses the theory aspects of the milestone. A complementary document summarizes the experimental efforts [1].

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modeling the Moisture Content and Dry Matter Loss in Dynamic Woody Biomass Storage Piles with Variable Extraction

The urgent need to mitigate climate change has spurred significant interest in renewable energy sources. This paper explores the storage and processing of woody biomass for biofuel production, considering the dynamic nature of biomass piles in real-world scenarios. A model has been developed to analyze moisture content changes and dry matter loss in woody biomass stored in piles prior to processing, taking into account varying extraction methods and environmental conditions. Case studies that examine the effects of different unpiling methods (FIFO, LIFO, and homogeneous) on moisture content and dry matter loss under various feedstock arrival rates and weather conditions are presented. Results indicate that unpiling methods significantly impact moisture content, with LIFO typically resulting in higher moisture content due to the utilization of fresher feedstock. Dry matter loss increases with pile size and time, emphasizing the importance of accurate modeling for assessing carbon emissions and feedstock quality. Furthermore, the model highlights the importance of process loading order and extraction methods in biomass storage, suggesting potential cost benefits associated with decreased moisture content. The difference between different extraction methods could vary the moisture content in the feedstock reaching the biofuel process by as much as 37.6%, however dry matter loss varies minimally for realistic pile changes. Overall, this research contributes to a better understanding of biomass storage dynamics and informs sustainable biofuel production practices.

Niska, Janel↗

Modeling the Structural Response from a Propagating High Explosive Using Smooth Particle Hydrodynamics

This report primarily concerns the use of two massively parallel finite element codes originally written and maintained at Lawrence Livermore National Laboratory. ALE3D is an explicit hydrodynamics code commonly employed to simulate wave propagation from high energy scenarios and the resulting interaction with nearby structures. This coupled response ensures that a structure is accurately applied with a blast loading varying both in space and time. Figure 1 illustrates the radial outward propagation of a pressure wave due to a center detonated spherical explosive originating from the lower left. The radial symmetry seen in this scenario is lost when instead a cylindrocal charge is detonated. Figure 2 indicates that a stronger, faster traveling pressure wave occurs in the direction of the normal axis to the cylinder. The ALE3D name is derived because of the use of arbitrary-Lagrange-Eulerian elements in which the mesh is allowed to advect; a process through which the mesh is modified to alleviate tanlging and general mesh distortion often cuased by high energy scenarios. The counterpart to an advecting element is a Lagrange element, whose mesh moves with the material. Ideally all structural components are kept Lagrange as long as possible to preserve accuracy of material variables and minimize advection related errors. Advection leads to mixed zoning, so using structural Lagrange elements also improves the visualization when post processing the results. A simplified representation of the advection process is shown in Figure 3. First the mesh is distorted due to material motion during the Lagrange step. The mesh is then shifted to an idealized and less distorted state to prevent irregular zones caused by the Lagrange motion. Lastly, the state variables are remapped to the elements of the newly constructed mesh. Note that Figure 3 represents a purely Eulerian mesh relaxation because the mesh is relocated back to the pre-Lagrange position. This is the case when the material flows through a still mesh. This is not typically done in an ALE3D analysis, especially if Lagrange elements exist. Deforming Lagrange elements would certainly tangle with a Eulerian mesh eventually. The best method in this case is to have an advecting mesh positioned as some relaxed version of the pre and post Lagrange step; this gives the best opportunity of modeling a high energy event with a combination of Lagrange and ALE elements. Dyne3D is another explicit dynamic analysis code, ParaDyn being the parallel version. ParaDyn is used for predicting the transient response of three dimensional structures using Lagrangian solid mechanics. Large deformation and mesh tangling is often resolved through the use of an element deletion scheme. This is useful to accommodate component failure, but if it is done purely as a means to preserve a useful mesh it can lead to problems because it does not maintain continuity of the material bulk response. Whatever medium exists between structural components is typically not modeled in ParaDyn. Instead, a structure either has a known loading profile applied or given initial conditions. The many included contact algorithms can calculate the loading response of materials if and when they collide. A recent implementation of an SPH module in which failed or deleted material nodes are converted to independent particles is currently being utilized for a variety of spall related problems and high velocity impact scenarios. Figure 4 shows an example of a projectile, given an initial velocity, and how it fails the first plate which generates SPH particles which then interact with and damage the second plate.

36 MATERIALS SCIENCE↗

Characterization of nanosized Al{sub 2}(WO{sub 4}){sub 3}

Graphical abstract: TEM method allows to detect small quantities of impurities not detectable by other methods. In our case impurities of W{sub 5}O{sub 14} are detected in Al{sub 2}(WO{sub 4}){sub 3} nanopowder. Highlights: {yields} Nanosized Al{sub 2}(WO{sub 4}){sub 3} by simple co-precipitation method. {yields} Spherical particles with mean size of 22 nm distributed between 10 and 40 nm at 630 {sup o}C. {yields} XRD, DTA and TEM confirm well defined products with perfect structure. {yields} TEM locality allows detection of impurities not detectable by XRD and DTA. -- Abstract: Nanosized aluminum tungstate Al{sub 2}(WO{sub 4}){sub 3} was prepared by co-precipitation reaction between Na{sub 2}WO{sub 4} and Al(NO{sub 3}){sub 3} aqueous solutions. The powder size and shape, as well as size distribution are estimated after different conditions of powder preparation. The purity of the final product was investigated by XRD and DTA analyses, using the single crystal powder as reference. Between the specimen and the reference no difference was detected. The crystal structure of Al{sub 2}(WO{sub 4}){sub 3} nanosized powder was confirmed by TEM (SAED, HRTEM). In additional, TEM locality allows to detect some W{sub 5}O{sub 14} impurities, which are not visible by conventional X-ray powder diffraction and thermal analyses.

36 MATERIALS SCIENCE↗

Multi-reference approach to the calculation of photoelectron spectra including spin-orbit coupling

X-ray photoelectron spectra provide a wealth of information on the electronic structure. The extraction of molecular details requires adequate theoretical methods, which in case of transition metal complexes has to account for effects due to the multi-configurational and spin-mixed nature of the many-electron wave function. Here, the restricted active space self-consistent field method including spin-orbit coupling is used to cope with this challenge and to calculate valence- and core-level photoelectron spectra. The intensities are estimated within the frameworks of the Dyson orbital formalism and the sudden approximation. Thereby, we utilize an efficient computational algorithm that is based on a biorthonormal basis transformation. The approach is applied to the valence photoionization of the gas phase water molecule and to the core ionization spectrum of the [Fe(H{sub 2}O){sub 6}]{sup 2+} complex. The results show good agreement with the experimental data obtained in this work, whereas the sudden approximation demonstrates distinct deviations from experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗