Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Vector Radiative Transfer Equation for Arbitrarily Shaped and Arbitrarily Oriented Particles: A Microphysical Derivation from Statistical Electromagnetics

We use the concepts of statistical electromagnetics to derive the general radiative transfer equation (RTE) describing multiple scattering of polarized light by sparse discrete random media consisting of arbitrarily shaped and arbitrarily oriented particles. The derivation starts from the volume integral and Lippmann-Schwinger equations for the electric field scattered by a fixed N-particle system and proceeds to the vector form of the Foldy-Lax equations and their approximate far-field version. We then assume that particle positions are completely random and derive the vector RTE by applying the Twersky approximation to the coherent electric field and the Twersky and ladder approximations to the coherency dyad of the diffuse field in the limit N -> infinity. The concluding section discusses the physical meaning of the quantities entering the general vector RTE and the assumptions made in its derivation.

Mishchenko, Michael I.↗
NASA NTRS2011

Design, Development and Test Challenges: Separation Mechanisms for the Orion Pad Abort-1 Flight Test

On May 6, 2010, NASA launched the first successful integrated flight test, Pad Abort-1, of the Orion Project from the White Sands Missile Range in Las Cruces, New Mexico. This test demonstrated the ability to perform an emergency pad abort of a full-scale 4.8 m diameter, 8200 kg crew capsule. During development of the critical separation mechanisms for this flight test, various challenges were overcome related to environments definition, installation complications, separation joint retraction speed, thruster ordnance development issues, load path validation and significant design loads increases. The Launch Abort System retention and release (LAS R&R) mechanism consisted of 6 discrete structural connections between the LAS and the crew module (CM) simulator, each of which had a preloaded tension tie, Superbolt torque-nut and frangible nut. During the flight test, the frangible nuts were pyrotechnically split, permitting the CM to separate from the LAS. The LAS separation event was the driving case in the shock environment for many co-located hardware items. During development testing, it was necessary to measure the source shock during the separation event so the predicted shock environment could be validated and used for certification testing of multiple hardware items. The Lockheed Martin test team measured the source separation shock due to the LAS R&R function, which dramatically decreased the predicted environment by 90% at 100 Hz. During development testing a hydraulic tensioner was used to preload the joint; however, the joint relaxation with the tensioner proved unsatisfactory so the design was modified to include a Superbolt torque-nut. The observed preload creep during lab testing was 4% after 30 days, with 2.5% occurring in the first 24 hours. The conversion of strain energy (preload) to kinetic energy (retraction) was measured to be 50-75%. Design features and careful monitoring of multiple strain gauges on each tension tie allowed a pure tensile load to be applied after stacking at the launch pad. Following installation, preload in each joint was monitored for 24 hours. Due to unforeseen complications and the influence of temperature on the portable data acquisition system, the team encountered difficulty in tracking the joint relaxation. In some cases, bond-line failure of the strain gauges occurred.

Dinsel, Alison↗
NASA NTRS2011

Direct Simulation of Multiple Scattering by Discrete Random Media Illuminated by Gaussian Beams

The conventional orientation-averaging procedure developed in the framework of the superposition T-matrix approach is generalized to include the case of illumination by a Gaussian beam (GB). The resulting computer code is parallelized and used to perform extensive numerically exact calculations of electromagnetic scattering by volumes of discrete random medium consisting of monodisperse spherical particles. The size parameters of the scattering volumes are 40, 50, and 60, while their packing density is fixed at 5%. We demonstrate that all scattering patterns observed in the far-field zone of a random multisphere target and their evolution with decreasing width of the incident GB can be interpreted in terms of idealized theoretical concepts such as forward-scattering interference, coherent backscattering (CB), and diffuse multiple scattering. It is shown that the increasing violation of electromagnetic reciprocity with decreasing GB width suppresses and eventually eradicates all observable manifestations of CB. This result supplements the previous demonstration of the effects of broken reciprocity in the case of magneto-optically active particles subjected to an external magnetic field.

Mackowski, Daniel W.↗
NASA NTRS2012

Radio Detections During Two State Transitions of the Intermediate-Mass Black Hole HLX-1

