Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shape Decomposition”

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 55 records · Page 3

Effect of Surface Nonequilibrium Thermochemistry in Simulation of Carbon Based Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver using finite-rate gas/surface interaction model provides time-accurate solutions for multidimensional ablation of carbon based charring ablators. The material thermal response code used in this study is the Two-dimensional Implicit Thermal Response and Ablation Program (TITAN), which predicts charring material thermal response and shape change on hypersonic space vehicles. Its governing equations include total energy balance, pyrolysis gas momentum conservation, and a three-component decomposition model. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation (DPLR) method. Loose coupling between material response and flow codes is performed by solving the surface mass balance in DPLR and the surface energy balance in TITAN. Thus, the material surface recession is predicted by finite-rate gas/surface interaction boundary conditions implemented in DPLR, and the surface temperature and pyrolysis gas injection rate are computed in TITAN. Two sets of gas/surface interaction chemistry between air and carbon surface developed by Park and Zhluktov, respectively, are studied. Coupled fluid-material response analyses of stagnation tests conducted in NASA Ames Research Center arc-jet facilities are considered. The ablating material used in these arc-jet tests was a Phenolic Impregnated Carbon Ablator (PICA). Computational predictions of in-depth material thermal response and surface recession are compared with the experimental measurements for stagnation cold wall heat flux ranging from 107 to 1100 Watts per square centimeter.

Chen, Yih-Kang↗

Effect of Non-Equilibrium Surface Thermochemistry in Simulation of Carbon Based Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver using non-equilibrium gas/surface interaction model provides time-accurate solutions for the multidimensional ablation of carbon based charring ablators. The material thermal response code used in this study is the Two-dimensional Implicit Thermal-response and AblatioN Program (TITAN), which predicts charring material thermal response and shape change on hypersonic space vehicles. Its governing equations include total energy balance, pyrolysis gas mass conservation, and a three-component decomposition model. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation (DPLR) method. Loose coupling between the material response and flow codes is performed by solving the surface mass balance in DPLR and the surface energy balance in TITAN. Thus, the material surface recession is predicted by finite-rate gas/surface interaction boundary conditions implemented in DPLR, and the surface temperature and pyrolysis gas injection rate are computed in TITAN. Two sets of nonequilibrium gas/surface interaction chemistry between air and the carbon surface developed by Park and Zhluktov, respectively, are studied. Coupled fluid-material response analyses of stagnation tests conducted in NASA Ames Research Center arc-jet facilities are considered. The ablating material used in these arc-jet tests was Phenolic Impregnated Carbon Ablator (PICA). Computational predictions of in-depth material thermal response and surface recession are compared with the experimental measurements for stagnation cold wall heat flux ranging from 107 to 1100 Watts per square centimeter.

Chen, Yih-Kanq↗

Atmospheric Energy Spectra in Global Kilometre-Scale Models

Eleven 40-day long integrations of five different global models with horizontal resolutions of less than 9 km are compared in terms of their global energy spectra. The method of normal-mode function decomposition is used to distinguish between balanced (Rossby wave; RW) and unbalanced (inertia-gravity wave; IGW) circulation. The simulations produce the expected canonical shape of the spectra, but their spectral slopes at mesoscales, and the zonal scale at which RW and IGW spectra intersect differ significantly. The partitioning of total wave energies into RWs an IGWs is most sensitive to the turbulence closure scheme and this partitioning is what determines the spectral crossing scale in the simulations, which differs by a factor of up to two. It implies that care must be taken when using simple spatial filtering to compare gravity wave phenomena in storm-resolving simulations, even when the model horizontal resolutions are similar. In contrast to the energy partitioning between the RWs and IGWs, changes in turbulence closure schemes do not seem to strongly affect spectral slopes, which only exhibit major differences at mesoscales. Despite their minor contribution to the global (horizontal kinetic plus potential available) energy, small scales are important for driving the global mean circulation. Our results support the conclusions of previous studies that the strength of convection is a relevant factor for explaining discrepancies in the energies at small scales. The models studied here produce the major large-scale features of tropical precipitation patterns. However, particularly at large horizontal wavenumbers, the spectra of upper tropospheric vertical velocity, which is a good indicator for the strength of deep convection, differ by factors of three or more in energy. High vertical kinetic energies at small scales are mostly found in those models that do not use any convective parameterisation.

