Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Regularization”

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 181 records · Page 10

srMO-BO-3GP: A sequential regularized multi-objective Bayesian optimization for constrained design applications using an uncertain Pareto classifier

Bayesian optimization (BO) is an efficient and flexible global optimization framework that is applicable to a very wide range of engineering applications. To leverage the capability of the classical BO, many extensions, including multi-objective, multi-fidelity, parallelization, and latent-variable modeling, have been proposed to address the limitations of the classical BO framework. In this work, we propose a novel multi-objective BO formalism, called srMO-BO-3GP, to solve multi-objective optimization problems in a sequential setting. Three different Gaussian processes (GPs) are stacked together, where each of the GPs is assigned with a different task. The first GP is used to approximate a single-objective computed from the multi-objective definition, the second GP is used to learn the unknown constraints, and the third one is used to learn the uncertain Pareto frontier. At each iteration, a multi-objective augmented Tchebycheff function is adopted to convert multi-objective to single-objective, where the regularization with a regularized ridge term is also introduced to smooth the single-objective function. Finally, we couple the third GP along with the classical BO framework to explore the convergence and diversity of the Pareto frontier by the acquisition function for exploitation and exploration. The proposed framework is demonstrated using several numerical benchmark functions, as well as a thermomechanical finite element model for flip-chip package design optimization.

42 ENGINEERING↗

WE-AB-204-01: Performance Characterization of Regularized-Reconstruction Algorithm for 90Y PET/CT Images

Purpose: {sup 90}Y PET/CT imaging and quantification have recently been suggested as an approach of treatment verification. However, due to low positron yield (32ppm), the {sup 90}Y-PET/CT images are very noisy. Iterative reconstruction techniques that employ regularization, e.g. block sequential regularized expectation maximization (BSREM) algorithm (recently implemented on GE scanners – QClear™), has the potential to increase quantitative accuracy with lower noise penalty compared to OSEM. Our aim is to investigate the performance of RR algorithms in {sup 90}Y PET/CT studies. Methods: A NEMA IEC phantom filled with 3GBq {sup 90}YCl{sub 2} (to simulate patient treatment) was imaged on GE-D690 for 1800s/bed. The sphere-to-background ratio of 7. The data were reconstructed using OSEM and BSREM with PSF modeling and TOF correction while varying the iterations (IT) from 1–6 with fixed 24subsets. For BSREM, the edge-preservation parameter (γ ) was 2 and the penalty-parameters (β) was varied 350–950. In all cases a post-reconstruction filter of 5.2mm (2pixel) transaxial and standard z-axis were used. Sphere average activity concentration (AC) and background standard deviation (SD) were then calculated from VOIs drawn in the spheres and background. Results: Increasing IT from 1to6, the %SD in OSEM increased from 30% to 80%, whereas %SD in BSREM images increased by <5% for all βs. BSREM with β=350 didn’t offer any improvement over OSEM (convergence of mean achieved at 2 IT, in this study). Increasing β from 350 to 950 reduced the AC accuracy of small spheres (<20mm) by 10% and noise from 40% to 20%, which resulted in CNR increase from 11 to 17. Conclusion: In count-limited studies such as {sup 90}Y PET/CT, BSREM can be used to suppress image noise and increase CNR at the expense of a relatively small decrease of quantitative accuracy. The BSREM parameters need to be optimized for each study depending on the radionuclides and count densities. Research reported in this publication was supported by the National Cancer Institute of the National Institutes of Health under Award Number R01CA13898 and was also partly supported by General Electric.

60 APPLIED LIFE SCIENCES↗

Temporal regularization of ultrasound-based liver motion estimation for image-guided radiation therapy

