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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Soft bootstrap and effective field theories

The soft bootstrap program aims to construct consistent effective field theories (EFT’s) by recursively imposing the desired soft limit on tree-level scattering amplitudes through on-shell recursion relations. A prime example is the leading two-derivative operator in the EFT of SU(N) × SU(N)/SU(N) nonlinear sigma model (NLSM), where O(p 2 ) amplitudes with an arbitrary multiplicity of external particles can be soft-bootstrapped. We extend the program to O(p 4 ) operators and introduce the “soft blocks,” which are the seeds for soft bootstrap. The number of soft blocks coincides with the number of independent operators at a given order in the derivative expansion and the incalculable Wilson coefficient emerges naturally. We also uncover a new soft-constructible EFT involving the “multi-trace” operator at the leading two-derivative order, which is matched to SO(N + 1)/SO(N) NLSM. In addition, we consider Wess-Zumino-Witten (WZW) terms, the existence of which, or the lack thereof, depends on the number of flavors in the EFT, after a novel application of Bose symmetry. Remarkably, we find agreements with grouptheoretic considerations on the existence of WZW terms in SU(N) NLSM for N ≥ 3 and the absence of WZW terms in SO(N) NLSM for N 6= 5.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Open associahedra and scattering forms

We continue the study of open associahedra associated with bi-color scattering amplitudes initiated in ref. [1]. We focus on the facet geometries of the open associahedra, uncovering many new phenomena such as fiber-product geometries. We then provide novel recursion procedures for calculating the canonical form of open associahedra, generalizing recursion relations for bounded polytopes to unbounded polytopes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

A Refined Global Segmentation of X-Ray CT Images for Multi-phase Geomaterials

X-ray computed tomography images of three-phase silica sand and glass bead specimens are analyzed and used to evaluate the segmentation performances of Otsu-, and recursion-based multilevel algorithms. A global image segmentation technique that combines iterative and recursive algorithms, namely a refined statistics-based global segmentation is proposed for segmenting multi-phase granular geomaterials. The performance of the proposed algorithm is tested by segmenting partially saturated silica sand and glass bead specimens. For the tested silica sand specimen, the refined statistics method estimated void ratio and degree of saturation were 0.67 and 39.35%. The estimates for the glass bead specimen yielded 0.64 and 43.49%, respectively. The true void ratio (0.66) and degree of saturation (37.71%) were determined with a user-controlled Image processing software package—Image-Pro. It was found that the proposed method estimated the void ratio and the degree of saturation with 1.52 and 4.35 percent errors for the silica sand and with 15.63 and 0.34 percent errors for the glass bead, respectively. The computational time of the proposed method was found to be shorter than other methods considered. Overall, it is concluded that the proposed technique performed better in segmenting three-phase granular geomaterials.

42 ENGINEERING↗
DOE OSTI2020

Decomposition Algorithms for Solving NP-hard Problems on a Quantum Annealer

NP-hard problems such as the maximum clique or minimum vertex cover problems, two of Karp’s 21 NP-hard problems, have several applications in computational chemistry, biochemistry and computer network security. Adiabatic quantum annealers can search for the optimum value of such NP-hard optimization problems, given the problem can be embedded on their hardware. However, this is often not possible due to certain limitations of the hardware connectivity structure of the annealer. This paper studies a general framework for a decomposition algorithm for NP-hard graph problems aiming to identify an optimal set of vertices. Our generic algorithm allows us to recursively divide an instance until the generated subproblems can be embedded on the quantum annealer hardware and subsequently solved. Furthermore, the framework is applied to the maximum clique and minimum vertex cover problems, and we propose several pruning and reduction techniques to speed up the recursive decomposition. The performance of both algorithms is assessed in a detailed simulation study.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Optimization of collision-free paths in a differential-drive robot by a smoothing piecewise-linear approach

