Search NASA⌕ Search

SEARCH · Search NASA

Results for “consistent discretization”

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 397 records · Page 22

Hybrid equations of motion for flexible multibody systems using quasi-coordinates

A variety of engineering systems, such as automobiles, aircraft, rotorcraft, robots, spacecraft, etc., can be modeled as flexible multibody systems. The individual flexible bodies are in general characterized by distributed parameters. In most earlier investigations they were approximated by some spatial discretization procedure, such as the classical Rayleigh-Ritz method or the finite element method. This paper presents a mathematical formulation for distributed-parameter multibody systems consisting of a set of hybrid (ordinary and partial) differential equations of motion in terms of quasi-coordinates. Moreover, the equations for the elastic motions include rotatory inertia and shear deformation effects. The hybrid set is cast in state form, thus making it suitable for control design.

Meirovitch, L.↗

Calculation of fully three-dimensional separated flow with an unsteady viscous-inviscid interaction method

Three previous papers have shown that the viscous inviscid numerical methods were capable of calculating separated turbulent flows. The 'Semi-inverse method' and the models of the first author allow the computation of massive separations, stalled flows, and shock wave-boundary layer interactions, in two-dimensional or quasi-three-dimensional conditions, from low speeds to supersonic speeds. The 'Semi-implicit' method for time-consistent coupling allows for the computation of time-accurate transonic separated flow and buffer computations in two-dimensions. The present paper shows that the viscous-inviscid interaction approach is also able to compute the fully three-dimensional flow separation phenomena. The method is based on a thin-layer approximation of the theory of 'Defect-Formulation' that provides the viscous-inviscid splitting of the Navier-Stokes equations. A parametric analytical modelling of the 3D-turbulent velocity profiles is involved. Numerically, the 3D-velocity profiles are discretized in the normal z-direction and driven by parametric integral 3D-equations in direct or inverse modes in the x-direction. The viscous-inviscid coupling is fully 3D and solved the time-consistency problem with an extension of the 'Semi-implicit' method previously suggested in two-dimensions. A 3D inviscid subroutine with TSP approximation is used. Results are obtained for transonic steady flows over wings with shock-induced transonic separation. The method provides results for highly three-dimensional flow separations, such as induced by a three-dimensional through at the wall of a flat plate. The 3D viscous-inviscid coupling and the 3D model of the velocity field provide three-dimensional instantaneous skin-friction lines whose patterns exhibit the same complex topology as Navier-Stokes solvers, with foci, nodes, and saddle-points.

Leballeur, J. C.↗

Calculations of particle-laden flows by second-moment closures

A differential second-moment closure is applied to calculate nonswirling turbulent particle-laden jet flow. The closure is effected by extending constant-density single-phase turbulence models to the continuous (carrier) phase of two-phase flow. The dispersed (discrete) phase is treated by stochastic-Lagrangian technique. To reduce the effects of numerical (false) diffusion on the predicted results, a higher order differencing scheme, namely the flux-spline, is employed. In addition, an experimental study is conducted to provide data of good quality, especially near the inlet, for model assessment. The experiment consists of a confined coaxial jet flow of air with glass beads. The beads' diameter is 105 microns and is injected from the inner jet with a mass loading ratio of 0.9. The computed results are in good agreement with experimental data, especially for the carrier phase. However, the predicted turbulence fluctuations associated with the discrete phase are lower than the data at the downstream region of the flow field.

Nikjooy, M.↗

DG-IMEX method for a two-moment model for radiation transport in the $\mathscr{O}$($v$/$c$) limit

Here, we consider neutral particle systems described by moments of a phase-space density and propose a realizability-preserving numerical method to evolve a spectral two-moment model for particles interacting with a background fluid moving with nonrelativistic velocities. The system of nonlinear moment equations, with special relativistic corrections to $\mathscr{O}$($v$/$c$), expresses a balance between phase-space advection and collisions and includes velocity-dependent terms that account for spatial advection, Doppler shift, and angular aberration. The model is conservative for the correct $\mathscr{O}$($v$/$c$) Eulerian-frame number density and is consistent, to $\mathscr{O}$($v$/$c$), with Eulerian-frame energy and momentum conservation. This model is closely related to the one promoted by Lowrie et al. and similar to models currently used to study transport phenomena in large-scale simulations of astrophysical environments. The proposed numerical method is designed to preserve moment realizability, which guarantees that the moments correspond to a nonnegative phase-space density. The realizability-preserving scheme consists of the following key components: (i) a strong stability-preserving implicit-explicit (IMEX) time-integration method; (ii) a discontinuous Galerkin (DG) phase-space discretization with carefully constructed numerical uxes; (iii) a realizability-preserving implicit collision update; and(iv) a realizability-enforcing limiter. In time integration, nonlinearity of the moment model necessitates solution of nonlinear equations, which we formulate as fixed-point problems and solve with tailored iterative solvers that preserve moment realizability with guaranteed global convergence. We also analyze the simultaneous Eulerian-frame number and energy conservation properties of the semi-discrete DG scheme and propose a "spectral redistribution" scheme that promotes Eulerian-frame energy conservation. Through numerical experiments, we demonstrate the accuracy and robustness of this DG-IMEX method and investigate its Eulerian-frame energy conservation properties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Mixed models and reduction method for dynamic analysis of anisotropic shells

