Search NASA⌕ Search

SEARCH · Search NASA

Results for “discrete surface”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Mars 2020's First Sample: The Fractured Rough Rock Unit on the Floor of Jezero Crater

A central goal of the Mars 2020 mission is to select and cache samples for future return to Earth. The first samples targeted for collection are from the crater-retaining, Crater Floor Fractured Rough (CF-Fr) unit of Stack et al., 2020.The CF-Fr unit is a topographically low unit in the current Jezero setting, likely overlain by morphologically discrete units, including possible Jezero delta deposits and eolian features. CF-Fr is an aerially extensive unit with lobate margins. Two distinct morphologies are observed: a locally-exposed lower expression, with flat relatively horizontal light-toned surfaces and polygonal fracturing, and an upper expression consisting of up to ~5vertical meters of massive, sometimes boulder-producing, material. The locally-exposed lower expression appears to represent a local ground level in which the upper material has been removed. Near the Octavia E. Butler landing site, the lower morphology is exposed as polygonally-fractured, light-toned bedrock that appears to grade continuously into higher-standing massive outcrops, often with no clearly exposed contact. Further south, a darker upper expression of the CF-Fr unit is more distinct and the unit exhibits some horizontal layering. The Perseverance rover will initially sample the flat-lying expression of this unit. Hypotheses for the origin of CF-Fr include fluvial, aeolian, or lacustrine sediment likely derived from the Jezero watershed, and/or pyroclastic material resulting from regional volcanism. Geochronology of the returned sample could be used to help constrain the timing of geological events in Jezero. It may also help constrain stratigraphic relationships with crater-retaining units outside of Jezero within the Nili Planum, which could be used to calibrate the cratering chronology of Mars. Paleomagnetic analyses of an oriented sample could establish the history of the martian dynamo and whether it persisted into the Hesperian. Finally, if present, secondary phases within the primary deposit would help constrain diagenetic conditions and inform post-depositional aqueous and potentially habitable environmental conditions.

Justin I Simon↗
NASA NTRSDate not supplied

Autonomous In-space Construction, Maintenance, and Reconfiguration Using Programmable Meta-Material

NASA ARC's Coded Structures Laboratory (CSL) is developing autonomous construction, maintenance, and reconfiguration technologies to meet long-duration and deep space infrastructure needs, in accordance with long-term NASA goals of "in-space reliance" and "mass-less exploration." We seek to achieve these capabilities by utilizing a "programmable meta-material" approach that integrates emerging advances in materials (mechanical meta-materials), manufacturing (cooperative mobile robotics), and autonomy (multi-agent planning algorithms). Through the ARMADAS project, we have shown assembly of high-performance engineered cellular materials using multiple cooperating mobile robotic assemblers. In this paper, we describe how such a programmable meta-material architecture may shift the paradigm of how we design, build, manufacture, and operate future space infrastructure and assets. The core of a programmable meta-material architecture consists of 3 main technology sub-areas: the structure, the assembly agents, and the assembly algorithms. We co-design these systems to ensure an adaptable system that can create and reconfigure structures from a base set of building block components. From this core technology, we can branch out and expand the capability of the system through additional secondary component types and robotic agents to perform activities such as inspections, maintenance, repair, payload installation, or perform power and communications interconnect. As these technologies mature, future designers will be able to utilize the system to rapidly integrate and operate assets in space or on planetary surfaces from a set of well-tested part library, or create their own modules to integrate into the system. A core trait to the development of this system is the automation approach. Because of the modular and functional discrete (pixel-like) nature of the structural system, a diverse set of powerful algorithms for analysis, planning, and simulation can be adapted and leveraged to optimize construction, maintenance, and dynamic reorganization (as hardware with programmable form and function). With an ability to free the design space from launch vehicle constraints and fundamentally shift how a mission is designed and conducted, we discuss the influence of a programmable meta-material architecture on mission design, build, and operations. For the "design phase", we discuss project lifecycle effects, costs, time, and performance. For the "build phase", we discuss reusability, ISRU, manufacturing, material logistics, and scalability. And for "operations", we discuss autonomy, maintenance and upgrades, reliability, and reconfiguration. Autonomy and modularity are the primary enabling traits of this system. Engineering systems that utilize a modular and reconfiguration building block approach such as digital communication and computation systems, currently lead all other areas of technology in size and complexity scalability. NASA is extending the benefits and flexibility of digital systems to hardware systems, to optimize materials lifecycle management and expand our space exploration mission capabilities.

in space assembly↗
NASA NTRS2023

Numerical Investigation of C L,max Prediction on the NASA High-Lift Common Research Model

