Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Performance of a Discrete Wavelet Transform for Compressing Plasma Count Data and its Application to the Fast Plasma Investigation on NASA's Magnetospheric Multiscale Mission

Plasma measurements in space are becoming increasingly faster, higher resolution, and distributed over multiple instruments. As raw data generation rates can exceed available data transfer bandwidth, data compression is becoming a critical design component. Data compression has been a staple of imaging instruments for years, but only recently have plasma measurement designers become interested in high performance data compression. Missions will often use a simple lossless compression technique yielding compression ratios of approximately 2:1, however future missions may require compression ratios upwards of 10:1. This study aims to explore how a Discrete Wavelet Transform combined with a Bit Plane Encoder (DWT/BPE), implemented via a CCSDS standard, can be used effectively to compress count information common to plasma measurements to high compression ratios while maintaining little or no compression error. The compression ASIC used for the Fast Plasma Investigation (FPI) on board the Magnetospheric Multiscale mission (MMS) is used for this study. Plasma count data from multiple sources is examined: resampled data from previous missions, randomly generated data from distribution functions, and simulations of expected regimes. These are run through the compression routines with various parameters to yield the greatest possible compression ratio while maintaining little or no error, the latter indicates that fully lossless compression is obtained. Finally, recommendations are made for future missions as to what can be achieved when compressing plasma count data and how best to do so.

Particles↗
NASA NTRS2015

Discrete Fourier Transform in a Complex Vector Space

An image-based phase retrieval technique has been developed that can be used on board a space based iterative transformation system. Image-based wavefront sensing is computationally demanding due to the floating-point nature of the process. The discrete Fourier transform (DFT) calculation is presented in "diagonal" form. By diagonal we mean that a transformation of basis is introduced by an application of the similarity transform of linear algebra. The current method exploits the diagonal structure of the DFT in a special way, particularly when parts of the calculation do not have to be repeated at each iteration to converge to an acceptable solution in order to focus an image.

Dean, Bruce H.↗
NASA NTRS2015

Calculated Low-Speed Steady and Time-Dependent Aerodynamic Derivatives for Some Airfoils Using a Discrete Vortex Method

This paper contains a collection of some results of four individual studies presenting calculated numerical values for airfoil aerodynamic stability derivatives in unseparated inviscid incompressible flow due separately to angle-of-attack, pitch rate, flap deflection, and airfoil camber using a discrete vortex method. Both steady conditions and oscillatory motion were considered. Variables include the number of vortices representing the airfoil, the pitch axis / moment center chordwise location, flap chord to airfoil chord ratio, and circular or parabolic arc camber. Comparisons with some experimental and other theoretical information are included. The calculated aerodynamic numerical results obtained using a limited number of vortices provided in each study compared favorably with thin airfoil theory predictions. Of particular interest are those aerodynamic results calculated herein (such as induced drag) that are not readily available elsewhere.

Riley, Donald R.↗
NASA NTRS2016

Calculated Low-Speed Steady and Time-Dependent Aerodynamic Derivatives for Several Different Wings Using a Discrete Vortex Method

Calculated numerical values for some aerodynamic terms and stability Derivatives for several different wings in unseparated inviscid incompressible flow were made using a discrete vortex method involving a limited number of horseshoe vortices. Both longitudinal and lateral-directional derivatives were calculated for steady conditions as well as for sinusoidal oscillatory motions. Variables included the number of vortices used and the rotation axis/moment center chordwise location. Frequencies considered were limited to the range of interest to vehicle dynamic stability (kb <.24 ). Comparisons of some calculated numerical results with experimental wind-tunnel measurements were in reasonable agreement in the low angle-of-attack range considering the differences existing between the mathematical representation and experimental wind-tunnel models tested. Of particular interest was the presence of induced drag for the oscillatory condition.

Riley, Donald R.↗
NASA NTRS2018

Discrete Assemblers Utilizing Conventional Motion Systems

An alternative to additive manufacturing is disclosed, introducing an end-to-end workflow in which discrete building blocks are reversibly joined to produce assemblies called digital materials. Described is the design of the bulk-material building blocks and the devices that are assembled from them. Detailed is the design and implementation of an automated assembler, which takes advantage of the digital material structure to avoid positioning errors within a large tolerance. To generate assembly sequences, a novel CAD/CAM workflow is described for designing, simulating, and assembling digital materials. The structures assembled using this process have been evaluated, showing that the joints perform well under varying conditions and that the assembled structures are functionally precise.

Langford, William Kai↗
NASA NTRS2017