DYAMOND↗

A general algorithm using finite element method for aerodynamic configurations at low speeds

A finite element algorithm for numerical simulation of two-dimensional, incompressible, viscous flows was developed. The Navier-Stokes equations are suitably modelled to facilitate direct solution for the essential flow parameters. A leap-frog time differencing and Galerkin minimization of these model equations yields the finite element algorithm. The finite elements are triangular with bicubic shape functions approximating the solution space. The finite element matrices are unsymmetrically banded to facilitate savings in storage. An unsymmetric L-U decomposition is performed on the finite element matrices to obtain the solution for the boundary value problem.

Balasubramanian, R.↗

Efficient computation of volume in flow predictions

An efficient method for calculating cell volumes for time-dependent three-dimensional flow predictions by finite volume calculations is presented. Eight arbitrary corner points are considered and the shape face is divided into two planar triangles. The volume is then dependent on the orientation of the partitioning. In the case of a hexahedron, it is noted that any open surface with a boundary that is a closed curve possesses a surface vector independent of the surface shape. Expressions are defined for the surface vector, which is independent of the partitioning surface diagonal used to quantify the volume. Using a decomposition of the cell volume involving two corners, with each the vertex of three diagonals and six corners which are vertices of one diagonal, gives portions which are tetrahedra. The resultant mesh is can be used for time-dependent finite volume calculations one requires less computer time than previous methods.

Vinokur, M.↗

A Parallel Cartesian Approach for External Aerodynamics of Vehicles with Complex Geometry

This workshop paper presents the current status in the development of a new approach for the solution of the Euler equations on Cartesian meshes with embedded boundaries in three dimensions on distributed and shared memory architectures. The approach uses adaptively refined Cartesian hexahedra to fill the computational domain. Where these cells intersect the geometry, they are cut by the boundary into arbitrarily shaped polyhedra which receive special treatment by the solver. The presentation documents a newly developed multilevel upwind solver based on a flexible domain-decomposition strategy. One novel aspect of the work is its use of space-filling curves (SFC) for memory efficient on-the-fly parallelization, dynamic re-partitioning and automatic coarse mesh generation. Within each subdomain the approach employs a variety reordering techniques so that relevant data are on the same page in memory permitting high-performance on cache-based processors. Details of the on-the-fly SFC based partitioning are presented as are construction rules for the automatic coarse mesh generation. After describing the approach, the paper uses model problems and 3- D configurations to both verify and validate the solver. The model problems demonstrate that second-order accuracy is maintained despite the presence of the irregular cut-cells in the mesh. In addition, it examines both parallel efficiency and convergence behavior. These investigations demonstrate a parallel speed-up in excess of 28 on 32 processors of an SGI Origin 2000 system and confirm that mesh partitioning has no effect on convergence behavior.

Aftosmis, M. J.↗

Evidence for exclusively inorganic formation of magnetite in Martian meteorite ALH84001

Magnetite crystals produced by terrestrial magnetotactic bacterium MV-1 are elongated on a [111] crystallographic axis, in a so-called truncated hexa-Octahedral shape. This morphology has been proposed to constitute a biomarker (i.e., formed only in biogenic processes). A subpopulation of magnetite crystals associated with carbonate globules in Martian meteorite ALH84001 is reported to have this morphology, and the observation has been taken as evidence for biological activity on Mars. In this study, we present evidence for the exclusively inorganic origin of [111]-elongated magnetite crystals in ALH84001. We report three-dimensional(3-D) morphologies for approx.1000 magnetite crystals extracted from: (1) thermal decomposition products of Fe-rich carbonate produced by inorganic hydrothermal precipitation in laboratory experiments; (2) carbonate globules in Martian meteoriteeALH84001; and (3) cells of magnetotactic bacterial strain MV-1. The 3-D morphologies were derived by fitting 3-D shape models to two-dimensional bright-field transmission-electron microscope (TEAM) images obtained at a series of viewing angles. The view down the {110} axes closest to the [111] elongation axis of magnetite crystals ([111]x{110) not equal to 0) provides a 2-D projection that uniquely discriminates among the three [111]-elongated magnetite morphologies found in these samples: [111]-elongated truncated hexaoctahedron ([111]-THO), [111]-elongated cubo-octahedron ([111]-ECO), and [111]-elongated simple octahedron ([111]-ESO). All [111] -elongated morphologies are present in the three types of sample, but in different proportions. In the ALH84001 Martian meteorite and in our inorganic laboratory products, the most common [111]-elongated magnetite crystal morphology is [111]-ECO. In contrast, the most common morphology for magnetotactic bacterial strain MV-1 is [111]-THO. These results show that: (1) the morphology of [111]-elongated magnetite crystals associated with the carbonate globules in Martian meteorite ALH84001 is replicated by an inorganic process; and (2) the most common crystal morphology for biogenic (MV-1) magnetite is distinctly different from that in both ALH84001 and our inorganic laboratory products. Therefore, [111]-elongated magnetite crystals in ALH84001 do not constitute, as previously claimed, a robust biosignature and, in fact, an exclusively inorganic origin for the magnetite is fully consistent with our results. Furthermore, the inorganic synthesis method, i.e., the thermal decomposition of hydrothermally precipitated Fe-rich carbonate, is a process analogue for formation of the magnetite on Mars. Namely, precipitation of carbonate globules from carbonate-rich hydrothermal solutions followed at some later time by a thermal pulse, perhaps in association with meteoritic impact or volcanic processes on the Martian surface.