Relativistic jets are streams of plasma moving at appreciable fractions of the speed of light. They have been observed from stellar-mass black holes (approx. 3 to 20 solar masses) as well as supermassive black holes (approx.. 10(exp 6) to 10(exp 9) Solar Mass) found in the centers of most galaxies. Jets should also be produced by intermediate-mass black holes (approx. 10(exp 2) to 10(exp 5) Solar Mass), although evidence for this third class of black hole has, until recently, been weak. We report the detection of transient radio emission at the location of the intermediate-mass black hole candidate ESO 243-49 HLX-1, which is consistent with a discrete jet ejection event. These observations also allow us to refine the mass estimate of the black hole to be between approx. 9 × 10(exp 3) Solar Mass and approx. 9 × 10(exp 4) Solar Mass.

transitions↗
NASA NTRS2018

Unstructured Grid Adaptation and Solver Technology for Turbulent Flows

Unstructured grid adaptation is a tool to control Computational Fluid Dynamics (CFD) discretization error. However, adaptive grid techniques have made limited impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Issues that prevent the use of adaptive grid methods are identified by applying unstructured grid adaptation methods to a series of benchmark cases. Once identified, these challenges to existing adaptive workflows can be addressed. Unstructured grid adaptation is evaluated for test cases described on the Turbulence Modeling Resource (TMR) web site, which documents uniform grid refinement of multiple schemes. The cases are turbulent flow over a Hemisphere Cylinder and an ONERA M6Wing. Adaptive grid force and moment trajectories are shown for three integrated grid adaptation processes with Mach interpolation control and output error based metrics. The integrated grid adaptation process with a finite element (FE) discretization produced results consistent with uniform grid refinement of fixed grids. The integrated grid adaptation processes with finite volume schemes were slower to converge to the reference solution than the FE method. Metric conformity is documented on grid/metric snapshots for five grid adaptation mechanics implementations. These tools produce anisotropic boundary conforming grids requested by the adaptation process.

Park, Michael A.↗
NASA NTRS2019

Adjoint-Based Mesh Adaptation and Shape Optimization for Simulations with Propulsion

We demonstrate a well-posed formulation of permeable boundary conditions and mass- flow-rate functionals for adjoint-based mesh refinement and shape optimization governed by the steady Euler equations. The boundary conditions are used to model propulsion- system effects of inlets and nozzles. A two-shock diffuser with an analytic solution is used to verify the implementation. Numerical examples show that the adjoint solution is smooth at the boundary, indicating that the discretization is adjoint consistent when exit pressure is specified at subsonic outflow, and stagnation temperature and pressure at subsonic inflow. The results focus on improving simulation techniques for low-boom aircraft analysis and design. By including mass-flow-rate outputs, we obtain reliable estimates of engine flow rates concurrently with nearfield pressure signatures without increasing simulation cost. We also demonstrate the importance of mass-flow-rate constraints in shape optimization by examining trade-offs between maximizing performance of a shrouded supersonic nozzle and minimizing shocks in its nearfield.

Nemec, Marian↗
NASA NTRSDate not supplied

Wall-Modeled Large Eddy Simulation Method for Unstructured-Grid Navier-Stokes Solvers

This paper reports on the implementation and assessment of a Wall-Modeled Large-Eddy Simulation (WMLES) methodology in an unstructured-grid, node-centered flow solver, FUN3D that is developed and supported at the NASA Langley Research Center. Finite-volume (FV) and finite-element (FE) discretization schemes considered in the study provide formal second-order spatial accuracy. Large-Eddy Simulations (LES) resolve large-scale turbulent-flow features and filter out small-scale effects using the Vreman subgrid-scale model. At solid-wall boundaries, a shear-stress model is employed to provide a proper boundary-flux closure. The nonlinear equations are integrated in time using either an optimized backward difference formula or an implicit multistage Runge-Kutta temporal scheme. The implicit equations at each time step are solved by strong nonlinear iteration schemes. WMLES demonstrations are shown for two high-lift configurations, namely, the McDonnell Douglas 30P30N multielement airfoil and a NASA High-Lift Common Research Model. Results show that the WMLES approaches implemented in the FV and FE discretization methods produce consistent solutions and are capable of capturing key aerodynamic characteristics and flow structures for high-lift configurations at a wide range of angles of attack including maximum-lift conditions. In the 30P30N example, correct trends in the variations of integrated aerodynamic forces and moments, surface pressure distributions, and boundary-layer profiles are captured as the Reynolds number is increased.

CFD; turbulence modeling; High-Lift flow simulatio↗
NASA NTRSDate not supplied