Iterative Discrete Ordinates Solution of the Equation for the Surface-Reflected Radiance

This paper presents a new method of numerical solution of the integral equation for the radiance reflected from an anisotropic surface. The equation relates the radiance at the surface level with BRDF and solutions of the standard radiative transfer problems for a slab with no reflection on its surfaces. It is also shown that the kernel of the equation satisfies the condition of the existence of a unique solution and the convergence of the successive approximations to that solution. The developed method features two basic steps: discretization on a 2D quadrature, and solving the resulting system of algebraic equations with successive over-relaxation method based on the Gauss-Seidel iterative process. Presented numerical examples show good coincidence between the surface-reflected radiance obtained with DISORT and the proposed method. Analysis of contributions of the direct and diffuse (but not yet reflected) parts of the downward radiance to the total solution is performed. Together, they represent a very good initial guess for the iterative process. This fact ensures fast convergence. The numerical evidence is given that the fastest convergence occurs with the relaxation parameter of 1 (no relaxation). An integral equation for BRDF is derived as inversion of the original equation. The potential of this new equation for BRDF retrievals is analyzed. The approach is found not viable as the BRDF equation appears to be an ill-posed problem, and it requires knowledge the surface-reflected radiance on the entire domain of both Sun and viewing zenith angles.

Alexander Radkevich↗
NASA NTRS2020

CoCoSim, a Code Generation Framework for Control/command Applications: An Overview of CoCoSim for Multi-Periodic Discrete Simulink Models

We present CoCoSim, a framework to support the design, code generation and analysis of discrete dataflow model expressed in Simulink. In this work, we specifically focus on the analysis and code generation of multi-periodic systems. For that CoCoSim provides two complementary approaches: the first amounts to encode the multiperiodic semantics in a pure-synchronous one – à la Lustre–, enabling the use of model checker for verifying properties. The second provides a faithful code generation into multiple communicating (mono)synchronous components – à la Prelude– that can be then simulated or embedded in the final platform with any real-time scheduler. These approaches have been experimented in various settings.

Bourbouh, Hamza↗
NASA NTRS2020

Analysis of Sample Acquisition Dynamics Using Discrete Element Method

The analysis presented in this paper is conducted in the framework of the Ocean Worlds Autonomy Testbed for Exploration Research and Simulation (OceanWATERS) project, currently under development at NASA Ames Research Center. OceanWATERS aims at designing a simulation environment which allows for testing autonomy of scientific lander missions to the icy moons of our solar system. Mainly focused on reproducing the end effector interaction with the inherent terrain, this paper introduces a novel discrete element method (DEM)-based approach to determine forces and torques acting on the lander’s scoop during the sample acquisition process. An accurate force feedback from the terrain on the scoop is required by fault-detection and autonomous decision-making algorithms to identify when the requested torque on the robotic arm’s joints exceeds the maximum available torque. Knowledge of the terrain force feedback significantly helps evaluating the arm’s links structural properties and properly selecting actuators for the joints. Models available in literature constitute a partial representation of the dynamics of the interaction. As an example, Balovnev derived an analytical expression of the vertical and horizontal force acting on a bucket while collecting a sample as a function of its geometry and velocity, soil parameters and reached depth. Although the model represents an adequate approximation of the two force components, it ignores the direction orthogonal to the scoop motion and neglects the torque. This work relies on DEM analysis to compensate for analytical models’ deficiencies and inaccuracies, i. e. provide force and torque 3D vectors, defined in the moving reference (body) frame attached to the scoop, at each instant of the sample collection process. Results from the first presented analysis relate to the specific OceanWATERS sampling strategy, which consists of collecting the sample through five consecutive passes with increasing depth, each pass following the same circularlinear- circular trajectory. Data is collected given a specific scoop design interacting with two types of bulk materials, which may characterize the surface of icy planetary bodies: snow and ice. Although specifically concerned with the OceanWATERS design, this first analysis provides the expected force trends for similar sampling strategies and allows to deduce phenomenological information about the general scooping process. In order to further instruct the community on the use of DEM tools as a solution to the sampling collection problem, two more analyses have been carried out, mainly focused on reducing the DEM computation time, which increases with a decrease in particle size. After running a set of identical simulations, where the only changing parameter is the size of the spherical particle, it is observed that the resulting force trajectories, starting from a given particle size, converge to the true trend. It is deducible that a further decrease in size yields negligible improvements in the accuracy, while it sensibly increases computation time. A final analysis aims at discussing limitations of approximating bulk material particles having a complex shape, e. g. ice fragments, with spheres, by comparing force trends resulting in the two cases for the same simulation scenario.

