Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variable Geometry”

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 1,441 records · Page 80

Optimally robust redundancy relations for failure detection in uncertain systems

All failure detection methods are based, either explicitly or implicitly, on the use of redundancy, i.e. on (possibly dynamic) relations among the measured variables. The robustness of the failure detection process consequently depends to a great degree on the reliability of the redundancy relations, which in turn is affected by the inevitable presence of model uncertainties. In this paper the problem of determining redundancy relations that are optimally robust is addressed in a sense that includes several major issues of importance in practical failure detection and that provides a significant amount of intuition concerning the geometry of robust failure detection. A procedure is given involving the construction of a single matrix and its singular value decomposition for the determination of a complete sequence of redundancy relations, ordered in terms of their level of robustness. This procedure also provides the basis for comparing levels of robustness in redundancy provided by different sets of sensors.

Lou, X.-C.↗

Quasi-perpendicular/quasi-parallel divisions of earth's bow shock

Computer-drawn diagrams of the boundaries between quasi-perpendicular and quasi-parallel areas of earth's bow shock are displayed for a few selected cone angles of static interplanetary magnetic field (IMF). The effect on the boundary of variable IMF in the foreshock is also discussed and shown for one nominal case. The boundaries demand caution in applying them to the realistic, dynamic conditions of the solar wind and in interpreting the effects of small cone angles on the distributions of structures at the shock. However, the calculated, first-order boundaries are helpful in defining areas of the shock where contributions from active structures inherent in quasi-parallel geometry may be distinguishable from those derived secondarily from upstream reflected ion dynamics. The boundaries are also compatible with known behavior of daytime ULF geomagnetic waves and pulsations according to models postulating that cone angle-controlled, time-dependent ULF activity around the subsolar point of the bow shock provides the source of geomagnetic excitation.

Greenstadt, E. W.↗

Kalman Filtering and RTS Smoothing for Arc-Jet Sample Edge Tracking

Accurate and temporally consistent measurements of test-article recession are required to characterize material response during arc-jet ablation experiments. However, image-based boundary measurements are often affected by segmentation noise, brightness variations, and frame-to-frame variability. This work extends arcjetCV [1] with filtering methods for tracking the evolving boundaries of hemispherical and ISO-Q test articles. Two boundary-tracking approaches based on Kalman filtering [2] were implemented. The first applies a point-wise Kalman filter followed by a Rauch–Tung–Striebel smoother [3] to individual boundary-point locations. The second applies the same filtering and smoothing framework to a reduced set of geometry-dependent shape parameters. The point-wise method reduces local frame-to-frame fluctuations while preserving spatial details along the detected boundary. The shape-parameter method provides a compact and geometrically constrained estimate of the sample contour. Figure 1 illustrates the two approaches for a hemispherical test article. Both methods improve temporal consistency and support more robust estimation of surface recession from arc-jet imagery. The two filtering approaches provide complementary representations of boundary evolution and improve the reliability of image-based recession measurements during arc-jet experiments.

recession measurement↗

Optimization of a High Pressure Turbine Blade and Sector-Based Annular Rig Design for Supercritical CO2 Power Cycle Representative Testing

Abstract As part of the ongoing research into the design of hardware for zero emission cycles, a first-stage high-pressure turbine (HPT) blade is optimized for a 300 MWe supercritical CO2 (sCO2) power cycle using the surrogate-assisted genetic algorithm optimizer in Numeca FINE/Design three-dimensional with objectives of increasing efficiency and decreasing heat load to the blade. Supercritical CO2 property tables are constructed from NIST REFPROP data for the condensable gas simulation in FINE/Turbo. A detailed mesh sensitivity study is performed for a baseline design to identify the proper-grid refinement and efficiently allocate resources for the optimization. Seventy design variables are selected for the initial population generation. Self-organizing maps are then used to focus the design variables on the most important ones affecting the objective functions. The optimization results in approximately 3000 three-dimensional Reynolds Averaged Navier Stokes simulations of different blade shapes with increases in efficiency of up to 0.85% and decreases in heat load of 14%. Families of blade shapes are identified for experimental testing in an annular rig at the Purdue Experimental Turbine Aerothermal Laboratory. A design to adapt the annular cascade for testing optimized geometries is introduced, which features eccentric radius sectors allowing for scaled-up geometries of sCO2 optimized blade profiles to be tested at design cycle representative conditions at high Reynolds numbers in dry air. Analysis into the effects of Reynolds number, working fluid, and geometric relations are presented to prove the efficacy of the test method.