An assessment of a Hybrid RANS/LES (HRLES) approach for C L,max prediction is presented for the NASA High-Lift Common Research Model (CRM-HL). Both free air and wind tunnel configurations of the CRM-HL are investigated and the results are compared to the QinetiQ wind tunnel experiments and two other numerical approaches: Reynolds Averaged Navier-Stokes (RANS) and Wall-Modeled Large Eddy Simulations (WMLES). For the free-air configuration, HRLES is shown to address some of the known shortcomings with RANS and prevent inboard and outboard flow separation particularly in the region of C L,max and post-stall. To achieve these improvements over RANS, LES-appropriate grids and numerical discretizations are required. It was also found that when applying HRLES to a RANS best practice grid and numerics that the HRLES method significantly underperformed RANS. For the in-tunnel configuration, HRLES showed good agreement with the loads, surface pressure and oil-flow photographs obtained in the experiment. HRLES was able to improve upon the RANS simulations, which showed a sharp loss of lift at the two highest angles-of-attack due to large scale inboard and outboard separation on the wing, by correctly predicting the corner flow separation and showing remarkably close agreement in the flow topologies with the experiment.

TTT↗
NASA NTRS1982

Hybrid state vector methods for structural dynamic and aeroelastic boundary value problems

A computational technique is developed that is suitable for performing preliminary design aeroelastic and structural dynamic analyses of large aspect ratio lifting surfaces. The method proves to be quite general and can be adapted to solving various two point boundary value problems. The solution method, which is applicable to both fixed and rotating wing configurations, is based upon a formulation of the structural equilibrium equations in terms of a hybrid state vector containing generalized force and displacement variables. A mixed variational formulation is presented that conveniently yields a useful form for these state vector differential equations. Solutions to these equations are obtained by employing an integrating matrix method. The application of an integrating matrix provides a discretization of the differential equations that only requires solutions of standard linear matrix systems. It is demonstrated that matrix partitioning can be used to reduce the order of the required solutions. Results are presented for several example problems in structural dynamics and aeroelasticity to verify the technique and to demonstrate its use. These problems examine various types of loading and boundary conditions and include aeroelastic analyses of lifting surfaces constructed from anisotropic composite materials.

Lehman, L. L.↗
DOE OSTI2022

Configuration Space Integration for Adsorbate Partition Functions: The Effect of Anharmonicity on the Thermophysical Properties of CO–Pt(111) and CH 3 OH–Cu(111)

A method for computing anharmonic thermophysical properties for adsorbates on metal surfaces has been extended to include libration, or frustrated rotation. Classical phase space integration is used with Monte Carlo sampling of the configuration space to obtain the partition function of CO on Pt(111) and CH 3 OH on Cu(111). A minima-preserving neural network potential energy surrogate is used within the integration routines. Direct state counting using discrete variable representation is used to benchmark the results. We find that the phase space integration approach is in excellent agreement with the direct state counting results. Comparison with standard models such as the harmonic oscillator indicates that anharmonicity contributes significantly to the thermodynamic properties of CH 3 OH on Cu(111). We find that there is also a considerable difference between the harmonic oscillator and phase space integration for CO on Pt(111), although the discrepancy can largely be attributed to the presence of multiple binding sites within the unit cell. We demonstrate that a multisite harmonic oscillator model might be sufficient for CO-Pt(111). A more thorough description of the potential energy surface, which can be achieved with phase space integration, is necessary for weakly bound adsorbates such as CH 3 OH. In conclusion, the thermophysical properties were used to calculate free energies of adsorption on the respective metals, and subsequently the equilibrium constants and Langmuir isotherms in relevant temperature ranges. The results show that the choice of model to obtain partition functions greatly affects the resulting surface coverages in kinetic models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. 2. Quantifying the Uncertainty

Uncertainty quantification of groundwater (GW) aquifer parameters is critical for efficient management and sustainable extraction of GW resources. These uncertainties are introduced by the data, model, and prior information on the parameters. We develop a Bayesian inversion framework that uses Interferometric Synthetic Aperture Radar (InSAR) surface deformation data to infer the laterally heterogeneous permeability of a transient linear poroelastic model of a confined GW aquifer. The Bayesian solution of this inverse problem takes the form of a posterior probability density of the permeability. Exploring this posterior using classical Markov chain Monte Carlo (MCMC) methods is computationally prohibitive due to the large dimension of the discretized permeability field and the expense of solving the poroelastic forward problem. However, in many partial differential equation (PDE)-based Bayesian inversion problems, the data are only informative in a few directions in parameter space. For the poroelasticity problem, we prove this property theoretically for a one-dimensional problem and demonstrate it numerically for a three-dimensional aquifer model. Here we design a generalized preconditioned Crank-Nicolson (gpCN) MCMC method that exploits this intrinsic low dimensionality by using a low-rank-based Laplace approximation of the posterior as a proposal, which we build scalably. The feasibility of our approach is demonstrated through a real GW aquifer test in Nevada. The inherently two-dimensional nature of InSAR surface deformation data informs a sufficient number of modes of the permeability field to allow detection of major structures within the aquifer, significantly reducing the uncertainty in the pressure and the displacement quantities of interest.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Numerical algorithms for water waves with background flow over obstacles and topography