Embedding Differential Dynamic Logic in PVS

Runtime assurance is a control framework where a complex controller operates under the observation of a monitor. If the monitor detects the controller exhibiting undesirable behavior, control is passed off to a trusted controller until a desirable state is regained. The runtime assurance architecture provides a layer of assurance to the system being controlled, but special care must be taken that the resulting overall system, consisting of the monitors and controllers, is behaving as intended. This talk aims to formally model and reason about runtime assurance-equipped systems as hybrid programs- which are models that consist of both discrete and continuous components. Using the verification tool Plaidypvs, safety properties of some examples involving RTA architectures is shown.

Formal Verification↗
NASA NTRSDate not supplied

Computational Modeling of Two Mars Powered Descent Vehicle Concepts Tested in the Langley Unitary Plan Wind Tunnel

Future human Mars missions will require powered descent starting at supersonic conditions, something which has never been done before at Mars. Computational powered descent flowfield simulations have been completed at full-scale Mars conditions, but the available ground test data are not suitable for calibrating computational uncertainties for aerodynamic interference on proposed Mars descent vehicles. Testing was conducted in the NASA Langley Unitary Plan Wind Tunnel in order to investigate the aerodynamic interference of sub-scale versions of two Mars powered descent vehicle concepts at supersonic Mach numbers (2.4 and 3.5): a model based on a blunt hypersonic inflatable aerodynamic decelerator (HIAD) and the second representing a more slender rigid vehicle with body flaps (CobraMRV). This paper covers computational flowfield predictions completed at wind tunnel conditions and comparisons to the test data. On the blowing HIAD models, the time-averaged pressure inboard of the nozzles was generally well-predicted, especially if the nozzles are canted outward, when the nozzles are located further from the nose. At intermediate CobraMRV thrust coefficients, CFD pressures are more accurately predicted than they are for the HIAD models, largely due to the nozzle locations and pointing directions. Overall, the CFD pressure coefficients were predicted within 0.2 of the steady pressure measurements for all blowing models, with smaller discrepancies at higher HIAD thrust, and larger discrepancies at lower HIAD CobraMRV thrust. All HIAD models were predicted to have a gradually decreasing axial force coefficient as the total thrust increases, in agreement with available pressure sensitive paint data. On models with canted nozzles or with nozzles further from the nose, the force coefficient was slightly higher for a given thrust. On the CobraMRV model, the CFD also shows consistent results between solvers and follows trends revealed in the data; the aerodynamic force coefficient remains near the non-blowing value at a tunnel Mach number of 2.4 regardless of thrust, and increases above that level at a Mach number of 3.5, consistent with the discrete pressure data. CFD analysis at tunnel and flight conditions will continue as flight system designs concepts mature.

Supersonic Retropropulsion↗
NASA NTRSDate not supplied

Computational Modeling of Two Mars Powered Descent Vehicle Concepts Tested in the Langley Unitary Plan Wind Tunnel

Future human Mars missions will require powered descent starting at supersonic conditions, something which has never been done before at Mars. Computational powered descent flowfield simulations have been completed at full-scale Mars conditions, but the available ground test data are not suitable for calibrating computational uncertainties for aerodynamic interference on proposed Mars descent vehicles. Testing was conducted in the NASA Langley Unitary Plan Wind Tunnel in order to investigate the aerodynamic interference of sub-scale versions of two Mars powered descent vehicle concepts at supersonic Mach numbers (2.4 and 3.5): a model based on a blunt hypersonic inflatable aerodynamic decelerator (HIAD) and the second representing a more slender rigid vehicle with body flaps (CobraMRV). This paper covers computational flowfield predictions completed at wind tunnel conditions and comparisons to the test data. On the blowing HIAD models, the time-averaged pressure inboard of the nozzles was generally well-predicted, especially if the nozzles are canted outward, when the nozzles are located further from the nose. At intermediate CobraMRV thrust coefficients, CFD pressures are more accurately predicted than they are for the HIAD models, largely due to the nozzle locations and pointing directions. Overall, the CFD pressure coefficients were predicted within 0.2 of the steady pressure measurements for all blowing models, with smaller discrepancies at higher HIAD thrust, and larger discrepancies at lower HIAD CobraMRV thrust. All HIAD models were predicted to have a gradually decreasing axial force coefficient as the total thrust increases, in agreement with available pressure sensitive paint data. On models with canted nozzles or with nozzles further from the nose, the force coefficient was slightly higher for a given thrust. On the CobraMRV model, the CFD also shows consistent results between solvers and follows trends revealed in the data; the aerodynamic force coefficient remains near the non-blowing value at a tunnel Mach number of 2.4 regardless of thrust, and increases above that level at a Mach number of 3.5, consistent with the discrete pressure data. CFD analysis at tunnel and flight conditions will continue as flight system designs concepts mature.

