Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Preconditioning”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

Multigrid Solution of the Navier-Stokes Equations at Low Speeds with Large Temperature Variations

Multigrid methods for the Navier-Stokes equations at low speeds and large temperature variations are investigated. The compressible equations with time-derivative preconditioning and preconditioned flux-difference splitting of the inviscid terms are used. Three implicit smoothers have been incorporated into a common multigrid procedure. Both full coarsening and semi-coarsening with directional fine-grid defect correction have been studied. The resulting methods have been tested on four 2D laminar problems over a range of Reynolds numbers on both uniform and highly stretched grids. Two of the three methods show efficient and robust performance over the entire range of conditions. In addition none of the methods have any difficulty with the large temperature variations.

Sockol, Peter M.

Algorithmic Enhancements to the VULCAN Navier-Stokes Solver

VULCAN (Viscous Upwind aLgorithm for Complex flow ANalysis) is a cell centered, finite volume code used to solve high speed flows related to hypersonic vehicles. Two algorithms are presented for expanding the range of applications of the current Navier-Stokes solver implemented in VULCAN. The first addition is a highly implicit approach that uses subiterations to enhance block to block connectivity between adjacent subdomains. The addition of this scheme allows more efficient solution of viscous flows on highly-stretched meshes. The second algorithm addresses the shortcomings associated with density-based schemes by the addition of a time-derivative preconditioning strategy. High speed, compressible flows are typically solved with density based schemes, which show a high level of degradation in accuracy and convergence at low Mach numbers (M less than or equal to 0.1). With the addition of preconditioning and associated modifications to the numerical discretization scheme, the eigenvalues will scale with the local velocity, and the above problems will be eliminated. With these additions, VULCAN now has improved convergence behavior for multi-block, highly-stretched meshes and also can solve the Navier-Stokes equations for very low Mach numbers.

Litton, D. K.

Increased Agrobacterium-mediated transformation and rooting efficiencies in canola (Brassica napus L.) from hypocotyl segment explants

An efficient protocol for the production of transgenic Brassica napus cv. Westar plants was developed by optimizing two important parameters: preconditioning time and co-cultivation time. Agrobacterium tumefaciens-mediated transformation was performed using hypocotyls as explant tissue. Two variants of a green fluorescent protein (GFP)-encoding gene--mGFP5-ER and eGFP--both under the constitutive expression of the cauliflower mosaic virus 35S promoter, were used for the experiments. Optimizing the preconditioning time to 72 h and co-cultivation time with Agrobacterium to 48 h provided the increase in the transformation efficiency from a baseline of 4% to 25%. With mGFP5-ER, the transformation rate was 17% and with eGFP it was 25%. Transgenic shoots were selected on 200 mg/l kanamycin. Rooting efficiency was 100% on half-strength Murashige and Skoog medium with 10 g/l sucrose and 0.5 mg/l indole butyric acid in the presence of kanamycin.

NASA Discipline Life Sciences Technologies

Chronic shear induces caveolae formation and alters ERK and Akt responses in endothelial cells

Caveolae are plasmalemmal domains enriched with cholesterol, caveolins, and signaling molecules. Endothelial cells in vivo are continuously exposed to shear conditions, and their caveolae density and location may be different from that of static cultured cells. Here, we show that chronic shear exposure regulates formation and localization of caveolae and caveolin-1 in bovine aortic endothelial cells (BAEC). Chronic exposure (1 or 3 days) of BAEC to laminar shear increased the total number of caveolae by 45-48% above static control. This increase was due to a rise in the luminal caveolae density without changing abluminal caveolae numbers or increasing caveolin-1 mRNA and protein levels. Whereas some caveolin-1 was found in the plasma membrane in static-cultured cells, it was predominantly localized in the Golgi. In contrast, chronic shear-exposed cells showed intense caveolin-1 staining in the luminal plasma membrane with minimum Golgi association. The preferential luminal localization of caveolae may play an important role in endothelial mechanosensing. Indeed, we found that chronic shear exposure (preconditioning) altered activation patterns of two well-known shear-sensitive signaling molecules (ERK and Akt) in response to a step increase in shear stress. ERK activation was blunted in shear preconditioned cells, whereas the Akt response was accelerated. These results suggest that chronic shear stimulates caveolae formation by translocating caveolin-1 from the Golgi to the luminal plasma membrane and alters cell signaling responses.

NASA Discipline Cell Biology

Thigmomorphogenesis: the role of ethylene in the response of Pinus taeda and Abies fraseri to mechanical perturbation

