Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Newton optimization”

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 271 records · Page 15

Design and Analysis of an X-Ray Mirror Assembly Using the Meta-Shell Approach

Lightweight and high resolution optics are needed for future space-based x-ray telescopes to achieve advances in high-energy astrophysics. Past missions such as Chandra and XMM-Newton have achieved excellent angular resolution using a full shell mirror approach. Other missions such as Suzaku and NuSTAR have achieved lightweight mirrors using a segmented approach. This paper describes a new approach, called meta-shells, which combines the fabrication advantages of segmented optics with the alignment advantages of full shell optics. Meta-shells are built by layering overlapping mirror segments onto a central structural shell. The resulting optic has the stiffness and rotational symmetry of a full shell, but with an order of magnitude greater collecting area. Several meta-shells so constructed can be integrated into a large x-ray mirror assembly by proven methods used for Chandra and XMM-Newton. The mirror segments are mounted to the meta-shell using a novel four point semi-kinematic mount. The four point mount deterministically locates the segment in its most performance sensitive degrees of freedom. Extensive analysis has been performed to demonstrate the feasibility of the four point mount and meta-shell approach. A mathematical model of a meta-shell constructed with mirror segments bonded at four points and subject to launch loads has been developed to determine the optimal design parameters, namely bond size, mirror segment span, and number of layers per meta-shell. The parameters of an example 1.3 m diameter mirror assembly are given including the predicted effective area. To verify the mathematical model and support opto-mechanical analysis, a detailed finite element model of a meta-shell was created. Finite element analysis predicts low gravity distortion and low sensitivity to thermal gradients.

module↗

Meta-shell Approach for Constructing Lightweight and High Resolution X-Ray Optics

Lightweight and high resolution optics are needed for future space-based x-ray telescopes to achieve advances in high-energy astrophysics. Past missions such as Chandra and XMM-Newton have achieved excellent angular resolution using a full shell mirror approach. Other missions such as Suzaku and NuSTAR have achieved lightweight mirrors using a segmented approach. This paper describes a new approach, called meta-shells, which combines the fabrication advantages of segmented optics with the alignment advantages of full shell optics. Meta-shells are built by layering overlapping mirror segments onto a central structural shell. The resulting optic has the stiffness and rotational symmetry of a full shell, but with an order of magnitude greater collecting area. Several meta-shells so constructed can be integrated into a large x-ray mirror assembly by proven methods used for Chandra and XMM-Newton. The mirror segments are mounted to the meta-shell using a novel four point semi-kinematic mount. The four point mount deterministically locates the segment in its most performance sensitive degrees of freedom. Extensive analysis has been performed to demonstrate the feasibility of the four point mount and meta-shell approach. A mathematical model of a meta-shell constructed with mirror segments bonded at four points and subject to launch loads has been developed to determine the optimal design parameters, namely bond size, mirror segment span, and number of layers per meta-shell. The parameters of an example 1.3 m diameter mirror assembly are given including the predicted effective area. To verify the mathematical model and support opto-mechanical analysis, a detailed finite element model of a meta-shell was created. Finite element analysis predicts low gravity distortion and low thermal distortion. Recent results are discussed including Structural Thermal Optical Performance (STOP) analysis as well as vibration and shock testing of prototype meta-shells.

Lightweight and high resolution optics↗

A self-consistent field formulation of excited state mean field theory

In this work, we show that, as in Hartree-Fock theory, the orbitals for excited state mean field theory can be optimized via a self-consistent one-electron equation in which electron-electron repulsion is accounted for through mean field operators. In addition to showing that this excited state ansatz is sufficiently close to a mean field product state to admit a one-electron formulation, this approach brings the orbital optimization speed to within roughly a factor of two of ground state mean field theory. The approach parallels Hartree Fock theory in multiple ways, including the presence of a commutator condition, a one-electron mean-field working equation, and acceleration via direct inversion in the iterative subspace. When combined with a configuration interaction singles Davidson solver for the excitation coefficients, the self-consistent field formulation dramatically reduces the cost of the theory compared to previous approaches based on quasi-Newton descent.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spectroscopic Determination of the Physical Conditions in Hot Optically Thin Sources

