Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical analysis”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

MPCV Exercise Operational Volume Analysis

In order to minimize the loss of bone and muscle mass during spaceflight, the Multi-purpose Crew Vehicle (MPCV) will include an exercise device and enough free space within the cabin for astronauts to use the device effectively. The NASA Digital Astronaut Project (DAP) has been tasked with using computational modeling to aid in determining whether or not the available operational volume is sufficient for in-flight exercise.Motion capture data was acquired using a 12-camera Smart DX system (BTS Bioengineering, Brooklyn, NY), while exercisers performed 9 resistive exercises without volume restrictions in a 1g environment. Data were collected from two male subjects, one being in the 99th percentile of height and the other in the 50th percentile of height, using between 25 and 60 motion capture markers. Motion capture data was also recorded as a third subject, also near the 50th percentile in height, performed aerobic rowing during a parabolic flight. A motion capture system and algorithms developed previously and presented at last years HRP-IWS were utilized to collect and process the data from the parabolic flight [1]. These motions were applied to a scaled version of a biomechanical model within the biomechanical modeling software OpenSim [2], and the volume sweeps of the motions were visually assessed against an imported CAD model of the operational volume. Further numerical analysis was performed using Matlab (Mathworks, Natick, MA) and the OpenSim API. This analysis determined the location of every marker in space over the duration of the exercise motion, and the distance of each marker to the nearest surface of the volume. Containment of the exercise motions within the operational volume was determined on a per-exercise and per-subject basis. The orientation of the exerciser and the angle of the footplate were two important factors upon which containment was dependent. Regions where the exercise motion exceeds the bounds of the operational volume have been identified by determining which markers from the motion capture exceed the operational volume and by how much. A credibility assessment of this analysis was performed in accordance with NASA-STD-7009 prior to delivery to the MPCV program.

mathematical models↗
NASA NTRS2008

A Modeling Technique and Representation of Failure in the Analysis of Triaxial Braided Carbon Fiber Composites

Quasi-static tests have been performed on triaxially braided carbon fiber composite materials with large unit cell sizes. The effects of different fibers and matrix materials on the failure mode were investigated. Simulations of the tests have been performed using the transient dynamic finite element code, LS-DYNA. However, the wide range of failure modes observed for the triaxial braided carbon fiber composites during tests could not be simulated using composite material models currently available within LS-DYNA. A macroscopic approach has been developed that provides better simulation of the material response in these materials. This approach uses full-field optical measurement techniques to measure local failures during quasi-static testing. Information from these experiments is then used along with the current material models available in LS-DYNA to simulate the influence of the braided architecture on the failure process. This method uses two-dimensional shell elements with integration points through the thickness of the elements to represent the different layers of braid along with a new analytical method for the import of material stiffness and failure data directly. The present method is being used to examine the effect of material properties on the failure process. The experimental approaches used to obtain the required data will be described, and preliminary results of the numerical analysis will be presented.

Littell, Justin D.↗
NASA NTRS2008

A Modeling Technique and Representation of Failure in the Analysis of Triaxial Braided Carbon Fiber Composites

Quasi-static tests have been performed on triaxially braided carbon fiber composite materials with large unit cell sizes. The effects of different fibers and matrix materials on the failure mode were investigated. Simulations of the tests have been performed using the transient dynamic finite element code, LS-DYNA. However, the wide range of failure modes observed for the triaxial braided carbon fiber composites during tests could not be simulated using composite material models currently available within LS-DYNA. A macroscopic approach has been developed that provides better simulation of the material response in these materials. This approach uses full-field optical measurement techniques to measure local failures during quasi-static testing. Information from these experiments is then used along with the current material models available in LS-DYNA to simulate the influence of the braided architecture on the failure process. This method uses two-dimensional shell elements with integration points through the thickness of the elements to represent the different layers of braid along with a new analytical method for the import of material stiffness and failure data directly. The present method is being used to examine the effect of material properties on the failure process. The experimental approaches used to obtain the required data will be described, and preliminary results of the numerical analysis will be presented.

Littell, Justin D.↗
NASA NTRS2008

Investigation, Analysis, and Testing of Self-contained Oxygen Generators

