Search NASA⌕ Search

SEARCH · Search NASA

Results for “Interior point methods”

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 163 records · Page 9

Effect of crack curvature on stress intensity factors for ASTM standard compact tension specimens

The stress intensity factors (SIF) are calculated using the method of lines for the compact tension specimen in tensile and shear loading for curved crack fronts. For the purely elastic case, it was found that as the crack front curvature increases, the SIF value at the center of the specimen decreases while increasing at the surface. For the higher values of crack front curvatures, the maximum value of the SIF occurs at an interior point located adjacent to the surface. A thickness average SIF was computed for parabolically applied shear loading. These results were used to assess the requirements of ASTM standards E399-71 and E399-81 on the shape of crack fronts. The SIF is assumed to reflect the average stress environment near the crack edge.

Alam, J.↗

Smart Refractory Sensors Development for Corrosion and Erosion Monitoring in High Temperature Systems

To optimize the operation and functioning of high temperature systems such as slagging gasifiers, coal boilers and glass/steel melters, it is important to monitor corrosion and erosion of refractory used in such systems. Corrosion test strategies are generally based on continuous gravimetric and chemical reactivity monitoring at operational temperatures (750°-1500°C). Both thermocouples and failure sensors and arrays would be useful to monitor the health of any refractory or coatings in these systems. Many of such type of sensors are installed into the systems through open access ports within the refractory; however, there are some disadvantages of this approach where corrosive/erosive gas and molten materials can penetrate and compromise the system. The current work presents the development and performance demonstration of smart refractory with embedded high temperature sensors such as thermocouples, thermistors, and various spallation/crack monitoring sensors, which may be used within a variety of refractory brick in different high temperature processes and applications. The main feature of this technology is that electroceramic based sensors are embedded into smart refractory without significantly impact to the intrinsic properties of the refractory. This technology circumvents the need to insert an isolated monolithic, stand-alone sensor into the refractory via an access port. This technological approach guarantees the integrity and the chemical stability of the materials used in the sensor fabrication within the harsh environment and does not introduce molten material (such as slag) penetration pathways within the refractory. One interesting and important aspect of this innovation is that these embedded sensors can be used to in situ monitoring processes such as chemical reactions and at the same time give information and a deeper understanding of the corrosion and erosion process of the refractory within the system. As stated above, the objective of our work is to develop high-temperature sensors composed of electroceramic materials that are chemically stable at high temperatures (750°-1500°C) and high pressures (up to 1000 psi) that can be used in monitoring corrosion and erosion process in refractory used in high energy systems. The high-temperature sensors investigated in this work were composed of various oxide composites directly embedded into the refractory oxides. The composites used for this work were synthesized by a mixed-oxide route. Metal oxides were inserted within a matrix material composed of refractory oxides (Al2O3, ZrO2, etc.). The physical and electrical properties were specifically manipulated by altering the level of percolation of the conductive species (metal oxides) within the refractory constituent (refractory oxide). Prior to the development of the high-temperature sensors, the oxides composites developed in this study were sintered up to 1600°C under oxidizing atmosphere in order to investigate densification, microstructural evolution, phase development, and their thermoelectrical performance as a function of the composition. The 4-point DC conductivity measurements were performed between 100°-1500°C. The sensors were fabricated from the composite materials by 3D-printing or screen-printing methods into the refractory brick during the consolidation process. An example of one of these embedded sensors consisted of an electroceramic-based thermocouple fabricated with two separate oxide composite compositions which were patterned to produce a couple within the interior of a refractory matrix. The thermocouple successfully displayed thermoelectric voltage trend (as a function of temperature), and the voltage was 220.0 mV around 1400 °C. Corrosion tests on the refractory embedded sensors were performed. To evaluate corrosion in the refractory brick an in-house glass composition was prepared and pressed into pellets and delivered into a pre-cut cavity in the brick. Corrosion experiments results showed the glass penetrated the brick over a 90 h period, and the penetration of the glass through the brick could be monitored by both an amperometric and voltametric based sensor. With this experiment, it was demonstrated that the embedded sensor could dynamically monitor the corrosion process.

20 FOSSIL-FUELED POWER PLANTS↗

Flow field characteristics of a 127-pin rod bundle with hexagonal spacer grids

