Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Recursion”

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 559 records · Page 31

Modeling occupancy-driven building loads for large and diversified building stocks through the use of parametric schedules

Building energy modeling provides a fundamental tool to assess the potential for energy efficiency to contribute to reducing world energy consumption and global emissions. Occupancy-related operations are a key source of uncertainty for building energy analysis, particularly for aggregated building stocks. At a district or city level, it is critical to estimate aggregated power load profiles for sizing power grid infrastructure, power plant capacity allocation, and energy efficiency measures. The stochastic nature of behavior-related operations complicates the creation of models that accurately capture building load profiles for entire building stocks. This research introduces a new methodology called parametric schedules to model occupancy-driven schedules for large and diverse building stocks. In contrast to computationally expensive methodologies proposed in the literature, our work does not use a recursive time-consuming step. Occupancy is estimated by the extrapolation of operation times directly from metered electric consumption data; occupancy-related schedules are stochastically assigned to each building model, guaranteeing diversity of operation times in the stock. Our procedure has been tested on a large, diversified data-set of 25,000 commercial buildings in Los Angeles, California. It proved to be able to adequately represent the stochastic schedules diversity of the stock and to refine the stock calibration process by 1%. This innovative approach represents a useful asset for utility companies, grid operators, urban planners, and balancing authorities, seeking to improve building stock modeling and better estimate the impact of energy conservation measures. – This work is part of a larger stock modeling tool called ComStock, which is under development by NREL.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Logical design of anti-prion agents using NAGARA

To accelerate the logical drug design procedure, we created the program “NAGARA,” a plugin for PyMOL, and applied it to the discovery of small compounds called medical chaperones (MCs) that stabilize the cellular form of a prion protein (PrP{sup C}). In NAGARA, we constructed a single platform to unify the docking simulation (DS), free energy calculation by molecular dynamics (MD) simulation, and interfragment interaction energy (IFIE) calculation by quantum chemistry (QC) calculation. NAGARA also enables large-scale parallel computing via a convenient graphical user interface. Here, we demonstrated its performance and its broad applicability from drug discovery to lead optimization with full compatibility with various experimental methods including Western blotting (WB) analysis, surface plasmon resonance (SPR), and nuclear magnetic resonance (NMR) measurements. Combining DS and WB, we discovered anti-prion activities for two compounds and tegobuvir (TGV), a non-nucleoside non-structural protein NS5B polymerase inhibitor showing activity against hepatitis C virus genotype 1. Binding profiles predicted by MD and QC are consistent with those obtained by SPR and NMR. Free energy analyses showed that these compounds stabilize the PrP{sup C} conformation by decreasing the conformational fluctuation of the PrP{sup C}. Because TGV has been already approved as a medicine, its extension to prion diseases is straightforward. Finally, we evaluated the affinities of the fragmented regions of TGV using QC and found a clue for its further optimization. By repeating WB, MD, and QC recursively, we were able to obtain the optimum lead structure. - Highlights: • NAGARA integrates docking simulation, molecular dynamics, and quantum chemistry. • We found many compounds, e.g., tegobuvir (TGV), that exhibit anti-prion activities. • We obtained insights into the action mechanism of TGV as a medical chaperone. • Using QC, we obtained useful information for optimization of the lead compound, TGV. • NAGARA is a convenient platform for drug discovery and lead optimization.

60 APPLIED LIFE SCIENCES↗

An efficient method to integrate polynomials over polytopes and curved solids

Here in this paper, we present an efficient approach to compute the integral of monomials and polynomials over polyhedra and regions defined by parametric curved boundary surfaces. We use Euler's theorem for homogeneous functions in combination with Stokes's theorem to reduce the integration of a monomial over a three-dimensional solid to its boundary. If the solid is a polytope, through a recursive application of these theorems, the integral is further reduced to just the evaluation of the monomial and its derivatives at the vertices of the polytope. The present approach is simpler than existing techniques that rely on repeated use of the divergence theorem, which require the antiderivative of the monomials and the projection of these functions onto hyperplanes. For convex and nonconvex polytopes, our approach does not introduce any approximation for the integration of monomials. For curved solid regions bounded by surfaces that admit a parameterization, the same approach yields simplified formulas to compute the integral of any homogeneous function, including monomials. For surfaces parameterized by polynomial surfaces (such as Bezier surface triangles and B-spline patches), the method yields machine-precision accuracy for the volumetric integration of monomials with an appropriate quadrature rule. Numerical examples over regions bounded by polynomial surfaces and rational surfaces are presented to establish the accuracy and efficiency of the method.