Self Contained Oxygen Generators (SCOGs) have widespread use in providing emergency breathing oxygen in a variety of environments including mines, submarines, spacecraft, and aircraft. These devices have definite advantages over storing of gaseous or liquid oxygen. The oxygen is not generated until a chemical briquette containing a chlorate or perchlorate oxidizer and a solid metallic fuel such as iron is ignited starting a thermal decomposition process allowing gaseous oxygen to be produced. These devices are typically very safe to store, easy to operate, and have primarily only a thermal hazard to the operator that can be controlled by barriers or furnaces. Tens of thousands of these devices are operated worldwide every year without major incident. This report examines the rare case of a SCOG whose behavior was both abnormal and lethal. This particular type of SCOG reviewed is nearly identical to a flight qualified version of SCOG slated for use on manned space vehicles. This Investigative Report is a compilation of a NASA effort in conjunction with other interested parties including military and aerospace to understand the causes of the particular SCOG accident and what preventative measures can be taken to ensure this incident is not repeated. This report details the incident and examines the root causes of the observed SCOG behavior from forensic evidence. A summary of chemical and numerical analysis is provided as a background to physical testing of identical SCOG devices. The results and findings of both small scale and full scale testing are documented on a test-by-test basis along with observations and summaries. Finally, conclusions are presented on the findings of this investigation, analysis, and testing along with suggestions on preventative measures for any entity interested in the safe use of these devices.

Keddy, Christopher P.↗
DOE OSTI2025

A detailed study of pre-heating effects in electron beam melting powder bed fusion process

Metal-based additive manufacturing processes, such as powder bed fusion with electron beam (PBF-EB) process, also referred to as electron beam melting (EBM), can produce high-density parts with minimal residual stresses due to the uniform and coherent preheating of the powder bed. However, understanding and controlling the multiple stages of preheating is required to enable the production of high-quality, consistent parts of various materials. This work presents a large-scale, multi-layer, three-dimensional numerical analysis focused on studying the preheating stages for predicting thermal history during the PBF-EB process. The model follows a continuous multi-stage cyclic process, that incorporates all the main stages of the PBF-EB process for 316 L stainless steel. This includes the gradual deposition of a new powder layer, the first and second preheating levels of the powder bed, and the energy deposition during melting (excluding the actual melt-pool behavior simulation). The model employs an adaptive time-scaling approach that automatically adjusts the energy deposition for each solution time-increment. This allows for localized changes in time-resolution over an otherwise computationally expensive multi-layer procedure. The material property variations are also taken into account, with an emphasis on the subtle irreversible changes in powder effective thermal conductivity after the two requisite preheating stages of the powder bed. This effect is studied using simplified conductivity models from the literature for partially sintered powder, validated by a dedicated experiment and numerical simulation. The large-scale model is then used to estimate the actual temperatures during first and second preheating levels for 316 L steel, which is not yet fully supported commercially for PBF-EB. Model predictions are corroborated by experiments, using and analyzing IR images, taken at the completion of each layer by the machine’s built-in infrared camera. The current model also incorporates a qualitative assessment for the effects of conductivity change during pre-heating, as well as evaluates the applicability of the time-scaling approach.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Hybrid Plasmonic Coupled‐Mode Metasurface Design for Shot‐Noise Limited Ppb‐Level Hydrogen Detection

Plasmonic hydrogen sensors have enabled hydrogen detection below parts-per-million (ppm) range by boosting the sensitivity using localized surface plasmonic resonant (LSPR) structures. However, the intrinsic optical losses of Palladium (Pd), the primary plasmonic metal used for hydrogen detection, result in a low quality (Q) factor LSPR, which fundamentally hinders further improvement. In this work, a hybrid plasmonic metasurface is proposed that couples Pd-based LSPR structure with an Au film supporting surface plasmon polariton mode (Au-SPP). The coupled near-perfect absorber resonance yields a spectrally narrow, high Q response that retains strong sensitivity to hydrogen while improving resonance localization. Numerical analysis shows that, under shot-noise-limited conditions, the limit of detection (LoD) can be improved by over threefold compared to the state-of-the-art designs. Furthermore, this hybrid plasmonic coupled-mode metasurface thus presents a promising pathway to achieve parts-per-billion-level (ppb-level) hydrogen detection with enhanced spectral precision and robustness.

coupled modes↗
DOE OSTI2024

On the Departure from Monin–Obukhov Surface Similarity and Transition to the Convective Mixed Layer

Large-eddy simulations are used to evaluate mean profile similarity in the convective boundary layer (CBL). Particular care is taken regarding the grid sensitivity of the profiles and the mitigation of inertial oscillations in the simulation spin-up. The nondimensional gradients Φ for wind speed and air temperature generally align with Monin–Obukhov similarity across cases but have a steeper slope than predicted within each profile. The same trend has been noted in several other recent studies. The Businger-Dyer relations are modified here with an exponential cutoff term to account for the decay in Φ to first-order approximation, yielding improved similarity from approximately 0.05z i to above 0.3z i , where z i is the CBL depth. The necessity for the exponential correction is attributed to an extended transition from surface scaling to zero gradient in the mixed layer, where the departure from Monin–Obukhov similarity may be negligible at the surface but becomes substantial well below the conventional surface layer height of 0.1 z i .

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Resonant Raman in armchair graphene nanoribbons from first-principles