A time-domain computational procedure is presented for predicting the dynamic response of laminated anisotropic shells. The two key elements of the procedure are: (1) use of mixed finite element models having independent interpolation (shape) functions for stress resultants and generalized displacements for the spatial discretization of the shell, with the stress resultants allowed to be discontinuous at interelement boundaries; and (2) use of a dynamic reduction method, with the global approximation vectors consisting of the static solution and an orthogonal set of Lanczos vectors. The dynamic reduction is accomplished by means of successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite element method is first used to generate the global approximation vectors. Then the Rayleigh-Ritz technique is used to generate a reduced system of ordinary differential equations in the amplitudes of these modes. The temporal integration of the reduced differential equations is performed by using an explicit half-station central difference scheme (Leap-frog method). The effectiveness of the proposed procedure is demonstrated by means of a numerical example and its advantages over reduction methods used with the displacement formulation are discussed.

Noor, A. K.↗

Euler solutions for blunt bodies using triangular meshes - Artificial viscosity forms and numerical boundary conditions

A finite volume method is used to calculate compressible inviscid flows over blunt bodies using, in general, unstructured grids. Artificial viscosity forms are derived based on a simplified least squares procedure. The extra second order terms are consistent with the governing equations, hence a systematic treatment of the numerical boundary conditions can be easily implemented. A special treatment of blunt bodies may be required. The discrete equations are linearized and the resulting system is solved by a relaxation method. Preliminary results indicate that the effect of the numerical dissipation is minimal. For subsonic flows over smooth bodies, the solution is practically vorticity-free and the total pressure loss is of the same order as the truncation error. Finally, some extensions of the present method are briefly discussed.

Winterstein, R.↗

Watertight Anisotropic Surface Meshing Using Quadrilateral Patches

This paper presents a simple technique for generating anisotropic surface triangulations using unstructured quadrilaterals when the CAD entity can be mapped to a logical rectangle. Watertightness and geometric quality measures are maintained and are consistent with the CAPRI default tessellator. These triangulations can match user specified criteria for chord-height tolerance, neighbor triangle dihedral angle, and maximum triangle side length. This discrete representation has hooks back to the owning geometry and therefore can be used in conjunction with these entities to allow for easy enhancement or modification of the tessellation suitable for grid generation or other downstream applications.

Haimes, Robert↗

Distribution and diagenesis of microfossils from the lower Proterozoic Duck Creek Dolomite, Western Australia

Two distinct generations of microfossils occur in silicified carbonates from a previously undescribed locality of the Lower Proterozoic Duck Creek Dolomite, Western Australia. The earlier generation occurs in discrete organic-rich clasts and clots characterized by microquartz anhedra; it contains a variety of filamentous and coccoidal fossils in varying states of preservation. Second generation microfossils consist almost exclusively of well-preserved Gunflintia minuta filaments that drape clasts or appear to float in clear chalcedony. These filaments appear to represent an ecologically distinct assemblage that colonized a substrate containing the partially degraded remains of the first generation community. The two assemblages differ significantly in taxonomic frequency distribution from previously described Duck Creek florules. Taken together, Duck Creek microfossils exhibit a range of assemblage variability comparable to that found in other Lower Proterozoic iron formations and ferruginous carbonates. With increasing severity of post-mortem alteration, Duck Creek microfossils appear to converge morphologically on assemblages of simple microstructures described from early Archean cherts. Two new species are described: Oscillatoriopsis majuscula and O. cuboides; the former is among the largest septate filamentous fossils described from any Proterozoic formation.

NASA Discipline Number 52-30↗

Discrete Adjoint-Based Design Optimization of Unsteady Turbulent Flows on Dynamic Unstructured Grids

An adjoint-based methodology for design optimization of unsteady turbulent flows on dynamic unstructured grids is described. The implementation relies on an existing unsteady three-dimensional unstructured grid solver capable of dynamic mesh simulations and discrete adjoint capabilities previously developed for steady flows. The discrete equations for the primal and adjoint systems are presented for the backward-difference family of time-integration schemes on both static and dynamic grids. The consistency of sensitivity derivatives is established via comparisons with complex-variable computations. The current work is believed to be the first verified implementation of an adjoint-based optimization methodology for the true time-dependent formulation of the Navier-Stokes equations in a practical computational code. Large-scale shape optimizations are demonstrated for turbulent flows over a tiltrotor geometry and a simulated aeroelastic motion of a fighter jet.