97 MATHEMATICS AND COMPUTING↗

The introspective may achieve more: Enhancing existing Geoscientific models with native-language emulated structural reflection

Geoscientific models manage myriad and increasingly complex data structures as trans -disciplinary models are integrated. They often incur significant redundancy with cross-cutting tasks. Reflection, the ability of a program to inspect and modify its structure and behavior at runtime, is known as a powerful tool to improve code reusability, abstraction, and separation of concerns. Reflection is rarely adopted in high-performance Geoscientific models, especially with Fortran, where it was previously deemed implausible. Practical constraints of language and legacy often limit us to feather-weight, native-language solutions. We demonstrate the usefulness of a structural-reflection-emulating, dynamically-linked metaObjects, gd . We show real-world examples including data structure self-assembly, effortless save/restart and upgrade to parallel I/O, recursive actions and batch operations. We share gd and a derived module that reproduces MATLAB-like structure in Fortran and C++. We suggest that both a gd representation and a Fortran-native representation are maintained to access the data, each for separate purposes. In conclusion, embracing emulated reflection allows generically-written codes that are highly re-usable across projects.

58 GEOSCIENCES↗

Computational assessment of smooth and rough parameter dependence of statistics in chaotic dynamical systems

An assumption of smooth response to small parameter changes, of statistics or long-time averages of a chaotic system, is generally made in the field of sensitivity analysis, and the parametric derivatives of statistical quantities are critically used in science and engineering. In this paper, we propose a numerical procedure to assess the differentiability of statistics with respect to parameters in chaotic systems. Here, we numerically show that the existence of the derivative depends on the Lebesgue integrability of a certain density gradient function, which we define as the derivative of logarithmic SRB density along the unstable manifold. We develop a recursive formula for the density gradient that can be efficiently computed along trajectories, and demonstrate its use in determining the differentiability of statistics. Our numerical procedure is illustrated on low-dimensional chaotic systems whose statistics exhibit both smooth and rough regions in parameter space.

42 ENGINEERING↗

Predicting the mechanical response of oligocrystals with deep learning

In this study, we employ data-driven homogenization approaches to predict the particular mechanical evolution of polycrystalline aggregates with tens of individual crystals. In these oligocrystals the differences in stress response due to microstructural variation is pronounced. Shell-like structures produced by metal-based additive manufacturing and the like make the prediction of the behavior of oligocrystals technologically relevant. The predictions of traditional homogenization theories based on grain volumes are not sensitive to variations in local grain neighborhoods. Direct simulation of the local response with crystal plasticity finite element methods is more detailed, but the computations are expensive. To represent the stress-strain response of a polycrystalline sample given its initial grain texture and morphology we have designed a novel neural network that incorporates a convolution component to observe and reduce the information in the crystal texture field and a recursive component to represent the causal nature of the history information. This model exhibits accuracy on par with crystal plasticity simulations at minimal computational cost per prediction.

36 MATERIALS SCIENCE↗

Parallel adaptive fluid-structure interaction simulation of explosions impacting on building structures

We pursue a level set approach to couple an Eulerian shock-capturing fluid solver with space-time refinement to an explicit solid dynamics solver for large deformations and fracture. The coupling algorithms considering recursively finer fluid time steps as well as overlapping solver updates are discussed in detail. Our ideas are implemented in the AMROC adaptive fluid solver framework and are used for effective fluid-structure coupling to the general purpose solid dynamics code DYNA3D. Beside simulations verifying the coupled fluid-structure solver and assessing its parallel scalability, the detailed structural analysis of a reinforced concrete column under blast loading and the simulation of a prototypical blast explosion in a realistic multistory building are presented.

Deiterding, Ralf↗

Branches, quivers, and ideals for knot complements

We generalize the F K invariant, i.e. Z^ for the complement of a knot K in the 3-sphere, the knots-quivers correspondence, and A-polynomials of knots, and find several interconnections between them. We associate an F K invariant to any branch of the A-polynomial of K and we work out explicit expressions for several simple knots. We show that these F K invariants can be written in the form of a quiver generating series, in analogy with the knots-quivers correspondence. We discuss various methods to obtain such quiver representations, among others using R-matrices. We generalize the quantum a-deformed A-polynomial to an ideal that contains the recursion relation in the group rank, i.e. in the parameter a, and describe its classical limit in terms of the Coulomb branch of a 3d-5d theory. We also provide t-deformed versions. Furthermore, we study how the quiver formulation for closed 3-manifolds obtained by surgery leads to the superpotential of 3d N = 2 theory T[M 3 ] and to the data of the associated modular tensor category MTC[M 3 ].