This paper presents a methodology for optimizing pre-calculated collision-free paths of differential-drive wheeled robots. The main advantage of this methodology is that optimization is done by considering the kinematics and mechanical constraints of the mobile robot. In accordance to this proposal, the optimized path is achieved by applying recursively a local smoothing on an initial path which is originally modeled as a one-dimensional piecewise linear function. By this recursive smoothing, it can be ensured that the original piecewise linear function can be transformed into a smooth one that fulfill the constraints established by the kinematic equations of the wheeled mobile in terms of a minimum radius of curvature. As a result of this, a trajectory which guarantees lower power consumption and lower mechanical wear, is obtained. To show the better performance of the proposed approach, numerical simulation results are contrasted to those obtained from other reported methods with regards to path length, minimum radius of curvature, cross track error, continuity and resulting acceleration.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Avoidance of a Landau pole by flat contributions in QED

We consider massless Quantum Electrodynamics in the momentum scheme and carry forward an approach based on Dyson–Schwinger equations to approximate both the β-function and the renormalized photon self-energy (Yeats, 2011). Starting from the Callan–Symanzik equation, we derive a renormalization group (RG) recursion identity which implies a non-linear ODE for the anomalous dimension and extract a sufficient but not necessary criterion for the existence of a Landau pole. This criterion implies a necessary condition for QED to have no such pole. Solving the differential equation exactly for a toy model case, we integrate the corresponding RG equation for the running coupling and find that even though the β-function entails a Landau pole it exhibits a flat contribution capable of decreasing its growth, in other cases possibly to the extent that such a pole is avoided altogether. Finally, by applying the recursion identity, we compute the photon propagator and investigate the effect of flat contributions on both spacelike and timelike photons. -- Highlights: •We present an approach to approximate both the β-function and the photon self-energy. •We find a sufficient criterion for the self-energy to entail the existence of a Landau pole. •We study non-perturbative ‘flat’ contributions that emerge within the context of our approach. •We discuss a toy model and how it is affected by flat contributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Green’s functions on a renormalized lattice: An improved method for the integer quantum Hall transition

Highlights: • Speed up of numerical calculations by exact renormalization group steps. • Practical approach to handle delicate scaling analysis. • Absence of evidence to prefer marginal corrections over power-law corrections. We introduce a performance-optimized method to simulate localization problems on bipartite tight-binding lattices. It combines an exact renormalization group step to reduce the sparseness of the original problem with the recursive Green’s function method. We apply this framework to investigate the critical behavior of the integer quantum Hall transition of a tight-binding Hamiltonian defined on a simple square lattice. In addition, we employ an improved scaling analysis that includes two irrelevant exponents to characterize the shift of the critical energy as well as the corrections to the dimensionless Lyapunov exponent. We compare our findings with the results of a conventional implementation of the recursive Green’s function method, and we put them into broader perspective in view of recent development in this field.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A systematic comparison of machine learning methods for modeling of dynamic processes applied to combustion emission rate modeling

Ten established, data-driven dynamic algorithms are surveyed and a practical guide for understanding these methods generated. Existing Python programming packages for implementing each algorithm are acknowledged, and the model equations necessary for prediction are presented. A case study on a coal-fired power plant’s NO x emission rates is performed, directly comparing each modeling method’s performance on a mutual system. Each model is evaluated by its root mean squared error (RMSE) on out-of-sample future horizon predictions. Optimal hyperparameters are identified using either an exhaustive search or genetic algorithm. The top five model structures of each method are used to recursively predict future NO x emission rates over a 60-step time horizon. The RMSE at each future timestep is determined, and the recursive output prediction trends compared against measurements in time. The GRU neural network is identified as the best candidate for representing the system, demonstrating accurate and stable predictions across the future horizon by all considered models, while satisfactory performance was observed in several of the ARX/NARX formulations. Finally, these efforts have contributed 1) a concise resource of multiple proven dynamic machine learning methods, 2) a practical guide explaining the use of these methods, effectively lowering the “barrier-to-entry” of deploying such models in control systems, 3) a comparison study evaluating each method’s performance on a mutual system, 4) demonstration of accurate multi-timestep emissions modeling suitable for systems-level control, and 5) generalizable results demonstrating the suitability of each method for prediction over a multi-step future horizon to other complex dynamic systems.