Supersonic Retropropulsion↗
NASA NTRS1992

A comprehensive X-ray study of the Small Magellanic Cloud

A comprehensive analysis of the Einstein data base on the SMC is presented. Twenty-four soures are classified as objects in the SMC, and 13 as Galactic stars, AGN, and clusters of galaxies. The detected discrete emitters account for less than half of the excess emission associated with the SMC. The extensive diffuse emission has an X-ray spectrum substantially softer than that of the discrete sources, which is consistent with the presence of X-ray-emitting coronal gas associated with the SMC. A detailed study of the discrete sources is presented. An X-ray source is found to vary in its X-ray luminosity by a factor larger than 10 over a period of about 1 yr and has a relatively hard X-ray spectrum which apparently changes with the X-ray luminosity. It is suggested that this source, associated with a B1-type star, is a new high-mass X-ray binary in the Small Cloud.

Wang, Qingde↗
NASA NTRS1996

A PDE Sensitivity Equation Method for Optimal Aerodynamic Design

The use of gradient based optimization algorithms in inverse design is well established as a practical approach to aerodynamic design. A typical procedure uses a simulation scheme to evaluate the objective function (from the approximate states) and its gradient, then passes this information to an optimization algorithm. Once the simulation scheme (CFD flow solver) has been selected and used to provide approximate function evaluations, there are several possible approaches to the problem of computing gradients. One popular method is to differentiate the simulation scheme and compute design sensitivities that are then used to obtain gradients. Although this black-box approach has many advantages in shape optimization problems, one must compute mesh sensitivities in order to compute the design sensitivity. In this paper, we present an alternative approach using the PDE sensitivity equation to develop algorithms for computing gradients. This approach has the advantage that mesh sensitivities need not be computed. Moreover, when it is possible to use the CFD scheme for both the forward problem and the sensitivity equation, then there are computational advantages. An apparent disadvantage of this approach is that it does not always produce consistent derivatives. However, for a proper combination of discretization schemes, one can show asymptotic consistency under mesh refinement, which is often sufficient to guarantee convergence of the optimal design algorithm. In particular, we show that when asymptotically consistent schemes are combined with a trust-region optimization algorithm, the resulting optimal design method converges. We denote this approach as the sensitivity equation method. The sensitivity equation method is presented, convergence results are given and the approach is illustrated on two optimal design problems involving shocks.

Borggaard, Jeff↗
NASA NTRS2009

A Critical Study of Agglomerated Multigrid Methods for Diffusion

Agglomerated multigrid techniques used in unstructured-grid methods are studied critically for a model problem representative of laminar diffusion in the incompressible limit. The studied target-grid discretizations and discretizations used on agglomerated grids are typical of current node-centered formulations. Agglomerated multigrid convergence rates are presented using a range of two- and three-dimensional randomly perturbed unstructured grids for simple geometries with isotropic and highly stretched grids. Two agglomeration techniques are used within an overall topology-preserving agglomeration framework. The results show that multigrid with an inconsistent coarse-grid scheme using only the edge terms (also referred to in the literature as a thin-layer formulation) provides considerable speedup over single-grid methods but its convergence deteriorates on finer grids. Multigrid with a Galerkin coarse-grid discretization using piecewise-constant prolongation and a heuristic correction factor is slower and also grid-dependent. In contrast, grid-independent convergence rates are demonstrated for multigrid with consistent coarse-grid discretizations. Actual cycle results are verified using quantitative analysis methods in which parts of the cycle are replaced by their idealized counterparts.

Thomas, James L.↗
NASA NTRS2011

A Critical Study of Agglomerated Multigrid Methods for Diffusion