The Astrophysics Plasma Emission Code and Database (APEC/APED), developed in part under this grant, have been upgraded to Version 1.3.0 and are now beginning to find widespread applications to X-ray spectral data from Chandra and XMM-Newton. These models represent the best theoretical data currently available and are reasonably complete below about 30 A. We are making significant progress in improving the spectral models between 30 and 90 A, by the addition of HULLAC calculations and NIST wavelengths for L-shell ions of Ne, Mg, Si, S, Ar, and Ca. While the output models are now public, the code itself is not ready for public release. We have made significant progress toward this goal through programming improvements, including optimized error codes, a major (elusive) bug fix, and additions to our testing protocols. Secondarily, we are beginning to investigate limited uses of APEC for X-ray photoionized plasma, as for opacity modeling. Stellar coronae are being used to benchmark the atomic data in APED as part of the Emission Line Project. The models appear to be in good agreement with the observations for most of the strong lines; however, we have identified significant discrepancies in the 3s/3d line ratios not only for Fe XVII, but also for Fe XVIII and XIX. The careful analysis of these high quality X-ray spectra has also shown that there are serious problems with some of the analysis techniques in widespread use. There are several manifestations of the problem, but the main problem is that line-to-continuum ratios may be seriously mismeasured. We continue to collaborate widely with atomic physicists in order to find the best atomic data, and to solicit calculations and measurements that are needed.

Brickhouse, Nancy↗

Modified Interior Distance Functions (Theory and Methods)