Nielsen, Eric J.↗

A numerical study of the temporal eigenvalue spectrum of the Blasius boundary layer

A numerical study is made of the temporal eigenvalue spectrum of the Orr-Sommerfeld equation for the Blasius boundary layer. Unlike channel flows, there is no mathematical proof that this flow has an infinite spectrum of discrete eigenvalues. The Orr-Sommerfeld equation is integrated numerically, and the eigenvalues located by tracing out the contour lines in the complex wave velocity plane on which the real and imaginary parts of the secular determinant are zero. The spectrum of plane Poiseuille flow is used as a guide to study the spectrum of an artificial two-wall flow which consists of two Blasius boundary layers. As the upper boundary of this flow moves to infinity, it is found that the portion of the spectrum with an infinite number of eigenvalues moves towards phase velocity equal to unity and the spacing between eigenvalues goes to zero. The original few eigenvalues found are the only discrete eigenvalues that exist for Blasius flow.

Mack, L. M.↗

SAM Theory Manual

The System Analysis Module (SAM) is an advanced and modern system analysis tool under development at Argonne National Laboratory for advanced non-LWR reactor safety analysis. It aims to provide fast-running, modest-fidelity, whole-plant transient analyses capabilities, which are essential for fast turnaround design scoping and engineering analyses of advanced reactor concepts. While SAM is being developed as a system-level modeling and simulation tool, advanced modeling techniques being implemented include a reduced-order three-dimensional module, pseudo 3-D conjugate heat transfer modeling in reactor core, flexible and multi-scale modeling of heat transfer between fluid and structures, in addition to the advances in software environments and design, and numerical methods. SAM aims to be a generic system-level safety analysis tool for advanced non-LWRs, including Liquid-Metal-cooled fast Reactors (LMR), Molten Salt Reactors (MSR), Fluoride-salt-cooled High-temperature Reactors (FHR), and High-Temperature Gas-cooled Reactors (HTGR). SAM takes advantage of advances in physical modeling, numerical methods, and software engineering to enhance its user experience and usability. It utilizes an object-oriented computational framework (MOOSE), and its underlying meshing and finite-element library and linear and non-linear solvers, to leverage the modern advanced software environments and numerical methods. This document provides the theoretical and technical basis of the code to help users understand the underlying physical models (such as governing equations, closure models, and component models), system modeling approaches, numerical discretization and solution methods, and the overall capabilities in SAM. As new code capabilities and features are added, the SAM Theory Manual will be updated periodically to keep it consistent with the state of the development.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Exploring the Regimes of Particle Behavior upon Impact via the Discrete Element Method

Discrete element method simulations are conducted to probe the various regimes of post-impact behavior of particles with solid surfaces. The impacting particles are described as spherical agglomerates consisting of smaller constituent (or primary) particles held together via surface adhesion. Under the influence of a wide range of impact velocities and particle surface energies, five distinct behavioral regimes—rebounding, vibration, fragmentation, pancaking, and shattering—are identified, and force transmission patterns are linked to post-impact behavior. In the rebounding regime, the coefficient of restitution decreases linearly as impact velocity increases and the particle agglomerate experiences compaction. In the fragmentation regime, rebound velocity generally decreases with increasing fragment size. The rebound velocity of fragments decreases with time except for the smallest fragments, which can increase in velocity due to collisions with other fragments of high velocity. Particle breakage in the pancaking regime does not follow common mechanistic models of breakage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unpowered Aerodynamic Characteristics of a 15-Percent Scale Model of a Twin-Engine Commuter Aircraft

An experimental investigation was conducted in the Ames 12-Foot Pressure Wind Tunnel to determine the unpowered aerodynamic characteristics of a 15-percent-scale model of a twin-engine commuter aircraft. Model longitudinal aerodynamic characteristics were examined at discrete flap deflections for various angle-of-attack and wind-tunnel-velocity ranges with the empennage on and off. Data are presented for the basic model configuration consisting of the fuselage, wing, basic wing leading edge, double slotted flaps, midengine nacelles, and empennage. Other configurations tested include a particle-span drooped leading edge (dropped outboard of the engine nacelles), a full-span drooped leading edge, low- and high-mounted engine nacelles, and a single-slotted flap. An evaluation was made of the model mounting system by comparing data obtained with the model mounted conventionally on the wind-tunnel model-support struts and the model inverted.