42 ENGINEERING↗
DOE OSTI2020

Time-series machine-learning error models for approximate solutions to parameterized dynamical systems

This work proposes a machine-learning framework for modeling the error incurred by approximate solutions to parameterized dynamical systems. In particular, we extend the machine-learning error models (MLEM) framework proposed in Ref. Freno and Carlberg (2019) to dynamical systems. The proposed Time-Series Machine-Learning Error Modeling (T-MLEM) method constructs a regression model that maps features – which comprise error indicators that are derived from standard a posteriori error-quantification techniques – to a random variable for the approximate-solution error at each time instance. The proposed framework considers a wide range of candidate features, regression methods, and additive noise models. We consider primarily recursive regression techniques developed for time-series modeling, including both classical time-series models (e.g., autoregressive models) and recurrent neural networks (RNNs), but also analyze standard non-recursive regression techniques (e.g., feed-forward neural networks) for comparative purposes. Finally, numerical experiments conducted on multiple benchmark problems illustrate that the long short-term memory (LSTM) neural network, which is a type of RNN, outperforms other methods and yields substantial improvements in error predictions over traditional approaches.

42 ENGINEERING↗
DOE OSTI2023

Moving-horizon false data injection attack design against cyber–physical systems

Systematic attack design is essential to understanding the vulnerabilities of cyber–physical systems (CPSs), to better design for resiliency. In particular, false data injection attacks (FDIAs) are well-known and have been shown to be capable of bypassing bad data detection (BDD) while causing targeted biases in resulting state estimates. However, their effectiveness against moving horizon estimators (MHE) is not well understood. In fact, this paper shows that conventional FDIAs are generally ineffective against MHE. One of the main reasons is that the moving window renders the static FDIA recursively infeasible. Here, this paper proposes a new attack methodology, moving-horizon FDIA (MH-FDIA), by considering both the performance of historical attacks and the current system’s status. Theoretical guarantees for successful attack generation and recursive feasibility are given. Numerical simulations on the IEEE-14 bus system further validate the theoretical claims and show that the proposed MH-FDIA outperforms state-of-the-art counterparts in both stealthiness and effectiveness. In addition, an experiment on a path-tracking control system of an autonomous vehicle shows the feasibility of the MH-FDIA in real-world nonlinear systems.

42 ENGINEERING↗
DOE OSTI2011

Extraction of Hydrological Proximity Measures from DEMs using Parallel Processing

Land surface topography is one of the most important terrain properties which impact hydrological, geomorphological, and ecological processes active on a landscape. In our previous efforts to develop a soil depth model based upon topographic and land cover variables, we extracted a set of hydrological proximity measures (HPMs) from a Digital Elevation Model (DEM) as potential explanatory variables for soil depth. These HPMs may also have other, more general modeling applicability in hydrology, geomorphology and ecology, and so are described here from a general perspective. The HPMs we derived are variations of the distance up to ridge points (cells with no incoming flow) and variations of the distance down to stream points (cells with a contributing area greater than a threshold), following the flow path. These HPMs were computed using the D-infinity flow model that apportions flow between adjacent neighbors based on the direction of steepest downward slope on the eight triangular facets constructed in a 3 x 3 grid cell window using the center cell and each pair of adjacent neighboring grid cells in turn. The D-infinity model typically results in multiple flow paths between 2 points on the topography, with the result that distances may be computed as the minimum, maximum or average of the individual flow paths. In addition, each of the HPMs, are calculated vertically, horizontally, and along the land surface. Previously, these HPMs were calculated using recursive serial algorithms which suffered from stack overflow problems when used to process large datasets, limiting the size of DEMs that could be analyzed using that method to approximately 7000 x 7000 cells. To overcome this limitation, we developed a message passing interface (MPI) parallel approach for calculating these HPMs. The parallel algorithms of the HPMs spatially partition the input grid into stripes which are each assigned to separate processes for computation. Each of those processes then uses a queue data structure to order the processing of cells so that each cell is visited only once and the cross-process communications that are a standard part of MPI are handled in an efficient manner. This parallel approach allows analysis of much larger DEMs as compared to the serial recursive algorithms. In this paper, we present the definitions of the HPMs, the serial and parallel algorithms used in their extraction and their potential applications in hydrology, geomorphology and ecology.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Fast approximation to supershell partition functions: Explicit forms of the coefficients