Abstract We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

Ambrose, David M.↗
DOE OSTI2023

Quantifying the Impact of Vertical Resolution on the Representation of Marine Boundary Layer Physics for Global-Scale Models

Abstract While GCM horizontal resolution has received the majority of scale improvements in recent years, ample evidence suggests that a model’s vertical resolution exerts a strong control on its ability to accurately simulate the physics of the marine boundary layer. Here we show that, regardless of parameter tuning, the ability of a single-column model (SCM) to simulate the subtropical marine boundary layer improves when its vertical resolution is improved. We introduce a novel objective tuning technique to optimize the parameters of an SCM against profiles of temperature and moisture and their turbulent fluxes, horizontal winds, cloud water, and rainwater from large-eddy simulations (LES). We use this method to identify optimal parameters for simulating marine stratocumulus and shallow cumulus. The novel tuning method utilizes an objective performance metric that accounts for the uncertainty in the LES output, including the covariability between model variables. Optimization is performed independently for different vertical grid spacings and value of time step, ranging from coarse scales often used in current global models (120 m, 180 s) to fine scales often used in parameterization development and large-eddy simulations (10 m, 15 s). Uncertainty-weighted disagreement between the SCM and LES decreases by a factor of ∼5 when vertical grid spacing is improved from 120 to 10 m, with time step reductions being of secondary importance. Model performance is shown to converge at a vertical grid spacing of 20 m, with further refinements to 10 m leading to little further improvement. Significance Statement In successive generations of computer models that simulate Earth’s atmosphere, improvements have been mainly accomplished by reducing the horizontal sizes of discretized grid boxes, while the vertical grid spacing has seen comparatively lesser refinements. Here we advocate for additional attention to be paid to the number of vertical layers in these models, especially in the model layers closest to Earth’s surface where climatologically important marine stratocumulus and shallow cumulus clouds reside. Our experiments show that the ability of a one-dimensional model to represent physical processes important to these clouds is strongly dependent on the model’s vertical grid spacing.

Meteorology & Atmospheric Sciences↗
NASA NTRS1973

Extravehicular mobility unit thermal simulator

The analytical methods, thermal model, and user's instructions for the SIM bay extravehicular mobility unit (EMU) routine are presented. This digital computer program was developed for detailed thermal performance predictions of the crewman performing a command module extravehicular activity during transearth coast. It accounts for conductive, convective, and radiative heat transfer as well as fluid flow and associated flow control components. The program is a derivative of the Apollo lunar surface EMU digital simulator. It has the operational flexibility to accept card or magnetic tape for both the input data and program logic. Output can be tabular and/or plotted and the mission simulation can be stopped and restarted at the discretion of the user. The program was developed for the NASA-JSC Univac 1108 computer system and several of the capabilities represent utilization of unique features of that system. Analytical methods used in the computer routine are based on finite difference approximations to differential heat and mass balance equations which account for temperature or time dependent thermo-physical properties.

Hixon, C. W.↗
NASA NTRS1983

Three-dimensional unsteady Euler equation solutions using flux vector splitting

A method for numerically solving the three dimensional unsteady Euler equations using flux vector splitting is developed. The equations are cast in curvilinear coordinates and a finite volume discretization is used. An explicit upwind second-order predictor-corrector scheme is used to solve the discretized equations. The scheme is stable for a CFL number of two and local time stepping is used to accelerate convergence for steady-state problems. Characteristic variable boundary conditions are developed and used in the far field and at surfaces. No additional dissipation terms are included in the scheme. Numerical results are compared with results from an existing three dimensional Euler code and experimental data.

Whitfield, D. L.↗
NASA NTRS1989

Precise solution for a finite set of spherical coefficients from equiangular gridded data

