Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial 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 469 records · Page 26

Pulse Response of the Short-Wave Infrared Detection System of the Ocean Color Instrument for the NASA PACE Mission

The Ocean Color Instrument (OCI) on NASA’s Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission is a hyperspectral Earth imager with a spatial resolution of 1 km x 1 km and a spectral resolution of 5 nm in 2.5 nm steps over 342-887 nm. In addition, OCI provides 7 discrete bands in the 940-2260 nm Short-Wave InfraRed (SWIR) range. The front-end optical imager is a rotating mirror-based system that images the ground scene onto a slit with an instantaneous field of view of 16 km x 1 km. For the SWIR bands, the slit-image is re-imaged onto a 16x1 micro-lens array that effectively acts as the focal plane since each lens element is fiber coupled to wavelength filtered InGaAs and HgCdTe Photo Diodes (PDs). The pulse response of the detection system is critical to OCI SWIR performance. We find that PDs introduce an inherent slow tail in the pulse response due to slow diffusion moving carriers in their n and p regions. We show that this introduces response errors ranging from 1 down to 0.01 % for up to tens of science pixels after the pulse depending on the PD design and materials. It is shown that the response is distinctly different for the InGaAs and HgCdTe PDs. We explain how the front-end design can further increase this error. Finally, we detail the cause of the slow pulse response tail, how to model it, its impact on OCI performance and how it is characterized and corrected to meet OCI requirements.

ocean color↗

Pulse Response of the Short-Wave Infrared Detection System of the Ocean Color Instrument for the NASA Pace Mission

The Ocean Color Instrument (OCI) on NASA’s Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission is a hyperspectral Earth imager with a spatial resolution of 1 km x 1 km and a spectral resolution of 5 nm in 2.5 nm steps over 342-887 nm. In addition, OCI provides 7 discrete bands in the 940-2260 nm Short-Wave InfraRed (SWIR) range. The front-end optical imager is a rotating mirror-based system that images the ground scene onto a slit with an instantaneous field of view of 16 km x 1 km. For the SWIR bands, the slit-image is re-imaged onto a 16x1 micro-lens array that effectively acts as the focal plane since each lens element is fiber coupled to wavelength filtered InGaAs and HgCdTe Photo Diodes (PDs). The pulse response of the detection system is critical to OCI SWIR performance. We find that PDs introduce an inherent slow tail in the pulse response due to slow diffusion moving carriers in their n and p regions. We show that this introduces response errors ranging from 1 down to 0.01 % for up to tens of science pixels after the pulse depending on the PD design and materials. It is shown that the response is distinctly different for the InGaAs and HgCdTe PDs. We explain how the front-end design can further increase this error. Finally, we detail the cause of the slow pulse response tail, how to model it, its impact on OCI performance and how it is characterized and corrected to meet OCI requirements.

ocean color↗

Perseverance’s Scanning Habitable Environments with Raman and Luminescence for Organics and Chemicals (SHERLOC) Investigation

The Scanning Habitable Environments with Raman and Luminescence for Organics and Chemicals (SHERLOC) is a robotic arm-mounted instrument on NASA’s Perseverance rover. SHERLOC has two primary boresights. The Spectroscopy boresight generates spatially resolved chemical maps using fluorescence and Raman spectroscopy coupled to microscopic images (10.1 μm/pixel). The second boresight is a Wide Angle Topographic Sensor for Operations and eNgineering (WATSON); a copy of the Mars Science Laboratory (MSL) Mars Hand Lens Imager (MAHLI) that obtains color images from microscopic scales (∼13 μm/pixel) to infinity. SHERLOC Spectroscopy focuses a 40 μs pulsed deep UV neon-copper laser (248.6 nm), to a ∼100 μm spot on a target at a working distance of ∼48 mm. Fluorescence emissions from organics, and Raman scattered photons from organics and minerals, are spectrally resolved with a single diffractive grating spectrograph with a spectral range of 250 to ∼370 nm. Because the fluorescence and Raman regions are naturally separated with deep UV excitation (<250 nm), the Raman region ∼ 800 – 4000 cm−1 (250 to 273 nm) and the fluorescence region (274 to ∼370 nm) are acquired simultaneously without time gating or additional mechanisms. SHERLOC science begins by using an Autofocus Context Imager (ACI) to obtain target focus and acquire 10.1 μm/pixel greyscale images. Chemical maps of organic and mineral signatures are acquired by the orchestration of an internal scanning mirror that moves the focused laser spot across discrete points on the target surface where spectra are captured on the spectrometer detector. ACI images and chemical maps (< 100 μm/mapping pixel) will enable the first Mars in situ view of the spatial distribution and interaction between organics, minerals, and chemicals important to the assessment of potential biogenicity (containing CHNOPS). Single robotic arm placement chemical maps can cover areas up to 7x7 mm in area and, with the < 10 min acquisition time per map, larger mosaics are possible with arm movements. This microscopic view of the organic geochemistry of a target at the Perseverance field site, when combined with the other instruments, such as Mastcam-Z, PIXL, and SuperCam, will enable unprecedented analysis of geological materials for both scientific research and determination of which samples to collect and cache for Mars sample return.