In a previous work, we derived a formula for supershell partition functions, which are the cornerstone of the Super-Transition-Array approach to radiative-opacity calculations. The new expression 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 (denoted Γ $k$ ) for which we obtained a recursion relation. In the present short paper, as a complement of the previous one, we give an explicit formula for the coefficients Γ $k$ . Another recursion relation is also provided, as well as an alternative expression involving Bell polynomials. In conclusion, the connections with cycle indexes of permutation groups and with elementary symmetric polynomials are outlined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Active Suspension Parameters Identification: An Algebraic Approach and its Application to Suspension Travel Control

The parameters of an active suspension need to be identified online, such that the suspension control system can be adapted to mechanical wear and load change. Recursive least squares and observer-based methods are frequently utilized to fulfill this purpose. However, they can yield slow parameter identification due to their asymptotic nature. We propose an algebraic identifier to estimate the parameters of an active suspension online, which does not maintain an asymptotic convergence phase. Simulation results demonstrate the effectiveness of the proposed algebraic approach.

Wang, Zejiang↗
DOE OSTI2012

Impact of Triple-Negative Phenotype on Prognosis of Patients With Breast Cancer Brain Metastases

Purpose: To elucidate survival times and identify potential prognostic factors in patients with triple-negative (TN) phenotype who harbored brain metastases arising from breast cancer and who underwent stereotactic radiosurgery (SRS). Methods and Materials: A total of 103 breast cancer patients with brain metastases were treated with SRS and then studied retrospectively. Twenty-four patients (23.3%) were TN. Survival times were estimated using the Kaplan-Meier method, with a log-rank test computing the survival time difference between groups. Univariate and multivariate analyses to predict potential prognostic factors were performed using a Cox proportional hazard regression model. Results: The presence of TN phenotype was associated with worse survival times, including overall survival after the diagnosis of primary breast cancer (43 months vs. 82 months), neurologic survival after the diagnosis of intracranial metastases, and radiosurgical survival after SRS, with median survival times being 13 months vs. 25 months and 6 months vs. 16 months, respectively (p < 0.002 in all three comparisons). On multivariate analysis, radiosurgical survival benefit was associated with non-TN status and lower recursive partitioning analysis class at the initial SRS. Conclusion: The TN phenotype represents a significant adverse prognostic factor with respect to overall survival, neurologic survival, and radiosurgical survival in breast cancer patients with intracranial metastasis. Recursive partitioning analysis class also served as an important and independent prognostic factor.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2012

A Phase I Study of Short-Course Accelerated Whole Brain Radiation Therapy for Multiple Brain Metastases