Purpose: Ultrasound-based motion estimation is an expanding subfield of image-guided radiation therapy. Although ultrasound can detect tissue motion that is a fraction of a millimeter, its accuracy is variable. For controlling linear accelerator tracking and gating, ultrasound motion estimates must remain highly accurate throughout the imaging sequence. This study presents a temporal regularization method for correlation-based template matching which aims to improve the accuracy of motion estimates. Methods: Liver ultrasound sequences (15–23 Hz imaging rate, 2.5–5.5 min length) from ten healthy volunteers under free breathing were used. Anatomical features (blood vessels) in each sequence were manually annotated for comparison with normalized cross-correlation based template matching. Five sequences from a Siemens Acuson™ scanner were used for algorithm development (training set). Results from incremental tracking (IT) were compared with a temporal regularization method, which included a highly specific similarity metric and state observer, known as the α–β filter/similarity threshold (ABST). A further five sequences from an Elekta Clarity™ system were used for validation, without alteration of the tracking algorithm (validation set). Results: Overall, the ABST method produced marked improvements in vessel tracking accuracy. For the training set, the mean and 95th percentile (95%) errors (defined as the difference from manual annotations) were 1.6 and 1.4 mm, respectively (compared to 6.2 and 9.1 mm, respectively, for IT). For each sequence, the use of the state observer leads to improvement in the 95% error. For the validation set, the mean and 95% errors for the ABST method were 0.8 and 1.5 mm, respectively. Conclusions: Ultrasound-based motion estimation has potential to monitor liver translation over long time periods with high accuracy. Nonrigid motion (strain) and the quality of the ultrasound data are likely to have an impact on tracking performance. A future study will investigate spatial uniformity of motion and its effect on the motion estimation errors.

61 RADIATION PROTECTION AND DOSIMETRY↗

Accelerated fast iterative shrinkage thresholding algorithms for sparsity-regularized cone-beam CT image reconstruction

Purpose: The development of iterative image reconstruction algorithms for cone-beam computed tomography (CBCT) remains an active and important research area. Even with hardware acceleration, the overwhelming majority of the available 3D iterative algorithms that implement nonsmooth regularizers remain computationally burdensome and have not been translated for routine use in time-sensitive applications such as image-guided radiation therapy (IGRT). In this work, two variants of the fast iterative shrinkage thresholding algorithm (FISTA) are proposed and investigated for accelerated iterative image reconstruction in CBCT. Methods: Algorithm acceleration was achieved by replacing the original gradient-descent step in the FISTAs by a subproblem that is solved by use of the ordered subset simultaneous algebraic reconstruction technique (OS-SART). Due to the preconditioning matrix adopted in the OS-SART method, two new weighted proximal problems were introduced and corresponding fast gradient projection-type algorithms were developed for solving them. We also provided efficient numerical implementations of the proposed algorithms that exploit the massive data parallelism of multiple graphics processing units. Results: The improved rates of convergence of the proposed algorithms were quantified in computer-simulation studies and by use of clinical projection data corresponding to an IGRT study. The accelerated FISTAs were shown to possess dramatically improved convergence properties as compared to the standard FISTAs. For example, the number of iterations to achieve a specified reconstruction error could be reduced by an order of magnitude. Volumetric images reconstructed from clinical data were produced in under 4 min. Conclusions: The FISTA achieves a quadratic convergence rate and can therefore potentially reduce the number of iterations required to produce an image of a specified image quality as compared to first-order methods. We have proposed and investigated accelerated FISTAs for use with two nonsmooth penalty functions that will lead to further reductions in image reconstruction times while preserving image quality. Moreover, with the help of a mixed sparsity-regularization, better preservation of soft-tissue structures can be potentially obtained. The algorithms were systematically evaluated by use of computer-simulated and clinical data sets.

60 APPLIED LIFE SCIENCES↗

Use of regularized principal component analysis to model anatomical changes during head and neck radiation therapy for treatment adaptation and response assessment

Purpose: To develop standard (SPCA) and regularized (RPCA) principal component analysis models of anatomical changes from daily cone beam CTs (CBCTs) of head and neck (H&N) patients and assess their potential use in adaptive radiation therapy, and for extracting quantitative information for treatment response assessment. Methods: Planning CT images of ten H&N patients were artificially deformed to create “digital phantom” images, which modeled systematic anatomical changes during radiation therapy. Artificial deformations closely mirrored patients’ actual deformations and were interpolated to generate 35 synthetic CBCTs, representing evolving anatomy over 35 fractions. Deformation vector fields (DVFs) were acquired between pCT and synthetic CBCTs (i.e., digital phantoms) and between pCT and clinical CBCTs. Patient-specific SPCA and RPCA models were built from these synthetic and clinical DVF sets. EigenDVFs (EDVFs) having the largest eigenvalues were hypothesized to capture the major anatomical deformations during treatment. Results: Principal component analysis (PCA) models achieve variable results, depending on the size and location of anatomical change. Random changes prevent or degrade PCA’s ability to detect underlying systematic change. RPCA is able to detect smaller systematic changes against the background of random fraction-to-fraction changes and is therefore more successful than SPCA at capturing systematic changes early in treatment. SPCA models were less successful at modeling systematic changes in clinical patient images, which contain a wider range of random motion than synthetic CBCTs, while the regularized approach was able to extract major modes of motion. Conclusions: Leading EDVFs from the both PCA approaches have the potential to capture systematic anatomical change during H&N radiotherapy when systematic changes are large enough with respect to random fraction-to-fraction changes. In all cases the RPCA approach appears to be more reliable at capturing systematic changes, enabling dosimetric consequences to be projected once trends are established early in a treatment course, or based on population models.

