Search NASA⌕ Search

SEARCH · Search NASA

Results for “ORTHOGONAL FUNCTION”

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 361 records · Page 20

ZERNIPAX: A fast and accurate Zernike polynomial calculator in Python

Zernike polynomials serve as an orthogonal basis on the unit disc, and have proven to be effective in optics simulations, astrophysics, and more recently in plasma simulations. Unlike Bessel functions, Zernike polynomials are inherently finite and smooth at the disc center (r=0), ensuring continuous differentiability along the axis. This property makes them particularly suitable for simulations, requiring no additional handling at the origin. We developed ZERNIPAX, an open-source Python package capable of utilizing CPU/GPUs, leveraging Google's JAX package and available on GitHub as well as the Python software repository PyPI. Furthermore, our implementation of the recursion relation between Jacobi polynomials significantly improves computation time compared to alternative methods by use of parallel computing while still performing more accurately for high-mode numbers.

Astrophysics↗

Improving the Quasi‐Biennial Oscillation via a Surrogate‐Accelerated Multi‐Objective Optimization

Accurate simulation of the quasi-biennial oscillation (QBO) is challenging due to uncertainties in representing convectively generated gravity waves. We develop an end-to-end uncertainty quantification workflow that calibrates these gravity wave processes in E3SM for a realistic QBO. Central to our approach is a domain knowledge-informed, compressed representation of high-dimensional spatio-temporal wind fields. By employing a parsimonious statistical model that learns the fundamental frequency from complex observations, we extract interpretable and physically meaningful quantities capturing key attributes. Building on this, we train a probabilistic surrogate model that approximates the fundamental characteristics of the QBO as functions of critical physics parameters governing gravity wave generation. Leveraging the Karhunen–Loève decomposition, our surrogate efficiently represents these characteristics as a set of orthogonal features, capturing cross-correlations among multiple physics quantities evaluated at different pressure levels and enabling rapid surrogate-based inference at a fraction of the computational cost of full-scale simulations. Finally, we analyze the inverse problem using a multi-objective approach. Our study reveals a tension between amplitude and period that constrains the QBO representation, precluding a single optimal solution. To navigate this, we quantify the bi-criteria trade-off and generate a set of Pareto optimal parameter values that balance the conflicting objectives. This integrated workflow improves the fidelity of QBO simulations and offers a versatile template for uncertainty quantification in complex geophysical models.

54 ENVIRONMENTAL SCIENCES↗

On the OLP prediction of the unstable modes of the flat plate turbulent boundary layer

The Orr et al. (OLP) stability analysis is reformulated to yield profiles from hot-wire velocity covariance data in the boundary layer of a flow over a flat plate. The proper orthogonal decomposition theorem (PODT) is used to maximize the mean square of the inner product of the velocity vector with a deterministic candidate vector function in Hilbert space. The OLP formulation couples the PODT for the eddy viscosity with a global extremum principle for the disturbance kinetic energy. Derivation of the governing equations is described analytically, as are the extraction of the mean velocity profile and the eddy viscosity. Finally, techniques for characterizing the dominant mode of the velocity field are introduced and applied to three-dimensional flow over a flat plate. Good agreement is found between the predictions, using laboratory data, and the measured energy and three-dimensional structures of the dominant eddies.

Wadia, A. R.↗

Spatial orientation and balance control changes induced by altered gravitoinertial force vectors