Golden, D. C.↗

Parametric binary dissection

Binary dissection is widely used to partition non-uniform domains over parallel computers. This algorithm does not consider the perimeter, surface area, or aspect ratio of the regions being generated and can yield decompositions that have poor communication to computation ratio. Parametric Binary Dissection (PBD) is a new algorithm in which each cut is chosen to minimize load + lambda x(shape). In a 2 (or 3) dimensional problem, load is the amount of computation to be performed in a subregion and shape could refer to the perimeter (respectively surface) of that subregion. Shape is a measure of communication overhead and the parameter permits us to trade off load imbalance against communication overhead. When A is zero, the algorithm reduces to plain binary dissection. This algorithm can be used to partition graphs embedded in 2 or 3-d. Load is the number of nodes in a subregion, shape the number of edges that leave that subregion, and lambda the ratio of time to communicate over an edge to the time to compute at a node. An algorithm is presented that finds the depth d parametric dissection of an embedded graph with n vertices and e edges in O(max(n log n, de)) time, which is an improvement over the O(dn log n) time of plain binary dissection. Parallel versions of this algorithm are also presented; the best of these requires O((n/p) log(sup 3)p) time on a p processor hypercube, assuming graphs of bounded degree. How PBD is applied to 3-d unstructured meshes and yields partitions that are better than those obtained by plain dissection is described. Its application to the color image quantization problem is also discussed, in which samples in a high-resolution color space are mapped onto a lower resolution space in a way that minimizes the color error.

Bokhari, Shahid H.↗

RIACS

Topics considered include: high-performance computing; cognitive and perceptual prostheses (computational aids designed to leverage human abilities); autonomous systems. Also included: development of a 3D unstructured grid code based on a finite volume formulation and applied to the Navier-stokes equations; Cartesian grid methods for complex geometry; multigrid methods for solving elliptic problems on unstructured grids; algebraic non-overlapping domain decomposition methods for compressible fluid flow problems on unstructured meshes; numerical methods for the compressible navier-stokes equations with application to aerodynamic flows; research in aerodynamic shape optimization; S-HARP: a parallel dynamic spectral partitioner; numerical schemes for the Hamilton-Jacobi and level set equations on triangulated domains; application of high-order shock capturing schemes to direct simulation of turbulence; multicast technology; network testbeds; supercomputer consolidation project.

Oliger, Joseph↗

Process modelling for materials preparation experiments

The main goals of the research under this grant consist of the development of mathematical tools and measurement techniques for transport properties necessary for high fidelity modelling of crystal growth from the melt and solution. Of the tasks described in detail in the original proposal, two remain to be worked on: development of a spectral code for moving boundary problems, and development of an expedient diffusivity measurement technique for concentrated and supersaturated solutions. We have focused on developing a code to solve for interface shape, heat and species transport during directional solidification. The work involved the computation of heat, mass and momentum transfer during Bridgman-Stockbarger solidification of compound semiconductors. Domain decomposition techniques and preconditioning methods were used in conjunction with Chebyshev spectral methods to accelerate convergence while retaining the high-order spectral accuracy. During the report period we have further improved our experimental setup. These improvements include: temperature control of the measurement cell to 0.1 C between 10 and 60 C; enclosure of the optical measurement path outside the ZYGO interferometer in a metal housing that is temperature controlled to the same temperature setting as the measurement cell; simultaneous dispensing and partial removal of the lower concentration (lighter) solution above the higher concentration (heavier) solution through independently motor-driven syringes; three-fold increase in data resolution by orientation of the interferometer with respect to diffusion direction; and increase of the optical path length in the solution cell to 12 mm.