Resonant Raman spectra of armchair graphene nanoribbons (AGNRs) are computed using Density Functional Theory (DFT) and third-order perturbation theory. Results are benchmarked against available experimental data and compared to previously used theoretical approaches based on the Placzek approximation. Comparable agreement with experiments is found for both previously and presently used methods. In addition, a numerical analysis is carried out to provide a justification for the resonant modeling method based on the use of the frequency-dependent dielectric tensor in the Placzek approximation. Finally, this work also provides additional predictions and references for wide AGNRs that might be investigated with Raman scattering experiments in the future.

42 ENGINEERING↗
DOE OSTI2025

Micromechanical response of SiC-OPyC layers in TRISO fuel particles

Tristructural isotropic (TRISO)–coated particle fuel is a proposed fuel for multiple advanced reactor concepts. The performance of the particle depends on whether the silicon carbide (SiC) layer remains intact to prevent the release of metallic and gaseous fission products. Mechanical fracture of the SiC layer is a potential failure mode under various fuel configurations and operating environments, including the potential transmission of matrix-originating cracks through TRISO particles. Furthermore, this study uses instrumented indentation techniques on cross-sectioned surrogate particles to examine the mechanical stability of the critical interface between SiC and the outer pyrolytic carbon (OPyC) layer. The observed behavior at the interface is rationalized by examining the radially dependent fracture behavior of the SiC layer and performing a numerical analysis to quantify the residual stresses that develop during the processing and cross-sectioning of the as-fabricated particle. Characterizing the SiC-OPyC interface of surrogate TRISO particles using nanoindentation provides unique insight into the interface's room-temperature residual stress and mechanical stability. The modeling efforts were used to investigate the experimental procedure further, and the results are presented herein to validate this fuel form's potential mechanical failure modes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Oscillating surge wave energy converter using a novel above-water power takeoff with belt-arc speed amplification

We investigate the performance of a novel power takeoff (PTO) featuring belt-arc speed amplification for oscillating surge wave energy converters (OSWECs), aiming to address the challenges of extremely low rotary speed and large torque under the low-frequency ocean wave excitations. The belt-arc design significantly increases the rotary speed of the generator, enabling generator downsizing and decreasing the powertrain friction losses. The design also allows for placing the generator above water, eliminating the need for high Ingress Protection ratings for the generator and powertrain, and potentially leading to substantial reductions in capital and maintenance costs. Using the linear potential wave theory, the dynamics of the integrated system are analyzed, and key parameters are identified. To validate the numerical analysis, a 1:10 scale model is designed, fabricated, and tested in a wave tank. Performance evaluations are conducted under regular and irregular wave conditions, with quantified effects of parameter tuning. The results reveal an optimal wave-to-electric efficiency of 48% under regular wave excitation and 20% under irregular excitation. Furthermore, these findings underscore the effectiveness of the proposed novel PTO design in addressing the challenges of low rotation speed and large torque inherent in OSWECs, demonstrating its ability to efficiently convert wave power into electricity.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Control-inspired design and power optimization of an active mechanical motion rectifier based power takeoff for wave energy converters

Ocean waves have high energy density and are persistent and predictable. Yet, converting wave energy to a useable form remains challenging. A significant hurdle is the oscillatory nature of waves resulting in the alternating loads, which necessitate the use of rectification at some stage of the energy conversion. This research effort presents a novel design of active mechanical motion rectifier (AMMR) for the power takeoff (PTO), which provides enhanced controllability and better power performance when compared to passive mechanical motion rectifiers (MMR). Inspired by transistors used in synchronous electrical rectifiers, the proposed design uses controllable electromagnetic clutches in the mechanical transmission to allow active engagement-disengagement control; thus, rectifying the oscillatory motion into a unidirectional rotation for high energy conversion efficiency and allowing the generator in unidirectional rotation to control the bidirectional wave capture structure for maximizing the power output. A semi-analytical computational approach is developed to efficiently evaluate the optimal power achieved using the proposed AMMR-based PTO and active control. It is found that the AMMR-based PTO design yields a higher optimal power than the previous passive MMR design across the wave spectrum. The influences of generator inertia and reactive power are discussed. Furthermore, the effects of control parameters on the power output and the optimal trajectories are analyzed. Wave tank tests with the AMMR prototype demonstrated the effectiveness of AMMR based PTO design and validated the numerical analysis.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Task-oriented machine learning surrogates for tipping points of agent-based models

