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 199 records · Page 11

Cascade Error Projection: A Learning Algorithm for Hardware Implementation

In this paper, we workout a detailed mathematical analysis for a new learning algorithm termed Cascade Error Projection (CEP) and a general learning frame work. This frame work can be used to obtain the cascade correlation learning algorithm by choosing a particular set of parameters. Furthermore, CEP learning algorithm is operated only on one layer, whereas the other set of weights can be calculated deterministically. In association with the dynamical stepsize change concept to convert the weight update from infinite space into a finite space, the relation between the current stepsize and the previous energy level is also given and the estimation procedure for optimal stepsize is used for validation of our proposed technique. The weight values of zero are used for starting the learning for every layer, and a single hidden unit is applied instead of using a pool of candidate hidden units similar to cascade correlation scheme. Therefore, simplicity in hardware implementation is also obtained. Furthermore, this analysis allows us to select from other methods (such as the conjugate gradient descent or the Newton's second order) one of which will be a good candidate for the learning technique. The choice of learning technique depends on the constraints of the problem (e.g., speed, performance, and hardware implementation); one technique may be more suitable than others. Moreover, for a discrete weight space, the theoretical analysis presents the capability of learning with limited weight quantization. Finally, 5- to 8-bit parity and chaotic time series prediction problems are investigated; the simulation results demonstrate that 4-bit or more weight quantization is sufficient for learning neural network using CEP. In addition, it is demonstrated that this technique is able to compensate for less bit weight resolution by incorporating additional hidden units. However, generation result may suffer somewhat with lower bit weight quantization.

Duong, Tuan A.↗

The Mathematics of Navigating the Solar System

In navigating spacecraft throughout the solar system, the space navigator relies on three academic disciplines - optimization, estimation, and control - that work on mathematical models of the real world. Thus, the navigator determines the flight path that will consume propellant and other resources in an efficient manner, determines where the craft is and predicts where it will go, and transfers it onto the optimal trajectory that meets operational and mission constraints. Mission requirements, for example, demand that observational measurements be made with sufficient precision that relativity must be modeled in collecting and fitting (the estimation process) the data, and propagating the trajectory. Thousands of parameters are now determined in near real-time to model the gravitational forces acting on a spacecraft in the vicinity of an irregularly shaped body. Completing these tasks requires mathematical models, analyses, and processing techniques. Newton, Gauss, Lambert, Legendre, and others are justly famous for their contributions to the mathematics of these tasks. More recently, graduate students participated in research to update the gravity model of the Saturnian system, including higher order gravity harmonics, tidal effects, and the influence of the rings. This investigation was conducted for the Cassini project to incorporate new trajectory modeling features in the navigation software. The resulting trajectory model will be used in navigating the 4-year tour of the Saturnian satellites. Also, undergraduate students are determining the ephemerides (locations versus time) of asteroids that will be used as reference objects in navigating the New Millennium's Deep Space 1 spacecraft autonomously.

Hintz, Gerald↗

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↗

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.↗

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↗

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.↗

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.↗

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.↗

Optimal reusable-tug and expendable-kickstage trajectories for high-energy planetary missions including correction for nodal precession

Equations are derived by using branched trajectory optimization techniques and the maximum principle to maximize the payload capability of a reusable tug/expendable kickstage vehicle configuration for planetary missions. The two stages and the payload are launched into a low earth orbit by a single space shuttle. The analysis includes correction for precession of the orbit. This correction is done by the tug. The tug propels the payload and the kickstage to an energy beyond earth escape and returns within a specified time to the precessed orbit. After separating from the tug, the kickstage accelerates the payload to the required injection conditions. Planetary injection conditions are specified by the mission energy and a fixed declination and right ascension of the outgoing asymptote. The multipoint boundary value problem resulting from the analysis is solved by a Newton-Raphson iteration technique. Partial derivatives of the boundary conditions are obtained by perturbing the initial conditions one at a time, integrating the trajectory and adjoint equations, and observing the changes in boundary conditions. Maximum payload capability is derived for two typical mission energies. In addition, the variations of several mission and stage parameters are also examined.

Borsody, J.↗

Multigrid Strategies for Viscous Flow Solvers on Anisotropic Unstructured Meshes

Unstructured multigrid techniques for relieving the stiffness associated with high-Reynolds number viscous flow simulations on extremely stretched grids are investigated. One approach consists of employing a semi-coarsening or directional-coarsening technique, based on the directions of strong coupling within the mesh, in order to construct more optimal coarse grid levels. An alternate approach is developed which employs directional implicit smoothing with regular fully coarsened multigrid levels. The directional implicit smoothing is obtained by constructing implicit lines in the unstructured mesh based on the directions of strong coupling. Both approaches yield large increases in convergence rates over the traditional explicit full-coarsening multigrid algorithm. However, maximum benefits are achieved by combining the two approaches in a coupled manner into a single algorithm. An order of magnitude increase in convergence rate over the traditional explicit full-coarsening algorithm is demonstrated, and convergence rates for high-Reynolds number viscous flows which are independent of the grid aspect ratio are obtained. Further acceleration is provided by incorporating low-Mach-number preconditioning techniques, and a Newton-GMRES strategy which employs the multigrid scheme as a preconditioner. The compounding effects of these various techniques on speed of convergence is documented through several example test cases.

Movriplis, Dimitri J.↗

Human Mars Entry, Descent, and Landing Architecture Study: Phase 3 Summary

Over the past four years, NASA has directed the Entry, Descent and Landing Architecture Study (EDLAS) to evaluate candidate technologies to deliver human-scale vehicles (carrying 20t payloads) to the surface of Mars. The study focused on two vehicles in Phase 3 (during 2018 and 2019), one low and one mid lift-to-drag (L/D) vehicle. This paper summarizes updates to the structural and aeroshell outer mold line design for the mid-L/D Mars entry lander vehicle that occurred during EDLAS Phase 3. The design efforts used parametric optimization to arrive at a solution for total vehicle primary and secondary structure mass, and vehicle center-of-gravity location. Additionally, the end of Phase 2 provided an updated surface payload manifest that required three 22 t landers. This paper also summarizes the results of the packaging feasibility study of the updated payload configurations for those three reference missions. The study verified an arrangement of the payload elements within the mid-L/D vehicle that satisfied the center of gravity location limitations for flight.

Calderon, Damien↗

Spectroscopic Determination of the Physical Conditions in Hot, Optically Thin Sources

This annual report covers the period from Oct 1, 2000 to Sep 30, 2001. The Astrophysics Plasma Emission Code and Database (APEC/APED), developed in part under this grant, have been upgraded to Version 1.1 and are now beginning to be used outside our research group in 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 angstroms. Stellar coronae are being used to benchmark the atomic data in APED as part of the Emission Line Project. Initial results suggest that the models for most of the strongest lines are in good agreement with the observations, including the H-like and He-like emission from dominant elements and the Fe L shell emission, both near 1 keV and in the extreme ultraviolet near 100 eV. At this point in time, we define 'good agreement' at the level of accuracy expected from the atomic physics, approximately 20-30%. In order to benchmark the spectral models beyond the canonical theoretical accuracy, we are working closely with the Chandra gratings calibration group to ensure that we are using the optimal calibration, primarily effective areas and line response functions.

Brickhouse, N.↗

Embedding impedance approximations in the analysis of SIS mixers

Future millimeter-wave radio astronomy instruments will use arrays of many SIS receivers, either as focal plane arrays on individual radio telescopes, or as individual receivers on the many antennas of radio interferometers. Such applications will require broadband integrated mixers without mechanical tuners. To produce such mixers, it will be necessary to improve present mixer design techniques, most of which use the three-frequency approximation to Tucker's quantum mixer theory. This paper examines the adequacy of three approximations to Tucker's theory: (1) the usual three-frequency approximation which assumes a sinusoidal LO voltage at the junction, and a short-circuit at all frequencies above the upper sideband; (2) a five-frequency approximation which allows two LO voltage harmonics and five small-signal sidebands; and (3) a quasi five-frequency approximation in which five small-signal sidebands are allowed, but the LO voltage is assumed sinusoidal. These are compared with a full harmonic-Newton solution of Tucker's equations, including eight LO harmonics and their corresponding sidebands, for realistic SIS mixer circuits. It is shown that the accuracy of the three approximations depends strongly on the value of omega R(sub N)C for the SIS junctions used. For large omega R(sub N)C, all three approximations approach the eight-harmonic solution. For omega R(sub N)C values in the range 0.5 to 10, the range of most practical interest, the quasi five-frequency approximation is a considerable improvement over the three-frequency approximation, and should be suitable for much design work. For the realistic SIS mixers considered here, the five-frequency approximation gives results very close to those of the eight-harmonic solution. Use of these approximations, where appropriate, considerably reduces the computational effort needed to analyze an SIS mixer, and allows the design and optimization of mixers using a personal computer.

Kerr, A. R.↗

Optimization of rod diameter in solid state lasers side pumped with multiple laser diode arrays

Results of a study to determine the optimum laser rod diameter for maximum output energy in a solid state neodymium laser transversely pumped with multiple laser diode arrays are reported here. Experiments were performed with 1.0 mm, 1.5 mm and 2.0 mm rod radii of both neodymium doped Y3Al5O12 (Nd:YAG) and La2Be2O5 (Nd:BeL) pumped with laser diode arrays having a maximum combined energy of 10.5 mJ. Equations were derived which predict the optimum rod radius and corresponding output mirror reflectivity for a given laser material and total pump energy. Predictions of the equations agreed well with the experiments for each of the laser materials which possessed significantly different laser properties from one another.

Sims, Newton, Jr.↗