Rosenberger, Franz↗

Influence of flow on interface shape stability in low gravity

The objectives are to: (1) Understand the influence in low gravity of flow on interface shape. For example, document and control the influence of axial flow on the Plateau-Rayleigh instability of a liquid bridge; and (2) Extend the ground-based density-matching technique of low gravity simulation to situations with flow; that is, develop Plateau chamber experiments for which flow can be controlled. Containerless containment of liquid by surface tension has broad importance in low gravity. For space vehicles, the behavior of liquid/gas interfaces is crucial to successful liquid management systems. In microgravity science, free interfaces are exploited in various applications. Examples include float-zone crystal growth, phase separation near the critical point of liquid mixtures (spinoidal decomposition) and quenching of miscibility gap molten metal alloys. In some cases, it is desired to stabilize the capillary instability while in others it is desired to induce capillary breakup. In all cases, understanding the stability of interface shape in the presence of liquid motion is central.

Steen, Paul H.↗

Subband And Transform Compression Of Video Signals

Class of hierarchical subband coders developed primarily for compression of image data at video rates. Offers good performance with limited computational complexity and with flexibility inherent in subband decomposition. Particular subband decomposition chosen for coders appears to hide large quantitative errors effectively, largely because decomposition occurs along two-dimensional spatial-frequency-domain boundaries resembling spatial-frequency-domain curves of constant sensitivity of human visual system. Curves found approximately diamond-shaped: thus, low-pass filtering for reduction of data ideally involves nonrectangular passbands.

Sauer, Ken↗

High-Speed Schlieren Analysis of Retropropulsion Jet in Mach 10 Flow

High-speed schlieren imaging of a 5-in.-diam, 70° sphere-cone model with a single centerline 15° half-angle nozzle was performed in a Mach 10 flow. Image sequences captured at 100 kHz were obtained with the nozzle plugged (fully blunted model), no nozzle flow, and with nozzle flow over a range of supply pressures. The freestream unit Reynolds number was kept constant for all runs at3.38×106m−1. Time-average information pertaining to the location, width, and stand-off distance of the bow shock, interface, terminal shock, and triple point as a function of jet pressure ratio are presented based on analysis of covariance images from each run. Measurement of jet plume boundary shape as a function of jet pressure ratio is also presented. Analyses of frequency content associated with the jet structure are shown, as are results of the proper orthogonal decomposition of these data. For all runs with nozzle flow, a dominant fundamental frequency of 2 kHz was observed. Finally, a discussion of bow shock unsteadiness resulting from interaction with freestream disturbances for the blunted and no nozzle flow is provided, as is a discussion on the formation of weak waves emanating from the bow shock resulting from this interaction.

Brett F Bathel↗

An efficient solution procedure for the thermoelastic analysis of truss space structures

A solution procedure is proposed for the thermal and thermoelastic analysis of truss space structures in periodic motion. In this method, the spatial domain is first descretized using a consistent finite element formulation. Then the resulting semi-discrete equations in time are solved analytically by using Fourier decomposition. Geometrical symmetry is taken advantage of completely. An algorithm is presented for the calculation of heat flux distribution. The method is demonstrated via a numerical example of a cylindrically shaped space structure.

Givoli, D.↗

Moisture and Temperature Influences on Nonlinear Vegetation Trends in Serengeti National Park