We present a machine learning framework bridging manifold learning, neural networks, Gaussian processes, and Equation-Free multiscale approach, for the construction of different types of effective reduced order models from detailed agent-based simulators and the systematic multiscale numerical analysis of their emergent dynamics. The specific tasks of interest here include the detection of tipping points, and the uncertainty quantification of rare events near them. Our illustrative examples are an event-driven, stochastic financial market model describing the mimetic behavior of traders, and a compartmental stochastic epidemic model on an Erdös-Rényi network. We contrast the pros and cons of the different types of surrogate models and the effort involved in learning them. Importantly, the proposed framework reveals that, around the tipping points, the emergent dynamics of both benchmark examples can be effectively described by a one-dimensional stochastic differential equation, thus revealing the intrinsic dimensionality of the normal form of the specific type of the tipping point. This allows a significant reduction in the computational cost of the tasks of interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Effect of Density Irregularities on Radio Frequency Wave Propagation in Ionospheric Plasmas

Density irregularities play a vital role in determining how radio frequency (RF) waves travel through plasmas. In the Earth’s ionosphere, these density irregularities also impact radio communication. In this study, we conduct a detailed numerical analysis of RF wave propagation in small-scale ionospheric density irregularities using the advanced Petra-M code. We focus specifically on high-frequency (HF) waves, ranging from 3 to 30 MHz, which are essential for military, amateur radio operators, and emergency communications. By introducing density structures, such as equatorial plasma bubbles derived from fluid simulations, we demonstrate that HF waves can scatter in multiple directions when they encounter these irregularities. Additionally, we observe significant mode conversion, where incoming electromagnetic waves transform into electrostatic modes within the density gradient layer. This shows that smaller density irregularities can greatly weaken signals or cause complete signal loss for receivers, emphasizing the need for increased awareness and innovative solutions in radio communication transmission.

