Search NASA⌕ Search

SEARCH · Search NASA

Results for “Conformable gridding”

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

Highly parallel sparse Cholesky factorization

Several fine grained parallel algorithms were developed and compared to compute the Cholesky factorization of a sparse matrix. The experimental implementations are on the Connection Machine, a distributed memory SIMD machine whose programming model conceptually supplies one processor per data element. In contrast to special purpose algorithms in which the matrix structure conforms to the connection structure of the machine, the focus is on matrices with arbitrary sparsity structure. The most promising algorithm is one whose inner loop performs several dense factorizations simultaneously on a 2-D grid of processors. Virtually any massively parallel dense factorization algorithm can be used as the key subroutine. The sparse code attains execution rates comparable to those of the dense subroutine. Although at present architectural limitations prevent the dense factorization from realizing its potential efficiency, it is concluded that a regular data parallel architecture can be used efficiently to solve arbitrarily structured sparse problems. A performance model is also presented and it is used to analyze the algorithms.

Gilbert, John R.↗

Prediction of unsteady aerodynamic loads in cascades using the linearized Euler equations on deforming grids

A linearized Euler solver for calculating unsteady flows in turbomachinery blade rows due to both incident gusts and blade motion is presented. Using the linearized Euler technique, one decomposes the flow into a mean (or steady) flow plus an unsteady, harmonically varying, small disturbance flow. Linear variable coefficient equations describe the small disturbance behavior of the flow, and are solved using a pseudo-time marching Lax-Wendroff scheme. For the blade motion problem, a harmonically deforming computational rid that conforms to the motion of vibrating blades eliminates large error producing mean flow gradient terms that would otherwise appear in the unsteady flow tangency boundary condition. The paper also presents a new, numerically exact, nonreflecting far-field boundary condition based on an eigenanalysis of the discretized equations. Computed flow solutions demonstrate the computational accuracy and efficiency of the present method. The solution of the linearized Euler equations requires one to two orders of magnitude less computer time than solution of the nonlinear Euler equations using traditional time-accurate time-marching techniques. In addition, the deformable grid significantly improves the accuracy of the solution.

Hall, Kenneth C.↗

Parallel 3D Mortar Element Method for Adaptive Nonconforming Meshes

High order methods are frequently used in computational simulation for their high accuracy. An efficient way to avoid unnecessary computation in smooth regions of the solution is to use adaptive meshes which employ fine grids only in areas where they are needed. Nonconforming spectral elements allow the grid to be flexibly adjusted to satisfy the computational accuracy requirements. The method is suitable for computational simulations of unsteady problems with very disparate length scales or unsteady moving features, such as heat transfer, fluid dynamics or flame combustion. In this work, we select the Mark Element Method (MEM) to handle the non-conforming interfaces between elements. A new technique is introduced to efficiently implement MEM in 3-D nonconforming meshes. By introducing an "intermediate mortar", the proposed method decomposes the projection between 3-D elements and mortars into two steps. In each step, projection matrices derived in 2-D are used. The two-step method avoids explicitly forming/deriving large projection matrices for 3-D meshes, and also helps to simplify the implementation. This new technique can be used for both h- and p-type adaptation. This method is applied to an unsteady 3-D moving heat source problem. With our new MEM implementation, mesh adaptation is able to efficiently refine the grid near the heat source and coarsen the grid once the heat source passes. The savings in computational work resulting from the dynamic mesh adaptation is demonstrated by the reduction of the the number of elements used and CPU time spent. MEM and mesh adaptation, respectively, bring irregularity and dynamics to the computer memory access pattern. Hence, they provide a good way to gauge the performance of computer systems when running scientific applications whose memory access patterns are irregular and unpredictable. We select a 3-D moving heat source problem as the Unstructured Adaptive (UA) grid benchmark, a new component of the NAS Parallel Benchmarks (NPB). In this paper, we present some interesting performance results of ow OpenMP parallel implementation on different architectures such as the SGI Origin2000, SGI Altix, and Cray MTA-2.

Feng, Huiyu↗

Boundary-conforming mapping applied to computations of highly deformed solidification interfaces