Rohit Bhartia↗

Investigation of Improved Wall Heat Flux Behavior Observed When Computing Hypersonic Flow Using a Face-averaged Nodal-gradient Approach On Tetrahedral Grids

We investigate the origin of a significant reduction in computed numerical “noise” in the spatial variation of computed wall heat flux observed when using a node-centered weighted least squares(WLSQ) gradient approach when solving the Reynolds-averaged Navier-Stokes (RANS) equations for hypersonic boundary layer flows on highly stretched tetrahedral grids. This investigation is conducted using a 2nd-order, cell-centered, finite-volume, discretization which we employ to contrast results obtained using two WLSQ gradient approaches: a conventional 3-D, cell-centered, node-neighbor, WLSQ gradient (NN-CCG) approach and a 3-D, face-averaged, node-centered gradient (F-ANG)approach. The F-ANG approach, is shown to produce significantly less “noise” in the spatial variation of wall heat flux, than the NN-CCG approach and we demonstrate the cause of this “noise” via a series of numerical experiments where we compute hypersonic laminar and turbulent boundary layer flows.We further investigate the level of grid resolution required for a tetrahedral grid to approach the fidelity of a grid resolved hexahedral grid solution obtained using an equivalent structured grid method.

Jeffery A. White↗

X-ray image of the Cygnus Loop

An X-ray picture of the Cygnus Loop, recorded with an imaging X-ray telescope, is presented. The results are highly suggestive of a limb-brightened shell of hot gas which results from the expansion of a blast wave into the interstellar medium. Spatial structure is clearly evident on scale sizes down to about 1/4 deg, a likely indication of the inhomogeneities in the interstellar medium. No evidence is found for any discrete X-ray-emitting remnant of the original supernova explosion.

Rappaport, S.↗

Space-time modeling using environmental constraints in a mobile robot system

Grid-based models of a robot's local environment have been used by many researchers building mobile robot control systems. The attraction of grid-based models is their clear parallel between the internal model and the external world. However, the discrete nature of such representations does not match well with the continuous nature of actions and usually serves to limit the abilities of the robot. This work describes a spatial modeling system that extracts information from a grid-based representation to form a symbolic representation of the robot's local environment. The approach makes a separation between the representation provided by the sensing system and the representation used by the action system. Separation allows asynchronous operation between sensing and action in a mobile robot, as well as the generation of a more continuous representation upon which to base actions.

Slack, Marc G.↗

Spatial simulation of boundary layer instability - Effects of surface roughness

The effects of an isolated, two-dimensional roughness element on the spatial development of instability waves in boundary layers are investigated by numerically integrating the two-dimensional, time-dependent, incompressible Navier-Stokes equations, using a finite difference/Chebyshev discretization. It is shown that (high) inviscid frequencies have higher growth rates than Tollmien-Schlichting frequencies, indicating that disturbances growing in the separation zone are controlled by the inviscid instability of the shear layer at the edge of the separation zone.

Danabasoglu, G.↗