A determination of nominal flow phenomena in liquid metal fast reactor (LMFR) fuel assemblies is critical toward generation-IV reactor development. Axially positioned spacer grids are used to maintain the geometry of hexagonal rod bundles and simultaneously introduce perturbations in the flow. Three-dimensional (3D) printed asymmetric honeycomb spacer grids were installed in a prototypical 127-pin LMFR fuel assembly model to study complex fluid dynamics interactions induced by the spacer grid and rods. To characterize flow dynamics in this intricate geometry, time-resolved particle image velocimetry (TR-PIV) using the matched-index-of-refraction method was employed to obtain non-intrusive velocity measurements for three axial planes (one near-wall and two interior planes) at a Reynolds number of 6000. The statistical TR-PIV results compared sub-channel-dependent normalized time-averaged velocity, velocity fluctuations, Reynolds stress, vorticity, and turbulence kinetic energy distributions. TR-PIV line profiles characterized downstream spacer grid flow dynamics. Two-point spatial and spatial–temporal cross-correlation fields revealed local coherent structures and quantified convection velocities of traveling vortices. Spatial–temporal decomposition using dynamic mode decomposition (DMD) applied to the near-wall vorticity fields extracted turbulent structures and flow instabilities in the wake region of the spacer grid, along with their decay and frequency rates. Reduced-order velocity fields from DMD reconstructions identified the most energy-containing coherent structures persistent in the near-wall region. This research provides experimental data sets and analyses of flow behavior in rod bundles with hexagonal spacer grids. In conclusion, the results are critical toward LMFR design and geometry optimization, crucial for the validation of computational fluid dynamics and reduced-order flow models.

42 ENGINEERING↗

An Accelerated Clip Algorithm for Unstructured Meshes: A Batch-Driven Approach

The clip technique is a popular method for visualizing complex structures and phenomena within 3D unstructured meshes. Meshes can be clipped by specifying a scalar isovalue to produce an output unstructured mesh with its external surface as the isovalue. Similar to isocontouring, the clipping process relies on scalar data associated with the mesh points, including scalar data generated by implicit functions such as planes, boxes, and spheres, which facilitates the visualization of results interior to the grid. In this paper, we introduce a novel batch-driven parallel algorithm based on a sequential clip algorithm designed for high-quality results in partial volume extraction. Our algorithm comprises five passes, each progressively processing data to generate the resulting clipped unstructured mesh. The novelty lies in the use of fixed-size batches of points and cells, which enable rapid workload trimming and parallel processing, leading to a significantly improved memory footprint and run-time performance compared to the original version. On a 32-core CPU, the proposed batch-driven parallel algorithm demonstrates a run-time speed-up of up to 32.6x and a memory footprint reduction of up to 4.37x compared to the existing sequential algorithm. The software is currently available under an open-source license in the VTK visualization system.

Tsalikis, Spiros↗

Numerical determination of the fundamental eigenvalue for the Laplace operator on a spherical domain

Methods for obtaining approximate solutions for the fundamental eigenvalue of the Laplace-Beltrami operator (i.e., the membrane eignevalue problem for the vibration equation) on the unit spherical surface are developed. Two types of spherical surface domains are considered: the interior of a spherical triangle, and the exterior of a great circle arc extending for less than pi radians (a spherical surface with a slit). In both cases, zero boundary conditions are imposed. In order to solve the resulting second-order elliptic partial differential equations in two independent variables, a finite difference approximation is employed. The fundamental eigenvalue is approximated by iteration utilizing the power method and point successive overrelaxation. Some numerical results are given and compared, in certain special cases, with analytical solutions to the eigenvalue problem. The significance of the numerical eigenvalue results is discussed in terms of the singularities in the solution of three-dimensional boundary-value problems near a polyhedral corner of the domain.

Walden, H.↗

Flow structure capturing on overset patched meshes

The present paper describes the application of a simple, robust stable implicit approach to solutions of the conservative equations of gasdynamics on either composite or overset meshes to capture flow structures such as shocks with high resolution in complex geometric domains. Without requiring special flux conservative operators, but rather, interpolating conservative variable data at mesh boundaries, the well posed upwind method provides accurate and oscillation free solutions, even when shocks cross the interior patch boundaries. In three problems with flow complexities that require grid refinement, the paper demonstrates the capability to conveniently carry out for gasdynamics the adaptive refined meshing strategy in overset patches proposed by Berger and Oliger (1984), and it extends this technique to rapidly convergent implicit methods for the Euler and Navier-Stokes equations. The numerical experiments show concretely, in a realistic aerodynamic problem, the savings in mesh points (about an order of magnitude here in two dimensions) for similar accuracy that flow structure aligned adaptive patched meshing affords compared to uniform grid refinement.

Venkatapathy, E.↗

Coupling finite element and integral equation solutions using decoupled boundary meshes