97 MATHEMATICS AND COMPUTING↗

A note on efficiently generating ionic configurations for opacity calculations

When calculating the spectral opacity of hot dense plasmas one often encounters the need to generate a list of detailed ionic configurations of bound states for each ion stage in the plasma. Here we present here a non-recursive algorithm for the efficient construction of such a list of states.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the computation of moments in the Super-Transition-Arrays model for radiative opacity calculations

In the Super-Transition-Array statistical method for the computation of radiative opacity of hot dense matter, the moments of the absorption or emission features involve partition functions with reduced degeneracies, occurring through the calculation of averages of products of subshell populations. Here, in the present work, we discuss several aspects of the computation of such peculiar partition functions, insisting on the precautions that must be taken in order to avoid numerical difficulties. In a previous work, we derived a formula for supershell partition functions, which takes the form of a functional of the distribution of energies within the supershell and allows for fast and accurate computations, truncating the number of terms in the expansion. The latter involves coefficients for which we obtained a recursion relation and an explicit formula. We show that such an expansion can be combined with the recurrence relation for shifted partition functions. We also propose, neglecting the effect of fine structure as a first step, a positive-definite formula for the Super-Transition-Array moments of any order, providing an insight into the asymmetry and sharpness of the latter. The corresponding formulas are free of alternating sums. Several ways to speed up the calculations are also presented.

74 ATOMIC AND MOLECULAR PHYSICS↗

An Error-Entropy Minimization Algorithm for Tracking Control of Nonlinear Stochastic Systems with Non-Gaussian Variables

This paper presents an error-entropy minimization tracking control algorithm for a class of dynamic stochastic system. The system is represented by a set of time-varying discrete nonlinear equations with non-Gaussian stochastic input, where the statistical properties of stochastic input are unknown. By using Parzen windowing with Gaussian kernel to estimate the probability densities of errors, recursive algorithms are then proposed to design the controller such that the tracking error can be minimized. The performance of the error-entropy minimization criterion is compared with the mean-square-error minimization in the simulation results.

Entropy↗

A method to identify weak points of interconnection of renewable energy resources

The increasing penetration of renewable energy resources (RERs) is challenging power grid planning and operation in maintaining grid reliability. Potential power grid reliability and stability issues may arise when large-scale RERs are connected to a weak power grid. These issues may become more significant due to the mutual interactions among RERs that are interconnected electrically close through a power network. Weakness analysis based on grid strength assessment is useful for identifying the potential weak grid issues. However, when taking into account the impact of the interactions among RERs, the weakness analysis becomes computationally challenging. Different combinations of points of interconnections (POIs) of RERs may have different impacts on grid strength at each POI. Due to the combination nature, such weakness analysis may be time-consuming when identifying the weakest combination of POIs from a large number of potential candidates in realistic power grids. To address this concern, the impact of the power network structure on grid strength is analyzed in this paper. It is found that the impact mainly depends on the impedance ratios at POIs. The findings can be summarized as two criteria to develop a recursive method for fast identification of the weakest combination in this paper. The efficacy of the proposed method is demonstrated through case studies on the IEEE 39-bus system and a realistic power grid.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Decentralized Data-Driven Estimation of Generator Rotor Speed and Inertia Constant Based on Adaptive Unscented Kalman Filter

This paper proposes an online decentralized and data-driven method for synchronous generator rotor speed and inertia estimation. First, an analytical relationship between the generator terminal voltage/current phasors and its internal rotor speed is derived using the Thevenin equivalent. The parameters of the latter are obtained via a recursive estimator, which allows us to estimate the generator's internal rotor speed from voltage/current phasors without the need for any generator information. Second, we reformulate the equivalent swing equation into the Kalman filter form, which enables us to use a priori-estimated internal rotor speed along with measured real and reactive power injections to obtain an estimate of the inertia for each online synchronous generator. This is achieved through an adaptive unscented Kalman filter. We demonstrate that the use of terminal bus frequency to approximate the rotor speed for inertia estimation by existing methods leads to large estimation biases/errors. Numerical simulations carried out on the IEEE 39-bus and Texas 2000-bus systems reveal that our proposed method is consistently more accurate than existing methods. Moreover, our proposed method is insensitive to both the type and location of system disturbances, and it is robust to different levels of measurement noise.