Analysis and Design of Variable Stiffness Composite Cylinders

An investigation of the possible performance improvements of thin circular cylindrical shells through the use of the variable stiffness concept is presented. The variable stiffness concept implies that the stiffness parameters change spatially throughout the structure. This situation is achieved mainly through the use of curvilinear fibers within a fiber-reinforced composite laminate, though the possibility of thickness variations and discrete stiffening elements is also allowed. These three mechanisms are incorporated into the constitutive laws for thin shells through the use of Classical Lamination Theory. The existence of stiffness variation within the structure warrants a formulation of the static equilibrium equations from the most basic principles. The governing equations include sufficient detail to correctly model several types of nonlinearity, including the formation of a nonlinear shell boundary layer as well as the Brazier effect due to nonlinear bending of long cylinders. Stress analysis and initial buckling estimates are formulated for a general variable stiffness cylinder. Results and comparisons for several simplifications of these highly complex governing equations are presented so that the ensuing numerical solutions are considered reliable and efficient enough for in-depth optimization studies. Four distinct cases of loading and stiffness variation are chosen to investigate possible areas of improvement that the variable stiffness concept may offer over traditional constant stiffness and/or stiffened structures. The initial investigation deals with the simplest solution for cylindrical shells in which all quantities are constant around the circumference of the cylinder. This axisymmetric case includes a stiffness variation exclusively in the axial direction, and the only pertinent loading scenarios include constant loads of axial compression, pressure, and torsion. The results for these cases indicate that little improvement over traditional laminates exists through the use of curvilinear fibers, mainly due to the presence of a weak link area within the stiffness variation that limits the ultimate load that the structure can withstand. Rigorous optimization studies reveal that even though slight increases in the critical loads can be produced for designs with an arbitrary variation of the fiber orientation angle, the improvements are not significant when compared to traditional design techniques that utilize ring stiffeners and frames. The second problem that is studied involves arbitrary loading of a cylinder with a stiffness variation that changes only in the circumferential direction. The end effects of the cylinder are ignored, so that the problem takes the form of an analysis of a cross-section for a short cylinder segment. Various load cases including axial compression, pressure, torsion, bending, and transverse shear forces are investigated. It is found that the most significant improvements in load-carrying capability exist for cases which involve loads that also vary around the circumference of the shell, namely bending and shear forces. The stiffness variation of the optimal designs contribute to the increased performance in two ways: lowering the stresses in the critical areas through redistribution of the stresses; and providing a relatively stiff region that alters the buckling behavior of the structure. These results lead to an in-depth optimization study involving weight optimization of a fuselage structure subjected to typical design constraints. Comparisons of the curvilinear fiber format to traditional stiffened structures constructed of isotropic and composite materials are included. It is found that standard variable stiffness designs are quite comparable in terms of weight and load-carrying capability yet offer the added advantage of tailorability of distinct regions of the structure that experience drastically different loading conditions. The last two problems presented in this work involve the nonlinear phenomenon of long tubes under bending. Though this scenario is not as applicable to fuselage structures as the previous problems, the mechanisms that produce the nonlinear effect are ideally suited to be controlled by the variable stiffness concept. This is due to the fact that the dominating influence for long cylinders under bending is the ovalization of the cross-section, which is governed mainly by the stiffness parameters of the cylindrical shell. Possible improvement of the critical buckling moments for these structures is investigated using either a circumferential or axial stiffness variation. For the circumferential case involving infinite length cylinders, it is found that slight improvements can be observed by designing structures that resist the cross-sectional deformation yet do not detract from the buckling resistance at the critical location. The results also indicate that buckling behavior is extremely dependent on cylinder length. This effect is most easily seen in the solution of finite length cylinders under bending that contain an axial stiffness variation. For these structures, the only mechanism that exhibits improved response are those that effectively shorten the length of the cylinder, thus reducing the cross-sectional deformation due to the forced restraint at the ends. It was found that the use of curvilinear fibers was not able to achieve this effect in sufficient degree to resist the deformation, but that ring stiffeners produced the desired response admirably. Thus, it is shown that the variable stiffness concept is most effective at improving the bending response of long cylinders through the use of a circumferential stiffness variation.