In this paper we introduced and developed the theory of Modified Interior Distance Functions (MIDF's). The MIDF is a Classical Lagrangian (CL) for a constrained optimization problem which is equivalent to the initial one and can be obtained from the latter by monotone transformation both the objective function and constraints. In contrast to the Interior Distance Functions (IDF's), which played a fundamental role in Interior Point Methods (IPM's), the MIDF's are defined on an extended feasible set and along with center, have two extra tools, which control the computational process: the barrier parameter and the vector of Lagrange multipliers. The extra tools allow to attach to the MEDF's very important properties of Augmented Lagrangeans. One can consider the MIDFs as Interior Augmented Lagrangeans. It makes MIDF's similar in spirit to Modified Barrier Functions (MBF's), although there is a fundamental difference between them both in theory and methods. Based on MIDF's theory, Modified Center Methods (MCM's) have been developed and analyzed. The MCM's find an unconstrained minimizer in primal space and update the Lagrange multipliers, while both the center and the barrier parameter can be fixed or updated at each step. The MCM's convergence was investigated, and their rate of convergence was estimated. The extension of the feasible set and the special role of the Lagrange multipliers allow to develop MCM's, which produce, in case of nondegenerate constrained optimization, a primal and dual sequences that converge to the primal-dual solutions with linear rate, even when both the center and the barrier parameter are fixed. Moreover, every Lagrange multipliers update shrinks the distance to the primal dual solution by a factor 0 less than gamma less than 1 which can be made as small as one wants by choosing a fixed interior point as a 'center' and a fixed but large enough barrier parameter. The numericai realization of MCM leads to the Newton MCM (NMCM). The approximation for the primal minimizer one finds by Newton Method followed by the Lagrange multipliers update. Due to the MCM convergence, when both the center and the barrier parameter are fixed, the condition of the MDF Hessism and the neighborhood of the primal ninimizer where Newton method is 'well' defined remains stable. It contributes to both the complexity and the numerical stability of the NMCM.

Polyak, Roman A.↗

Preburner of Staged Combustion Rocket Engine

A regeneratively cooled LOX/hydrogen staged combustion assembly system with a 400:1 expansion area ratio nozzle utilizing an 89,000 Newton (20,000 pound) thrust regeneratively cooled thrust chamber and 175:1 tubular nozzle was analyzed, assembled, and tested. The components for this assembly include two spark/torch oxygen-hydrogen igniters, two servo-controlled LOX valves, a preburner injector, a preburner combustor, a main propellant injector, a regeneratively cooled combustion chamber, a regeneratively cooled tubular nozzle with an expansion area ratio of 175:1, an uncooled heavy-wall steel nozzle with an expansion area ratio of 400:1, and interconnecting ducting. The analytical effort was performed to optimize the thermal and structural characteristics of each of the new components and the ducting, and to reverify the capabilities of the previously fabricated components. The testing effort provided a demonstration of the preburner/combustor chamber operation, chamber combustion efficiency and stability, and chamber and nozzle heat transfer.

Yost, M. C.↗

Enhancing durability and activity toward oxygen evolution reaction using single-site Re-doped NiFeO x catalysts at ampere-level

NiFeO x materials are known as among the most active catalysts toward oxygen evolution reaction (OER) for hydrogen generation in alkaline media. Nevertheless, the long-term durability of NiFeO x catalysts for OER is still too far to the industrial application. Herein, we prepared a NiFeReO x catalyst with single-site Re dopants and observed that the single-site Re dopants could significantly enhance the durability without compromising the activity. A cell voltage of 1.82 V without iR correction is noted at the current density of 3000 mA cm –2 in anion-exchange membrane water electrolyzer (AEMWE) with NiFeReO x catalyst, and a very small degradation is observed under 2000 and 1000 mA cm –2 , which remarkably outperforms the pristine NiFeO x . Additionally, the overpotential of 305 mV at 10 mA cm –2 is achieved with the NiFeReO x catalyst, which is lower than 50 mV compared with the pristine NiFeO x catalyst, together with a smaller Tafel slope of 54.3 mV dec -1 . The boosted OER durability and activity of the NiFeReO x catalyst could be attributed to the strong electron-withdrawing property of Re 7+ single atoms leading to the electronic structure optimization and stabilization of Ni/Fe active sites. Our insights propose a new path for designing NiFeO x catalysts with high durability and activity toward OER.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimization of Time-Dependent Particle Tracing Using Tetrahedral Decomposition

An efficient algorithm is presented for computing particle paths, streak lines and time lines in time-dependent flows with moving curvilinear grids. The integration, velocity interpolation and step-size control are all performed in physical space which avoids the need to transform the velocity field into computational space. This leads to higher accuracy because there are no Jacobian matrix approximations or expensive matrix inversions. Integration accuracy is maintained using an adaptive step-size control scheme which is regulated by the path line curvature. The problem of cell-searching, point location and interpolation in physical space is simplified by decomposing hexahedral cells into tetrahedral cells. This enables the point location to be done analytically and substantially faster than with a Newton-Raphson iterative method. Results presented show this algorithm is up to six times faster than particle tracers which operate on hexahedral cells yet produces almost identical particle trajectories.

Kenwright, David↗

A Crewed Mission to Apophis Using a Hybrid Bimodal Nuclear Thermal Electric Propulsion (BNTEP) System

A BNTEP system is a dual propellant, hybrid propulsion concept that utilizes Bimodal Nuclear Thermal Rocket (BNTR) propulsion during high thrust operations, providing 10's of kilo-Newtons of thrust per engine at a high specific impulse (Isp) of 900 s, and an Electric Propulsion (EP) system during low thrust operations at even higher Isp of around 3000 s. Electrical power for the EP system is provided by the BNTR engines in combination with a Brayton Power Conversion (BPC) closed loop system, which can provide electrical power on the order of 100's of kWe. High thrust BNTR operation uses liquid hydrogen (LH2) as reactor coolant propellant expelled out a nozzle, while low thrust EP uses high pressure xenon expelled by an electric grid. By utilizing an optimized combination of low and high thrust propulsion, significant mass savings over a conventional NTR vehicle can be realized. Low thrust mission events, such as midcourse corrections (MCC), tank settling burns, some reaction control system (RCS) burns, and even a small portion at the end of the departure burn can be performed with EP. Crewed and robotic deep space missions to a near Earth asteroid (NEA) are best suited for this hybrid propulsion approach. For these mission scenarios, the Earth return V is typically small enough that EP alone is sufficient. A crewed mission to the NEA Apophis in the year 2028 with an expendable BNTEP transfer vehicle is presented. Assembly operations, launch element masses, and other key characteristics of the vehicle are described. A comparison with a conventional NTR vehicle performing the same mission is also provided. Finally, reusability of the BNTEP transfer vehicle is explored.

Electric Propulsion↗

Extrusion‐Spheronization of Mock Energetic Materials

The primary method for producing plastic bonded explosive (PBX) granules, or “prills”, has remained relatively unchanged for 70 years despite the complex nature of the process. In this work, we demonstrate the feasibility of using an extrusion‐spheronization technique to produce prills for PBX applications. We begin by detailing an inert formulation with similar properties of PBXs and then demonstrate the extrusion‐spheronization processing of these materials. A study is then performed where the spheronization process time of the extruded materials is varied and the resultant prills are morphologically characterized. Further, these prills are then pressed into high‐density articles and subject to compression testing to elucidate trends in process, properties, and performance. It was found that, for our formulation (95 wt.% melamine/5 wt.% polymer binder) and process, a spheronization time of 60 s yielded relatively uniform particles that exhibited improved poured, tapped, and pressed densities. Mechanical strength did not have a strong trend with process time as all spheronized materials had similar peak compression stress at failure. After further optimization, extrusion‐spheronization may be a promising path for future PBX formulation.

extrusion spheronization↗

Robust Event Simulation Variants Endowed to Ns-3 for General Exploration

Sandia's ns-3 contributions are modifications to the open source ns-3 network simulator that enable simulation speedup. One improvement, for example, removes the need for simulating the transmission and receipt of packets between nodes that are too distant to actually be able to communicate. In large-scale simulations this optimization has shown significant gains in performance. Additional contributions will pave the way for parallel discrete simulation (PDES) in ns-3. SAND2020-12456 O Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Dickson, Joseph↗

Local practically safe extremum seeking with assignable rate of attractivity to the safe set

We present Assignably Safe Extremum Seeking (ASfES), an algorithm designed to minimize a measured, static objective function while maintaining a measured, static metric of safety (a control barrier function or CBF) to be positive in a practical sense. We ensure that for trajectories with safe initial conditions, the violation of safety can be made arbitrarily small through appropriately chosen design constants. We also guarantee an assignable “attractivity” rate: from unsafe initial conditions, the trajectories approach the safe set, in the sense of the measured CBF, at a rate no slower than a user-assigned rate. Similarly, from safe initial conditions, the trajectories approach the unsafe set, in the sense of the CBF, no faster than the assigned attractivity rate. The feature of assignable attractivity is not present in the semiglobal version of safe extremum seeking, where the semiglobality of convergence is achieved by slowing the adaptation. We also demonstrate local convergence of the parameter to a neighborhood of the minimum of a quadratic objective function constrained to the safe set with a linear CBF. The ASfES algorithm and analysis are multivariable, but we also extend the algorithm to a Newton-Based ASfES scheme which we show is only useful in the scalar case. The proven properties of the designs are illustrated through simulation examples.

42 ENGINEERING↗

The SWIFT Gamma-Ray Burst X-Ray Telescope

The Swift Gamma-Ray Burst Explorer is designed to make prompt multi-wavelength observations of Gamma-Ray Bursts and GRB afterglows. The X-ray Telescope enables Swift to determine GRB positions with a few arcseconds accuracy within 100 seconds of the burst onset. The XRT utilizes a mirror set built for JET-X and an XMM-Newton/ EPIC MOS CCD detector to provide a sensitive broad-band (0.2-10 keV) X-ray imager with an effective area of more than 120 sq cm at 1.5 keV, a field of view of 23.6 x 23.6 arcminutes, and an angular resolution of 18 arcseconds (HPD). The detection sensitivity is 2x10(exp 14) erg/sq cm/s in 10(exp 4) seconds. The instrument provides automated source detection and position reporting within 5 seconds of target acquisition. It can also measure the redshifts of GRBs with Iron line emission or other spectral features. The XRT operates in an auto-exposure mode, adjusting the CCD readout mode automatically to optimize the science return as the source intensity fades. The XRT measures spectra and lightcurves of the GRB afterglow beginning about a minute after the burst and follows each burst for days or weeks. We provide an overview of the X-ray Telescope scientific background from which the systems engineering requirements were derived, with specific emphasis on the design and qualification aspects from conception through to launch. We describe the impact on cleanliness and vacuum requirements for the instrument low energy response and to maintain the high sensitivity to the fading signal of the Gamma-ray Bursts.

Hill, J. E.↗

J/Psi Photoproduction Near Threshold With CLAS12

The structure of the proton is comprised of quarks and a sea of gluons. A mechanism that can extract the characteristics of the hidden-color correlations of the nuclear wavefunction is the production of charm near threshold. Due to the fact that momentum transfer is large near threshold in the production of J/ψ , all three valence quarks must act coherently to ex- change energy for the reaction to occur. Models have been developed to predict the nature of J/ψ photoproduction at these specific energies. These include production mechanisms with the two-gluon and three-gluon exchanges. The transferred momentum dependence of the differential cross sections are sensitive to the gluonic form factors, which describe the distribution of color charge in the proton. The CLAS12 detector is capable of measuring J/ψ photoproduction at the energy range close to the threshold. Work showcased in this dissertation encompasses the preparation of the CLAS12 experiments, including the optimization of tracking reconstruction through the study of the Torus magnetic field. In terms of software, contributions were made to the CLAS12 Event Builder, a key stage of reconstruction where event-by-event information is summarized for efficient data analysis. After the run periods were successfully completed, analysis of the RG-A data commenced and an analysis framework was developed to measure the differential and total cross sections of J/ψ photoproduction in Hall B.

Newton, Joseph↗

A Residuals Approach to Filtering, Smoothing and Identification for Static Distributed Systems

An approach for state estimation and identification of spatially distributed parameters embedded in static distributed (elliptic) system models is advanced. The method of maximum likelihood is used to find parameter values that maximize a likelihood functional for the system model, or equivalently, that minimize the negative logarithm of this functional. To find the minimum, a Newton-Raphson search is conducted that from an initial estimate generates a convergent sequence of parameter estimates. For simplicity, a Gauss-Markov approach is used to approximate the Hessian in terms of products of first derivatives. The gradient and approximate Hessian are computed by first arranging the negative log likelihood functional into a form based on the square root factorization of the predicted covariance of the measurement process. The resulting data processing approach, referred to here by the new term of predicted data covariance square root filtering, makes the gradient and approximate Hessian calculations very simple. A closely related set of state estimates is also produced by the maximum likelihood method: smoothed estimates that are optimal in a conditional mean sense and filtered estimates that emerge from the predicted data covariance square root filter.

Rodriguez, G.↗

A flexible linear diffusion acceleration to k-eigenvalue neutron transport with SN discontinuous finite element method

In this paper, we derive a flexible linear diffusion acceleration (LDA) for k-eigenvalue neutron transport discretized with discontinuous finite element method (DFEM) and discrete ordinates(SN). This LDA is based on our two pieces of previous works: the flexible non linear diffusion acceleration (NDA) for DFEM-SN and LDA for k-eigenvalue neutron transport using pre-conditioned Jacobian-free Newton-Krylov with self-adjoint angular flux (SAAF), continuous finite element method(CFEM), and SN. We point out the differences between LDA and NDA for DFEM-SN and the difference between DFEM-SN and SAAF-CFEM-SN for LDA. Numerical tests are presented to compare the convergence behaviour of NDA and LDA. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Scalable Techniques for Stochastic Power Flow Problems (Final Report)

The proposed research focuses on developing scalable algorithms for two-stage security-constrained OPF problems with AC power flow constraints, a class of problems complicated by (i) scale arising from a scenario representation; and (ii) the presence of nonlinearity, nonconvexity, and possibly second-stage discreteness or complementarity. Unfortunately, most existing solvers cannot contend with both challenges simultaneously; accordingly, the proposed research focuses on developing solution techniques that can both scale with the number of scenarios and contend with nonconvexity and second-stage complementarity. We consider three avenues for addressing such problems: (i) Variable sample-size SQP (VS-SQP) methods that combine sparse Quasi-Newton updates with a scalable variance-reduced stochastic gradient scheme for stochastic QP subproblems, allowing for contending with second-stage complementarity via regularization; (ii) Variable sample-size stochastic Interior-point (VS-sIP) schemes that propose a sampling-based regularized (to allow for contending with complementarity) interior-point schemes in which a Schur-complement technique is employed for decomposing the Newton direction computation step; (iii) Variable sample-size tractable ADMM (VS-tADMM) schemes combine variable sample-sizes with carefully designed techniques for resolving each of the nonconvex updates (by leveraging the QCQP structures). We intend to compare the three schemes using performance profiles in terms of solution quality, scalability, etc. and then select one scheme which will then be developed and further refined in Python for purposes of the GO competition.

42 ENGINEERING↗

eddy Users Manual

eddy is a collection of tools - nonlinear solvers, meshing, post-processing, visualization, optimization, etc. - for performing scale-resolving simulations of multi-physics applications. The framework is designed to enable advanced R&D on a variety of topics by leveraging a mature capability for scale resolving simulations, and simultaneously be an appropriate tool for application analysis and support. Currently, eddy is at a relatively low technical readiness level (TRL), and users and developers should maintain appropriate expectations. The technical details behind eddy are outlined in several publications which can be consulted for more information [1–10]. The solvers are built around an unstructured high-order capability, and heavily utilize the tensor product sum-factorization approach for efficiency. The unsteady formulation utilizes a fully implicit space-time approach with a matrix-free Newton- Krylov method. A primitive steady-state solver is available for testing purposes, but is not expected to converge for all but simple verification cases. The Navier-Stokes fluid solvers do not support either RANS or hybrid-RANS capability, only LES and wall-modeled LES approaches. All of the solvers within eddy support three modes of operation: a primal solve of the full nonlinear problem, and two linearization approaches of the primal solve - the ad joint and the tangent solution. Details on how to select and use these three modes are outlined in Sec. 3.

Murman, Scott M.↗