Engineering↗

Investigation of the Effect of Framework Flexibility on Adsorption in SIFSIX-3-Cu using a Machine-Learned Force Field

Metal-organic frameworks (MOFs) are a promising class of adsorbents. The performance of MOF sorbents relies on high selectivity and low regeneration energy. This work focuses on the use of machine learned force fields (MLFFs) to model adsorption in a flexible MOF, SIFSIX-3-Cu. A DeePMD-based MLFF was trained to reproduce DFT (PBE+D3) energies, forces, and stresses, using an iterative sampling scheme combining sampling based on molecular dynamics, Monte Carlo, and geometry optimization to capture both attractive and repulsive regions of the potential energy surface. Flexibility of the MOF was explicitly included in this model. Hybrid Monte Carlo/molecular dynamics (MC/MD) simulations using the MLFF predicted adsorption isotherms in good agreement with experimental data for a range of pressures (40 Pa – 104 Pa) in contrast to rigid models, which overpredict CO2 adsorption at low pressures. The improvement was the result of a description of the variability of fluorine-fluorine diagonal distances at adsorption sites. This detailed description of flexibility afforded by the MLFF resulted in more accurate predictions adsorption isotherms when compared to the experimentally measured values. These results underscore the importance of including framework flexibility when modeling adsorption phenomena in MOFs, particularly for low pressure applications and provide a robust procedure for training MLFF models for MOFs.

Atomistic Simulation↗

Thermal Comptonization and Disk Thermal Reprocessing in NGC3516

We present an application of the thermal Comptonization/disk reprocessing model recently proposed by Zdziarski, Lubiniski, and Smith. We show that the absence of strong optical variations in the presence of strong concurrent X-ray variations, similar to those found by Hubble Space Telescope (HST)/Rossi X-Ray Timing Explorer (RXTE) monitoring observations of NGC 3516, can be explained by changing the geometry of the Comptonizing plasma rather than the accretion disk itself. The total X-ray luminosity of the Comptonizing plasma must decrease as its spatial extent increases. In contrast, the disk inner radius must be roughly fixed in order not to produce optical/ultraviolet color variations stronger than observed. By including emission due to internal viscous dissipation in the disk, we can roughly match the optical and X-ray flux levels and variability amplitudes seen from NGC 3516 during the HST/RXTE campaign.

Chiang, James↗

Design Study of Wafer Seals for Future Hypersonic Vehicles

Future hypersonic vehicles require high temperature, dynamic seals in advanced hypersonic engines and on the vehicle airframe to seal the perimeters of movable panels, flaps, and doors. Current seals do not meet the demanding requirements of these applications, so NASA Glenn Research Center is developing improved designs to overcome these shortfalls. An advanced ceramic wafer seal design has shown promise in meeting these needs. Results from a design of experiments study performed on this seal revealed that several installation variables played a role in determining the amount of leakage past the seals. Lower leakage rates were achieved by using a tighter groove width around the seals, a higher seal preload, a tighter wafer height tolerance, and a looser groove length. During flow testing, a seal activating pressure acting behind the wafers combined with simulated vibrations to seat the seals more effectively against the sealing surface and produce lower leakage rates. A seal geometry study revealed comparable leakage for full-scale wafers with 0.125 and 0.25 in. thicknesses. For applications in which lower part counts are desired, fewer 0.25-in.-thick wafers may be able to be used in place of 0.125-in.-thick wafers while achieving similar performance. Tests performed on wafers with a rounded edge (0.5 in. radius) in contact with the sealing surface resulted in flow rates twice as high as those for wafers with a flat edge. Half-size wafers had leakage rates approximately three times higher than those for full-size wafers.

Dunlap, Patrick H.↗

A computer code for multiphase all-speed transient flows in complex geometries. MAST version 1.0

The operation of the MAST code, which computes transient solutions to the multiphase flow equations applicable to all-speed flows, is described. Two-phase flows are formulated based on the Eulerian-Lagrange scheme in which the continuous phase is described by the Navier-Stokes equation (or Reynolds equations for turbulent flows). Dispersed phase is formulated by a Lagrangian tracking scheme. The numerical solution algorithms utilized for fluid flows is a newly developed pressure-implicit algorithm based on the operator-splitting technique in generalized nonorthogonal coordinates. This operator split allows separate operation on each of the variable fields to handle pressure-velocity coupling. The obtained pressure correction equation has the hyperbolic nature and is effective for Mach numbers ranging from the incompressible limit to supersonic flow regimes. The present code adopts a nonstaggered grid arrangement; thus, the velocity components and other dependent variables are collocated at the same grid. A sequence of benchmark-quality problems, including incompressible, subsonic, transonic, supersonic, gas-droplet two-phase flows, as well as spray-combustion problems, were performed to demonstrate the robustness and accuracy of the present code.