Ethylene production was monitored for 48 h in two half-sibs of Pinus taeda L. grown in the greenhouse and given mechanical perturbation (MP) by flexing; and for 22 h in Abies fraseri (Pursh) Poir. grown in the field and exposed to wind-mediated MP. Both species produced a peak of ethylene 18 h after MP. Seedlings of P. taeda exposed to MP for the duration of the growing season (preconditioned) produced less ethylene compared to non-MP controls, with a peak production at 8 h. One half-sib which responded to MP by an increase in radial growth produced 16 times more ethylene than another half-sib which had no significant change in radial growth. Preconditioned A. fraseri produced no significant quantities of ethylene after MP. The production of wound ethylene appears to be different from MP-induced ethylene. When an ethylene-generating solution was applied to P. taeda seedlings, it mimicked many of the morphological and mechanical characteristics of MP seedlings. The putative role of ethylene in the thigmomorphogenetic response is addressed.

Non-NASA Center

An Upwind Solver for the National Combustion Code

An upwind solver is presented for the unstructured grid National Combustion Code (NCC). The compressible Navier-Stokes equations with time-derivative preconditioning and preconditioned flux-difference splitting of the inviscid terms are used. First order derivatives are computed on cell faces and used to evaluate the shear stresses and heat fluxes. A new flux limiter uses these same first order derivatives in the evaluation of left and right states used in the flux-difference splitting. The k-epsilon turbulence equations are solved with the same second-order method. The new solver has been installed in a recent version of NCC and the resulting code has been tested successfully in 2D on two laminar cases with known solutions and one turbulent case with experimental data.

Sockol, Peter M.

Method to Enhance the Operation of an Optical Inspection Instrument Using Spatial Light Modulators

For many aspheric and freeform optical components, existing interferometric solutions require a custom computer-generated hologram (CGH) to characterize the part. The overall objective of this research is to develop hardware and a procedure to produce a combined, dynamic, Hartmann/ Digital Holographic interferometry inspection system for a wide range of advanced optical components, including aspheric and freeform optics. This new instrument would have greater versatility and dynamic range than currently available measurement systems. The method uses a spatial light modulator to pre-condition wavefronts for imaging, interferometry, and data processing to improve the resolution and versatility of an optical inspection instrument. Existing interferometers and Hartmann inspection systems have either too small a dynamic range or insufficient resolution to characterize conveniently unusual optical surfaces like aspherical and freeform optics. For interferometers, a specially produced, computer-generated holographic optical element is needed to transform the wavefront to within the range of the interferometer. A new hybrid wavefront sensor employs newly available spatial light modulators (SLMs) as programmable holographic optical elements (HOEs). The HOE is programmed to enable the same instrument to inspect an optical element in stages, first by a Hartmann measurement, which has a very large dynamic range but less resolution. The first measurement provides the information required to precondition a reference wave that avails the measurement process to the more precise phase shifting interferometry. The SLM preconditions a wavefront before it is used to inspect an optical component. This adds important features to an optical inspection system, enabling not just wavefront conditioning for null testing and dynamic range extension, but also the creation of hybrid measurement procedures. This, for example, allows the combination of dynamic digital holography and Hartmann sensing procedures to cover a virtually unlimited dynamic range with high resolution. Digital holography technology brings all of the power and benefits of digital holographic interferometry to the requirement, while Hartmann-type wavefront sensors bring deflectometry technologies to the solution. The SLM can be used to generate arbitrary wavefronts in one leg of the interferometer, thereby greatly simplifying its use and extending its range. The SLM can also be used to modify the system into a dynamic Shack-Hartmann system, which is useful for optical components with large amounts of slope. By integrating these capabilities into a single instrument, the system will have tremendous flexibility to measure a variety of optical shapes accurately.

Trolinger, James

Simulations of Continuous Descent Operations with Arrival-management Automation and Mixed Flight-deck Interval Management Equipage

Air traffic management simulations conducted in the Airspace Operations Laboratory at NASA Ames Research Center have addressed the integration of trajectory-based arrival-management automation, controller tools, and Flight-Deck Interval Management avionics to enable Continuous Descent Operations (CDOs) during periods of sustained high traffic demand. The simulations are devoted to maturing the integrated system for field demonstration, and refining the controller tools, clearance phraseology, and procedures specified in the associated concept of operations. The results indicate a variety of factors impact the concept's safety and viability from a controller's perspective, including en-route preconditioning of arrival flows, useable clearance phraseology, and the characteristics of airspace, routes, and traffic-management methods in use at a particular site. Clear understanding of automation behavior and required shifts in roles and responsibilities is important for controller acceptance and realizing potential benefits. This paper discusses the simulations, drawing parallels with results from related European efforts. The most recent study found en-route controllers can effectively precondition arrival flows, which significantly improved route conformance during CDOs. Controllers found the tools acceptable, in line with previous studies.