A new boundary-conforming mapping is developed for the calculation of highly deformed cellular solidification interfaces in a model of directional solidification of a binary alloy. The mapping is derived through a variational formulation that is designed so that the grid penetrates the grooves between cells along the interface without causing a loss of ellipticity of the mapping equations. A finite element/Newton method is presented for simultaneous solution of the free boundary problem described by the solutal model of directional solidification and the mapping equations. Results are compared to previous calculations and demonstrate the importance of accurate representation of the interface shape for understanding the solution structure.

Tsiveriotis, K.↗

Numerical procedure for three-dimensional hypersonic viscous flow over aerobrake configuration

A numerical method, which is simpler and more efficient than others currently in use, is proposed for the computation of the full viscous flow over an aerobrake body in hypersonic stream at high altitude. It treats the shock layer surrounding the blunt forebody and the near wake behind the base simultaneously by formulating the Navier-Stokes equations in conformal and azimuthal-angle coordinates. The computational domain is confined by the body wall, outflow surface and the shock, which is adjusted along the coordinate normal to the wall in the course of iterations. Because of the optimal grid and a well developed alternating direction implicit factorization technique for the governing equations, reasonably accurate results can be obtained with a 28 x 36 x 7 grid and 400 time-marching iterations. Excellent agreement of shock location is found between the present result and the schlieren photograph. Details of the base flow and shear layer impingement on the cylindrical aft body are presented for an adiabatic wall case.

Li, C. P.↗

Computational methods for hypersonic viscous flow over finite ellipsoid-cones at incidence

A numerical method, which is simpler than others currently in use, is proposed for determining the full viscous flow over a finite body in hypersonic stream at high altitude. It treats the shock layer surrounding the blunt foebody and the near wake behind the base simultaneously by formulating the Navier-Stokes equations in conformal and azimuthal-angle coordinates. The computational domain is confined to the body wall, outflow surface and the bow shock, which is adjusted along the coordinate normal to the wall in the course of iterations. Because of the optimal grid and a well developed alternating direction implicit factorization technique for the governing equations, reasonably accurate results can be obtained on a 30 by 36 by 6 grid with 400 time-marching iterations. Results for body shapes belonging to the ellipsoid-cone family are compared with the experimental data for the Apollo command module and the Viking aeroshell. Validation of the method based on self-consistency is also discussed.

Li, C. P.↗

Harnessing the power of ab initio calculations, distributed computing and machine learning to efficiently locate extreme molecules for use in carbon-based solar cells (Final Technical Report)

The use of high-throughput virtual screening (HTVS) tools is a powerful tool to expedite the materials discovery of commercially relevant materials. In previous years, our group has developed a molecular discovery platform to generate libraries in order to obtain suitable candidates for different applications, starting from the Harvard Clean Energy Project. This platform is suitable to test in-silico on traditional supercomputing clusters and shared resources, for example, in the IBM World Community Grid. In this project, we used the molecular discovery platform to create and screen a library of candidates of organic photovoltaic (OPVs) molecules. Based on a set of candidates created with combinations of molecular moieties, we were able to filter, by conformation stability, the energy of electronic orbitals and approximated power conversion efficiencies (PCE). To improve the predictions of orbital energies calculated and the PCEs, we used Gaussian Process regression and two sets of molecules. These sets correspond to electronic structure calculations of a higher level of theory and experimental PCE values, respectively. Finally, we selected a subset of the best candidates (molecules with a PCE higher than 10%) to understand its absorbance properties with TD-DFT. This project has demonstrated the capabilities of our molecular discovery platform for HTVS. Finally, machine learning can help us to introduce more complex effects included in bulk conditions and computational intensive calculations on models.

14 SOLAR ENERGY↗

Towards an Entropy Stable Spectral Element Framework for Computational Fluid Dynamics