To better understand the mechanisms of human adaptation to rotating environments, we exposed 19 healthy subjects and 8 vestibular-deficient subjects ("abnormal"; four bilateral and four unilateral lesions) to an interaural centripetal acceleration of 1 g (resultant 45 degrees roll-tilt of 1.4 g) on a 0.8-m-radius centrifuge for periods of 90 min. The subjects sat upright (body z-axis parallel to centrifuge rotation axis) in the dark with head stationary, except during 4 min of every 10 min, when they performed head saccades toward visual targets switched on at 3- to 5-s intervals at random locations (within +/- 30 degrees) in the earth-horizontal plane. Eight of the normal subjects also performed the head saccade protocol in a stationary chair adjusted to a static roll-tilt angle of 45 degrees for 90 min (reproducing the change in orientation but not the magnitude of the gravitoinertial force on the centrifuge). Eye movements, including voluntary saccades directed along perceived earth- and head-referenced planes, were recorded before, during, and immediately after centrifugation. Postural center of pressure (COP) and multisegment body kinematics were also gathered before and within 10 min after centrifugation. Normal subjects overestimated roll-tilt during centrifugation and revealed errors in perception of head-vertical provided by directed saccades. Errors in this perceptual response tended to increase with time and became significant after approximately 30 min. Motion-sickness symptoms caused approximately 25% of normal subjects to limit their head movements during centrifugation and led three normal subjects to stop the test early. Immediately after centrifugation, subjects reported feeling tilted 10 degrees in the opposite direction, which was in agreement with the direction of their earth-referenced directed saccades. Postural COP, segmental body motion amplitude, and hip-sway frequency increased significantly after centrifugation. These postural effects were short-lived, however, with a recovery time of several postural test trials (minutes). There were also asymmetries in the direction of postcentrifugation COP and head tilt which depended on the subject's orientation during the centrifugation adaptation period (left ear or right ear out). The amount of total head movements during centrifugation correlated poorly or inversely with postcentrifugation postural stability, and the most unstable subject made no head movements. There was no decrease in postural stability after static tilt, although these subjects also reported a perceived tilt briefly after return to upright, and they also had COP asymmetries. Abnormal subjects underestimated roll-tilt during centrifugation, and their directed saccades revealed permanent spatial distortions. Bilateral abnormal subjects started out with poor postural control, but showed no postural decrements after centrifugation, while unilateral abnormal subjects had varying degrees of postural decrement, both in their everyday function and as a result of experiencing the centrifugation. In addition, three unilateral, abnormal subjects, who rode twice in opposite orientations, revealed a consistent orthogonal pattern of COP offsets after centrifugation. These results suggest that both orientation and magnitude of the gravitoinertial vector are used by the central nervous system for calibration of multiple orientation systems. A change in the background gravitoinertial force (otolith input) can rapidly initiate postural and perceptual adaptation in several sensorimotor systems, independent of a structured visual surround.

Clinical Trial↗

An Algorithm for Atom-Centered Lossy Compression of the Atomic Orbital Basis in Density Functional Theory Calculations

Large atomic-orbital (AO) basis sets of at least triple and preferably quadruple-ζ (QZ) size are required to adequately converge Kohn–Sham density functional theory (DFT) calculations toward the complete basis set limit. However, incrementing the cardinal number by one nearly doubles the AO basis dimension, and the computational cost scales as the cube of the AO dimension, so this is very computationally demanding. Here, in this work, we develop and test a threshold-based natural atomic orbital (NAO) scheme in which ϵ-NAOs are obtained as eigenfunctions of atomic blocks of the density matrix in a one-center orthogonalized representation. This enables compression of the AO basis that is optimal for a given threshold, 10 –ϵ , by discarding NAOs with occupation numbers below that threshold. Extensive pilot test calculations using the Hartree–Fock functional and taking the converged density matrix as input suggest that a threshold of 10 –5 can yield a compression factor (ratio of AO to compressed ϵ-NAO dimension) between 2.5 and 4.5 for the QZ pc-3 basis. The errors in relative energies are typically less than 0.1 kcal/mol when the compressed basis is used instead of the uncompressed basis. Between 10 and 100 times smaller errors (i.e., usually less than 0.01 kcal/mol) can be obtained with a threshold 10 –7 , while the compression factor is typically between 2 and 2.5.

basis sets↗

A Concept for a Mobile Remote Manipulator System

A conceptual design for a Mobile Remote Manipulator System (MRMS) is presented. This concept does not require continuous rails for mobility (only guide pins at truss hardpoints) and is very compact, being only one bay square. The MRMS proposed is highly maneuverable and is able to move in any direction along the orthogonal guide pin array under complete control at all times. The proposed concept would greatly enhance the safety and operational capabilities of astronauts performing EVA functions such as structural assembly, payload transport and attachment, space station maintenance, repair or modification, and future spacecraft construction or servicing. The MRMS drive system conceptual design presented is a reasonably simple mechanical device which can be designed to exhibit high reliability. Developmentally, all components of the proposed MRMS either exist or are considered to be completely state of the art designs requiring minimal development, features which should enhance reliability and minimize costs.