arrival management

Thermodynamics in the Suppressed Phase of the Madden-Julian Oscillation Using a Multiplatform Strategy

The Madden-Julian Oscillation (MJO) represents a prominent mode of intraseasonal tropical variability. It is manifest by coherent large-scale changes in atmospheric circulation, convection, and thermodynamic processes. Preconditioning of the environment prior to the active phase of the MJO has been noted, but the balance of theorized mechanisms to accomplish this process remains unresolved. Further, there is a lack of consensus on the means by which primary initiation of an MJO event occurs. Observational and modeling efforts have recently been undertaken to advance our understanding of the physical underpinnings governing MJO development. However these intensive studies are often limited in space and/or time and are potentially subject to model deficiencies. Satellite observations, especially those providing vertical resolution of temperature and moisture, provide an opportunity to expand our knowledge of processes critical to MJO initiation and preconditioning. This work will provide an analysis of suppressed phase thermodynamics with an emphasis on the use of a complementary suite of satellite observations including AIRS/AMSU-A profiles, CERES radiative fluxes, and cloud properties observed by MODIS. Emphasis of this work will regard the distribution of cloud regimes, their radiative-convective effects, and their relationship to moist static energy during the recharge and suppressed stages of MJO initiation and eastward propagation. The analyses will make use of cloud regimes from MODIS observations to provide a compositing technique that enables the identification of systematic connections between different cloud regimes and the larger scale environment. Within these cloud regimes, the relationship between the associated cloud-radiative effects observed by CERES, vertically-resolved and vertically-integrated thermodynamics using AIRS/AMSU-A observations, and atmospheric boundary layer fluxes will be demonstrated.

Roberts, J. Brent

The Scintillation Prediction Observations Research Task (SPORT): an International Science Mission Using a Cubesat

The Scintillation Prediction Observations Research Task (SPORT) is a 6U CubeSat mission to address the compelling but difficult problem of understanding the preconditions leading to equatorial plasma bubbles. The scientific literature describes the preconditions in both the plasma drifts and the density profiles related to bubble formations that occur several hours later in the evening. Most of the scientific discovery has resulted from observations at a single site, within a single longitude sector, from Jicamarca, Peru. SPORT will provide a systematic study of the state of the pre-bubble conditions at all longitudes sectors to enhance understanding between geography and magnetic geometry. SPORT is an international partnership between National Aeronautics and Space Administration (NASA), the Brazilian National Institute for Space Research (INPE), and the Technical Aeronautics Institute under the Brazilian Air Force Command Department (DCTA/ITA), and encouraged by U.S. Southern Command. This talk will present an overview of the SPORT mission, observation strategy, and science objectives to improve predictions of ionospheric disturbances that affect radio propagation of telecommunication signals. The science goals will be accomplished by a unique combination of satellite observations from a nearly circular middle inclination orbit and the extensive operation of ground based observations from South America near the magnetic equator.

Scintillation Prediction Observations

Effect of Environments on Parametric Degradation in Polymer Tantalum Capacitors

Chip polymer tantalum capacitors (CPTCs) are more sensitive to environments compared to conventional, MnO2 cathode capacitors. Contrary to MnO2 capacitors, CPTCs can degrade with time of storage even in dry conditions, that might limit their use for space applications. Performance of CPTCs depends on the amount of moisture in the slug, and varies with time of storage after manufacturing. For this reason, characteristics of the parts and qualification testing should be carried out after preconditioning to stabilize moisture content. This work analyzes kinetics of moisture sorption and desorption in different types of CPTCs from three manufacturers. Variations of capacitance, dissipation factor, ESR, and leakage currents caused by exposure to humid, dry, and radiation environments are discussed. Examples of the effects of long-term exposure to vacuum and humid environments at 85%RH are presented. Temperature dependencies of the characteristic times of moisture diffusion allow for assessments of the bake-out times that are necessary to avoid pop-corning failures after soldering and to precondition parts before qualification testing.

Teverovsky, Alexander

Development of Segregated Thermal-Hydraulics Solvers in MOOSE