An important goal of geodesy is to determine the anomalous potential and its derivatives outside of the earth. Representing the surface anomalies by a series of spherical harmonics is useful since it is then possible to do a term by term solution of Laplace's equation and upward continuation. The problem of finding such a spherical harmonic series for anomaly values given on an equiangular surface grid is addressed. (This is a first step toward the more complicated problem of finding a function such that locally averaged values fit a grid of mean anomalies.) Three approaches to this fitting problem are discussed and compared: the discrete Fourier technique, the discrete integral technique, and a new approach. The peculiar nature of the equiangular grid, with its increasing density of (noisy) data toward the poles, causes each method to exhibit a different type of difficulty. The new method is shown to be practical as well as precise since the numerical conditioning problems which appear can be successfully handled by such well-known techniques as a (simple) Kalman filter.

Zucker, Paul A.↗
NASA NTRS1991

Heat transfer in the turbulent boundary layer with a step change in surface roughness

Measurements of Stanton numbers, velocity profiles, temperature profiles, and turbulence intensity profiles are reported for turbulent flat plate boundary layer flows with a step change in surface roughness. The first 0.9 m length of the test surface is roughened with 1.27 mm diameter hemispheres spaced 2 base diameters apart in a staggered array. The remaining 1.5 m length is smooth. The experiments show that the step change from a rough to a smooth surface has a dramatic effect on the convective heat transfer. In many cases, the Stanton number drops below the smooth-wall correlation immediately downstream of the change in roughness. The Stanton number measurements are compared with predictions using the discrete element method with excellent results.

Taylor, Robert P.↗
NASA NTRS1992

A discrete vortex model for predicting wing rock of slender wings

The fluid mechanism responsible for generating wing rock of slender sharp-edged delta wings was investigated using an unsteady discrete vortex model developed for that purpose, which is based on results of experimental investigations. Combined experimental and computational results indicate that wing rock is sustained by a lag in the position of the leading edge vortices normal to the surface. Results of comutations also indicate that certain complex aerodynamic problems may be governed primarily by unsteady inviscid phenomena.

Arena, Andrew S., Jr.↗
NASA NTRS2005

A Method for Aircraft Concept Selection Using Multicriteria Interactive Genetic Algorithms

The problem of aircraft concept selection has become increasingly difficult in recent years as a result of a change from performance as the primary evaluation criteria of aircraft concepts to the current situation in which environmental effects, economics, and aesthetics must also be evaluated and considered in the earliest stages of the decision-making process. This has prompted a shift from design using historical data regression techniques for metric prediction to the use of physics-based analysis tools that are capable of analyzing designs outside of the historical database. The use of optimization methods with these physics-based tools, however, has proven difficult because of the tendency of optimizers to exploit assumptions present in the models and drive the design towards a solution which, while promising to the computer, may be infeasible due to factors not considered by the computer codes. In addition to this difficulty, the number of discrete options available at this stage may be unmanageable due to the combinatorial nature of the concept selection problem, leading the analyst to arbitrarily choose a sub-optimum baseline vehicle. These concept decisions such as the type of control surface scheme to use, though extremely important, are frequently made without sufficient understanding of their impact on the important system metrics because of a lack of computational resources or analysis tools. This paper describes a hybrid subjective/quantitative optimization method and its application to the concept selection of a Small Supersonic Transport. The method uses Genetic Algorithms to operate on a population of designs and promote improvement by varying more than sixty parameters governing the vehicle geometry, mission, and requirements. In addition to using computer codes for evaluation of quantitative criteria such as gross weight, expert input is also considered to account for criteria such as aeroelasticity or manufacturability which may be impossible or too computationally expensive to consider explicitly in the analysis. Results indicate that concepts resulting from the use of this method represent designs which are promising to both the computer and the analyst, and that a mapping between concepts and requirements that would not otherwise be apparent is revealed.

Buonanno, Michael↗
NASA NTRS2010

Stability of Hypersonic Boundary Layers on a Cone at an Angle of Attack

The stability and receptivity of a three-dimensional hypersonic boundary layer over a 7deg half-angle straight cone at an angle of attack of 6deg is numerically investigated at a freestream Mach number of 6.0 and a Reynolds number of 10.4x10(exp 6)/m. The generation and evolution of stationary crossflow vortices are also investigated by performing simulations with three-dimensional roughness elements located on the surface of the cone. The flow fields with and without the roughness elements are obtained by solving the full Navier- Stokes equations in cylindrical coordinates using a fifth-order accurate weighted essentially non-oscillatory (WENO) scheme for spatial discretization and a third-order total-variation-diminishing (TVD) Runge-Kutta scheme for temporal integration. Stability computations produced azimuthal wavenumbers in the range of m approx. 20-50 for the most amplified traveling disturbances and in the range of m approx.30-70 for the stationary disturbances. The frequency of the unstable second-mode ranges from 400 kHz to 900 kHz along the windward ray. The N-Factor computations predicted transition would occur more forward on the sides of the cone as compared to the transition fronts near the windward and the leeward rays. The simulations also show the crossflow vortices originating from the nose region propagate towards the leeward ray. No perturbations were observed toward the windward half of the cone.