Tatting, Brian F.↗

Variational Asymptotic Homogenization of Finitely Deformed Viscoelastic-Viscoplastic Composites

The objective of this paper is to develop a constitutive model for finitely deformed viscoelastic-viscoplastic materials and a micromechanics approach to homogenizing composites consisting of such materials. The development of the constitutive model involves establishing a thermodynamic framework based on finite strain theory, developing a viscoelasticity and a viscoplasticity model based on the thermodynamic framework, developing a radial return algorithm based on a classic framework, and deriving a closed-from incremental constitutive relation in the spatial configuration. The development of the micromechanics approach involves pulling-back the above constitutive relation to the material configuration, formulating a variational statement with the resulting constitutive relation, discretizing variational statement in a finite-dimensional space, and solving the discretized variational statement using an Euler–Newton predictor–corrector method. The constitutive model is calibrated via monotonic uniaxial compression tests on a polymer, and the calibrated model is validated by comparing its predictions with the cyclic test data. It is shown capable of characterizing viscoelasticity, viscoplasticity, and complex loading paths. The micromechanics approach’s capabilities are demonstrated through homogenizing a unidirectional fiber-reinforced composite, subjected to uniaxial, biaxial, and shear loading, at different strain rates. It is demonstrated to be capable of handling rate dependence and complex loading paths. The present framework can be further improved by implementing more sophisticated viscoelasticity and viscoplasticity models in future work.

Finite element analysis↗

Automation of measurement of spatial damping parameters

A traveling wave unit was built to process measurement data on spatial damping parameters and to simplify statistical treatment of the results. Mean values of the loss coefficient, lengths of flexual waves, and wave propagation velocities for a given set of discrete functions were determined. Results are given in a table.

Akolzin, A. N.↗

A space-time discretization procedure for wave propagation problems

Higher order compact algorithms are developed for the numerical simulation of wave propagation by using the concept of a discrete dispersion relation. The dispersion relation is the imprint of any linear operator in space-time. The discrete dispersion relation is derived from the continuous dispersion relation by examining the process by which locally plane waves propagate through a chosen grid. The exponential structure of the discrete dispersion relation suggests an efficient splitting of convective and diffusive terms for dissipative waves. Fourth- and eighth-order convection schemes are examined that involve only three or five spatial grid points. These algorithms are subject to the same restrictions that govern the use of dispersion relations in the constructions of asymptotic expansions to nonlinear evolution equations. A new eighth-order scheme is developed that is exact for Courant numbers of 1, 2, 3, and 4. Examples are given of a pulse and step wave with a small amount of physical diffusion.

Davis, Sanford↗

Implementing system simulation of C3 systems using autonomous objects

The basis of all conflict recognition in simulation is a common frame of reference. Synchronous discrete-event simulation relies on the fixed points in time as the basic frame of reference. Asynchronous discrete-event simulation relies on fixed-points in the model space as the basic frame of reference. Neither approach provides sufficient support for autonomous objects. The use of a spatial template as a frame of reference is proposed to address these insufficiencies. The concept of a spatial template is defined and an implementation approach offered. Discussed are the uses of this approach to analyze the integration of sensor data associated with Command, Control, and Communication systems.

Rogers, Ralph V.↗

An optimized finite-difference scheme for wave propagation problems

Two fully-discrete finite-difference schemes for wave propagation problems are presented, a maximum-order scheme and an optimized (or spectral-like) scheme. Both combine a seven-point spatial operator and an explicit six-stage time-march method. The maximum-order operator is fifth-order in space and is sixth-order in time for a linear problem with periodic boundary conditions. The phase and amplitude errors of the schemes obtained using Fourier analysis are given and compared with a second-order and a fourth-order method. Numerical experiments are presented which demonstrate the usefulness of the schemes for a range of problems. For some problems, the optimized scheme leads to a reduction in global error compared to the maximum-order scheme with no additional computational expense.

Zingg, D. W.↗

Semi-discrete Galerkin solution of the compressible boundary-layer equations with viscous-inviscid interaction