Purpose: To define the maximum tolerated dose (MTD) of a SHort-course Accelerated whole brain RadiatiON therapy (SHARON) in the treatment of patients with multiple brain metastases. Methods and Materials: A phase 1 trial in 4 dose-escalation steps was designed: 12 Gy (3 Gy per fraction), 14 Gy (3.5 Gy per fraction), 16 Gy (4 Gy per fraction), and 18 Gy (4.5 Gy per fraction). Eligibility criteria included patients with unfavorable recursive partitioning analysis (RPA) class > or =2 with at least 3 brain metastases or metastatic disease in more than 3 organ systems, and Eastern Cooperative Oncology Group (ECOG) performance status {<=}3. Treatment was delivered in 2 days with twice-daily fractionation. Patients were treated in cohorts of 6-12 to define the MTD. The dose-limiting toxicity (DLT) was defined as any acute toxicity {>=}grade 3, according to the Radiation Therapy Oncology Group scale. Information on the status of the main neurologic symptoms and quality of life were recorded. Results: Characteristics of the 49 enrolled patients were as follows: male/female, 30/19; median age, 66 years (range, 23-83 years). ECOG performance status was <3 in 46 patients (94%). Fourteen patients (29%) were considered to be in recursive partitioning analysis (RPA) class 3. Grade 1-2 acute neurologic (26.4%) and skin (18.3%) toxicities were recorded. Only 1 patient experienced DLT (neurologic grade 3 acute toxicity). With a median follow-up time of 5 months (range, 1-23 months), no late toxicities have been observed. Three weeks after treatment, 16 of 21 symptomatic patients showed an improvement or resolution of presenting symptoms (overall symptom response rate, 76.2%; confidence interval 0.95: 60.3-95.9%). Conclusions: Short-course accelerated radiation therapy in twice-daily fractions for 2 consecutive days is tolerated up to a total dose of 18 Gy. A phase 2 study has been planned to evaluate the efficacy on overall survival, symptom control, and quality of life indices.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Is Androgen Deprivation Therapy Necessary in All Intermediate-Risk Prostate Cancer Patients Treated in the Dose Escalation Era?

Purpose: The benefit of adding androgen deprivation therapy (ADT) to dose-escalated radiation therapy (RT) for men with intermediate-risk prostate cancer is unclear; therefore, we assessed the impact of adding ADT to dose-escalated RT on freedom from failure (FFF). Methods: Three groups of men treated with intensity modulated RT or 3-dimensional conformal RT (75.6-78 Gy) from 1993-2008 for prostate cancer were categorized as (1) 326 intermediate-risk patients treated with RT alone, (2) 218 intermediate-risk patients treated with RT and ≤6 months of ADT, and (3) 274 low-risk patients treated with definitive RT. Median follow-up was 58 months. Recursive partitioning analysis based on FFF using Gleason score (GS), T stage, and pretreatment PSA concentration was applied to the intermediate-risk patients treated with RT alone. The Kaplan-Meier method was used to estimate 5-year FFF. Results: Based on recursive partitioning analysis, intermediate-risk patients treated with RT alone were divided into 3 prognostic groups: (1) 188 favorable patients: GS 6, ≤T2b or GS 3+4, ≤T1c; (2) 71 marginal patients: GS 3+4, T2a-b; and (3) 68 unfavorable patients: GS 4+3 or T2c disease. Hazard ratios (HR) for recurrence in each group were 1.0, 2.1, and 4.6, respectively. When intermediate-risk patients treated with RT alone were compared to intermediate-risk patients treated with RT and ADT, the greatest benefit from ADT was seen for the unfavorable intermediate-risk patients (FFF, 74% vs 94%, respectively; P=.005). Favorable intermediate-risk patients had no significant benefit from the addition of ADT to RT (FFF, 94% vs 95%, respectively; P=.85), and FFF for favorable intermediate-risk patients treated with RT alone approached that of low-risk patients treated with RT alone (98%). Conclusions: Patients with favorable intermediate-risk prostate cancer did not benefit from the addition of ADT to dose-escalated RT, and their FFF was nearly as good as patients with low-risk disease. In patients with GS 4+3 or T2c disease, the addition of ADT to dose-escalated RT did improve FFF.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Dosimetric Predictors of Duodenal Toxicity After Intensity Modulated Radiation Therapy for Treatment of the Para-aortic Nodes in Gynecologic Cancer