A method is outlined for calculating scattered fields from inhomogeneous penetrable objects using a coupled finite element-integral equation solution. The finite element equation can efficiently model fields in penetrable and inhomogeneous regions, while the integral equation exactly models fields on the finite element mesh boundary and in the exterior region. By decoupling the interior finite element and exterior integral equation meshes, considerable flexibility is found in both the number of field expansion points as well as their density. Only the nonmetal portions of the object need be modeled using a finite element expansion; exterior perfect conducting surfaces are modeled using an integral equation with a single unknown field since E(tan) is identically zero on these surfaces. Numerical convergence, accuracy, and stability at interior resonant frequencies are studied in detail.

Cwik, Tom↗

An improved higher order panel method for linearized supersonic flow

An improved higher order panel method for linearized supersonic flow is described. Each panel, defined by four points on the surface, is divided into eight subpanels in such a way that all subpanel and panel edges are contiguous. By prescribing a quadratic distribution of the doublet on each subpanel, the doublet strength is made strictly continuous on the paneled surface. A linear source distribution is also used. Numerical results are smoother and in better agreement with experiment than the previous method with less strict continuity. A brief discussion of superinclined panels used to eliminate interior interference in nacelles is included.

Ehlers, F. E.↗

Three-stage sorption type cryogenic refrigeration systems and methods employing heat regeneration

A three-stage sorption type cryogenic refrigeration system, each stage containing a fluid having a respectively different boiling point, is presented. Each stage includes a compressor in which a respective fluid is heated to be placed in a high pressure gaseous state. The compressor for that fluid which is heated to the highest temperature is enclosed by the other two compressors to permit heat to be transferred from the inner compressor to the surrounding compressors. The system may include two sets of compressors, each having the structure described above, with the interior compressors of the two sets coupled together to permit selective heat transfer therebetween, resulting in more efficient utilization of input power.

Bard, Steven↗

Improved Maneuver Reconstructions for the GRAIL Orbiters

Maneuver reconstructions for the Gravity Recovery and Interior Laboratory (GRAIL) A and B lunar orbiters were improved through updates to the orbit determination filter and dynamic models. Consistent reconstructions of the 27 GRAIL A and B maneuvers from the Trans-Lunar Cruise phase in the fall of 2011 through the Transition to Science Formation phase in February 2012 were performed. Standard methods of orbit determination were applied incorporating the latest dynamic models and filter strategies developed by the GRAIL Navigation and Science Teams, including a high resolution, 420 x 420 degree and order lunar spherical harmonic gravity field model. For Trans-Lunar Cruise for GRAIL-A (TLC-A), all maneuvers executed with delta V errors below 5.50 +/- 0.50 mm/s and pointing errors below 0.25 degrees. GRAIL-A lunar orbit maneuvers had delta V errors below 30.0 mm/s and pointing errors below 0.51 degrees. For TLC-B, all maneuvers executed with delta V errors below 8.60 +/- 1.41 mm/s and pointing errors below 0.300 degrees. GRAIL-B maneuvers in lunar orbit executed with maximum delta V errors of 25.0 mm/s and pointing error of 0.43 degrees. These maneuver reconstructions will enable the GRAIL Navigation Team to better characterize the main engine performance of each spacecraft. This will help the Navigation Team to navigate low (greater than 8 km) altitude orbits during the extended mission phase in the fall of 2012.

lunar mapping↗

Shuttle Imaging Radar B (SIR-B) Weddell Sea ice observations - A comparison of SIR-B and scanning multichannel microwave radiometer ice concentrations

Ice concentrations over the Weddell Sea were studied using SIR-B data obtained during the October 1984 mission, with special attention given to the effect of ocean waves on the radar return at the ice edge. Sea ice concentrations were derived from the SIR-B data using two image processing methods: the classification scheme at JPL and the manual classification method at Scott Polar Research Institute (SPRI), England. The SIR ice concentrations were compared with coincident concentrations from the Nimbus-7 SMMR. For concentrations greater than 40 percent, which was the smallest concentration observed jointly by SIR-B and the SMMR, the mean difference between the two data sets for 12 points was 2 percent. A comparison between the JPL and the SPRI SIR-B algorithms showed that the algorithms agree to within 1 percent in the interior ice pack, but the JPL algorithm gives slightly greater concentrations at the ice edge (due to the fact that the algorithm is affected by the wind waves in these areas).

Martin, Seelye↗

Convection effects in protein crystal growth