60 APPLIED LIFE SCIENCES↗

Charged fermions tunneling from regular black holes

We study Hawking radiation of charged fermions as a tunneling process from charged regular black holes, i.e., the Bardeen and ABGB black holes. For this purpose, we apply the semiclassical WKB approximation to the general covariant Dirac equation for charged particles and evaluate the tunneling probabilities. We recover the Hawking temperature corresponding to these charged regular black holes. Further, we consider the back-reaction effects of the emitted spin particles from black holes and calculate their corresponding quantum corrections to the radiation spectrum. We find that this radiation spectrum is not purely thermal due to the energy and charge conservation but has some corrections. In the absence of charge, e = 0, our results are consistent with those already present in the literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Electromagnetic wave propagation with negative phase velocity in regular black holes

We discuss the propagation of electromagnetic plane waves with negative phase velocity in regular black holes. For this purpose, we consider the Bardeen model as a nonlinear magnetic monopole and the Bardeen model coupled to nonlinear electrodynamics with a cosmological constant. It turns out that the region outside the event horizon of each regular black hole does not support negative phase velocity propagation, while its possibility in the region inside the event horizon is discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Propagation of excitation in long 1D chains: Transition from regular quantum dynamics to stochastic dynamics

The quantum dynamics problem for a 1D chain consisting of 2N + 1 sites (N Much-Greater-Than 1) with the interaction of nearest neighbors and an impurity site at the middle differing in energy and in coupling constant from the sites of the remaining chain is solved analytically. The initial excitation of the impurity is accompanied by the propagation of excitation over the chain sites and with the emergence of Loschmidt echo (partial restoration of the impurity site population) in the recurrence cycles with a period proportional to N. The echo consists of the main (most intense) component modulated by damped oscillations. The intensity of oscillations increases with increasing cycle number and matrix element C of the interaction of the impurity site n = 0 with sites n = {+-}1 (0 < C {<=} 1; for the remaining neighboring sites, the matrix element is equal to unity). Mixing of the components of echo from neighboring cycles induces a transition from the regular to stochastic evolution. In the regular evolution region, the wave packet propagates over the chain at a nearly constant group velocity, embracing a number of sites varying periodically with time. In the stochastic regime, the excitation is distributed over a number of sites close to 2N, with the populations varying irregularly with time. The model explains qualitatively the experimental data on ballistic propagation of the vibrational energy in linear chains of CH{sub 2} fragments and predicts the possibility of a nondissipative energy transfer between reaction centers associated with such chains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Electrostatic self-force in the field of an (n + 1)-dimensional black hole: Dimensional regularization

The self-energy of a classical charged particle localized at a relatively large distance outside the event horizon of an (n + 1)-dimensional Schwarzschild–Tangherlini black hole for an arbitrary n ≥ 3 is calculated. An expression for the electrostatic Green function is derived in the first two orders of the perturbation theory. Dimensional regularization is proposed to be used to regularize the corresponding formally divergent expression for the self-energy. The derived expression for the renormalized self-energy is compared with the results of other authors.

79 ASTRONOMY AND ASTROPHYSICS↗

NSVZ Relation in Supersymmetric Theories Regularized by Higher Derivatives

We investigate the NSVZ equation in N=1 supersymmetric gauge theories, which relates the β ‑function to the anomalous dimension of the matter superfields. In particular, we argue that it is closely connected with the non-renormalization theorem for the three-point gauge-ghost vertices (in which one external leg corresponds to the quantum gauge superfield). Using finiteness of these vertices the exact NSVZ β -function can be equivalently presented in the form of a relation between the β -function and the anomalous dimensions of the quantum gauge superfield, of the Faddeev–Popov ghosts, and of the matter superfields. This equation allows explaining how the NSVZ equation appears in the perturbation theory for theories regularized by higher covariant derivatives and constructing the NSVZ scheme in all orders in the case of using this regularization. The results are verified by an explicit three-loop calculation of the terms quartic in the Yukawa couplings.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Posterior Regularized Bayesian Neural Network