Purpose: To determine the incidence of duodenal toxicity in patients receiving intensity modulated radiation therapy (IMRT) for treatment of para-aortic nodes and to identify dosimetric parameters predictive of late duodenal toxicity. Methods and Materials: We identified 105 eligible patients with gynecologic malignancies who were treated with IMRT for gross metastatic disease in the para-aortic nodes from January 1, 2005, through December 31, 2009. Patients were treated to a nodal clinical target volume to 45 to 50.4 Gy with a boost to 60 to 66 Gy. The duodenum was contoured, and dosimetric data were exported for analysis. Duodenal toxicity was scored according to Radiation Therapy Oncology Group criteria. Univariate Cox proportional hazards analysis and recursive partitioning analysis were used to determine associations between dosimetric variables and time to toxicity and to identify the optimal threshold that separated patients according to risk of toxicity. Results: Nine of the 105 patients experienced grade 2 to grade 5 duodenal toxicity, confirmed by endoscopy in all cases. The 3-year actuarial rate of any duodenal toxicity was 11.7%. A larger volume of the duodenum receiving 55 Gy (V55) was associated with higher rates of duodenal toxicity. The 3-year actuarial rates of duodenal toxicity with V55 above and below 15 cm{sup 3} were 48.6% and 7.4%, respectively (P<.01). In Cox univariate analysis of dosimetric variables, V55 was associated with duodenal toxicity (P=.029). In recursive partitioning analysis, V55 less than 13.94% segregated all patients with duodenal toxicity. Conclusions: Dose-escalated IMRT can safely and effectively treat para-aortic nodal disease in gynecologic malignancies, provided that care is taken to limit the dose to the duodenum to reduce the risk of late duodenal toxicity. Limiting V55 to below 15 cm{sup 3} may reduce the risk of duodenal complications. In cases where the treatment cannot be delivered within these constraints, consideration should be given to other treatment approaches such as resection or initial chemotherapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

Predictive Parameters of Symptomatic Hematochezia Following 5-Fraction Gantry-Based SABR in Prostate Cancer

Purpose: This study identified predictors of high-grade late hematochezia (HH) following 5-fraction gantry-based stereotactic ablative radiation therapy (SABR). Methods and Materials: Hematochezia data for 258 patients who received 35 to 40 Gy SABR in 5-fractions as part of sequential phase 2 prospective trials was retrieved. Grade 2 or higher late rectal bleeding was labeled HH. Hematochezia needing steroid suppositories, 4% formalin, or 1 to 2 sessions of argon plasma coagulation (APC) was labeled grade 2. More than 2 sessions of APC, blood transfusion, or a course of hyperbaric oxygen was grade 3 and development of visceral fistula, grade 4. Various dosimetric and clinical factors were analyzed using univariate and multivariate analyses. Receiver operating characteristic (ROC) curve analysis and recursive partitioning analysis were used to determine clinically valid cut-off points and identify risk groups, respectively. Results: HH was observed in 19.4%, grade ≥3 toxicity in 3.1%. Median follow-up was 29.7 months (interquartile range [IQR]: 20.6-61.7) Median time to develop HH was 11.7 months (IQR: 9.0-15.2) from the start of radiation. At 2 years, cumulative HH was 4.9%, 27.2%, and 42.1% in patients who received 35 Gy to prostate (4-mm planning target volume [PTV] margin), 40 Gy to prostate (5-mm PTV margin), and 40 Gy to prostate/seminal vesicles (5-mm PTV margin), respectively (P<.0001). In the ROC analysis, volume of rectum receiving radiation dose of 38 Gy (V38) was a strong predictor of HH with an area under the curve of 0.65. In multivariate analysis, rectal V38 (≥2.0 cm{sup 3}; odds ratio [OR]: 4.7); use of anticoagulants in the follow-up period (OR: 6.5) and presence of hemorrhoids (OR: 2.7) were the strongest predictors. Recursive partitioning analysis showed rectal V38 < 2.0 cm{sup 3}, and use of anticoagulants or rectal V38 ≥ 2.0 cm{sup 3} plus 1 other risk factor resulted in an HH risk of >30%. Conclusions: Rectal V38 and 2 clinical factors were strong predictors of HH following 5-fraction SABR. Planning constraints should keep rectal V38 below 2.0 cm{sup 3}.

62 RADIOLOGY AND NUCLEAR MEDICINE↗