bus frequency↗

Extended real-time voltage instability identification method based on synchronized phasor measurements

This paper presents an extended adaptive approach designed to accurately estimate the Thévenin equivalent parameters using phasor measurements at a given bus for measurements lying in any of four quadrants of the PQ-plane. The improvement is achieved by using a new condition to properly update the estimated parameters after an initial guess. Based on an adaptive philosophy, the proposed approach can correctly account for the intrinsic nonlinearities of a power system, can provide a real-time estimation of Thévenin parameters, and does not require network topology knowledge. The method is validated using the Kundur 2-area system, showing estimation improvements compared to the current adaptive approach and the classical recursive least-squares method. The proposed approach is able to estimate both sides of the system with respect to the measurement bus. In addition, a data-driven voltage stability index is developed. To illustrate the performance of the proposed approach in a larger power system, a voltage stability assessment is carried out on the IEEE 39-bus system, considering the action of overexcitation limiters of generators and nonlinear loads. The proposed approach is suitable for applications that require an accurate Thévenin equivalent estimation in real-time, such as for voltage stability assessment. The new approach provides a reliable tool for the system operators to make proper and timely decisions.

42 ENGINEERING↗

Integrating Epidermal Growth Factor Receptor Assay With Clinical Parameters Improves Risk Classification for Relapse and Survival in Head-and-Neck Squamous Cell Carcinoma

Purpose: Epidermal growth factor receptor (EGFR) overexpression has been consistently found to be an independent predictor of local-regional relapse (LRR) after radiotherapy. We assessed the extent by which it can refine risk classification for overall survival (OS) and LRR in patients with head-and-neck squamous cell carcinoma (HNSCC). Methods and Materials: EGFR expression in locally advanced HNSCC was measured by immunohistochemistry in a series of patients randomized to receive accelerated or conventional radiation regimens in a Phase III trial. Subsequently, data of the two series were pooled (N = 533) for conducting a recursive partitioning analysis that incorporated clinical parameters (e.g., performance status, primary site, T and N categories) and four molecular markers (EGFR, p53, Ki-67, and microvessel density). Results: This study confirmed that patients with higher than median levels of tumor EGFR expression had a lower OS (relative risk [RR]: 1.90, p = 0.0010) and a higher LRR (RR: 1.91, p = 0.0163). Of the four markers analyzed, only EGFR was found to contribute to refining classification of patients into three risk classes with distinct OS and LRR outcomes. The addition of EGFR to three clinical parameters could identify patients having up to a fivefold difference in the risk of LRR. Conclusions: Adding pretreatment EGFR expression data to known robust clinical prognostic variables improved the estimation of the probability for OS and LRR after radiotherapy. Its use for stratifying or selecting patients with defined tumor feature and pattern of relapse for enrollment into clinical trials testing specific therapeutic strategy warrants further investigation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Shorter-Course Whole-Brain Radiotherapy for Brain Metastases in Elderly Patients

Purpose: Many patients with brain metastases receive whole-brain radiotherapy (WBRT) alone. Using 10 Multiplication-Sign 3 Gy in 2 weeks is the standard regimen in most centers. Regarding the extraordinarily poor survival prognosis of elderly patients with multiple brain metastases, a shorter WBRT regimen would be preferable. This study compared 10 Multiplication-Sign 3 Gy with 5 Multiplication-Sign 4 Gy in elderly patients ({>=}65 years). Methods and Materials: Data from 455 elderly patients who received WBRT alone for brain metastases were retrospectively analyzed. Survival and local (= intracerebral) control of 293 patients receiving 10 Multiplication-Sign 3 Gy were compared with 162 patients receiving 5 Multiplication-Sign 4 Gy. Eight additional potential prognostic factors were investigated including age, gender, Karnofsky performance score (KPS), primary tumor, number of brain metastases, interval from tumor diagnosis to WBRT, extracerebral metastases, and recursive partitioning analysis (RPA) class. Results: The 6-month overall survival rates were 29% after 5 Multiplication-Sign 4 Gy and 21% after 10 Multiplication-Sign 3 Gy (p = 0.020). The 6-month local control rates were 12% and 10%, respectively (p = 0.32). On multivariate analysis, improved overall survival was associated with KPS {>=} 70 (p < 0.001), only one to three brain metastases (p = 0.029), no extracerebral metastasis (p = 0.012), and lower RPA class (p < 0.001). Improved local control was associated with KPS {>=} 70 (p < 0.001), breast cancer (p = 0.029), and lower RPA class (p < 0.001). Conclusions: Shorter-course WBRT with 5 Multiplication-Sign 4 Gy was not inferior to 10 Multiplication-Sign 3 Gy with respect to overall survival or local control in elderly patients. 5 Multiplication-Sign 4 Gy appears preferable for the majority of these patients.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Effect of Tumor Subtype on Survival and the Graded Prognostic Assessment for Patients With Breast Cancer and Brain Metastases