Mikulus, M. M., Jr.↗

Updating the Finite Element Model of the Aerostructures Test Wing using Ground Vibration Test Data

Improved and/or accelerated decision making is a crucial step during flutter certification processes. Unfortunately, most finite element structural dynamics models have uncertainties associated with model validity. Tuning the finite element model using measured data to minimize the model uncertainties is a challenging task in the area of structural dynamics. The model tuning process requires not only satisfactory correlations between analytical and experimental results, but also the retention of the mass and stiffness properties of the structures. Minimizing the difference between analytical and experimental results is a type of optimization problem. By utilizing the multidisciplinary design, analysis, and optimization (MDAO) tool in order to optimize the objective function and constraints; the mass properties, the natural frequencies, and the mode shapes can be matched to the target data to retain the mass matrix orthogonality. This approach has been applied to minimize the model uncertainties for the structural dynamics model of the Aerostructures Test Wing (ATW), which was designed and tested at the National Aeronautics and Space Administration (NASA) Dryden Flight Research Center (DFRC) (Edwards, California). This study has shown that natural frequencies and corresponding mode shapes from the updated finite element model have excellent agreement with corresponding measured data.

Lung, Shun-fat↗

Updating the Finite Element Model of the Aerostructures Test Wing Using Ground Vibration Test Data

Improved and/or accelerated decision making is a crucial step during flutter certification processes. Unfortunately, most finite element structural dynamics models have uncertainties associated with model validity. Tuning the finite element model using measured data to minimize the model uncertainties is a challenging task in the area of structural dynamics. The model tuning process requires not only satisfactory correlations between analytical and experimental results, but also the retention of the mass and stiffness properties of the structures. Minimizing the difference between analytical and experimental results is a type of optimization problem. By utilizing the multidisciplinary design, analysis, and optimization (MDAO) tool in order to optimize the objective function and constraints; the mass properties, the natural frequencies, and the mode shapes can be matched to the target data to retain the mass matrix orthogonality. This approach has been applied to minimize the model uncertainties for the structural dynamics model of the aerostructures test wing (ATW), which was designed and tested at the National Aeronautics and Space Administration Dryden Flight Research Center (Edwards, California). This study has shown that natural frequencies and corresponding mode shapes from the updated finite element model have excellent agreement with corresponding measured data.

Lung, Shun-Fat↗

High-Speed Schlieren Analysis of Retropropulsion Jet in Mach 10 Flow

High-speed schlieren imaging of a 5-in.-diam, 70° sphere-cone model with a single centerline 15° half-angle nozzle was performed in a Mach 10 flow. Image sequences captured at 100 kHz were obtained with the nozzle plugged (fully blunted model), no nozzle flow, and with nozzle flow over a range of supply pressures. The freestream unit Reynolds number was kept constant for all runs at3.38×106m−1. Time-average information pertaining to the location, width, and stand-off distance of the bow shock, interface, terminal shock, and triple point as a function of jet pressure ratio are presented based on analysis of covariance images from each run. Measurement of jet plume boundary shape as a function of jet pressure ratio is also presented. Analyses of frequency content associated with the jet structure are shown, as are results of the proper orthogonal decomposition of these data. For all runs with nozzle flow, a dominant fundamental frequency of 2 kHz was observed. Finally, a discussion of bow shock unsteadiness resulting from interaction with freestream disturbances for the blunted and no nozzle flow is provided, as is a discussion on the formation of weak waves emanating from the bow shock resulting from this interaction.

Brett F Bathel↗

Polarization‐Engineered Near‐Field Generation Using a Hybrid Tip–Antenna System