Traditional NNs often lack the ability for uncertainty quantification. Bayesian NNs(BNNs) could help measure the confidence level by using distributions in NNs modeling. Besides, knowledge is commonly available and could improve the performance of BNNs if it can be properly incorporated. In this work, we propose a novel Posterior-Regularized BNN(PR-BNN) model by incorporating soft and hard constraints as a posterior regularization term. We also propose an augmented Lagrangian method and stochastic optimization algorithm for efficient updating via Monte Carlo sampling. The simulations and case studies for solar PV plants have shown the performance improvement of the proposed model over traditional BNNs.

97 MATHEMATICS AND COMPUTING↗

Georgia Tech Accelerated, Compressed, and Regularized Compute of Kinetic-based PDEs (Final Report)

This report summarizes the collaborative effort between Lawrence Livermore National Laboratory and Georgia Tech to enhance the BoBa library for tensor train computation in PDE solvers, with a target on kinetic equations and their continuum limits. We aimed to reduce computational cost and memory usage by replacing traditional array-based computations with tensor trains. We examined the compressibility of time-evolving solutions to the Euler equations with discontinuities. We also explored using the first invsicid and linear regularization of the compressible flow equations via the information geometric regularization (IGR). We explored this in a tensor train formulation. To identify that inverse terms in the IGR equations pose problems for tensor train formulations and investigate efficient methods for batched inversion of tensor trains.

97 MATHEMATICS AND COMPUTING↗

Regularizing Invertible Neural Networks for Airfoil Design Through Dimension Reduction

This work examines how dimension reduction can improve the performance of invertible neural networks (INN) for airfoil design. Design workflows are typically expensive, relying on many evaluations of high fidelity computational fluid dynamics (CFD) models. Furthermore, the inverse design problem is typically ill-posed. That is, multiple valid solutions exist that satisfy the design criteria. Regularization can reduce this inverse design space and simplify the problem. We study the use of subspace-based input dimension reduction to act as a regularizer for the INN model and improve the recovery of new airfoil shapes with desired performance characteristics. We find that the dimension reduction identifies two dominant modes, relating to airfoil thickness and camber, that optimally determine the airfoil's aerodynamics. We demonstrate the capability of the proposed INN model to generate 100 airfoils that satisfy the specific aerodynamic and structural characteristics.

airfoil design↗

SYSTEMATIC REGULARITY OF HEMISPHERIC SUNSPOT AREAS OVER THE PAST 140 YEARS

Solar magnetic activity varies with time in the two hemispheres in different ways. The hemispheric interconnection of solar activity phenomena provides an important clue to understanding the dynamical behavior of solar dynamo actions. In this paper, several analysis approaches are proposed to analyze the systematic regularity of hemispheric asynchronism and amplitude asymmetry of long-term sunspot areas during solar cycles 9–24. It is found that, (1) both the hemispheric asynchronism and the amplitude asymmetry of sunspot areas are prevalent behaviors and are not anomalous, but the hemispheric asynchronism exhibits a much more regular behavior than the amplitude asymmetry; (2) the phase-leading hemisphere returns back to the identical hemisphere every 8 solar cycles, and the secular periodic pattern of hemispheric phase differences follows 3 (south leading) + 5 (north leading) solar cycles, which probably corresponds to the Gleissberg cycle; and (3) the pronounced periodicities of (absolute and normalized) asymmetry indices and lines of synchronization (LOSs) are not identical: the significant periodic oscillations are 80.65 ± 6.31, 20.91 ± 0.40, and 13.45 ± 0.16 years for the LOS values, and 51.34 ± 2.48, 8.83/8.69 ± 0.07, and 3.77 ± 0.02 years for the (absolute and normalized) asymmetry indices. The analysis results improve our knowledge on the hemispheric interrelation of solar magnetic activity and may provide valuable constraints for solar dynamo models.

79 ASTRONOMY AND ASTROPHYSICS↗

Unexpected series of regular frequency spacing of δ Scuti stars in the non-asymptotic regime. II. Sample-Echelle diagrams and rotation

