SEARCH · Engineering Papers

Results for “interior point method”

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.

SourceResearch connectionFrom yearThrough year

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

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Experimental verification of the asymtotic modal analysis method as applied to a rectangular acoustic cavity excited by structural vibration

An experiment was performed on a rigid wall rectangular acoustic cavity driven by a flexible plate mounted in a quarter of one end wall and excited by white noise. The experiment was designed so that the assumptions of Asymptotic Modal Analysis (AMA) were satisfied for certain bandwidths and center frequencies. Measurements of sound pressure levels at points along the boundaries and incrementally into tbe interior were taken. These were compared with the theoretical results predicted with AMA, and found to be in good agreement, particularly for moderate (1/3 octave) bandwidths and sufficiently high center frequencies. Sound pressure level measurements were also taken well into the cavity interior at various points along the 5 totally rigid walls. The AMA theory, including boundary intensification effects, was shown to be accurate provided the assumption of large number of acoustic modes is satisfied, and variables such as power spectra of the wall acceleration, frequency, and damping are slowly varying in the frequency of bandwidth.

Peretti, L. F.↗

NASA NTRS1972

Subsonic, transonic, and supersonic nozzle flow by the inverse technique.

The inverse technique is used to obtain a mathematically and physically consistent solution of the flowfield in a nozzle from the mass generation surface through the supersonic region. The inverse method employs an assumed centerline function which is of the Cauchy type in that the values and the derivatives of the function are known. Since the Cauchy boundary conditions can give rise to numerical instabilities, the governing gasdynamic equations for rotational steady flow were transformed into a form which puts the geometry into a rectangular shape, and which spaces the network of interior points more finely in regions of the greatest gradients of the dependent variables. For arbitrarily specified centerline data, the solution of the governing flow equations may not exist, and if it does it may not depend continuously on the data.

Norton, D. J.↗

NASA NTRS1988

An upwind differencing scheme for the time-accurate incompressible Navier-Stokes equations

The two-dimensional incompressible Navier-Stokes equations are solved in a time-accurate manner in using the method of pseudocompressibility. Using this method, subiterations in pseudotime are required to satisfy the continuity equation at each time step. An upwind differencing scheme based on flux-difference splitting is used to compute the convective terms. The upwind differencing is biased based on the sign of the local eigenvalue of the Jacobian matrix. Third-order or fifth-order spatial accuracy is maintained throughout the interior grid points. The equations are solved using an implicit line-relaxation scheme. This solution scheme is stable and is capable of running at large time steps in pseudotime, leading to fast convergence for each physical time step. A variety of computed results are presented to validate the present scheme. Results for the flow over an oscillating plate are compared with the exact analytic solution, good agreement is seen. Excellent comparison is obtained between the computed solution and the analytical results for inviscid channel flow with an oscillating back pressure. Flow solutions over a circular cylinder with vortex shedding are also presented. Finally, the flow past an airfoil at -90 deg angle-of-attack is also computed.

Rogers, Stuart E.↗

NASA NTRS1983

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.↗

DOE OSTI2023

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↗

DOE OSTI2023

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↗

DOE OSTI2024

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↗

NASA NTRS1977

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.↗

NASA NTRS1985

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.↗

NASA NTRS1992

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↗

NASA NTRS1978

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.↗

NASA NTRS1992

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↗

NASA NTRS2012

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↗

NASA NTRS1987

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↗

NASA NTRS1988

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.↗

NASA NTRS1977

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.↗

NASA NTRS1985

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.↗

NASA NTRS2016

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.↗