Precise control of light polarization at the nanoscale is critical for accessing chiral optical responses and manipulating spin–photon interactions in advanced materials. Yet, conventional scattering-type near-field probes predominantly generate out-of-plane linear polarization and offer little control over phase or polarization state. Here, in this study, we introduce a polarization-engineered near-field methodology based on a combined metallic tip and planar dipole nanoantenna system. Using full-wave electromagnetic simulations, we show that the tip acts as a vertically oriented plasmonic resonator, while the antenna supports an in-plane dipolar mode. By tuning the tip–antenna geometry and tip height, the two orthogonal field components attain comparable amplitudes and a controllable ∼90° phase offset, producing circularly polarized nano-light in the antenna gap. The proposed system effectively functions as a nanoscale quarter-wave plate, converting linearly polarized illumination into circularly polarized hotspots without external polarization optics. This method establishes an experimentally accessible route toward polarization-programmable near-field nanoscopy, enabling chiral spectroscopy, selective excitation of spin/valley degrees of freedom, and quantum optical investigations at the nanoscale.

36 MATERIALS SCIENCE↗

Improving NOvA's Sterile Neutrino Search with the Booster Neutrino Beam

The NOvA experiment’s most recent search for eV-scale sterile neutrinos is systematically limited in the region of parameter space where $\Delta m^2_{41} \gtrsim 1~\mathrm{eV}^2$. This region of parameter space is preferred by sterile neutrino interpretations of current experimental anomalies; improving sensitivity here is high-priority. When added directly into the fit, additional data samples which are subject to orthogonal systematic uncertainties act as in-situ constraints, breaking the degeneracy between systematic uncertainties and sterile-induced oscillations. The NOvA experiment consists of two functionally identical detectors, 14.6 mrad off-axis of the NuMI beam, with the Near Detector (Far Detector) 1 km (810 km) from the beam source. The Near Detector’s position on-site at Fermilab means that it is also able to observe neutrinos from a second neutrino beam, the BNB, 160 mrad off-axis. NOvA has been taking BNB data since 2015, but has not yet analysed these data. The BNB and NuMI are subject to different beam-related uncertainties, allowing us to leverage this sample as an in-situ constraint. This poster will present the current status, preliminary simulations, and potential additional uses of this unique experimental setup

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Rewritable Surface‐Grafted Polymer Brushes with Dynamic Covalent Linkages

Abstract Surface grafting of polymer brushes drastically modifies surface properties, including wettability, compatibility, responsiveness, etc. A broad variety of functionalities can be introduced to the surface via different types of polymers. Bringing together properties of two or more types of polymer brushes to one surface opens up even more possibilities in brush‐modified materials. However, while it is generally feasible to introduce several chemical compositions along the brushes via copolymerization, it is challenging to vary the types of polymer brushes along a surface. Although previous studies have demonstrated binary brushes via orthogonal polymerization techniques or partial deactivation/regrafting, they commonly limit the number of polymer types to two. Here, we propose a strategy to introduce dynamic covalent diketoenamine linkages at the root of polymer brushes. The grafting density could be precisely tuned during surface functionalization. The surface‐anchored polymer brushes were cleaved by the addition of small molecule amines. New polymer brushes can be regrafted from the surface following refunctionalization of exposed sites. The maneuverability allows tuning of the types and densities of the polymer brushes, pointing the way to the preparation of a new generation of well‐defined brush‐modified materials with mixed grafts, with potential applications in the design of smart materials and surfaces.

Cui, Feichen↗

On-lattice kinetic Monte Carlo approaches for modeling molecular anisotropy in resveratrol crystallization

Stilbenes are a class of organic compounds with broad-ranging pharmaceutical and agricultural applications, which are typically isolated and purified through recrystallization. We are motivated by reducing experimental waste and optimizing yield via developing predictive simulations for processing-dependent crystal morphologies. Using resveratrol as a model stilbene system, we have developed an approach for simulating crystallization with molecular resolution using on-lattice kinetic Monte Carlo. In this work, we highlight modifications to the Stochastic Parallel PARticle Kinetic Simulator (SPPARKS) software package, which were essential to this application. Key enhancements include the incorporation of non-orthogonal cell shapes and monomer anisotropy approximations using bound hard spheres. This new SPPARKS application has been applied to resveratrol with attachment energy libraries obtained from density functional theory, resulting in excellent agreement with experimental morphology prediction.