Morgan, D. G.↗

Atmospheric emissions photometric imaging experiment

The atmospheric emissions photometric imaging experiment was flown on Spacelab 1 to study faint natural and artificial atmospheric emission phenomena. The instrument imaged optical emission in the region 2000 to 7500 angstroms with a television system consisting of two optical channels, one wide-angle and one telephoto. A third optical channel imaged onto the photocathode of a microchannel plate photomultiplier tube that has 100 discrete anodes. A hand-held image intensifier camera with an objective grating permitted spectral analysis of the earth's airglow and the Shuttle glow. Preliminary data show magnesium ion emission features in the lower ionosphere as well as the spacecraft glow spectrum.

Mende, S. B.↗

Space station integrated propulsion and fluid systems study. Space station program fluid management systems databook

Commonality and integration of propulsion and fluid systems associated with the Space Station elements are being evaluated. The Space Station elements consist of the core station, which includes habitation and laboratory modules, nodes, airlocks, and trusswork; and associated vehicles, platforms, experiments, and payloads. The program is being performed as two discrete tasks. Task 1 investigated the components of the Space Station architecture to determine the feasibility and practicality of commonality and integration among the various propulsion elements. This task was completed. Task 2 is examining integration and commonality among fluid systems which were identified by the Phase B Space Station contractors as being part of the initial operating capability (IOC) and growth Space Station architectures. Requirements and descriptions for reference fluid systems were compiled from Space Station documentation and other sources. The fluid systems being examined are: an experiment gas supply system, an oxygen/hydrogen supply system, an integrated water system, the integrated nitrogen system, and the integrated waste fluids system. Definitions and descriptions of alternate systems were developed, along with analyses and discussions of their benefits and detriments. This databook includes fluid systems descriptions, requirements, schematic diagrams, component lists, and discussions of the fluid systems. In addition, cost comparison are used in some cases to determine the optimum system for a specific task.

Bicknell, B.↗

Buckling and postbuckling of composite panels with cutouts subjected to combined edge shear and temperature change

The results of a detailed study of the buckling and postbuckling responses of composite panels with central circular cutouts are presented. The panels are subjected to combined edge shear and temperature change. The panels are discretized by using a two-field degenerate solid element with each of the displacement components having a linear variation throughout the thickness of the panel. The fundamental unknowns consist of the average mechanical strains through the thickness and the displacement components. The effects of geometric nonlinearities and laminated anisotropic material behavior are included. The stability boundary, postbuckling response and the hierarchical sensitivity coefficients are evaluated. The hierarchical sensitivity coefficients measure the sensitivity of the buckling and postbuckling responses to variations in the panel stiffnesses, and the material properties of both the individual layers and the constituents (fibers and matrix). Numerical results are presented for composite panels with central circular cutouts subjected to combined edge shear and temperature change, showing the effects of variations in the hole diameter, laminate stacking sequence and fiber orientation, on the stability boundary and postbuckling response and their sensitivity to changes in the various panel parameters.

Noor, Ahmed K.↗

Chandra Observation of the X-ray Source Population of NGC 6946

We present the results of a study of discrete X-ray sources in NGC 6946 using a deep Chandra ACIS observation. Based on the slope of the log N-log S distribution and the general correlation of sources with the spiral arms, we infer that the overall discrete source sample in NGC 6946 is dominated by high mass X-ray binaries, in contrast to the source distributions in M31 and the Milky Way. This is consistent with the higher star formation rate in NGC 6946 than in those galaxies. We find that the strong X-ray sources in the region of the galactic center do not correlate in detail with images of the region in the near-IR, although one of them may be coincident with the galactic center. The non-central ultra-luminous X-ray source in NGC 6946, previously identified with a supernova remnant, has an X-ray spectrum and luminosity that is inconsistent with either a traditional pulsar wind nebula or a blast wave remnant.

Holt, S. S.↗

Catalyst comprising Ir or Ir and Ru for hydrazine decomposition

A catalyst for hydrazine decomposition consisting essentially of a carrier having a pore volume of at least 0.1 cubic centimeters per gram and a specific surface area, measured in square meters per gram, equal to 195 (C.sub.p + 0.013 + 0.736 V.sub.p) where C.sub.p is the specific heat capacity of the carrier at about 25.degree. C in calories per gram per degree and V.sub.p is the pore volume of the carrier in cubic centimeters per gram and metal of the group consisting of iridium, and mixtures consisting of iridium and ruthenium deposited on said carrier in an amount between 20% and about 40% by weight of the catalyst and distributed through the pores thereof in discrete particles sufficiently separated from each other so that they do not sinter or fuse together when the catalyst is at hydrazine decomposition temperature.

Armstrong, Warren E.↗