A semi-discrete Galerkin (SDG) method is under development to model attached, turbulent, and compressible boundary layers for transonic airfoil analysis problems. For the boundary-layer formulation the method models the spatial variable normal to the surface with linear finite elements and the time-like variable with finite differences. A Dorodnitsyn transformed system of equations is used to bound the infinite spatial domain thereby providing high resolution near the wall and permitting the use of a uniform finite element grid which automatically follows boundary-layer growth. The second-order accurate Crank-Nicholson scheme is applied along with a linearization method to take advantage of the parabolic nature of the boundary-layer equations and generate a non-iterative marching routine. The SDG code can be applied to any smoothly-connected airfoil shape without modification and can be coupled to any inviscid flow solver. In this analysis, a direct viscous-inviscid interaction is accomplished between the Euler and boundary-layer codes through the application of a transpiration velocity boundary condition. Results are presented for compressible turbulent flow past RAE 2822 and NACA 0012 airfoils at various freestream Mach numbers, Reynolds numbers, and angles of attack.

Day, Brad A.↗

Aeroacoustic computation of a gust-plate interaction via MacCormack schemes

We apply second, fourth, and sixth order spatially accurate variations of the MacCormack scheme to calculate the noise radiated by the interaction of a flat plate with an oncoming gust. For the given gust wavelength and numerical discretization, the fourth and sixth order schemes each are effective in capturing the occurring acoustic waves. However, for the case of the sixth order scheme, the high order of extrapolation applied at the far field boundaries necessitates an extension of the computational boundaries from their prescribed location.

Martin, James E.↗

Multigrid Method for Modeling Multi-Dimensional Combustion with Detailed Chemistry

A highly accurate and efficient numerical method is developed for modeling 3-D reacting flows with detailed chemistry. A contravariant velocity-based governing system is developed for general curvilinear coordinates to maintain simplicity of the continuity equation and compactness of the discretization stencil. A fully-implicit backward Euler technique and a third-order monotone upwind-biased scheme on a staggered grid are used for the respective temporal and spatial terms. An efficient semi-coarsening multigrid method based on line-distributive relaxation is used as the flow solver. The species equations are solved in a fully coupled way and the chemical reaction source terms are treated implicitly. Example results are shown for a 3-D gas turbine combustor with strong swirling inflows.

Zheng, Xiaoqing↗

Clouds on Neptune: Motions, Evolution, and Structure

The aims of our original proposal were these: (1) improving measurements of Neptune's circulation, (2) understanding the spatial distribution of cloud features, (3) discovery of new cloud features and understanding their evolutionary process, (4) understanding the vertical structure of zonal cloud patterns, (5) defining the structure of discrete cloud features, and (6) defining the near IR albedo and light curve of Triton. Towards these aims we proposed analysis of existing 1996 groundbased NSFCAM/IRTF observations and nearly simultaneous WFPC2 observations from the Hubble Space Telescope. We also proposed to acquire new observations from both HST and the IRTF.

Sromovsky, Larry A.↗

Computational Methods for Jet Noise Simulation

The purpose of our project is to develop, analyze, and test novel numerical technologies central to the long term goal of direct simulations of subsonic jet noise. Our current focus is on two issues: accurate, near-field domain truncations and high-order, single-step discretizations of the governing equations. The Direct Numerical Simulation (DNS) of jet noise poses a number of extreme challenges to computational technique. In particular, the problem involves multiple temporal and spatial scales as well as flow instabilities and is posed on an unbounded spatial domain. Moreover, the basic phenomenon of interest, the radiation of acoustic waves to the far field, involves only a minuscule fraction of the total energy. The best current simulations of jet noise are at low Reynolds number. It is likely that an increase of one to two orders of magnitude will be necessary to reach a regime where the separation between the energy-containing and dissipation scales is sufficient to make the radiated noise essentially independent of the Reynolds number. Such an increase in resolution cannot be obtained in the near future solely through increases in computing power. Therefore, new numerical methodologies of maximal efficiency and accuracy are required.

Goodrich, John W.↗