Entropy stable (SS) discontinuous spectral collocation formulations of any order are developed for the compressible Navier-Stokes equations on hexahedral elements. Recent progress on two complementary efforts is presented. The first effort is a generalization of previous SS spectral collocation work to extend the applicable set of points from tensor product, Legendre-Gauss-Lobatto (LGL) to tensor product Legendre-Gauss (LG) points. The LG and LGL point formulations are compared on a series of test problems. Although being more costly to implement, it is shown that the LG operators are significantly more accurate on comparable grids. Both the LGL and LG operators are of comparable efficiency and robustness, as is demonstrated using test problems for which conventional FEM techniques suffer instability. The second effort generalizes previous SS work to include the possibility of p-refinement at non-conforming interfaces. A generalization of existing entropy stability machinery is developed to accommodate the nuances of fully multi-dimensional summation-by-parts (SBP) operators. The entropy stability of the compressible Euler equations on non-conforming interfaces is demonstrated using the newly developed LG operators and multi-dimensional interface interpolation operators.

Carpenter, Mark H.↗

UAV Inspection of Electrical Transmission Infrastructure with Path Conformance Autonomy and Lidar-Based Geofences NASA Report on UTM Reference Mission Flights at Southern Company Flights November 2016

Flights at low altitudes in close proximity to electrical transmission infrastructure present serious navigational challenges: GPS and radio communication quality is variable and yet tight position control is needed to measure defects while avoiding collisions with ground structures. To advance unmanned aerial vehicle (UAV) navigation technology while accomplishing a task with economic and societal benefit, a high voltage electrical infrastructure inspection reference mission was designed. An integrated air-ground platform was developed for this mission and tested in two days of experimental flights to determine whether navigational augmentation was needed to successfully conduct a controlled inspection experiment. The airborne component of the platform was a multirotor UAV built from commercial off-the-shelf hardware and software, and the ground component was a commercial laptop running open source software. A compact ultraviolet sensor mounted on the UAV can locate 'hot spots' (potential failure points in the electric grid), so long as the UAV flight path adequately samples the airspace near the power grid structures. To improve navigation, the platform was supplemented with two navigation technologies: lidar-to-polyhedron preflight processing for obstacle demarcation and inspection distance planning, and trajectory management software to enforce inspection standoff distance. Both navigation technologies were essential to obtaining useful results from the hot spot sensor in this obstacle-rich, low-altitude airspace. Because the electrical grid extends into crowded airspaces, the UAV position was tracked with NASA unmanned aerial system traffic management (UTM) technology. The following results were obtained: (1) Inspection of high-voltage electrical transmission infrastructure to locate 'hot spots' of ultraviolet emission requires navigation methods that are not broadly available and are not needed at higher altitude flights above ground structures. (2) The sensing capability of a novel airborne UV detector was verified with a standard ground-based instrument. Flights with this sensor showed that UAV measurement operations and recording methods are viable. With improved sensor range, UAVs equipped with compact UV sensors could serve as the detection elements in a self-diagnosing power grid. (3) Simplification of rich lidar maps to polyhedral obstacle maps reduces data volume by orders of magnitude, so that computation with the resultant maps in real time is possible. This enables real-time obstacle avoidance autonomy. Stable navigation may be feasible in the GPS-deprived environment near transmission lines by a UAV that senses ground structures and compares them to these simplified maps. (4) A new, formally verified path conformance software system that runs onboard a UAV was demonstrated in flight for the first time. It successfully maneuvered the aircraft after a sudden lateral perturbation that models a gust of wind, and processed lidar-derived polyhedral obstacle maps in real time. (5) Tracking of the UAV in the national airspace using the NASA UTM technology was a key safety component of this reference mission, since the flights were conducted beneath the landing approach to a heavily used runway. Comparison to autopilot tracking showed that UTM tracking accurately records the UAV position throughout the flight path.

Moore, Andrew J.↗

Emerging Pathways to Upgrade the US Housing Stock: A Review of the Home Energy Upgrade Literature