crystallization↗

Numerical construction of the Hill functions.

As an aid in the numerical construction of Hill functions and their derivatives, an algorithm using local coordinates and an expansion in Legendre polynomials is proposed. The algorithm is shown to possess sufficient stability, and the orthogonality of the Legendre polynomials simplifies the computation when the Ritz-Galerkin technique is used.

Segethova, J.↗

Design of a walking robot

Carnegie Mellon University's Autonomous Planetary Exploration Program (APEX) is currently building the Daedalus robot; a system capable of performing extended autonomous planetary exploration missions. Extended autonomy is an important capability because the continued exploration of the Moon, Mars and other solid bodies within the solar system will probably be carried out by autonomous robotic systems. There are a number of reasons for this - the most important of which are the high cost of placing a man in space, the high risk associated with human exploration and communication delays that make teleoperation infeasible. The Daedalus robot represents an evolutionary approach to robot mechanism design and software system architecture. Daedalus incorporates key features from a number of predecessor systems. Using previously proven technologies, the Apex project endeavors to encompass all of the capabilities necessary for robust planetary exploration. The Ambler, a six-legged walking machine was developed by CMU for demonstration of technologies required for planetary exploration. In its five years of life, the Ambler project brought major breakthroughs in various areas of robotic technology. Significant progress was made in: mechanism and control, by introducing a novel gait pattern (circulating gait) and use of orthogonal legs; perception, by developing sophisticated algorithms for map building; and planning, by developing and implementing the Task Control Architecture to coordinate tasks and control complex system functions. The APEX project is the successor of the Ambler project.

Whittaker, William↗

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↗

Trigonometric Polynomials For Estimation Of Spectra

Orthogonal sets of trigonometric polynomials used as suboptimal substitutes for discrete prolate-spheroidal "windows" of Thomson method of estimation of spectra. As used here, "windows" denotes weighting functions used in sampling time series to obtain their power spectra within specified frequency bands. Simplified windows designed to require less computation than do discrete prolate-spheroidal windows, albeit at price of some loss of accuracy.

Greenhall, Charles A.↗

NASA supporting studies for microgravity research on eye movements

The purpose of the work on this project was to provide support for ground-based studies on the effects of gravity on eye movements. The effects of microgravity on the optokinetic eye movements of humans are investigated. OKN was induced by having subjects watch 3.3 deg stripes moving at 35 deg/s for 45 s in a binocular, head-fixed apparatus. The field (hor., 88 deg; vert., 72 deg), was rotated about axes that were upright or tilted 45 deg or 90 deg. The head was upright or tilted 45 deg on the body. Head-horizontal (yaw axis) and head-vertical (pitch axis) components of OKN were recorded with electro-oculography (EOG). Slow phase velocity vectors were determined relative to gravity. With the head upright, the axis of eye rotation during yaw axis OKN was coincident with the stimulus axis and the spatial vertical. With the head tilted 45 deg on the body, a persistent vertical component of eye velocity developed during yaw axis stimulation, and there was an average shift of the axis of eye rotation toward the spatial vertical of approximately 18 deg in six subjects. During oblique optokinetic stimulation with the head upright, the axis of eye rotation shifted 12 deg toward the spatial vertical. When the head was tilted, the axis of eye rotation rotated to the other side of the spatial vertical by 5.4 deg during the same oblique stimulation. This counter-rotation of the axis of eye rotation is similar to the 'Muller (E) effect', in which the perception of the upright counter-rotates to the opposite side of the spatial vertical when subjects are tilted in darkness. The data were simulated by a model of OKN. Despite the short OKAN time constants, strong horizontal to vertical cross-coupling was produced if the horizontal and vertical time constants were in proper ratio, and there was no suppression of nystagmus orthogonal to the stimulus direction. This shows that the spatial orientation of OKN can be due to a restructuring of the system matrix of velocity storage as a function of gravity. It is concluded that although human OKAN is weak, velocity storage orients the slow phase velocity of OKN towards the spatial vertical.

Cohen, Bernard↗