Protein crystals for X-ray diffraction study are usually grown resting on the bottom of a hanging drop of a saturated protein solution, with slow evaporation to the air in a small enclosed cell. The evaporation rate is controlled by hanging the drop above a reservoir of water, with its saturation vapor pressure decreased by a low concentration of a passive solute. The drop has a lower solute concentration, and its volume shrinks by evaporation until the molecular concentrations match. Protein crystals can also be grown from a seed crystal suspended or supported in the interior of a supersaturated solution. The main analysis of this report concerns this case because it is less complicated than hanging-drop growth. Convection effects have been suggested as the reason for the apparent cessation of growth at a certain rather small crystal size. It seeems that as the crystal grows, the number of dislocations increases to a point where further growth is hindered. Growth in the microgravity environment of an orbiting space vehicle has been proposed as a method for obtaining larger crystals. Experimental observations of convection effects during the growth of protein crystals have been reported.

Roberts, Glyn O.↗

Nanoindentation of Micrograins in Polycrystals under Multiaxial Stress: Control of Abrupt & Stochastic Plastic Events

Deforming objects in the world of the small has proven a much more complex task than just bending a common-sized object like a spoon. Over the last decade, it has become clear that the application of stress on micron-sized objects, such as micropillars, leads to abrupt and stochastic plastic strain jumps that do not satisfy Weibull statistics (as expected for uncorrelated defect yielding) but instead seem to follow highly wide, power-law distributions, pointing towards cooperative dislocation plasticity. Are power-law distributed plastic events the rule or the exception in crystals? If the rule, can we design protocols to control and suppress them in order to achieve deterministic and smooth microforming processes? We propose to combine experiments and theory to unravel the mysteries of abrupt microplastic deformation in a natural dislocation environment: micron-sized grains (\micrograins") of polycrystals. We propose a novel experimental design that utilizes nanoindentation and 4-point bending techniques on the substrate of polycrystalline thin films of two typical FCC metals (Al and Cu), and make concrete theoretical conjectures, based on intuition built from previous theoretical and experimental efforts in micropillar compression and nanoindentation. Based on a wide set of observations over the last decade, we put forward two distinct hypotheses that we propose to investigate in both theory and experiments: Our first conjecture is that crystal plasticity in the central region of micron-sized grains has micropillar-like stochastic features: power-law plastic strain jumps at stresses above the macroscale yield strength, displaying strong size and rate effects. We will demonstrate this similarity in spherical nanoindentation experiments in the central-region of multi-orientation grains of Al and Cu thin films. Al and Cu are selected for their disparity in stacking fault energy and consequently, their different propensity for cross-slip relaxation; however, our protocol may be applied to any crystalline material of technological interest. We will test our hypothesis by performing two-dimensional extended discrete dislocation plasticity (2.5D-DDP) simulations, inspired by comparable successful studies of micropillar plasticity, and extended to model nanoindentation with two slip systems and on a surface of a bicrystal with microhard grain-boundaries (GB), ie. not allowing dislocation absorption (116). Our model will include effective three dimensional dislocation mechanisms such as dislocation multiplication and double-cross-slip assisted glide. Second, we hypothesize that the application of uniform in-plane stress on our thin films below the yield point can smoothen the micrograin plastic response during nanoindentation. On a technological level, we propose that indentation-induced localization (sink-in) may be diminished by the application of concurrent in-plane stress that activates slow, viscoplastic relaxations near GBs. We aim to demonstrate our hypothesis through a 4-point bending apparatus that introduces a homogeneous in-plane stress state before indentation. We will perform simulations of a coarse-grained mesoscale theory of plasticity in a system with microhard GBs, benchmarked through 2.5D-DDP simulations, which will capture the effects of slow viscoplastic relaxations. Our nanoindentation abrupt event statistics studies, either in the micrograin interior or near GBs, will be correlated to direct observations of induced surface dislocation microstructures. For we will use TEM-based methods (PACOM/TKD) in combination with EBSD. Indents placed in the center of grains will be analyzed based on the orientation of the grain, while those near to GBs will be analyzed based on the orientation of each grain, the distance from the GB and the 5 degrees of freedom of the GB.

36 MATERIALS SCIENCE↗

Method for fabricating a mass spectrometer inlet leak

A method of fabricating an inlet leak for measuring chemically reactive gases is proposed. The method involves the steps of preparing a selected internal surface of a high melting point metal tubing for sealing to a glass tubing. A portion of the glass tubing is inserted into the prepared portion of the metal tubing, the glass tubing having an outer diameter substantially the same as the inner diameter of the metal tubing. The inserted glass tubing is sealed to the internal surface of the metal tubing; the glass tubing is drawn under a magnifying device to reduce the entire interior diameter thereof to substantially two microns; and an end portion of the drawing glass tubing is broken away.

Harris, R. F.↗

Triangular Berstein-bezier Patches, a Survey and New Results

Triangular Bernstein-Bezier patches are an alternative to the standard rectangular ones. They are maps of triangular domains into R cubed, where the patches are described by control nets similar to the rectangular case. One difference in the definition of surface normals at the patch corners: rectangular patches have inherent problems with so-called degenerate corners, i.e., edges collapsing into one point. Triangular patches do not have these problems. Since triangular shapes arise rather frequently in the design of complicated surfaces (e.g., interior car body panels), it seems that for such surfaces triangular patches are a promising addition to the usually employed rectangular patches. Another major application is the formulation of interpolants to solve the scattered data problem: here Bernstein-Bezier patches provide a very powerful tool; they can be used to discover new schemes as well as to discover new properties of existing methods.

Farin, G.↗

Experimental Investigation of the Distribution of Shock Effects in Regolith Impact Ejecta Using an Ejecta Recovery Chamber

Because the mass-flux of solar system meteoroids is concentrated in the approx. 200 microns size range, small-scale impacts play a key role in driving the space weathering of regoliths on airless bodies. Quantifying this role requires improved data linking the mass, density and velocity of the incoming impactors to the nature of the shock effects produced, with particular emphasis on effects, such as production of impact melt and vapor, that drive the optical changes seen in space weathered regoliths. Of particular importance with regard to space weathering is understanding not only the composition of the shock melt created in small-scale impacts, but also how it is partitioned volumetrically between the local impact site and more widely distributed ejecta. To improve the ability of hypervelocity impact experiments to obtain this type of information, we have developed an enclosed sample target chamber with multiple-geometry interior capture cells for in-situ retention of ejecta from granular targets. A key design objective was to select and test capture cell materials that could meet three requirements: 1) Capture ejecta fragments traveling at various trajectories and velocities away from the impact point, while inducing minimal additional damage relative to the primary shock effects; 2) facilitate follow-up characterization of the ejecta either on or in the cell material by analytical SEM, or ex-situ by microprobe, TEM and other methods; and 3) enable the trajectories of the captured and characterized ejecta to be reconstructed relative to the target.