The residential buildings sector is responsible for about 20% of total US energy use. In order to achieve climate goals, we need ways to reduce carbon emissions and energy use in this sector. In addition, resiliency, electric grid stability, emergency survivability and other energy and building-related issues are becoming increasingly important challenges. New homes in most of the US meet various energy codes and are reasonably energy efficient. However, the vast majority of energy use is from existing homes that were not required to conform to energy performance requirements. It is becoming imperative to reach as many of these existing homes as possible and find ways to improve their energy-related performance. This must be done in such a way that it meets the needs and desires of homeowners and building occupants, as well as those of the contractors and design professionals engaged in doing the upgrades themselves. Energy retrofits of homes started in the 1970’s in response to the energy crisis, however, these retrofits were very limited in scope and relatively few homes were upgraded. Those homes that have been upgraded generally still have much scope for improvement. A huge effort is needed to get to scale to address the energy use in housing. The target population is effectively every home in the country, whether a large suburban single-family home, or a small downtown apartment. In order to provide a framework for analysis and the basis for plans to get to large-scale retrofits of homes, this literature review summarizes the state-of-the art in the US buildings industry. It identifies where more research, engineering, or technology is needed, as well as relevant industry trends, such as electrification, one-stop shop program design and others. It also examines other key topics, such as availability of financing, minimizing household disruption, and engaging home owners and occupants. This literature review builds on a similar review from several years ago (Less and Walker, 2014). The current review focuses on efforts in the intervening years. This literature review is part of a larger DOE study of deep energy upgrades that includes industry surveys and development of cost-stack analyses. For this review, we gathered data not just from the published literature, but also from practitioners in conjunction with other aspects of the larger DOE study. In some cases, we refer to comments from specific individuals or companies, or refer to specific products by name. This is not intended as an endorsement, but rather to provide clarity on sources of information and examples of relevant technologies.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

The use of back corrections in multistep methods of numerical integration

Generalized multistep methods for the numerical solution of nonlinear systems of ordinary differential equations are introduced which allow the correction of previously computed values of the solution at internal points of the grid. These methods are shown to possess enhanced numerical stability. Preliminary numerical results indicate that for some satellite orbit problems these methods also possess greater overall efficiency. A uniformly converging theory of error propagation is presented which is valid for nonasymptotic values of the step size. Experimental results are seen to conform with theory.

Beaudet, P.↗

Algebraic grid generation with boundary orthogonality control

This paper presents a new approach in applying the blending function method to the generation of structured finite-difference meshes. It is shown that the surface equation developed by Coons (1967) can be easily adapted to generate two-dimensional meshes, and its extension to hyperspace allows the construction of three-dimensional meshes. The robust method in its basic form blends among the boundary curves/surfaces conforming to their intrinsic slopes. When parameter-controlled boundary-slope corrections are introduced, it is possible to establish varying degrees of orthogonality along the boundaries. Additional controls are made possible by moving points in the parameter space, which affects the movement of points in the physical space. The combination of these controls provides a simple and yet powerful means for exploiting the blending-function method to satisfy the needs of practical CFD applications. Results from a validation program for two-dimensions are shown for a variety of geometries.

Luh, Raymond Ching-Chung↗

Analysis of Fibrous Felts for Flexible Ablators Using Synchrotron Hard X-Ray Micro-Tomography

We analyzed the material properties of low-density felts that are used as substrates for new-generation flexible and conformal carbon/phenolic ablators, and compared them with those of a rigid carbon fiber preform that is used to manufacture rigid carbon/phenolic ablators. Micro-tomography measurements were obtained using synchrotron X-rays, allowing the characterization of the materials microstructure at the scale of the fibers. Using the tomography voxels as computational grids, we computed tortuosity and room temperature conductivity. In addition we performed micro-scale simulations of the oxidation of carbon fibers using a random walk model for oxygen diffusion and a probability law to model surface reactions.

Ablation↗

A numerical analysis of tonal acoustics in rotor-stator interactions

In this study, the unsteady, thin-layer, Navier-Stokes equations are solved using a system of patched grids for a rotor-stator configuration of an axial turbine. The study focuses on the plurality of spinning modes that are present in such an interaction. The propagation of these modes in the upstream and downstream regions is analyzed and compared with numerical results. It was found that the numerically calculated tonal acoustics could be affected by the type of numerical boundary conditions employed at the inlet and exit of the computational boundaries and the grid spacing in the upstream and downstream regions. Results in the form of surface pressure amplitudes and the spectra of turbine tones and their far field behavior are presented. Numerical results and experimental data are compared wherever possible. The 'mode-content' for different harmonics of blade-passage frequency is shown to conform with that predicted by a kinematical analysis.

Rangwalla, A. A.↗

Computation of three-dimensional flow about aerobrake configurations