Balakumar, P.↗
DOE OSTI2026

Spatially Accelerated Winding Numbers for Curved Geometry

The generalized winding number (GWN) is a scalar field that supports robust containment queries on curved geometry, including non-watertight, overlapping, and nested boundary representations. While queries can be easily parallelized over samples, direct evaluation on parametric curves and surfaces remains costly for large and complex models. Fast, state-of-the-art GWN approaches leverage a spatial index to approximate the GWN, typically coupled with a Taylor expansion which approximates the GWN contribution for far clusters of geometric primitives. However, such methods operate only on discrete inputs such as triangle meshes and point clouds, and would introduce containment errors near boundaries if applied to curved input. We extend support for fast GWN evaluation over arbitrary collections of NURBS curves in 2D and trimmed NURBS patches in 3D via a Bounding Volume Hierarchy that stores efficiently precomputed moment data in the hierarchy nodes. When querying the hierarchy, approximations for far clusters are used alongside direct evaluation for nearby NURBS primitives, achieving sub-linear complexity while preserving the geometric features in the vicinity of the query point. Central to our performance improvements is an adaptive subdivision strategy for NURBS primitives during a preprocessing phase, creating better spatial partitions while retaining the same accuracy for containment decisions as a direct evaluation. We demonstrate the performance and accuracy of our approach across a large collection of 2D and 3D datasets.

Computer science↗
NASA NTRSDate not supplied

Recent BOLT Discrete-Roughness Trip Results from the 20-Inch Mach 6 Tunnel

- Air Force Office of Scientific Research (AFOSR) primary sponsor of BOLT (for Boundary Layer Transition) sounding rocket flight project as continuation of previous HIFiRE flights (1 & 5) - Sounding rockets provide cost-effective flight research tool - Study hypersonic BoLT on increasingly complex geometries - BOLT shape includes concave surfaces and swept leading edges - Two missions, two flights: - BOLT mission at Esrange, Sweden - Johns Hopkins University Applied Physics Laboratory as lead - Focus on transition onset measurements - Flew June 2021 (staging anomaly) - BOLT II mission at Wallops Flight Facility, Virginia - Texas A&M and CUBRC provide primary leadership - Focus on turbulence measurements - Flew March 2022 - Flight named in memory of Mike Holden - NASA’s support to both BOLT missions provided through the Hypersonic Technologies Project (HTP)

Hypersonic↗
NASA NTRSDate not supplied

Precision in Assembled Discrete Lattice Space Structures for Next-Generation ISAM Applications

Robust autonomous robotic assembly of large-scale space structures has been a long-term and challenging goal to enable higher quality in-space communication and science instrumentation. For optical observatory support structures and antenna structures, the challenge is the strict dimensional precision requirements (generally RMS surface error) driven by the operational radiation wavelength. Early theoretical work in the area linked RMS surface error of a reflector support truss plate to the natural passive dynamic mode frequencies of the structure, as well as to the error distribution in dimensions of constituent structural elements. Prior robotic truss assembly demonstrations focused on designing ultra-high precision structural elements, joints, and robotic actuators. Recently, an alternative space structure assembly strategy based on a programmable matter approach (NASA ARMADAS) was demonstrated. This approach uses lattice building blocks (voxels) that are reversibly, mechanically joined into a bulk lattice structure by robots that locomote in and on the structure itself. Such a system can achieve high-level autonomy with low computation, robustly assemble utilizing inexpensive and imprecise robots, and efficiently build structures several orders of magnitude larger than the assembly robots. However, for instrumentation support structure applications, the resulting precision of these building block-based lattice structures is not well studied. Since they are demonstrated with many more assembly units than prior art trusses, it is unclear whether the same precision design approaches apply. In this work, we study in simulation the effects of voxel geometry, error in distribution, assembly resolution (module size), and assembly geometry on the error of both beam and plate lattice structures. While average RMS error of a plate increases with plate size, increasing plate thickness quickly collapses RMS error towards a limit that is on the order of the error of the constituent parts. We validate our models against previously published precision measurements of built systems. Results from this study will guide manufacturing precision requirements, as well as designs, for future robotically assembled structural applications and establish feasibility for different applications.

Christine E Gregg↗