Christoffersen, R.↗

Radiation/convection coupling in rocket motors and plumes

The three commonly used propellant systems - H2/O2, RP-1/O2, and solid propellants - primarily radiate as molecular emitters, non-scattering small particles, and scattering larger particles, respectively. Present technology has accepted the uncoupling of the radiation analysis from that of the flowfield. This approximation becomes increasingly inaccurate as one considers plumes, interior rocket chambers, and nuclear rocket propulsion devices. This study will develop a hierarchy of methods which will address radiation/convection coupling in all of the aforementioned propulsion systems. The nature of the radiation/convection coupled problem is that the divergence of the radiative heat flux must be included in the energy equation and that the local, volume-averaged intensity of the radiation must be determined by a solution of the radiative transfer equation (RTE). The intensity is approximated by solving the RTE along several lines of sight (LOS) for each point in the flowfield. Such a procedure is extremely costly; therefore, further approximations are needed. Modified differential approximations are being developed for this purpose. It is not obvious which order of approximations are required for a given rocket motor analysis. Therefore, LOS calculations have been made for typical rocket motor operating conditions in order to select the type approximations required. The results of these radiation calculations, and the interpretation of these intensity predictions are presented herein.

Farmer, R. C.↗

Effect of the V{sub As}V{sub Ga} complex defect doping on properties of the semi-insulating GaAs

The different position V{sub As}V{sub Ga} cluster defect doping in semi-insulating (SI) GaAs has been studied by first-principles calculation based on hybrid density functional theory. Our calculated results show that EL6 level is formed due to the V{sub As}V{sub Ga} complex defect, which is very close to the experimental result. It provides the explanation of the absorption of laser with the wavelength beyond in semi-insulating GaAs. The formation energy of V{sub As}V{sub Ga} complex defect is found to decrease from surface to interior gradually. The conduction band minima and valence band maxima of GaAs (001) surface with the V{sub As}V{sub Ga} complex defect are all located at Γ point, and some defect levels are produced in the forbidden band. In contrast, the conduction band minima and valence band maxima of GaAs with the interior V{sub As}V{sub Ga} complex defect are not located at the same k-point, so it might involve the change of momentum in the electron transition process. The research will help strengthen the understanding of photoelectronic properties and effectively guide the preparation of the SI-GaAs materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