Chen, C. P.↗

Simulation of shading and algal growth in experimental raceways

Many algae research groups use elevated experimental raceways to characterize algal biomass productivity. Side walls and a relatively large paddlewheel shade the culture and lower productivity, particularly in winter and with low culture depth. This paper presents a methodology to calculate shading for a given raceway shape: (1) develop a mesh of the inner surface geometry of the raceway; (2) offset the mesh horizontally in the x and y directions on the water surface as a function of wall height, solar zenith angle, and solar azimuth angle; (3) calculate the shaded area for each projected mesh with the shoelace algorithm; (4) use regression to develop a response surface for shaded area as a function of horizontal x and y offset. For each time step, the shaded area is calculated with the regression equation based on the x and y offset at that time step. The shading model was added to the Huesemann Algae Biomass Growth (HABG) model, which predicts algae growth rate as a function of water temperature and light intensity in algae culture layers. The shading model was evaluated with biomass concentrations from three Regional Algal Feedstock Testbed (RAFT) experiments with three different species of algae. The shading algorithm reduced the HABG estimate of productivity and improved the agreement with observed productivities. In the model, productivity was sensitive to the assumed light distribution below the water surface. Averaging light in each layer resulted in higher calculated growth rate than calculating growth separately in shaded and unshaded areas. Shaded area was sensitive to raceway orientation, water depth, and time of year. A north-south raceway orientation had less average shading in winter but more average shading in summer than an east-west orientation.

59 BASIC BIOLOGICAL SCIENCES↗

A Four‐Coordinate Pr 4+ Imidophosphorane Complex

Abstract The imidophosphorane ligand, [NP t Bu 3 ] − ( t Bu= tert ‐butyl), enables isolation of a pseudo‐tetrahedral, tetravalent praseodymium complex, [Pr 4+ (NP t Bu 3 ) 4 ] ( 1‐Pr ), which is characterized by a suite of physical characterization methods including single‐crystal X‐ray diffraction, electron paramagnetic resonance, and L 3 ‐edge X‐ray near‐edge spectroscopies. Variable‐temperature direct‐current magnetic susceptibility data, supported by multiconfigurational quantum chemical calculations, demonstrate that the electronic structure diverges from the isoelectronic Ce 3+ analogue, driven by increased crystal field. The four‐coordinate environment around Pr 4+ in 1‐Pr , which is unparalleled in reported extended solid systems, provides a unique opportunity to study the interplay between crystal field splitting and spin‐orbit coupling in a molecular tetravalent lanthanide within a pseudo‐tetrahedral coordination geometry.

Boggiano, Andrew C.↗

A Four‐Coordinate Pr 4+ Imidophosphorane Complex

Abstract The imidophosphorane ligand, [NP t Bu 3 ] − ( t Bu= tert ‐butyl), enables isolation of a pseudo‐tetrahedral, tetravalent praseodymium complex, [Pr 4+ (NP t Bu 3 ) 4 ] ( 1‐Pr ), which is characterized by a suite of physical characterization methods including single‐crystal X‐ray diffraction, electron paramagnetic resonance, and L 3 ‐edge X‐ray near‐edge spectroscopies. Variable‐temperature direct‐current magnetic susceptibility data, supported by multiconfigurational quantum chemical calculations, demonstrate that the electronic structure diverges from the isoelectronic Ce 3+ analogue, driven by increased crystal field. The four‐coordinate environment around Pr 4+ in 1‐Pr , which is unparalleled in reported extended solid systems, provides a unique opportunity to study the interplay between crystal field splitting and spin‐orbit coupling in a molecular tetravalent lanthanide within a pseudo‐tetrahedral coordination geometry.

Boggiano, Andrew C.↗

Fully implicit crystal plasticity models representing orientations with modified Rodrigues parameters