Purpose: The diagnosis-specific Graded Prognostic Assessment (GPA) was published to clarify prognosis for patients with brain metastases. This study refines the existing Breast-GPA by analyzing a larger cohort and tumor subtype. Methods and Materials: A multi-institutional retrospective database of 400 breast cancer patients treated for newly diagnosed brain metastases was generated. Prognostic factors significant for survival were analyzed by multivariate Cox regression and recursive partitioning analysis (RPA). Factors were weighted by the magnitude of their regression coefficients to define the GPA index. Results: Significant prognostic factors by multivariate Cox regression and RPA were Karnofsky performance status (KPS), HER2, ER/PR status, and the interaction between ER/PR and HER2. RPA showed age was significant for patients with KPS 60 to 80. The median survival time (MST) overall was 13.8 months, and for GPA scores of 0 to 1.0, 1.5 to 2.0, 2.5 to 3.0, and 3.5 to 4.0 were 3.4 (n = 23), 7.7 (n = 104), 15.1 (n = 140), and 25.3 (n = 133) months, respectively (p < 0.0001). Among HER2-negative patients, being ER/PR positive improved MST from 6.4 to 9.7 months, whereas in HER2-positive patients, being ER/PR positive improved MST from 17.9 to 20.7 months. The log-rank statistic (predictive power) was 110 for the Breast-GPA vs. 55 for tumor subtype. Conclusions: The Breast-GPA documents wide variation in prognosis and shows clear separation between subgroups of patients with breast cancer and brain metastases. This tool will aid clinical decision making and stratification in clinical trials. These data confirm the effect of tumor subtype on survival and show the Breast-GPA offers significantly more predictive power than the tumor subtype alone.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Gamma Knife Radiosurgery for Treatment of Cerebral Metastases From Non-Small-Cell Lung Cancer

Purpose: To evaluate clinical and physico-dosimetric variables affecting clinical outcome of patients treated with Gamma Knife radiosurgery (GKRS) for brain metastases from non-small cell lung cancer (NSCLC). Methods and Materials: Between 2001 and 2006, 373 patients (298 men and 75 women, median age 65 years) with brain metastases from NSCLC underwent GKRS. All of them had KPS {>=} 60%, eight or fewer brain metastases, confirmed histopathological diagnosis and recent work-up (<3 months). Thirty-five patients belonged to recursive partitioning analysis (RPA) Class I, 307 patients were in RPA Class II, 7 patients were in RPA Class III. Median tumor volume was 3.6 cm{sup 3}. Median marginal dose was 22.5 Gy at 50% isodose.; median 10 Gy and 12 Gy isodose volumes were 30.8 cm{sup 3} and 15.8 cm{sup 3}, respectively. Follow-up with MRI was performed every 3 months. Overall survival data were collected from internal database, telephone interviews, and identifying registries. Results: Mean follow-up after GKRS was 51 months (range, 6 to 96 months); mean overall survival was 14.2 months. Of 373 patients, 29 were alive at time of writing, 104 had died of cerebral progression, and 176 had died of systemic progression. In 64 cases it was not possible to ascertain the cause. Univariate and multivariate analysis were adjusted for the following: RPA class, surgery, WBRT, age, gender, number of lesions, median tumor volume, median peripheral dose, and 10 Gy and 12 Gy volumes. Identified RPA class and overall tumor volume >5 cc were the only two covariates independently predictive of overall survival in patients who died of cerebral progression. Conclusions: Global volume of brain disease should be the main parameter to consider for performing GKRS, which is a first-line therapy for patient in good general condition and controlled systemic disease.

62 RADIOLOGY AND NUCLEAR MEDICINE↗