Ellipsoid, cone and cylinder aerobrake configurations are analyzed to provide comparison data between experimental and model predictions. An analytical model was devised to account for the shock layer ahead of the body and in the near-wake region in terms of the Navier-Stokes equations expressed in conformal polar and azimuthal-angle coordinates. Using polar coordinates simplified the equations by mapping the body onto a sphere, a procedure which also reduced the magnitude of the discretization errors. The equations are then solved using an alternating direction implicit (ADI) factorization technique. Computations were carried out for Mach 3-10 at various grid resolutions and compared with available wind tunnel data. The model generated pressure distributions, heat transfer coefficients and velocity profile data that agreed relatively well with experimental data at a reduced computational cost. Further work is necessary to identify the location of shocks and to model flows about asymmetric configurations.

Li, C. P.↗

Critical role of structural order in bipolar redox-active molecules for organic redox flow batteries

Bipolar redox-active molecules (BRMs) have been suggested as a means to address crossover related issues in all-organic redox flow batteries (RFBs). In such species, electron donors (anolytes) and electron acceptors (catholytes) are linked by a chain, and the same chemical composition is used in both anode and cathode compartments to reduce chemical gradients and unbalanced ion transport. The resultant RFBs resemble the aqueous vanadium RFBs, but offer greater design flexibility and potentially more favorable electrochemical and physicochemical properties. Yet the complex trade-offs in these properties have been a complication. Here we developed a fundamental rubric to uncover the likely origins in the performance metrics of these BRMs that are salient for their use in RFBs. Methylene linked phenothiazine and phthalimide moieties were employed as model BRMs and variations in their properties and molecular conformations were evaluated by systematically varying the methylene linker. Our results revealed that even minimalistic changes in the linker length resulted in dramatic oscillations in the solubility and stability. Using crystallography, quantum chemistry, and molecular dynamics, we clarify that the seemingly fluctuating behavior is due to (1) the inter- and intra-molecular charge transfer between the donor and acceptor through bond or/and through space; (2) the formation of distinctive packing/clustering motifs. Both interactions are strongly dependent on the molecular conformation. As these structural factors modify the electronic structure of both electroactive functional groups and do not change monotonically, the structure–property response in BRMs is often complex. Finally, accordingly, these behaviors need to be taken into account in developing BRMs and knowledge of these structural factors will allow the more rational design of BRMs for grid-scale energy storage.

25 ENERGY STORAGE↗

Control theory based airfoil design using the Euler equations

This paper describes the implementation of optimization techniques based on control theory for airfoil design. In our previous work it was shown that control theory could be employed to devise effective optimization procedures for two-dimensional profiles by using the potential flow equation with either a conformal mapping or a general coordinate system. The goal of our present work is to extend the development to treat the Euler equations in two-dimensions by procedures that can readily be generalized to treat complex shapes in three-dimensions. Therefore, we have developed methods which can address airfoil design through either an analytic mapping or an arbitrary grid perturbation method applied to a finite volume discretization of the Euler equations. Here the control law serves to provide computationally inexpensive gradient information to a standard numerical optimization method. Results are presented for both the inverse problem and drag minimization problem.

Jameson, Antony↗

High-precision early mission narrow angle sciene with the Space Interferometry Mission

We have developed a technique that allows SIM to measure relative stellar positions with an accuracy of 1 micro-arcsecond at any time during its 5-yr mission. Unlike SIM's standard narrow-angle approach, Gridless Narrow Angle Astrometry (GNAA) does not rely on the global reference frame of grid stars that reaches full accuracy after 5 years. GNAA is simply the application of traditional single-telescope narrow angle techniques to SIM's narrow angle optical path delay measurements. In GNAA, a set of reference stars and a target star are observed at several baseline orientations. A linearized model uses delay measurements to solve for star positions and baseline orientations. A conformal transformation maps observations at different epochs to a common reference frame. The technique works on short period signals (P=days to months), allowing it to be applied to many of the known extra-solar planets, intriguing radio/X- ray binaries, and other periodic sources. The technique's accuracy is limited in the long-term by false acceleration due to a combination of reference star and target star proper motion. The science capability 1 micro-arcsecond astrometric precision - is unique to SIM.

Space Interferometry Mission↗