Here, this work describes a crystal plasticity formulation combining several mathematical, numerical, and implementation choices to produce a highly efficient model. Specifically, the key choices in the implementation are (1) representing orientations with modified Rodrigues parameters, (2) implementing a fully coupled implicit time integration for the elastic stretch, the crystal orientations, and the model internal variables, (3) implementing the model in the NEML2 constitutive modeling framework, based on PyTorch, to vectorize the calculations and port the computation to GPUs and other hardware accelerators, and (4) an exact implementation of the consistent tangent matrix, even for arbitrary coupling to other field variables beyond the displacements, like temperature, neutron fluence, etc. The first two features of the model are, to our knowledge, novel. The paper considers each of these choices individually as well as the final model as a whole. This includes a full description of modified Rodrigues parameters, their advantages over other representations of orientations, the mathematical formulae and tools required to implement a model with modified Rodrigues parameters, and a detailed description of the geometry of the space of modified Rodrigues parameters (in an appendix). It also includes a description of a fully implicit time integration scheme for the orientations and the advantages in representing orientations with modified Rodrigues parameters in implementing such a model. The work then assess, via numerical examples, the advantages of fully coupled implicit time integration versus more common decoupled and explicit time integration schemes. These studies demonstrate the computational advantages of fully coupled integration versus other time integration algorithms, though the performance of the competing models depends on the complexity of the underlying single crystal model. The study concludes by demonstrating that the choice of time integration method affects the sharpness of the predicted texture, with explicit methods for integrating the orientations overestimating texture sharpness and implicit methods underestimating texture sharpness.

Crystal plasticity↗

A comparison of single-cycle versus multiple-cycle proof testing strategies

An evaluation of single-cycle and multiple-cycle proof testing (MCPT) strategies for SSME components is described. Data for initial sizes and shapes of actual SSME hardware defects are analyzed statistically. Closed-form estimates of the J-integral for surface flaws are derived with a modified reference stress method. The results of load- and displacement-controlled stable crack growth tests on thin IN-718 plates with deep surface flaws are summarized. A J-resistance curve for the surface-cracked configuration is developed and compared with data from thick compact tension specimens. The potential for further crack growth during large unload/reload cycles is discussed, highlighting conflicting data in the literature. A simple model for ductile crack growth during MCPT based on the J-resistance curve is used to study the potential effects of key variables. The projected changes in the crack size distribution during MCPT depend on the interactions between several key parameters, including the number of proof cycles, the nature of the resistance curve, the initial crack size distribution, the component boundary conditions (load vs. displacement control), and the magnitude of the applied load or displacement. The relative advantages of single-cycle and multiple-cycle proof testing appear to be specific, therefore, to individual component geometry, material, and loading.

Hudak, S. J., Jr.↗

A new form of the Saint-Venant equations for variable topography

The solution stability of river models using the one-dimensional (1D) Saint-Venant equations can be easily undermined when source terms in the discrete equations do not satisfy the Lipschitz smoothness condition for partial differential equations. Although instability issues have been previously noted, they are typically treated as model implementation issues rather than as underlying problems associated with the form of the governing equations. This study proposes a new reference slope form of the Saint-Venant equations to ensure smooth slope source terms and eliminate one source of potential numerical oscillations. It is shown that a simple algebraic transformation of channel geometry provides a smooth reference slope while preserving the correct cross-section flow area and the total Piezometric pressure gradient that drives the flow. The reference slope method ensures the slope source term in the governing equations is Lipschitz continuous while maintaining all the underlying complexity of the real-world geometry. The validity of the mathematical concept is demonstrated with the open-source Simulation Program for River Networks (SPRNT) model in a series of artificial test cases and a simulation of a small urban creek. Validation comparisons are made with analytical solutions and the Hydrologic Engineering Center's River Analysis System (HEC-RAS) model. The new method reduces numerical oscillations and instabilities without requiring ad hoc smoothing algorithms.

54 ENVIRONMENTAL SCIENCES↗

Saturnian kilometric radiation - Statistical properties and beam geometry

An adequate description of the average statistical properties of Saturn's radio emissions is needed for a study of these emissions. Data regarding these properties are presented, and the implications for source location, beaming, and magnetic surface anomalies are discussed. A description is presented of the average properties of the Saturn kilometric radiation (SKR) as observed from two locations by the Voyager 1 planetary radio astronomy (PRA) instrument for a 2-month period centered on the November 12, 1980 encounter. An analysis of the occurrence pattern of SKR as a function of Saturn's rotation phase has shown that SKR occurrence is not continuous, but variable and roughly periodic. The statistical SKR properties obtained strongly constrain possible source locations. Several source locations are possible, but most intriguing is the region at high latitudes near the noon meridian.

Kaiser, M. L.↗

A new class of variable-radii diffraction optics for high-resolution x-ray spectroscopy at the National Ignition Facility (invited)