The simulation of fluid flows is an essential part of the design and analysis of nuclear systems. Algorithms able to simulate flows at different fidelity levels are available in the Multiphysics Object-Oriented Simulation Environment (MOOSE) and MOOSE-based applications such as Pronghorn \cite{novak2018pronghorn}, Pronghorn-Subchannel, RELAP-7, and SAM. Currently, significant effort is being invested in the development of coarse-mesh Computational Fluid Dynamics (CFD) capabilities within MOOSE and Pronghorn for the simulation of Generation IV nuclear reactors. Traditionally, the solution algorithms in MOOSE have relied on Newton or quasi-Newton methods (such as the preconditioned Jacobian-free Newton-Krylov method) where residuals and Jacobians (or approximations thereof) are constructed. Both Newton and quasi-Newton methods require the solution of a linear system at each nonlinear Newton iteration with the Jacobian as the system matrix. The Jacobian contains blocks originating from all variables in the problem (i.e., for thermal-hydraulics at least pressure, velocities, and temperature). Due to the formulation of the problem in a general multiphysics setting on unstructured mesh, creating a good preconditioner for the linear system can be challenging, thus many fluid applications have utilized direct solver-based methods such as LU factorization. However, with increasing system size and complexity in multi-dimensional problems, the direct solution of linear systems becomes computationally expensive both in execution time and and memory. For this reason, recent effort has focused on adapting segregated solution algorithms for CFD problems in MOOSE. These algorithms use fixed-point iteration between segregated systems whose assembly and preconditioning are easier those of the monolithic system. Initial results show that the segregated solution algorithm outperforms the monolithic approach in terms of memory usage and for large 3D problems in terms of CPU time as well.

42 ENGINEERING

Development of Segregated Thermal-Hydraulics Solvers in MOOSE

The simulation of fluid flows is an essential part of the design and analysis of nuclear systems. Algorithms able to simulate flows at different fidelity levels are available in the Multiphysics Object-Oriented Simulation Environment (MOOSE) and MOOSE-based applications such as Pronghorn \cite{novak2018pronghorn}, Pronghorn-Subchannel, RELAP-7, and SAM. Currently, significant effort is being invested in the development of coarse-mesh Computational Fluid Dynamics (CFD) capabilities within MOOSE and Pronghorn for the simulation of Generation IV nuclear reactors. Traditionally, the solution algorithms in MOOSE have relied on Newton or quasi-Newton methods (such as the preconditioned Jacobian-free Newton-Krylov method) where residuals and Jacobians (or approximations thereof) are constructed. Both Newton and quasi-Newton methods require the solution of a linear system at each nonlinear Newton iteration with the Jacobian as the system matrix. The Jacobian contains blocks originating from all variables in the problem (i.e., for thermal-hydraulics at least pressure, velocities, and temperature). Due to the formulation of the problem in a general multiphysics setting on unstructured mesh, creating a good preconditioner for the linear system can be challenging, thus many fluid applications have utilized direct solver-based methods such as LU factorization. However, with increasing system size and complexity in multi-dimensional problems, the direct solution of linear systems becomes computationally expensive both in execution time and and memory. For this reason, recent effort has focused on adapting segregated solution algorithms for CFD problems in MOOSE. These algorithms use fixed-point iteration between segregated systems whose assembly and preconditioning are easier those of the monolithic system. Initial results show that the segregated solution algorithm outperforms the monolithic approach in terms of memory usage and for large 3D problems in terms of CPU time as well.

42 ENGINEERING

A fully implicit, asymptotic-preserving, semi-Lagrangian algorithm for the time dependent anisotropic heat transport equation

In this paper, we extend the operator-split asymptotic-preserving, semi-Lagrangian algorithm for time dependent anisotropic heat transport equation proposed in Chacón et al. (2014) [18] to use a fully implicit time integration with backward differentiation formulas. The proposed implicit method can deal with arbitrary heat-transport anisotropy ratios $\mathcal{X}$∥ /$ \mathcal{X}$⟂ $\ggg$ 1 (with $\mathcal{X}$∥, $ \mathcal{X}$⟂ the parallel and perpendicular heat diffusivities, respectively) in complicated magnetic field topologies in an accurate and efficient manner. Further, the implicit algorithm is second-order accurate temporally and demonstrates an accurate treatment at boundary layers (e.g., island separatrices), which was not ensured by the operator-split implementation. The condition number of the resulting algebraic system is independent of the anisotropy ratio, and is inverted with preconditioned GMRES. We propose a simple preconditioner that renders the finite-dimensional linear operator compact, resulting in mesh-independent convergence rates for topologically simple magnetic fields, and convergence rates scaling as ~ (NΔt) 1/4 (with N the total mesh size and Δt the timestep) in topologically complex magnetic-field configurations. We demonstrate the accuracy and performance of the approach with test problems of varying complexity, including an analytically tractable boundary-layer problem in a straight magnetic field, and a topologically complex magnetic field featuring magnetic islands with extreme anisotropy ratios $\mathcal{X}$∥ /$ \mathcal{X}$⟂ = 10 10 ) .

