Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “field 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 415 records · Page 23

Removing cosmic-ray hits from multiorbit HST Wide Field Camera images

We present an optimized algorithm that removes cosmic rays ('CRs') from multiorbit Hubble Space Telescope (HST) Wide Field/Planetary Camera ('WF/PC') images. It computes the image noise in every iteration from the WF/PC CCD equation. This includes all known sources of random and systematic calibration errors. We test this algorithm on WF/PC stacks of 2-12 orbits as a function of the number of available orbits and the formal Poissonian sigma-clipping level. We find that the algorithm needs greater than or equal 4 WF/PC exposures to locate the minimal sky signal (which is noticeably affected by CRs), with an optimal clipping level at 2-2.5 x sigma(sub Poisson). We analyze the CR flux detected on multiorbit 'CR stacks,' which are constructed by subtracting the best CR filtered images from the unfiltered 8-12 orbit average. We use an automated object finder to determine the surface density of CRS as a function of the apparent magnitude (or ADU flux) they would have generated in the images had they not been removed. The power law slope of the CR 'counts' (gamma approximately = 0.6 for N(m) m(exp gamma)) is steeper than that of the faint galaxy counts down to V approximately = 28 mag. The CR counts show a drop off between 28 less than or approximately V less than or approximately 30 mag (the latter is our formal 2 sigma point source sensitivity without spherical aberration). This prevents the CR sky integral from diverging, and is likely due to a real cutoff in the CR energy distribution below approximately 11 ADU per orbit. The integral CR surface density is less than or approximately 10(exp 8)/sq. deg, and their sky signal is V approximately = 25.5-27.0 mag/sq. arcsec, or 3%-13% of our NEP sky background (V = 23.3 mag/sq. arcsec), and well above the EBL integral of the deepest galaxy counts (B(sub J) approximately = 28.0 mag/sq. arcsec). We conclude that faint CRs will always contribute to the sky signal in the deepest WF/PC images. Since WFPC2 has approximately 2.7x lower read noise and a thicker CCD, this will result in more CR detections than in WF/PC, potentially affecting approximately 10%-20% of the pixels in multiorbit WFPC2 data cubes.

Windhorst, Rogier A.↗

Three-dimensional magnetohydrodynamic modeling of auto-magnetizing liner implosions on the Z accelerator

Auto-magnetizing (AutoMag) liners are cylindrical tubes that employ helical current flow to produce strong internal axial magnetic fields prior to radial implosion on ~100 ns timescales. AutoMag liners have demonstrated strong uncompressed axial magnetic field production (>100 T) and remarkable implosion uniformity during experiments on the 20 MA Z accelerator. However, both axial field production and implosion morphology require further optimization to support the use of AutoMag targets in magnetized liner inertial fusion (MagLIF) experiments. Data from experiments studying the initiation and evolution of dielectric flashover in AutoMag targets on the Mykonos accelerator have enabled the advancement of magnetohydrodynamic (MHD) modeling protocols used to simulate AutoMag liner implosions. Implementing these protocols using ALEGRA has improved the comparison of simulations to radiographic data. Specifically, both the liner in-flight aspect ratio and the observed width of the encapsulant-filled helical gaps during implosion in ALEGRA simulations agree more closely with radiography data compared to previous GORGON simulations. Although simulations fail to precisely reproduce the measured internal axial magnetic field production, improved agreement with radiography data inspired the evaluation of potential design improvements with newly developed modeling protocols. Three-dimensional MHD simulation studies focused on improving AutoMag target designs, specifically seeking to optimize the axial magnetic field production and enhance the cylindrical implosion uniformity for MagLIF. Importantly, by eliminating the driver current prepulse and reducing the initial inter-helix gap widths in AutoMag liners, simulations indicate that the optimal 30–50 T range of precompressed axial magnetic field for MagLIF on Z can be accomplished concurrently with improved cylindrical implosion uniformity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optimization of the N ∥ Upshift in the DIII-D high field side lower hybrid current drive experiment

High field side lower hybrid current drive (LHCD) is one potential candidate for efficient non-inductive current drive in tokamak power plants, and the first test of this technology will occur on the DIII-D tokamak during the 2024 campaign. Previous LFS launch experiments operated in the multi-pass regime and relied on scrape-off layer interactions to close the spectral gap. In the DIII-D experiment, single-pass damping is achievable via an upshift in the parallel refractive index N ∥ caused by mode converting twice (slow → fast → slow). This mode conversion affects the ray trajectories and can lead to enhanced N ∥ upshift depending on where mode conversion occurs. Compared to multi-pass absorption experiments, the optimization of launched N ∥ and plasma parameters can be counter-intuitive: increased density may increase efficiency and smaller N ∥,launch tend to damp closer to the separatrix. A hard x-ray camera installed to measure the bremsstrahlung (50–250 keV) radiation from LHCD-generated fast electrons is capable of verifying the trends reporting in this paper through comparison to the ray-tracing/Fokker–Planck codes GENRAY/CQL3D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Self-Consistent Field Methods for Excited States in Strong Magnetic Fields: A Comparison Between Energy - and Variance-based Approaches

An efficient implementation of geometrical derivatives at the Hartree–Fock (HF) and current-density functional theory (CDFT) levels is presented for the study of molecular structure in strong magnetic fields. The required integral derivatives are constructed using a hybrid McMurchie–Davidson and Rys quadrature approach, which combines the amenability of the former to the evaluation of derivative integrals with the efficiency of the latter for basis sets with high angular momentum. In addition to its application to evaluating derivatives of four-center integrals, this approach is also applied to gradients using the resolution-of-the-identity approximation, enabling efficient optimization of molecular structure for many-electron systems under a strong magnetic field. The CDFT contributions have been implemented for a wide range of density functionals up to and including the meta-GGA level with current-density dependent contributions and (range-separated) hybrids for the first time. Illustrative applications are presented to the OH and benzene molecules, revealing the rich and complex chemistry induced by the presence of an external magnetic field. Challenges for geometry optimization in strong fields are highlighted, along with the requirement for careful analysis of the resulting electronic structure at each stationary point. The importance of correlation effects is examined by comparison of results at the HF and CDFT levels. The present implementation of molecular gradients at the CDFT level provides a cost-effective approach to the study of molecular structure under strong magnetic fields, opening up many new possibilities for the study of chemistry in this regime.

Comparison↗

Mathmatical modeling for diffractive optics

We consider a 'diffractive optic' to be a biperiodic surface separating two half-spaces, each having constant constitutive parameters; within a unit cell of the periodic surface and across the transition zone between the two half-spaces, the constitutive parameters can be a continuous, complex-valued function. Mathematical models for diffractive optics have been developed, and implemented as numerical codes, both for the 'direct' problem and for the 'inverse' problem. In problems of the 'direct' class, the diffractive optic is specified, and the full set of Maxwell's equations is cast in a variational form and solved numerically by a finite element approach. This approach is well-posed in the sense that existence and uniqueness of the solution can be proved and specific convergence conditions can be derived. An example of a metallic grating at a Wood anomaly is presented as a case where other approaches are known to have convergence problems. In problems of the 'inverse' class, some information about the diffracted field (e.g., the far-field intensity) is given, and the problem is to find the periodic structure in some optimal sense. Two approaches are described: phase reconstruction in the far-field approximation; and relaxed optimal design based on the Helmholtz equation. Practical examples are discussed for each approach to the inverse problem, including array generators in the far-field case and antireflective structures for the relaxed optimal design.

Dobson, David↗

A Generalized Variational Principle with Applications to Excited State Mean Field Theory

We present a generalization of the variational principle that is compatible with any Hamiltonian eigenstate that can be specified uniquely by a list of properties. This variational principle appears to be compatible with a wide range of electronic structure methods, including mean field theory, density functional theory, multireference theory, and quantum Monte Carlo. Like the standard variational principle, this generalized variational principle amounts to the optimization of a nonlinear function that, in the limit of an arbitrarily flexible wave function, has the desired Hamiltonian eigenstate as its global minimum. Unlike the standard variational principle, it can target excited states and select individual states in cases of degeneracy or near-degeneracy. As an initial demonstration of how this approach can be useful in practice, we employ it to improve the optimization efficiency of excited state mean field theory by an order of magnitude. With this improved optimization, we are able to demonstrate that the accuracy of the corresponding second-order perturbation theory rivals that of singles-and-doubles equation-of-motion coupled cluster in a substantially broader set of molecules than could be explored by our previous optimization methodology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multiconfiguration self-consistent-field calculation of the dipole moment function of CO/X 1 sigma +/

Using the optimized valence configurations (OVC) multiconfiguration self-consistent-field (MCSCF) method, the dipole moment function for the ground state of CO in the vicinity of the equilibrium internuclear distance has been calculated. The OVC MCSCF calculation results are compared with existing Hartree-Fock and configuration interaction treatments of this molecule at single points and also the dipole moment function deduced from experimental infrared intensities. A general prescription for constructing OVC wavefunctions for diatomic molecules is also presented.

Billingsley, F. P., II↗

Propulsive trajectory optimization to minimize surface contamination

MOTIVATION: We present an optimization technique for propulsive vehicles that autonomously minimizes contamination during surface approach and landing. In addition to short-range hoppers, the optimization technique is also fully applicable to traditional orbit-to-surface landers. This study addresses scenarios where surface alterations from propulsion events are counterproductive or hazardous to the mission objectives. This is of immediate interest for landers (whether human or robotic), that may rely on pristine soils collected in the immediate vicinity of landing sites to accomplish science investigations, mining, or ISRU surface operations. Such missions are averse to various surface-plume interactions such as thermal scoring, physical agitation, and contamination. The capability can be applied with minimal impact to the baseline mission concept. METHODS: Optimization algorithms have been developed to calculate descent trajectories and maneuvers, thrust magnitude, and attitude for various mission cases. These parameters are determined as an optimal solution when minimizing either fuel consumption, contamination deposited at the landing site, or some weighted combination of both. Among constraints imposed on the solution, we examined pitch rate, vertical takeoff and vertical landing (VTVL) requirements, size of the contamination zone, and minimum ground clearance during flight. This tool provides unique, non-intuitive solutions and can be a valuable resource for mission planners. RESULTS: A variety of agile trajectory solutions were obtained, each yielding different reductions in landing site contamination and corresponding to only modest increases in fuel consumption. Several optimal trajectories were obtained by varying the contamination weight in the fitness function. As expected, when the contamination weight is zero, the trajectory appears close to parabolic since the optimization scheme only attempts to minimize for fuel utilization, yielding essentially, the expected ballistic trajectory. Notably for contamination weights greater than zero, trajectory inflections are observed in the descent phase, which manifests as hovering or additional, mini “pseudo hops” before the final touchdown. A trajectory inflection is characterized by arresting the majority of the spacecraft vertical velocity component at a coordinate outside of the landing target, and without violating ground clearance constraints. FUTURE WORK: Our optimization technique is ready for laboratory or field demonstrations to validate the sophisticated maneuvering solutions obtained for fuel optimization and surface preservation. An appropriate testbed would validate the optimal guidance algorithms, the navigation system, and sensor suite by emulating vehicle flight in closed loop robotic tests. Critically, these algorithms could then be ported to flight software for implementation.

surface contamination↗

Recent advances in efficient long-life, eye-safe solid state and CO2 lasers for laser radar applications

The key problems in the development of eye-safe solid-state lasers are discussed, taking into account the energy transfer mechanisms between the complicated energy level manifolds of the Tm, Ho, Er ion dopants in hosts with decreasing crystal fields such as YAG or YLF. Optimization of energy transfer for efficient lasing through choice of dopant concentration, power density, crystal field and temperature is addressed. The tailoring of energy transfer times to provide efficient energy extraction for short pulses used in DIAL and Doppler lidar is considered. Recent advances in Pt/SnO2 oxide catalysts and other noble metal/metal oxide combinations for CO2 lasers are discussed. Emphasis is given to the dramatic effects of small quantities of H2O vapor for increasing the activity and lifetime of Pt/SnO2 catalysts and to increased lifetime operation with rare isotope (C-12)(O-18)2 lasing mixtures.

Hess, R. V.↗

Interferometers adaptations to lidars

To perform daytime measurements of the density and temperature by Rayleigh lidar, it is necessary to select the wavelength with a very narrow spectral system. This filter is composed by an interference filter and a Fabry Perot etalon. The Fabry Perot etalon is the more performent compound, and it is necessary to build a specific optic around it. The image of the entrance pupil or the field diaphragm is at the infinite and the other diaphragm is on the etalon. The optical quality of the optical system is linked to the spectral resolution of the system to optimize the reduction of the field of view. The resolution is given by the formula: R = 8(xD/Fd)exp 2 where R = lambda/delta(lambda), x = diameter of the field diaphragm, D = diameter of the reception mirror, F = focal length of the telescope, and d = useful diameter of the etalon. In the Doppler Rayleigh lidars, the PF interferometer is the main part of the experiment and the exact spectral adaptation is the most critical problem. In the spectral adaptation of interferometers, the transmittance of the system will be acceptable if the etalon is exactly adjusted to the wavelength of the laser. It is necessary to work with a monomode laser, and adjust the shift to the bandpass of the interferometer. We are working with an interferometer built with molecular optical contact. This interferometer is put in a special pressure closed chamber.

Porteneuve, J.↗

Optimal two-impulse transfer between specified terminal states of Keplerian orbits

The integrals of motion for an inverse-square force field and a necessary condition for optimality are used in a transfer orbit co-ordinate system to formulate the solution of the optimal two-impulse transfer between fixed position and velocity vectors on Keplerian orbits. In this co-ordinate system, the equations reveal two asymptotes that are useful in analysing the solution. It is shown that there are only two real extremals (both minimums) which are separated by an asymptote. An approximate analytic solution is also obtained by ignoring a quadratic term whose coefficient is approximately zero for a large class of orbit transfer problems.

Jezewski, D. J.↗

Optimization of Low-Thrust Spiral Trajectories by Collocation

As NASA examines potential missions in the post space shuttle era, there has been a renewed interest in low-thrust electric propulsion for both crewed and uncrewed missions. While much progress has been made in the field of software for the optimization of low-thrust trajectories, many of the tools utilize higher-fidelity methods which, while excellent, result in extremely high run-times and poor convergence when dealing with planetocentric spiraling trajectories deep within a gravity well. Conversely, faster tools like SEPSPOT provide a reasonable solution but typically fail to account for other forces such as third-body gravitation, aerodynamic drag, solar radiation pressure. SEPSPOT is further constrained by its solution method, which may require a very good guess to yield a converged optimal solution. Here the authors have developed an approach using collocation intended to provide solution times comparable to those given by SEPSPOT while allowing for greater robustness and extensible force models.

Falck, Robert D.↗

Searching for optimal THz generation through calculations of the asymmetry of photoelectron momentum distributions by an improved strong-field-approximation method

From the solution of the time-dependent Schrödinger equation (TDSE), it has been shown that in a two-color laser pulse the relative phase that optimizes terahertz (THz) radiation also optimizes the asymmetry of the angular distribution of photoelectrons. Here, in this work, we show that a second-order strong-field approximation can accurately calculate the asymmetry of photoelectrons and thus can locate the phase delay that optimizes THz generation, doing so four orders of magnitude faster than TDSE calculations. We further trace that this is possible because THz emission originates from a free-free radiative transition between the rescattered electron and the target ion, similar to other electron-ion collisions in the laser field, e.g., high-harmonic generation and nonsequential double ionization. Our results pave the way to locate the optimal phase delay for any typical laser pulses for maximal THz generation in the laboratory.

74 ATOMIC AND MOLECULAR PHYSICS↗

A Reduced Order Model of the Linearized Incompressible Navier-Strokes Equations for the Sensor/Actuator Placement Problem

A reduced order modeling approach of the Navier-Stokes equations is presented for the design of a distributed optimal feedback kernel. This approach is based oil a Krylov subspace method where significant modes of the flow are captured in the model This model is then used in all optimal feedback control design where sensing and actuation is performed oil tile entire flow field. This control design approach yields all optimal feedback kernel which provides insight into the placement of sensors and actuators in the flow field. As all evaluation of this approach, a two-dimensional shear layer and driven cavity flow are investigated.

Allan, Brian G.↗

Selecting a spatial resolution for estimation of per-field green leaf area index

For any application of multispectral scanner (MSS) data, a user is faced with a number of choices concerning the characteristics of the data; one of these is their spatial resolution. A pilot study was undertaken to determine the spatial resolution that would be optimal for the per-field estimation of green leaf area index (GLAI) in grassland. By reference to empirically-derived data from three areas of grassland, the suitable spatial resolution was hypothesized to lie in the lower portion of a 2-18 m range. To estimate per-field GLAI, airborne MSS data were collected at spatial resolutions of 2 m, 5 m and 10 m. The highest accuracies of per-field GLAI estimation were achieved using MSS data with spatial resolutions of 2 m and 5 m.

Curran, Paul J.↗

Predicting System Response at Unmeasured Locations

Traditional techniques to derive dynamic specification for components have a great deal of uncertainty. One of the major sources of uncertainty is that the number of response measurements in the operational system environment is insufficient to determine the component motion. This inadequacy is due to logistical limitations for data recording in field testing and space limitations for accelerometers, strain gages and associated wiring. Available measurements are often some distance from the component and therefore do not represent component motion. Typical straight-line envelopes of these unrepresentative measurements guarantee an increase in the uncertainty. In this paper multiple methods are attempted to expand a sparse set of field test measurements on a system to responses of interest that cannot be measured in the field due to the limitations. Proof of concept is demonstrated on the Modal Analysis Test Vehicle (MATV). The responses of interest, known as “truth responses”, are measured in a system vibration environment along with an optimized sparse set of 30 field responses. Methods to expand the field responses to the truth responses are demonstrated by comparing the acceleration spectral density of the expanded response to the measured response. Two methods utilize a validated finite element model of the MATV. One is developed from purely experiment based frequency response functions of a laboratory pre-test. Finally, these approaches are designed to drastically reduce the uncertainty of the component in-service motion as a basis for developing specifications that are guaranteed to be conservative with a known (instead of unknown) conservatism.

42 ENGINEERING↗

Spatial Statistical Data Fusion (SSDF)

As remote sensing for scientific purposes has transitioned from an experimental technology to an operational one, the selection of instruments has become more coordinated, so that the scientific community can exploit complementary measurements. However, tech nological and scientific heterogeneity across devices means that the statistical characteristics of the data they collect are different. The challenge addressed here is how to combine heterogeneous remote sensing data sets in a way that yields optimal statistical estimates of the underlying geophysical field, and provides rigorous uncertainty measures for those estimates. Different remote sensing data sets may have different spatial resolutions, different measurement error biases and variances, and other disparate characteristics. A state-of-the-art spatial statistical model was used to relate the true, but not directly observed, geophysical field to noisy, spatial aggregates observed by remote sensing instruments. The spatial covariances of the true field and the covariances of the true field with the observations were modeled. The observations are spatial averages of the true field values, over pixels, with different measurement noise superimposed. A kriging framework is used to infer optimal (minimum mean squared error and unbiased) estimates of the true field at point locations from pixel-level, noisy observations. A key feature of the spatial statistical model is the spatial mixed effects model that underlies it. The approach models the spatial covariance function of the underlying field using linear combinations of basis functions of fixed size. Approaches based on kriging require the inversion of very large spatial covariance matrices, and this is usually done by making simplifying assumptions about spatial covariance structure that simply do not hold for geophysical variables. In contrast, this method does not require these assumptions, and is also computationally much faster. This method is fundamentally different than other approaches to data fusion for remote sensing data because it is inferential rather than merely descriptive. All approaches combine data in a way that minimizes some specified loss function. Most of these are more or less ad hoc criteria based on what looks good to the eye, or some criteria that relate only to the data at hand.

Braverman, Amy J.↗