Catanoso, Damiana↗
NASA NTRS2019

Real-time Exponential Curve Fits Using Discrete Calculus

This paper presents an improved solution for curve fitting data to an exponential equation (Y = AeBt + C). This improvement is in four areas ? speed, stability, determinant processing time, and the removal of limits. The solution presented in this paper avoids iterative techniques and their stability errors by using three mathematical ideas ? discrete calculus, a special relationship (between exponential curves and the Mean Value Theorem for Derivatives), and a simple linear curve fit algorithm. This method can also be applied to fitting data to the general power law equation Y = AxB + C and the general geometric growth equation Y = AkBt + C.

Rowe, Geoffrey K.↗
NASA NTRS2019

Stabilization of a Swept-Wing Boundary Layer by Discrete Roughness Elements at High Reynolds Numbers

Direct numerical simulations (DNS) are performed to study potential stabilizing ef- fect of spanwise periodic discrete roughness elements (DREs) on cross ow instabilities in a spatially developing three-dimensional boundary layer over an in nite-swept natural- laminar- ow wing at a freestream Mach number of 0:75 and a chord Reynolds number of approximately 25 million. In the DNS, both the spanwise periodic DREs and distributed roughness in the leading-edge region are implemented to simulate a typical experimen- tal scenario in which multiple steady cross ow modes including the most unstable mode (i.e., the \target" mode) emerge because of the presence of naturally distributed surface roughness in the leading edge region and spanwise periodic control cylinders of subcritical wavelength are used to force small-wavelength disturbances (i.e., the control mode) for damping the target mode. The DNS results show that the e ectiveness of DRE control is sensitive to roughness diameter, height, and chordwise placement. For the DRE parame- ters considered in this study, the stabilizing e ect on the target mode is small within the computational domain that ended at about 35% of the chord.

Nicholson, Gary↗
NASA NTRSDate not supplied

Recent Advances in Discrete Crack Modeling Applied to Laminated Composites with Emphasis on: Floating Node Method, VCCT and Cohesive Zone Modeling

The present talk will provide an overview of the work performed during the Advanced Composites Project (ACP) on the development, and verification and validation of the Floating Node Method (FNM) as well as the Virtual Crack Closure Technique (VCCT) and cohesive zone modeling (CZM). The FNM is a finite element-based technique to represent crack networks. The complex nature of matrix crack-delamination interactions observed in unidirectional (UD) tape laminates suggests that a methodology such as the FNM method may be required to be able to accurately simulate damage progression in these laminates. Simulating crack onset and growth within the context of the FNM relies on techniques such as VCCT and CZM. The talk is organized to provide, via select examples, an overview of the breadth of the Verification & Validation (V&V) exercises performed during the ACP, and how these challenged the state-of-the art and guided further developments in discrete crack modeling, while helping to establish confidence in the progress made and map the challenges ahead. The performance of the VCCT and CZM individually, and in combination with the FNM can be assessed through verification exercises. These exercises typically consist of a comparison of simulation results to known numerical or analytical solutions. Verification is key to identify implementation issues and limitations that, otherwise, may remain undetected and cloud any subsequent validation efforts. Indeed, a subset of these numerical exercises led to further developments of the VCCT and the FNM method as will be illustrated. Before embarking on the subsequent validation of the framework, it is critical to have adequate characterization data. However, the testing campaign conducted revealed material responses that challenged the state-of-the-art and required further developments. The developments in CZM technology associated with the modeling of the responses of hybrid interfaces (fabric/UD) will be given as an example. Finally, the talk will conclude with a summary of the validation exercises performed under quasi-static and fatigue loadings, highlighting some of the key achievements, outstanding challenges and lessons learned.

finite elements↗
NASA NTRSDate not supplied

Implementation of Charged Particle Behavior in Discrete Element Method (DEM) Simulations

To understand the behavior of charged lunar regolith when perturbed by lunar landers, it is important to couple the grain dynamics with mechanical and electrical particle interactions. To accomplish this, improvements have been made to a discrete element method (DEM) software package to include both short- and long-range interactions between spherical particles. Short-range interactions rely on contact between the particles, such as electrical conduction and triboelectric charge transfer driven by work functions. Long-range interactions act at a distance between every pairing of particles, such as electrostatic forces and gravitational forces. Results from simulations between a few particles are compared with theory to verify these added behaviors prior to scaling up to more complex scenarios. The radii, initial charges, electrical conductivities, work functions, and separation of the particles are varied and the resultant charges as well as the time required to reach the final state are determined.

Electrostatics↗
NASA NTRS1993

Discrete Recurrent Neural Networks for Grammatical Inference

Recurrent neural networks have recently been shown to have the ability to learn regular and context-free grammars from examples. We show that while conventional analog recurrent networks try to form clusters in activation space to represent discrete states of the grammars during learning, and can be successful in doing so, the clusters so formed tend to become unstable as longer and longer test input strings are presented to the network.

Neural↗
NASA NTRSDate not supplied

Implementation of Charged Particle Behavior in Discrete Element Method (DEM) Simulations

Lunar landers will agitate the surface of the Moon with an exhaust plume during descent which will, due to the particulate nature of the lunar regolith, loosen and eject grains from the surface. This ejection is not only coupled with the charged plume gas, but also results in significant particle-particle interactions. Settling of these grains after plume effects have subsided takes much longer than expected in a ballistic trajectory. The prevailing hypothesis attributes this behavior to the accumulated charge on the particles. We are thus developing a discrete element method (DEM) approach to explore these charged particle interactions on the lunar surface. The Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) Improved for General Granular and Granular Heat Transfer Simulations (LIGGGHTS) software package provides a DEM modeling framework for granular interactions. It includes many complexities such as non-spherical particle shapes, cohesion and frictional forces, and heat transfer, but has no provision for inter-particle electrostatic forces and charge transfer that are important to examine in the lunar environment. In this work, a standard Coulomb potential and a Yukawa potential are integrated into the LIGGGHTS framework to provide a basis for particle-particle electrostatic interactions, as well as a gravitational potential to enable inter-grain gravitational interactions. A preliminary approach to charge transfer between grains incorporating properties such as work function and electrical conductivity to the library of available material characteristics will be presented. Several scenarios have been simulated that include charged particle interactions within a diffuse granular gas, settling of charged grains into a regolith bed, sliding of granular material along an incline, and vibration of settled grains to produce a behavior similar to fluidization. There are numerous challenges to incorporate realistic interactions between complex lunar particles. Currently, grains are modeled to behave as if the entirety of the charge acts at the center of mass, such as conductors with spherical symmetry and insulators with homogeneously distributed charge. We are developing improvements that will include the use of non-spherical particle geometries, as well as reasonable approximations of insulating/dielectric materials that have non-uniform charge distributions. The cases simulated thus far will be examined in a relevant environment within a vacuum chamber to validate the simulations. These simulations will be bounded by experiments utilizing high-speed camera observations of the motion for validation. The grains in the experiment will exchange charge during their motion and this can be quantified by collection within a charge measurement device such as a Faraday cup. Such a device may be modeled within the software by defining an integration region and computing the contained charge as a function of simulation time, allowing for side-by-side comparison of simulated and measured bulk charging results. Any differences will be reconciled by updating the mathematical mechanisms described within the simulation suite. Successfully combining results from experiments within a relevant environment into the LIGGGHTS framework will improve modeling of the charged grain dynamics experienced on the Moon to provide insights into dust behavior for future lunar exploration missions.