While long-term vegetation greening trends have appeared across large land areas over the late 20th century, uncertainty remains in identifying and attributing finer-scale vegetation changes and trends, particularly across protected areas. Serengeti National Park (SNP) is a critical East African protected area, where seasonal vegetation cycles support vast populations of grazing herbivores and a host of ecosystem dynamics. Previous work has shown how non-climate drivers (e.g. land use) shape the SNP ecosystem, but it is still unclear to what extent changing climate conditions influence SNP vegetation, particularly at finer spatial and temporal scales. We fill this research gap by evaluating long-term (1982–2016) changes in SNP leaf area index (LAI) in relation to both temperature and moisture availability using Ensemble Empirical Mode Decomposition and Principal Component Analysis with regression techniques. We find that SNP LAI trends are nonlinear, display high sub-seasonal variation, and are influenced by lagged changes in both moisture and temperature variables and their interactions. LAI during the long rains (e.g. March) exhibits a greening-to-browning trend reversal starting in the early 2000s, partly due to antecedent precipitation declines. In contrast, LAI during the short rains (e.g. November, December) displays browning-to-greening alongside increasing moisture availability. Rising temperature trends also have important, secondary interactions with moisture variables to shape these SNP vegetation trends. Our findings show complex vegetation-climate interactions occurring at important temporal and spatial scales of the SNP, and our rigorous statistical approaches detect these complex climate-vegetation trends and interactions, while guarding against spurious vegetation signals.

Moisture↗

Carbon Nanotubes/Nanofibers by Plasma Enhanced Chemical Vapour Deposition

Plasma enhanced chemical vapour deposition (PECVD) has been recently used for the production of vertically aligned carbon nanotubedfibers (CN) directly on substrates. These structures are potentially important technologically as electron field emitters (e.g. microguns, microwave amplifiers, displays), nanoelectrodes for sensors, filter media, superhydrophobic surfaces and thermal interface materials for microelectronics. A parametric study on the growth of CN grown by glow discharge dc-PECVD is presented. In this technique, a substrate containing thin film Ni catalyst is exposed to C2H2 and NH3 gases at 700 C. Without plasma, this process is essentially thermal CVD which produces curly spaghetti-like CN as seen in Fig. 1 (a). With the plasma generated by biasing the substrate at -6OOV, we observed that the CN align vertically during growth as shown in Fig. l(b), and that the magnitude of the applied substrate bias affects the degree of alignment. The thickness of the thin film Ni catalyst was found to determine the average diameter and inversely the length of the CN. The yield and density of the CN were controlled by the use of different diffusion barrier materials under the Ni catalyst. Patterned CN growth [Fig. l(c)], with la variation in CN diameter of 4.1% and 6.3% respectively, is achieved by lithographically defining the Ni thin film prior to growth. The shape of the structures could be varied from very straight nanotube-like to conical tip-like nanofibers by increasing the ratio of C2H2 in the gas flow. Due to the plasma decomposition of C2H2, amorphous carbon (a-C) is an undesirable byproduct which could coat the substrate during CN growth. Using a combination of depth profiled Auger electron spectroscopy to study the substrate and in-situ mass spectroscopy to examine gas phase neutrals and ions, the optimal conditions for a-C free growth of CN is determined.

Teo, K. B. K.↗

A hybrid perturbation-Galerkin technique for partial differential equations

A two-step hybrid perturbation-Galerkin technique for improving the usefulness of perturbation solutions to partial differential equations which contain a parameter is presented and discussed. In the first step of the method, the leading terms in the asymptotic expansion(s) of the solution about one or more values of the perturbation parameter are obtained using standard perturbation methods. In the second step, the perturbation functions obtained in the first step are used as trial functions in a Bubnov-Galerkin approximation. This semi-analytical, semi-numerical hybrid technique appears to overcome some of the drawbacks of the perturbation and Galerkin methods when they are applied by themselves, while combining some of the good features of each. The technique is illustrated first by a simple example. It is then applied to the problem of determining the flow of a slightly compressible fluid past a circular cylinder and to the problem of determining the shape of a free surface due to a sink above the surface. Solutions obtained by the hybrid method are compared with other approximate solutions, and its possible application to certain problems associated with domain decomposition is discussed.

Geer, James F.↗

Experimental verification of a large flexible manipulator

A large experimental lightweight manipulator would be useful for material handling, for welding, or for ultrasonic inspection of a large structure, such as an airframe. The flexible parallel link mechanism is designed for high rigidity without increasing weight. This constrained system is analyzed by singular value decomposition of the constraint Jacobian matrix. A verification of the modeling using the assumed mode method is presented. Eigenvalues and eigenvectors of the linearized model are compared to the measured system natural frequencies and their associated mode shapes. The modeling results for large motions are compared to the time response data from the experiments. The hydraulic actuator is verified.

Lee, Jac Won↗