A sequence search method was developed for searching for regular frequency spacing in δ Scuti stars by visual inspection (VI) and algorithmic search. The sample contains 90 δ Scuti stars observed by CoRoT. An example is given to represent the VI. The algorithm (SSA) is described in detail. The data treatment of the CoRoT light curves, the criteria for frequency filtering, and the spacings derived by two methods (i.e., three approaches: VI, SSA, and FT) are given for each target. Echelle diagrams are presented for 77 targets for which at least one sequence of regular spacing was identified. Comparing the spacing and the shifts between pairs of echelle ridges revealed that at least one pair of echelle ridges is shifted to midway between the spacing for 22 stars. The estimated rotational frequencies compared to the shifts revealed rotationally split doublets, triplets, and multiplets not only for single frequencies, but for the complete echelle ridges in 31 δ Scuti stars. Furthermore, using several possible assumptions for the origin of the spacings, we derived the large separation (${\rm{\Delta }}\nu $) that are distributed along the mean density versus large separations relation derived from stellar models.

79 ASTRONOMY AND ASTROPHYSICS↗

Fuel Stratification Effects on Gasoline Compression Ignition with a Regular-Grade Gasoline on a Single-Cylinder Medium-Duty Diesel Engine at Low Load

Prior research studies have investigated a wide variety of gasoline compression ignition (GCI) injection strategies and the resulting fuel stratification levels to maintain control over the combustion phasing, duration, and heat release rate. Previous GCI research at the US Department of Energy’s Oak Ridge National Laboratory has shown that for a combustion mode with a low degree of fuel stratification, called “partial fuel stratification” (PFS), gasoline range fuels with anti-knock index values in the range of regular-grade gasoline (~87 anti-knock index or higher) provides very little controllability over the timing of combustion without significant boost pressures. On the contrary, heavy fuel stratification (HFS) provides control over combustion phasing but has challenges achieving low temperature combustion operation, which has the benefits of low NOX and soot emissions, because of the air handling burdens associated with the required high exhaust gas recirculation rates. Furthermore, this work investigates HFS and PFS combustion, efficiency, and emissions performance on a single-cylinder, medium-duty engine with a regular-grade gasoline (91 research octane number) at 1,200 rpm, 4.3 bar, and 3.0 nominal gross indicated mean effective pressure operating points with boost levels similar to those in a medium-duty diesel application. Authority of combustion phasing with main injection timing sweeps for HFS and second injection timing sweeps and fuel split sweeps for PFS are shown. In addition, this work is discussed in the context of previous findings with a light-duty diesel platform, and next steps and future direction for this work are presented.

33 ADVANCED PROPULSION SYSTEMS↗

Deep Neural Networks are Adaptive to Function Regularity and Data Distribution in Approximation and Estimation

Deep learning has exhibited remarkable results across diverse areas. To understand its success, substantial research has been directed towards its theoretical foundations. Nev- ertheless, the majority of these studies examine how well deep neural networks can model functions with uniform regularities. In this paper, we explore a different angle: how deep neural networks can adapt to varying degrees of smoothness in functions and nonuni- form data distributions across different locations and scales. More precisely, we focus on a broad class of functions defined by nonlinear tree-based approximation methods. This class encompasses a range of function types, such as functions with uniform regularities and discontinuous functions. We develop nonparametric approximation and estimation theories for this class using deep ReLU networks. Our results show that deep neural networks are adaptive to the nonuniform smoothness of functions and nonuniform data distributions at different locations and scales. We apply our results to several function classes, and derive the corresponding approximation and generalization errors. The validity of our results is demonstrated through numerical experiments.

97 MATHEMATICS AND COMPUTING↗

Physically regularized machine learning emulators of aerosol activation

Abstract. The activation of aerosol into cloud droplets is an important step in the formation of clouds and strongly influences the radiative budget of the Earth. Explicitly simulating aerosol activation in Earth system models is challenging due to the computational complexity required to resolve the necessary chemical and physical processes and their interactions. As such, various parameterizations have been developed to approximate these details at reduced computational cost and accuracy. Here, we explore how machine learning emulators can be used to bridge this gap in computational cost and parameterization accuracy. We evaluate a set of emulators of a detailed cloud parcel model using physically regularized machine learning regression techniques. We find that the emulators can reproduce the parcel model at higher accuracy than many existing parameterizations. Furthermore, physical regularization tends to improve emulator accuracy, most significantly when emulating very low activation fractions. This work demonstrates the value of physical constraints in machine learning model development and enables the implementation of improved hybrid physical and machine learning models of aerosol activation into next-generation Earth system models.

58 GEOSCIENCES↗