Electrostatics↗
NASA NTRS2021

Prompt Response of the Dayside Magnetosphere to Discrete Structures within the Sheath Region of a Coronal Mass Ejection

A sequence of discrete solar wind structures within the sheath region of an interplanetary coronal mass ejection (ICME) on Nov 06, 2015 caused a series of compressions and releases of the dayside magnetosphere. Each compression resulted in a brief adiabatic enhancement of ions (electrons) driving bursts of electromagnetic ion cyclotron (whistler mode chorus) wave growth across the dayside magnetosphere. Fine-structured rising tones were observed in the EMIC wave bursts, resulting in nonlinear scattering of relativistic electrons in the outer radiation belt. Multipoint observations allow us to study the spatial structure and evolution of these sheath structures as they propagate Earthward from L1 as well as the spatio-temporal characteristics of the magnetospheric response. This event highlights the importance of fine-scale solar wind structure, in particular within complex sheath regions, in driving dayside phenomena within the inner magnetosphere.

L. W. Blum↗
NASA NTRSDate not supplied

Investigation of a Smooth Local Correlation-based Transition Model in a Discrete-Adjoint Aerodynamic Shape Optimization Algorithm

A smooth local correlation-based transition model is fully coupled to a RANS-based Newton-Krylov flow solver and discrete-adjoint gradient-based optimization algorithm. The free-transition optimization framework is evaluated using lift-constrained drag minimizations of airfoils at design conditions ranging from light to single-aisle aircraft and an infinite swept wing at design conditions representative of a transonic strut-braced wing aircraft. The impact of the streamwise grid resolution on the ability of the optimization algorithm to delay boundary-layer transition is investigated, with the results demonstrating that streamwise grid resolution requirements increase as the transition length decreases with increasing Reynolds number. The optimization problem at the light aircraft design conditions is demonstrated to be multi-modal, with the optimization algorithm producing two distinct designs: one with a thin, reflexed trailing edge and steep pressure recovery regions, the other with increased aft loading, with the latter design outperforming the former. A drag minimization of an airfoil at transonic design conditions demonstrates that the optimization algorithm successfully trades a decrease in viscous drag by delaying boundary-layer transition with an increase in wave drag, while the drag minimization of an infinite swept wing demonstrates the capability of the optimizational gorithm to delay both Tollmien-Schlichting and stationary crossflow instabilities.