Kim, Eun-Hwa [Princeton Plasma Physics Laboratory ↗
DOE OSTI2026

Initial electron cyclotron heating/current drive scoping study for the SMall Aspect Ratio Tokamak (SMART)

In this work an initial electron cyclotron heating/current drive (EC H&CD) scoping study is presented for the SMall Aspect Ratio Tokamak (SMART), constructed and operated by the Plasma Science and Fusion Technology (PSFT) Laboratory of the University of Seville. We consider two planned phases of operation with different magnitude of magnetic field (B 0.4, and 1 T). This EC H&CD scoping study is carried out by using the ray tracing code TRAVIS, which allows to quickly perform a series of scans on the main SMART plasma scenarios. A scan in the poloidal and toroidal launching angles is performed for the B 0.4 T plasma scenario assuming three launching locations (mid-plane, off-mid-plane, and top-launcher) and the extraordinary (X-) mode polarization. The absorption occurs at the second harmonic and, overall, the highest power absorption is obtained launching the EC beam from the top of the machine. A numerical analysis of the fundamental X-mode EC H&CD start-up regime for the B 1.0 T plasma scenario with two different wave frequencies (28 and 35 GHz) is presented showing a very high current drive efficiency.

Bertelli, Nicola [Princeton Plasma Physics Laborat↗
DOE OSTI2024

X-ray sources for in situ wavelength calibration of x-ray imaging crystal spectrometers

X-ray sources for a range of wavelengths are being considered for in situ calibration of X-ray Imaging Crystal Spectrometers (XICSs) and for monitoring line shifts due to changes in the crystal temperature, which can vary during experimental operation over a day [A. Ince-Cushman et al., Rev. Sci. Instrum. 79, 10E302 (2008), L. Delgado-Aparicio et al., Plasma Phys. Control. Fusion 55, 125011 (2013)]. Such crystal temperature dependent shifts, if not accounted for, could be erroneously interpreted as Doppler shifts leading to errors in plasma flow-velocity measurements. The x-ray sources encompass characteristic x-ray lines falling within the wavelength range of 0.9–4.0 Å, relevant for the XICSs on present and future fusion devices. Several technological challenges associated with the development of x-ray sources for in situ calibration are identified and are being addressed in the design of multiple x-ray tubes, which will be installed inside the spectrometer housing of the XICS for the JT-60SA tokamak. These x-ray sources will be especially useful for in situ calibration between plasma discharges. Here, in this paper, laboratory experiments are described that were conducted with a Cu x-ray source, a heated quartz (102) crystal, and a pixelated Pilatus detector to measure the temperature dependent shifts of the Cu Kα 1 and Kα 2 lines at 1.5405 and 1.5443 Å, respectively, and to evaluate the 2d-lattice constant for the Bragg reflecting crystal planes as a function of temperature, which, in the case of in situ wavelength calibration, would have to be used for numerical analysis of the x-ray spectra from the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Time Correlations from Steady-State Expectation Values

Recovering properties of correlation functions is typically challenging. On the one hand, experimentally, it requires measurements with a temporal resolution finer than the system’s dynamics. On the other hand, analytical or numerical analysis requires solving the system evolution. Here, we use recent results of quantum metrology with continuous measurements to derive general lower bounds on the relaxation and second-order correlation times that are both easy to calculate and measure. These bounds are based solely on steady-state expectation values and their derivatives with respect to a system parameter, and can be readily extended to the autocorrelation of arbitrary observables. We validate our method on two examples of critical quantum systems: a critical driven-dissipative resonator, where the bound matches analytical results for the dynamics, and the infinite-range Ising model, where only the steady state is solvable, and thus the bound provides information beyond the reach of existing analytical approaches. Our results can be applied to the experimental characterization of ultrafast systems and to the theoretical analysis of many-body models whose dynamics are hard to compute.

Górecki, Wojciech [INFN, Pavia] (ORCID:00000001991↗
DOE OSTI2026

Tensor-polarized parton density in the 𝑁 → Δ transition from the large-𝑁 𝑐 light-cone wave function

The tensor-polarized parton density is defined by the forward matrix element of a partonic operator in the 𝑁 → Δ transition. In this work, we investigate it by employing the large-𝑁 𝑐 light-cone wave function derived from the mean-field approach. The mean-field picture is based on low-energy effective dynamics in the large-𝑁 𝑐 limit, where the baryon wave function is formulated in the rest frame. By exploiting the covariance of the mean-field solution, we derive the corresponding large-𝑁 𝑐 light-cone wave function—decomposed unambiguously into 3⁢𝑄 , 5⁢𝑄 , 7⁢𝑄 , and higher Fock components—in the infinite momentum frame. Evaluating the overlap of these wave functions, we derive an overlap representation of the tensor-polarized parton density in the 𝑁 → Δ transition and find that the leading contribution arises from the 5⁢𝑄 Fock sector. This indicates that the tensor-polarized parton density directly probes the genuine 5⁢𝑄 component and is governed by chiral dynamics. Our numerical analysis shows that the 𝑁 → Δ tensor-polarized parton density is suppressed, consistent with standard large-𝑁 𝑐 expectations. Finally, we establish connections among the tensor-polarized parton density, the generalized parton distribution 𝐻 𝑋 , and the energy-momentum tensor form factor 𝐹 4 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Higher-order Van Hove singularities in kagome topological bands

Motivated by the growing interest in band structures featuring higher-order Van Hove singularities (HOVHS), we investigate a spinless fermion kagome system characterized by nearest-neighbor (NN) and next-nearest-neighbor (NNN) hopping amplitudes. While NN hopping preserves time-reversal symmetry, NNN hopping, akin to chiral hopping on the Haldane lattice, breaks time-reversal symmetry and leads to the formation of topological bands with Chern numbers ranging from 𝐶 = ±1 to ±4. We perform analytical and numerical analysis of the energy bands near the high-symmetry points Γ, ±𝐊, and 𝐌 𝑖 (𝑖 = 1, 2, and 3), which uncover a rich and complex landscape of HOVHS, controlled by the magnitude and phase of the NNN hopping. We observe power-law divergences in the density of states (DOS), 𝜌⁡(𝜀)∼|𝜀| −𝜈 , with exponents 𝜈 = 1/2, 1/3, 1/4, which can significantly affect the anomalous Hall response at low temperatures when the Fermi level crosses the HOVHS. Additionally, the NNN hopping induces the formation of higher Chern number bands 𝐶 = ±2, ±4 in the middle of the spectrum obeying a sublattice interference whereupon electronic states are maximally localized in each of the sublattices when the momentum approaches the three high-symmetry points 𝐌 𝑖 (𝑖 = 1, 2, and 3) on the Brillouin zone boundary. Finally, this classification of HOVHS in kagome systems provides a platform to explore unconventional electronic orders induced by electronic correlations.

Chern insulators↗