A new class of crystal shapes has been developed for x-ray spectroscopy of point-like or small (a few mm) emission sources. These optics allow for dramatic improvement in both achievable energy resolution and total throughput of the spectrometer as compared with traditional designs. This class of crystal shapes, collectively referred to as the Variable-Radii Spiral (VR-Spiral), utilize crystal shapes in which both the major and minor radii are variable. A crystal using this novel VR-Spiral shape has now been fabricated for high-resolution Extended X-ray Absorption Fine Structure (EXAFS) experiments targeting the Pb-L 3 (13.0 keV) absorption edge at the National Ignition Facility. The performance of this crystal has been characterized in the laboratory using a microfocus x-ray source, showing that high-resolution high-throughput EXAFS spectra can be acquired using this geometry. Importantly, these successful tests show that the complex three-dimensional crystal shape is manufacturable with the required precision needed to realize the expected performance of better than 5 eV energy resolution while using a 30 mm high crystal. An improved generalized mathematical form for VR-Spiral shapes is also presented allowing improved optimization as compared to the first sinusoidal-spiral based design. Finally, this new formulation allows VR-Spiral spectrometers to be designed at any magnification with optimized energy resolution at all energies within the spectrometer bandwidth.

47 OTHER INSTRUMENTATION↗

3-D Unsteady Model for Be-Steam Reaction in Water-Cooled Ceramic Breeder Blanket

Here, the design of the water-cooled ceramic breeder (WCCB) blanket includes beryllium multiplier layers. At elevated temperature, beryllium reacts with steam in an exothermic reaction producing beryllium oxide and hydrogen. Such situation may occur in WCCB in the case of the rupture of one of the cooling pipes in the blanket module. This process occurs locally in a complex 3-D geometry of the blanket containing several different granular levels and a network of cooling pipes and structural supports. The process is also inherently unsteady since reaction rate depends on concentration of steam and pure beryllium which changes in time. In order to perform the detailed analysis of the process, the model of the reacting flowthrough porous media was developed and introduced into 3-D computational fluid dynamics code. In this model, granular beds are introduced as porous solids simplifying the model geometry and reducing typical mesh size to manageable amount of tens of millions of elements. A reaction rate between solid beryllium and steam is obtained from experimental results, and depends on temperature and concentration of the reactants. The differential equation for beryllium oxide fraction is introduced, allowing obtaining distributions of beryllium oxide in space and time. Multicomponent flow consisting of a homogenous mixture of steam and hydrogen is considered flowing through the porous solid with variable properties. Sink and source terms for steam and hydrogen fractions are determined by local beryllium oxide mass fraction source according to the molar ratios of beryllium steam reaction. The conjugated heat transfer approach is applied to calculate heat transfer in support structures as well as coolant flow, simultaneously with the porous medium steam flow in a blanket's granular beds. The model is validated using experimental data on beryllium steam reaction for granular bed samples.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An Adaptively-Refined, Cartesian, Cell-Based Scheme for the Euler and Navier-Stokes Equations

A Cartesian, cell-based scheme for solving the Euler and Navier-Stokes equations in two dimensions is developed and tested. Grids about geometrically complicated bodies are generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, polygonal 'cut' cells are created. The geometry of the cut cells is computed using polygon-clipping algorithms. The grid is stored in a binary-tree data structure which provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite-volume formulation. The convective terms are upwinded, with a limited linear reconstruction of the primitive variables used to provide input states to an approximate Riemann solver for computing the fluxes between neighboring cells. A multi-stage time-stepping scheme is used to reach a steady-state solution. Validation of the Euler solver with benchmark numerical and exact solutions is presented. An assessment of the accuracy of the approach is made by uniform and adaptive grid refinements for a steady, transonic, exact solution to the Euler equations. The error of the approach is directly compared to a structured solver formulation. A non smooth flow is also assessed for grid convergence, comparing uniform and adaptively refined results. Several formulations of the viscous terms are assessed analytically, both for accuracy and positivity. The two best formulations are used to compute adaptively refined solutions of the Navier-Stokes equations. These solutions are compared to each other, to experimental results and/or theory for a series of low and moderate Reynolds numbers flow fields. The most suitable viscous discretization is demonstrated for geometrically-complicated internal flows. For flows at high Reynolds numbers, both an altered grid-generation procedure and a different formulation of the viscous terms are shown to be necessary. A hybrid Cartesian/body-fitted grid generation approach is demonstrated. In addition, a grid-generation procedure based on body-aligned cell cutting coupled with a viscous stensil-construction procedure based on quadratic programming is presented.

Coirier, William John↗