97 MATHEMATICS AND COMPUTING

A scalable multidimensional fully implicit solver for Hall magnetohydrodynamics

We propose an optimally performant fully implicit algorithm for the Hall magnetohydrodynamics (HMHD) equations based on multigrid-preconditioned Jacobian-free Newton-Krylov methods. HMHD is a challenging system to solve numerically because it supports stiff fast dispersive waves. The preconditioner is formulated using an operator-split approximate block factorization (Schur complement), informed by physics insight. We use a vector-potential formulation (instead of a magnetic field one) to allow a clean segregation of the problematic $\nabla$ x $\nabla$ x operator in the electron Ohm's law subsystem. This segregation allows the formulation of an effective damped block-Jacobi smoother for multigrid. We demonstrate by analysis that our proposed block-Jacobi iteration is convergent and has the smoothing property. The resulting HMHD solver is verified linearly with wave propagation examples, and nonlinearly with the GEM challenge reconnection problem by comparison against another HMHD code. We demonstrate the excellent algorithmic and parallel performance of the algorithm up to 16384 MPI tasks in two dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC

Unorthodox parallelization for Bayesian quantum state estimation

Quantum state tomography (QST) allows for the reconstruction of quantum states through measurements and some inference technique under the assumption of repeated state preparations. Bayesian inference provides a promising platform to achieve both efficient QST and accurate uncertainty quantification, yet is generally plagued by the computational limitations associated with long Markov chains. In this work, we present a novel Bayesian QST approach that leverages modern distributed parallel computer architectures to efficiently sample a D-dimensional Hilbert space. Using a parallelized preconditioned Crank–Nicholson Metropolis–Hastings algorithm, we demonstrate our approach on simulated data and experimental results from IBM Quantum systems up to four qubits, showing significant speedups through parallelization. Although highly unorthodox in pooling independent Markov chains, our method proves remarkably practical, with validation ex post facto via diagnostics like the intrachain autocorrelation time. We conclude by discussing scalability to higher-dimensional systems, offering a path toward efficient and accurate Bayesian characterization of large quantum systems.

Bayesian inference

A two-level GPU-accelerated incomplete LU preconditioner for general sparse linear systems

This paper presents a parallel preconditioning approach based on incomplete LU (ILU) factorizations in the framework of Domain Decomposition (DD) for general sparse linear systems. We focus on distributed memory parallel architectures, specifically, those that are equipped with graphic processing units (GPUs). In addition to block-Jacobi, we present general purpose two-level ILU Schur complement-based approaches, where different strategies are presented to solve the coarse-level reduced system. These strategies are combined with modified ILU methods in the construction of the coarse-level operator, in order to effectively remove smooth errors by targeting an algebraically smooth vector. We leverage available GPU-based sparse matrix kernels to accelerate the setup and the solve phases of the proposed ILU preconditioner. We evaluate the efficiency of the proposed methods as a smoother for algebraic multigrid (AMG) and as a preconditioner for Krylov subspace methods on challenging anisotropic diffusion problems and a collection of general sparse matrices.

97 MATHEMATICS AND COMPUTING

Speedup of UEDGE Parameter Scans Using Machine-Learning Optimized OpenMP Parallelization and a Continuation Solver

This article presents the OpenMP parallelization of the preconditioning Jacobian assembly and right‐hand side residual evaluation in UEDGE. A continuation algorithm, utilizing the internal NKSOL implicit Jacobian‐Free Newton‐Krylov solver to efficiently scan physical parameters, is also presented. The implemented parallelization reduces the computational time for a benchmark scan run on 32 threads by compared to the serial version when using trained random forest regression models to identify the optimal decomposition of the system of equations. Random forest regression models applied to the UEDGE time‐dependent and continuation solver algorithms did not yield meaningful improvement in computational performance. A benchmark DIII‐D gas injection rate scan in the 0.35–0.75 kA interval, performed on a test cluster using the parallelized code and continuation solver, produced 1066 steady‐state solutions with a 22 s average wall‐clock computational time per steady‐state solution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY