SEARCH · Search NASA
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.
Report discusses applicability of conjugate-gradient algorithms to computation of forward dynamics of robotic manipulators. Rapid computation of forward dynamics essential to teleoperation and other advanced robotic applications. Part of continuing effort to find algorithms meeting requirements for increased computational efficiency and speed. Method used for iterative solution of systems of linear equations.
Discretization of the Stokes equations produces a symmetric indefinite system of linear equations. For stable discretizations, a variety of numerical methods have been proposed that have rates of convergence independent of the mesh size used in the discretization. In this paper, we compare the performance of four such methods: variants of the Uzawa, preconditioned conjugate gradient, preconditioned conjugate residual, and multigrid methods, for solving several two-dimensional model problems. The results indicate that where it is applicable, multigrid with smoothing based on incomplete factorization is more efficient than the other methods, but typically by no more than a factor of two. The conjugate residual method has the advantage of being both independent of iteration parameters and widely applicable.
The development of efficient iterative solution methods for the numerical solution of two- and three-dimensional compressible Navier-Stokes equations is discussed. Iterative time marching methods have several advantages over classical multi-step explicit time marching schemes, and non-iterative implicit time marching schemes. Iterative schemes have better stability characteristics than non-iterative explicit and implicit schemes. In this work, another approach based on the classical conjugate gradient method, known as the Generalized Minimum Residual (GMRES) algorithm is investigated. The GMRES algorithm has been used in the past by a number of researchers for solving steady viscous and inviscid flow problems. Here, we investigate the suitability of this algorithm for solving the system of non-linear equations that arise in unsteady Navier-Stokes solvers at each time step.
The Advanced Wavefront Sensing and Control Testbed (AWCT) is built as a versatile facility for developing and demonstrating, in hardware, the future technologies of wave front sensing and control algorithms for active optical systems. The testbed includes a source projector for a broadband point-source and a suite of extended scene targets, a dispersed fringe sensor, a Shack-Hartmann camera, and an imaging camera capable of phase retrieval wavefront sensing. The testbed also provides two easily accessible conjugated pupil planes which can accommodate the active optical devices such as fast steering mirror, deformable mirror, and segmented mirrors. In this paper, we describe the testbed optical design, testbed configurations and capabilities, as well as the initial results from the testbed hardware integrations and tests.
The development of the Modular High-Temperature Gas-Cooled Reactor is a significant milestone in advanced nuclear reactor technology. One of the concerns for the reactor’s safe operation is the effects of a loss-of-flow accident (LOFA) where the coolant circulators are tripped, and forced coolant flow through the core is lost. Depending on the steam generator placement, loop or intracore natural circulation develops to help transfer heat from the core to the reactor cavity, cooling system. This paper investigates the fundamental physical phenomena associated with intracore coolant natural circulation flow in a one-sixth Computational Fluid Dynamics (CFD) model of the Oregon State University High Temperature Test Facility (OSU HTTF) following a loss-of-flow accident transient. This study employs conjugate heat transfer and steady-state flow along with an SST k-ω turbulence model to characterize the phenomenon of core channel-to-channel natural convection. Previous studies have revealed the importance of complex flow distribution in the inlet and outlet plenums with the potential to generate hot coolant jets. For this reason, complete upper and lower plenum volumes are included in the analyzed computational domain. CFD results also include parametric studies performed for a mesh sensitivity analysis, generated using the STAR-CCM+ software. The resulting channel axial velocities and flow directions support the test facility scaling analysis and similarity group distortions calculation.
The theoretical and computational aspects related to the application of the Conjugate Gradient FFT (CGFFT) method in computational electromagnetics are examined. The advantages of applying the CGFFT method to a class of large scale scattering and radiation problems are outlined. The main advantages of the method stem from its iterative nature which eliminates a need to form the system matrix (thus reducing the computer memory allocation requirements) and guarantees convergence to the true solution in a finite number of steps. Results are presented for various radiators and scatterers including thin cylindrical dipole antennas, thin conductive and resistive strips and plates, as well as dielectric cylinders. Solutions of integral equations derived on the basis of generalized impedance boundary conditions (GIBC) are also examined. The boundary conditions can be used to replace the profile of a material coating by an impedance sheet or insert, thus, eliminating the need to introduce unknown polarization currents within the volume of the layer. A general full wave analysis of 2-D and 3-D rectangular grooves and cavities is presented which will also serve as a reference for future work.
The Holographic Enhanced Remote Sensing System (HERSS) consists of three primary subsystems: (1) an Image Acquisition System (IAS); (2) a Digital Image Processing System (DIPS); and (3) a Holographic Generation System (HGS) which multiply exposes a thermoplastic recording medium with sequential 2-D depth slices that are displayed on a Spatial Light Modulator (SLM). Full-parallax holograms were successfully generated by superimposing SLM images onto the thermoplastic and photopolymer. An improved HGS configuration utilizes the phase conjugate recording configuration, the 3-SLM-stacking technique, and the photopolymer. The holographic volume size is currently limited to the physical size of the SLM. A larger-format SLM is necessary to meet the desired 6 inch holographic volume. A photopolymer with an increased photospeed is required to ultimately meet a display update rate of less than 30 seconds. It is projected that the latter two technology developments will occur in the near future. While the IAS and DIPS subsystems were unable to meet NASA goals, an alternative technology is now available to perform the IAS/DIPS functions. Specifically, a laser range scanner can be utilized to build the HGS numerical database of the objects at the remote work site.
The dynamics of an externally pumped phase-conjugate resonator that uses four-wave mixing in BaTi03 is investigated experimentally. The emergence of spatiotemporal instabilities as the degree of transverse confinement is varied by changing the Fresnel number is described in detail. Local intensity time series show that relaxing the transverse confinement leads the system from a stationary state to periodic, then to quasi-periodic, motions and finally to chaotic behavior. In some regions of parameter space two frequencies are identified in the power spectra of the time series, indicating a route to chaos following the Ruelle-Takens-Newhouse scheme. Wave-front topological defects are identified by interferometry. As the system's confinement is varied, the phase-defect density and the spatial correlation index are found to follow similar trends, indicating that the observed spatiotemporal dynamics may indeed be an example of defect-mediated turbulence.
Direct solvers currently dominate commercial finite element structural software, but do not scale well in the fine granularity regime targeted by emerging parallel processors. Substructure based iterative solvers--often called also domain decomposition algorithms--lend themselves better to parallel processing, but must overcome several obstacles before earning their place in general purpose structural analysis programs. One such obstacle is the solution of systems with many or repeated right hand sides. Such systems arise, for example, in multiple load static analyses and in implicit linear dynamics computations. Direct solvers are well-suited for these problems because after the system matrix has been factored, the multiple or repeated solutions can be obtained through relatively inexpensive forward and backward substitutions. On the other hand, iterative solvers in general are ill-suited for these problems because they often must restart from scratch for every different right hand side. In this paper, we present a methodology for extending the range of applications of domain decomposition methods to problems with multiple or repeated right hand sides. Basically, we formulate the overall problem as a series of minimization problems over K-orthogonal and supplementary subspaces, and tailor the preconditioned conjugate gradient algorithm to solve them efficiently. The resulting solution method is scalable, whereas direct factorization schemes and forward and backward substitution algorithms are not. We illustrate the proposed methodology with the solution of static and dynamic structural problems, and highlight its potential to outperform forward and backward substitutions on parallel computers. As an example, we show that for a linear structural dynamics problem with 11640 degrees of freedom, every time-step beyond time-step 15 is solved in a single iteration and consumes 1.0 second on a 32 processor iPSC-860 system; for the same problem and the same parallel processor, a pair of forward/backward substitutions at each step consumes 15.0 seconds.
The Goddard Earth Observing System (GEOS-S) Data Assimilation System is now configurable to run strong constraint 4DVAR. The system combines the 4DVAR-capable Grid-point Statistical Interpolation analysis with the GEOS-S general circulation model and an early version of the tangent linear and adjoint models of the finite-volume hydrodynamics with simplified physics. The GMAO 4DVAR uses a Lanczos-based conjugate gradient algorithm in a nested resolution inner-loop setting. A range of tests and experiments will have been carried out by the time of the Symposium. Results on balance-related issues, resolution configuration, and choice of assimilation time window will be shown. A comparison with the Global Modeling and Assimilation Office (GMAO) 3DVAR and a discussion of future directions of development will also be presented.
Radiation encountered on an extended duration space mission (estimates of 3 Sieverts for a mission to Mars) poses a threat not only to human health, but also to the quality, nutritional value, and palatability of the food system. Free radicals generated by radiation interaction with foods may initiate many unwanted reactions including: 1) autoxidation in lipids that alters flavor, odor, and concentrations of essential fatty acids, and 2) depletion of antioxidants food products and dietary supplements. Studies have shown that antioxidants may provide long term health protection from oxidative stress caused by radiation exposure; therefore, consumption of antioxidants will be important. Stability of essential fatty acids is also important for astronauts long-term health status. The objectives of this study were to characterize the effects of low dose gamma-radiation on lipids and antioxidants by monitoring oxidation and reducing power, respectively, in model systems. Select oils and antioxidants were exposed to levels of gamma-radiation ranging from 0 to 1000 Gy (1 Gy = 1 Sv) using a Gammacell 220 and stored at ambient or elevated temperatures (65 C) for up to 3 months prior to analysis. A Fricke dosimeter was used to verify differences between the radiation doses administered. Primary and secondary products of lipid oxidation in soybean and peanut oils were monitored using conjugated diene and 2-thiobarbituric acid (TBARs) assays. Changes in fatty acid composition and formation and vitamin E levels were also measured. The reducing power of antioxidant compounds, including vitamins C and E and beta-carotene, was determined using the ferric reducing antioxidant power (FRAP) assay. Significant differences (alpha =0.05) were present between all radiation doses tested using the Fricke dosimeter. Increasing radiation doses above 3 Sv resulted in significantly (alpha =0.05) elevated levels of oxidation and free fatty acids in soybean and peanut oils. Decreases in concentrations of essential fatty acids upon increasing radiation doses were also observed. Increasing radiation doses caused significant (alpha =0.05) decreases in reducing power and hence the effectiveness of vitamins C and E as well as beta-carotene. This work establishes a need for quantifying the effects of space relevant radiation doses in the development of a food system for an extended duration mission and for identifying threshold radiation levels that will impact the useful shelf-life of the variety of foods that will be sent. Eventual rancidity of lipids and the loss of antioxidant bioprotective effects are major concerns for the acceptability and nutritional profile of a food system.
Adaptive optics (AO) improves the spatial resolution of microscopy by correcting optical aberrations. While its application has been well established in microscopy modalities utilizing a circular pupil, its adaptation to systems with non-circular pupils, such as Bessel-focus two-photon fluorescence microscopy (2PFM) with an annular pupil, remains relatively uncharted. Herein, we present a modal focal AO (MFAO) method for Bessel-focus 2PFM. Measuring and correcting aberration using a spatial light modulator placed in conjugation with the focal plane of the microscope objective, MFAO employs Zernike annular polynomials — a first in AO implementation — to achieve performance on par with a previous zonal AO method, but with a notably simplified optical configuration. We validated the performance of MFAO in correcting artificial and sample-induced aberrations, as well as in in vivo imaging of zebrafish larvae and mouse brains. By expanding the application of modal AO to annular pupils as well as aberration measurement and correction to a wavefront modulator at the objective focal plane, MFAO represents a notable advancement in the implementation of AO in microscopy.
Tight-binding calculations of the length dependence of the third-order molecular hyperpolarizability for polyenes and polyynes are reported. The pi-electron wave functions were determined by exploiting the limited translational symmetry of the molecules. Perturbation theory was used to calculate the longitudinal component of the electronic nonresonant hyperpolarizability. This is the first two-'band' calculation of third-order hyperpolarizabilities on finite pi-electron systems of varying length. In contrast to the results of the one-'band' models, the hyperpolarizability densities increase rapidly and then, after about 10-15 repeating units, approach an asymptotic value.
Non-thermal advanced fuel fusion trades the requirement of a large amount of recirculating tritium in the system for that of large recirculating power. Phase space engineering technologies utilizing externally injected electromagnetic fields can be applied to meet the challenge of maintaining non-thermal particle distributions at a reasonable cost. The physical processes of the phase space engineering are studied from a theoretical and algorithmic perspective. It is emphasized that the operational space of phase space engineering is limited by the underpinning symplectic dynamics of charged particles. The phase space incompressibility according to the Liouville theorem is just one of many constraints, and Gromov's non-squeezing theorem determines the minimum footprint of the charged particles on every conjugate phase space plane. In this sense and level of sophistication, the mathematical abstraction of phase space engineering is symplectic topology. To simulate the processes of phase space engineering, such as the Maxwell demon and electromagnetic energy extraction, and to accurately calculate the minimum footprints of charged particles, recently developed structure-preserving geometric algorithms can be used. The family of algorithms conserves exactly, on discretized spacetime, symplecticity and thus incompressibility, non-squeezability, and symplectic capacities. The algorithms apply to the dynamics of charged particles under the influence of external electromagnetic fields as well as the charged particle–electromagnetic field system governed by the Vlasov–Maxwell equations.
In this paper we have investigated the condensation of water vapor in a short tube. A numerical model of condensation heat transfer was incorporated in a flow network code. The flow network code that we have used in this paper is Generalized Fluid System Simulation Program (GFSSP). GFSSP is a finite volume based flow network code. Four different condensation models were presented in the paper. Soliman's correlation has been found to be the most stable in low flow rates which is of particular interest in this application. Another highlight of this investigation is conjugate or coupled heat transfer between solid or fluid. This work was done in support of NASA's International Space Station program.
The Wide-Field Infrared Survey Telescope (WFIRST) is a 2.4-meter-diameter space telescope NASA program. The payload will include a coronagraph instrument (CGI). The CGI designs under development use deformable mirrors (DM) to create a point spread function (PSF) with a dark region around the obscured star object. Electric field conjugation (EFC) is an iterative nonlinear optimization procedure that uses measurements of the electric field at the image to determine the DM displacements to modify the PSF to create the region of high contrast in the image. EFC requires a numerical model of the coronagraph to calculate the Jacobian of the system, which is used, along with regularization, to solve for the DM displacements for each iteration of the nonlinear optimization. Ideally, the coronagraph is aligned and calibrated, and the calibration data are used in the numerical model for calculating the Jacobian. However, calibration and alignment measurements always contain uncertainty resulting in calibration error. Therefore, the Jacobian calculated from the numerical model is not an exact representation of the physical coronagraph, and the resulting DM solution for an EFC iteration does not have the exact impact on the electric field of the coronagraph as predicted by the EFC. The result is slow convergence, and, as will be shown, drives the use of more restrictive regularization. Using Monte Carlo trials, we investigate the effect of calibration error on EFC convergence and regulatization. Comparison to results from the High Contrast Imaging Testbed Hybrid Lyot Coronagraph are also presented.
Approaches for estimating the composition of the matrix phase of alloys from the melt composition are reviewed. The first method is based on assigning essentially fixed stoichiometry to precipitating phases and is typified by PHACOMP. The second method uses analytical geometry to interpret phase diagrams and is applicable to a two-phase region of a six-component Ni-base system. The geometric method is also applicable to commercial Ni-base superalloys.
Continuous counter‐current microfluidic liquid–liquid extraction performs separations by flowing immiscible liquids in opposing directions within a single flow channel. In principle, this flow arrangement enables a large number of theoretical separation units in a small footprint, without using interstage valving, pumping, and phase separation. Despite its potential for excellent separation performance, this microfluidic scheme rarely appears in literature due to the requirement for capillary forces to be greater than hydrodynamic forces for stable flow. We present a novel microfluidic device and flow approaches that overcome this force‐balance challenge, enabling stable, long‐duration continuous counter‐current flow. Additionally, we cover a suite of methodologies for quantifying the performance of the microfluidic device, revealing the number of theoretical equilibrium stages achieved. The enabling technologies include a woven mesh screen‐based microfluidic device architecture that is easily fabricated outside of a clean room, surface functionalization strategies to promote conjugate (organic/aqueous) wettability, flow approaches to eliminate bubbles and carryover, and computer‐aided flow automation with optical measurement of extraction performance. The reported experiments lasted for over 36 h, terminated only at experiment conclusion, where the device still exhibited good performance. Automated Raman spectroscopy was used for solute quantitation of the ternary system tert‐butanol in a toluene/water matrix, a ternary system that was specifically chosen to analyze the device's performance with a small solute partition ratio and to enable in‐line Raman measurements of solute concentrations in both phases. The microfluidic device possessed a 55 mm contact length and a 38.5 µL internal volume. During counter‐current flow, we observed approximately 37 equilibrium stages (37 ± 13) based on a best‐fit of the solute fraction remaining in the aqueous phase using a Kremser Group Method analysis.