Agglomerated multigrid techniques used in unstructured-grid methods are studied critically for a model problem representative of laminar diffusion in the incompressible limit. The studied target-grid discretizations and discretizations used on agglomerated grids are typical of current node-centered formulations. Agglomerated multigrid convergence rates are presented using a range of two- and three-dimensional randomly perturbed unstructured grids for simple geometries with isotropic and stretched grids. Two agglomeration techniques are used within an overall topology-preserving agglomeration framework. The results show that multigrid with an inconsistent coarse-grid scheme using only the edge terms (also referred to in the literature as a thin-layer formulation) provides considerable speedup over single-grid methods but its convergence deteriorates on finer grids. Multigrid with a Galerkin coarse-grid discretization using piecewise-constant prolongation and a heuristic correction factor is slower and also grid-dependent. In contrast, grid-independent convergence rates are demonstrated for multigrid with consistent coarse-grid discretizations. Convergence rates of multigrid cycles are verified with quantitative analysis methods in which parts of the two-grid cycle are replaced by their idealized counterparts.

Nishikawa, Hiroaki↗
NASA NTRS1991

Discretized Newton-relaxation solution of high resolution flux-difference split schemes

A method is presented for the implicit solution of high resolution flux-difference split schemes for which it is impractical to obtain the Jacobian matrix analytically. The method is used frequently in numerical schemes for solving nonlinear systems of equations, but apparently has not been used to date for approximating the Jacobian of the Roe numerical flux vectors. The method consists of obtaining a discretized Jacobian by using simple finite-differences of the flux vectors, and then using this discretized Jacobian in a Newton-relaxation solution method. The method is applied to both compressible and incompressible flows. Numerical results are presented that demonstrate the method works for both inviscid and viscous flows.

Whitfield, David L.↗
NASA NTRS1981

Scattering of electromagnetic waves from random media with strong permittivity fluctuations

By taking into account the singularity of the dyadic Green's function in the renormalization method, a theory is derived for vector electromagnetic wave propagation in a random medium with large permittivity fluctuations and with anisotropic correlation function. The strong fluctuation theory is then applied to a discrete scatterer problem in which the permittivity can assume only two values. The results are found to be consistent with those derived from discrete scatterer theory for all values of dielectric constants of the scatterers.

Tsang, L.↗
NASA NTRS1987

Particle simulation of plasmas on the massively parallel processor

Particle simulations, in which collective phenomena in plasmas are studied by following the self consistent motions of many discrete particles, involve several highly repetitive sets of calculations that are readily adaptable to SIMD parallel processing. A fully electromagnetic, relativistic plasma simulation for the massively parallel processor is described. The particle motions are followed in 2 1/2 dimensions on a 128 x 128 grid, with periodic boundary conditions. The two dimensional simulation space is mapped directly onto the processor network; a Fast Fourier Transform is used to solve the field equations. Particle data are stored according to an Eulerian scheme, i.e., the information associated with each particle is moved from one local memory to another as the particle moves across the spatial grid. The method is applied to the study of the nonlinear development of the whistler instability in a magnetospheric plasma model, with an anisotropic electron temperature. The wave distribution function is included as a new diagnostic to allow simulation results to be compared with satellite observations.

Gledhill, I. M. A.↗
NASA NTRS1994

The development of a solid-state hydrogen sensor for rocket engine leakage detection

Hydrogen propellant leakage poses significant operational problems in the rocket propulsion industry as well as for space exploratory applications. Vigorous efforts have been devoted to minimizing hydrogen leakage in assembly, test, and launch operations related to hydrogen propellant. The objective has been to reduce the operational cost of assembling and maintaining hydrogen delivery systems. Specifically, efforts have been made to develop a hydrogen leak detection system for point-contact measurement. Under the auspices of Lewis Research Center, the Electronics Design Center at Case Western Reserve University, Cleveland, Ohio, has undertaken the development of a point-contact hydrogen gas sensor with potential applications to the hydrogen propellant industry. We envision a sensor array consisting of numbers of discrete hydrogen sensors that can be located in potential leak sites. Silicon-based microfabrication and micromachining techniques are used in the fabrication of these sensor prototypes. Evaluations of the sensor are carried out in-house at Case Western Reserve University as well as at Lewis Research Center and GenCorp Aerojet, Sacramento, California. The hydrogen gas sensor is not only applicable in a hydrogen propulsion system, but also usable in many other civilian and industrial settings. This includes vehicles or facility use, or in the production of hydrogen gas. Dual space and commercial uses of these point-contacted hydrogen sensors are feasible and will directly meet the needs and objectives of NASA as well as various industrial segments.

Liu, Chung-Chiun↗