AATT↗
NASA NTRS2021

Discrete anisotropic radiative transfer modelling of solar-induced chlorophyll fluorescence: Structural impacts in geometrically explicit vegetation canopies

Solar-induced fluorescence (SIF) is a subtle but informative optical signal of vegetation photosynthesis. Remotely sensed SIF integrates environmental, physiological and structural changes that alter photosynthesis at leaf, plant and canopy scales. Radiative transfer models are ideally suited to investigate the complex sources of variability in the SIF signal to guide the interpretation of SIF retrievals from airborne and space-borne platforms. Here, we coupled the Fluspect-Cx model of leaf optical properties and chlorophyll-a fluorescence with the Discrete Anisotropic Radiative Transfer (DART) model to upscale SIF from individual leaves to three-dimensional (3D) structurally explicit canopies. For one-dimensional homogeneous (turbid-like) canopies, DART-SIF was nearly identical to SIF simulated in two existing models, SCOPE and mSCOPE (RMSE <0.221 W./sq.m.μm.sr). DART simulations in geometrically explicit 3D canopies offered four important insights regarding the influence of vegetation structure on the multi-angular top-of-canopy SIF signal. First, changes in the 3D canopy architecture of maize crops, represented by leaf density (leaf area index), and plant clumping (canopy closure) had a larger impact on SIF than the modelled photosynthetic efficiency distinction between sun-adapted and shade-adapted foliage. Second, clumping of leaves at the crop and stand levels was identified as one of the key driving factors of multi-angular anisotropy of red and far-red SIF (686 and 740 nm) for both maize and eucalyptus canopies. Third, non-photosynthetic woody material had a significant impact on top-of-canopy SIF in modelled 3D forest stands. Wood shadowing decreased the photosynthetically active radiation absorbed by green leaves, and consequently the SIF emissions, by 10% in sparse and 17% in dense eucalyptus stands. The wood obstruction (blocking) effect, quantified as a relative difference of SIF escape probabilities from canopies with and without wood in the nadir viewing direction, decreased far-red SIF by 4–6% but it had a smaller and sometimes positive influence (by less than 2%) on red SIF. Fourth, DART 3D radiative budget profiles revealed that the majority of the SIF signal from a dense eucalyptus stand originated from the top 25% of the simulated canopy. Interestingly, the introduction of bark-covered woody elements did not alter the simulated balance and omnidirectional escape factor of red SIF in this upper canopy part but did raise significantly both of them in case of far-red SIF. These results demonstrate the importance of 3D radiative transfer and radiative budget simulations for investigating SIF interactions in structurally complex plant canopies and for a better understanding of spatiotemporal and multi-angular remote sensing SIF observations.

Zbyněk Malenovský↗
NASA NTRS2020

Entropy Stable Nonconforming Discretizations with the Summation-By-Parts Property for Curvilinear Coordinates

The entropy conservative/stable algorithm of Friedrichet al.(2018) for hyperbolic conservation laws on nonconforming p-refined/coarsened Cartesian grids, is extended to curvilinear grids for the compressible Euler equations. The primary focus is on constructing appropriate coupling procedures across the curvilinear nonconforming interfaces. A simple and flexible approach is proposed that uses interpolation operators from one element to the other. On the element faces,the analytic metrics are used to construct coupling terms, while metric terms in the volume are approximated to satisfy a discretization of the geometric conservation laws. The resulting scheme is entropy conservative/stable, elementwise conservative, and freestream preserving. The accuracy and stability properties of the resulting numerical algorithm are shown to be comparable to those ofthe original conforming scheme (∼p+ 1 convergence) in the context of the isentropic Euler vortex and the inviscid Taylor–Green vortex problems on manufactured high order grids.

David Del Rey Fernandez↗