Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “criticality simulations”

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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Measurements of the Ionization Coefficient of Simulated Iron Micrometeoroids

The interpretation of meteor radar observations has remained an open problem for decades. One of the most critical parameters to establish the size of an incoming meteoroid from radar echoes is the ionization coefficient, beta, which still remains poorly known. Here we report on new experiments to simulate micrometeoroid ablation in laboratory conditions to measure beta for iron particles impacting N2, air,CO2, and He gases. This new data set is compared to previous laboratory data where we find agreement except for He and air impacts greater than 30 kms. We calibrate the Jones model of beta(v) and provide fit parameters to these gases and find agreement with all gases except CO2 and high-speed air impacts where we observe beta(sub air) greater than 1 for velocities greater than 70 kms. These data therefore demonstrate potential problems with using the Jones model for CO2 atmospheres as well as for high-speed meteors on Earth.

Thomas, Evan↗
NASA NTRSDate not supplied

Validation of Fitness for Duty Standards Using Pre- and Post-Flight Capsule Egress and Suited Functional Performance Tasks in Simulated Reduced Gravity (Pilot Egress Fitness)

As NASA prepares for exploration missions, it will be critical to understand and build upon known crew capabilities. Specifically, we need to know the functional capacities of crew at different stages of a mission and characterize these profiles so that the appropriate vehicle requirements and mission concepts of operations (ConOps) can be developed. One of the most complex phases of any missions is when there is a transition between gravity environments. For instance, both physiological adaptation to microgravity and subsequent re-entry into a gravity environment result in reduced functional capacity, even with rigorous adherence to inflight countermeasures. Quantification of the astronauts’ post-landing functional performance is necessary to design ConOps for exploration missions. Specifically, these two high-risk tasks may have to be performed soon after gravity transitions: •Unassisted capsule egress task after return to Earth •Planetary extravehicular activity (EVA) soon after landing on Mars or the Moon This study has been broken down into two phases. A pilot phase to assess the overall feasibility and demonstrate the capability to do these tasks shortly after landing and the full Egress Fitness study, which will be part of the CIPHER complement. This study uses a functional task approach to characterize performance of these high-risk tasks in long-duration ISS crewmembers before flight and shortly after return to Earth. Prior to any testing, each astronaut subject completes a suit fit check to ensure adequate sizing and mobility to be able to complete the EVA tasks. Pilot Egress Fitness pre-flight and post-flight testing includes an Earth based emergency egress out of a functional capsule mockup, and a short Mars gravity EVA simulation including suit donning, hatch egress, ladder descent, task board cable operations, baggage transfer over sand/rocky regolith, alignment with a rear entry port, and suit egress. Pre-flight testing can occur almost anytime pre-flight, but post-flight scheduling is much more critical with the capsule egress test occurring 1-4 hours after landing and the planetary EVA approximately 18-36 hours after landing. Data collected for both tasks includes task completion time, photo, and video. The EVA portion also includes collection of metabolic and heart rate. Pilot Egress Fitness study has completed baseline pre-flight testing on 3 astronaut subjects. Post-flight testing of the first two subjects will be completed in November 2021 and the final subject around March 2022. Results and lessons learned from this pilot study will inform the full Egress Fitness study, which will incorporate additional pre-flight sessions, longer EVA tasks, and post-flight testing on R+1, 4, and 7 to characterize the timeframe of recovery

J R Norcross↗
NASA NTRS2011

The Plasma Wake Downstream of Lunar Topographic Obstacles: Preliminary Results from 2D Particle Simulations

Anticipating the plasma and electrical environments in permanently shadowed regions (PSRs) of the moon is critical in understanding local processes of space weathering, surface charging, surface chemistry, volatile production and trapping, exo-ion sputtering, and charged dust transport. In the present study, we have employed the open-source XOOPIC code [I] to investigate the effects of solar wind conditions and plasma-surface interactions on the electrical environment in PSRs through fully two-dimensional pattic1e-in-cell simulations. By direct analogy with current understanding of the global lunar wake (e.g., references) deep, near-terminator, shadowed craters are expected to produce plasma "mini-wakes" just leeward of the crater wall. The present results (e.g., Figure I) are in agreement with previous claims that hot electrons rush into the crater void ahead of the heavier ions, fanning a negative cloud of charge. Charge separation along the initial plasma-vacuum interface gives rise to an ambipolar electric field that subsequently accelerates ions into the void. However, the situation is complicated by the presence of the dynamic lunar surface, which develops an electric potential in response to local plasma currents (e.g., Figure Ia). In some regimes, wake structure is clearly affected by the presence of the charged crater floor as it seeks to achieve current balance (i.e. zero net current to the surface).

Zimmerman, Michael I.↗
DOE OSTI2020

Validation and application of a multiphase CFD model for hydrodynamics, temperature field and RTD simulation in a pilot-scale biomass pyrolysis vapor phase upgrading reactor

Accurate prediction of transport phenomena is critical for VPU reactor design, optimization, and scale-up. The current study focused on the validation and application of a multiphase CFD model within an open-source code MFiX for hydrodynamics, temperature field, and residence time distribution (RTD) simulation in a non-reacting circulating fluidized bed riser for biomass pyrolysis vapor phase upgrading (VPU). First, an Eulerian-Eulerian approach three-dimensional CFD model was employed to simulate the pilot-scale VPU riser on the supercomputer Joule. Excellent quantitative agreement between experimental and simulated results was achieved for pressure drops and temperature field in a range of operating conditions. Then the validated multiphase CFD model was applied to predict gas and solid residence time distributions (RTDs) since prediction and analysis of RTD is an important tool to study the complex multiphase flow behavior and mixing inside chemical reactors. The predictions show that solid mean residence time is 3.5 times the gas residence time; the solid RTD is more sensitive to the process gas flow rate than the solids circulation rate.

09 BIOMASS FUELS↗
DOE OSTI2021

Adaptive critic design-based reinforcement learning approach in controlling virtual inertia-based grid-connected inverters

In this report, an adaptive critic design (ACD) approach is proposed to control the phase and voltage of a grid-connected virtual synchronous generator (VSG). The penetration of fast responding inertia-less power converters significantly affect the stability of the power system, especially weak systems such as micro grids. The concept of virtual inertia addresses this concern by virtually emulating the behavior of a synchronous generator. However, the conventional VSG is designed based on two conditions: (i) fixed operating point and (ii) inductive grid connections. The performance of VSGs in low-voltage semi-resistive microgrids is far from optimal. To overcome the aforementioned concerns, a heuristic dynamic programing (HDP) approach is proposed to optimally control grid-connected VSGs. The neural-network-based inherence of the HDP enables the proposed technique to adapt to any impedance angle. The HDP controller includes two subnetworks: (i) the action network that controls the system optimally and (ii) the critic network, which evaluates the effectiveness of the action network. The simulation and experimental results are provided to evaluate the effectiveness of the proposed technique. As shown, the HDP-based approach illustrates a better performance in comparison with the conventional PI-based VSG in various operating conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Hybrid simulations of sub-cyclotron compressional and global Alfvén eigenmode stability in spherical tokamaks

A comprehensive numerical study has been conducted in order to investigate the stability of beam-driven, sub-cyclotron-frequency compressional Alfvén eigenmodes (CAEs) and global Alfvén eigenmodes (GAEs) in low-aspect-ratio plasmas for a wide range of beam parameters. The presence of CAEs and GAEs has previously been linked to anomalous electron temperature profile flattening at high beam powers in NSTX experiments, prompting a further examination of the conditions necessary for their excitation. Linear simulations have been performed with the hybrid MHD–kinetic initial value code HYM in order to capture the general Doppler-shifted cyclotron resonance that drives the modes. Additionally, three distinct types of modes were found in the simulations—co-CAEs, cntr-GAEs, and co-GAEs—with differing spectral and stability properties. The simulations revealed that unstable GAEs are more ubiquitous than unstable CAEs, which is consistent with experimental observations, as they are excited at lower beam energies and generally have larger growth rates. Local analytic theory is used to explain key features of the simulation results, including the preferential excitation of different modes based on beam injection geometry and the growth rate dependence on the beam injection velocity, critical velocity, and degree of velocity space anisotropy. The background damping rate is inferred from the simulations and analytically estimated for relevant sources absent from the simulation model, indicating that co-CAEs are closer to marginal stability than modes driven by the cyclotron resonances.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Theory and simulation of gap distance transitions in planar magnetically insulated transmission lines

Magnetically insulated transmission lines (MITLs) are critical components for efficiently transporting power from pulsed power sources to a diverse range of loads, including magneto-inertial confinement fusion and radiation production. This study investigates the behavior of planar MITLs, taking into account the source voltage and source and load impedances, using two-dimensional theory and simulations. Self-limited behavior, where the MITL operating impedance (voltage divided by current) is insensitive to the value of load impedance, is extended to small-to-large gap distance MITL transitions. Load-limited behavior, in which the load impedance is below the self-limited impedance, causes the MITL operating impedance to adjust to match the load impedance and can be extended to describe large-to-small transitions. In both cases, the operating impedance of the larger MITL section changes to match the self-limited operating impedance of the smaller section. This behavior extends to multiple transitions and to both low and high impedance loads. Assuming equilibration occurs (pulse duration larger than MITL transit time), the entire MITL operating impedance becomes the self-limited impedance of the smallest-gap MITL section or the load impedance, whichever is lower.

Darr, Adam M. [Sandia National Laboratories (SNL-N↗
DOE OSTI2022

Dispersion-enhanced sequential batch sampling for adaptive contour estimation

In computer simulation and optimal design, sequential batch sampling offers an appealing way to iteratively stipulate optimal sampling points based upon existing selections and efficiently construct surrogate modeling. Nonetheless, the issue of near duplicates poses tremendous quandary for sequential learning. It refers to the situation that selected critical points cluster together in each sampling batch, which are individually but not collectively informative towards the optimal design. Near duplicates severely diminish the computational efficiency as they barely contribute extra information towards update of the surrogate. To address this issue, we impose a dispersion criterion on concurrent selection of sampling points, which essentially forces a sparse distribution of critical points in each batch, and demonstrate the effectiveness of this approach in adaptive contour estimation. Specifically, we adopt Gaussian process surrogate to emulate the simulator, acquire variance reduction of the critical region from new sampling points as a dispersion criterion, and combine it with the modified expected improvement (EI) function for critical batch selection. The critical region here is the proximity of the contour of interest. This proposed approach is vindicated in numerical examples of a two-dimensional four-branch function, a four-dimensional function with a disjoint contour of interest and a time-delay dynamic system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Simulation of Particulate Transport for Delivery of Solid Amendments into the Subsurface: FY24 Status Report

For particulate-based amendments to be viable for field-scale remediation at the Hanford Site (e.g., 200 DV-1 Operable Unit), particles need to be delivered a sufficient radial distance from an injection well and retained at concentrations high enough for effective treatment. An accurate description of the particle radius of influence (ROI) is critical for developing an overall remediation strategy. However, field-scale particle simulations are currently limited due to insufficient simulation capabilities and a lack of experimental data to validate and parameterize particle transport models. To help build toward field-scale deployment, this fiscal year (FY) we have (1) developed a pre screening tool to estimate particle transport, (2) implemented particle transport models within PFLOTRAN, and (3) conducted preliminary estimations of particle ROI. While field-scale numerical simulations will ultimately be necessary before remedy design and field implementation, we have developed a pre-screening tool that offers valuable estimations of expected particle injectability and ROI in a 1-D system. The advantage of the tool is that it does not require extensive laboratory experiments and instead makes predictions based solely on routine laboratory measurements. This tool can assist in down-selection and decision-making by identifying which particle amendment systems are worth pursuing in future laboratory experiments, such as 1-D column tests and beyond. With any system, scaling up from the lab to the field presents challenges. Currently, there is no field data available for model calibration or validation. However, the theoretical particle models being developed herein are the best tools available to guide progress toward field deployment. To help bridge this gap and verify model predictions, larger-scale lab experiments are being proposed. To advance simulation capabilities, six particle transport models are being integrated into the reactive transport simulator PFLOTRAN. These include colloid filtration theory (CFT) and five additional particle transport models (M1-M5). Each model, from M1 to M5, progressively incorporates additional particle transport and retention processes. Ultimately, the simplest model capable of accurately describing 1-D column data will be selected and parameterized. During FY24, the CFT and M1 model have been fully implemented within PFLOTRAN. Using an existing 1 D column experiment, the two currently implemented particle transport models (CFT and M1), and associated parameters, were fit to this experiment. While simpler model formulations are helpful for estimations, these formulations could not fully describe particle transport and retention behavior in the previous 1-D column experiment. Thus, additional complexities will need to be considered, which will be accounted for in the M2-M5 model formulations. Additionally, because a viscous, shear thinning fluid was required to keep particles in suspension, considerations for flow will also need to also be accounted for. Therefore, a new immiscible two-phase flow mode is currently being implemented in PFLOTRAN. With some modifications, this new flow module could also support simulation of non-Newtonian liquid amendments, foams, and emulsions. We also estimated the expected ROI of solid amendments using 1-D simulations. The average predicted ROI was approximately 15 ft for micron-sized zero valent iron (mZVI) suspended in xanthan gum (XG). Using the pre screening tool and ROI estimates, additional amendment-delivery laboratory characterization and experiments are proposed. The results from additional experiments can be used to validate and parametrize particulate transport model formulations, which will ultimately provide predictive capabilities for field amendment-delivery systems.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Critical Clearing Time for Droop-Controlled Grid-Forming Inverters With Circular Current Limiting and Virtual Impedance Current Limiting

The critical clearing time is a crucial parameter for assessing the large-signal stability of grid-forming (GFM) inverters. Variations in current limiters can significantly impact an inverter's transient behavior during large disturbances. Here, this letter introduces a method based on the power-angle relationship to evaluate the fault recovery of a droop-controlled GFM inverter under a three-phase bolted fault scenario, considering both the circular current limiter and the virtual impedance method. High-fidelity simulations validate the accuracy of the proposed formulas in estimating the critical clearing time of the GFM inverter.

42 ENGINEERING↗
DOE OSTI2024

EQSIM: Exascale Predictions of Earthquake Effects on Critical Infrastructure

The great “San Francisco” earthquake of 1906 is one of the most recognized, and sobering, demonstrations of the havoc that can be caused by the sudden and violent movement of Earth’s tectonic plates. The estimated 7.9-magnitude quake and subsequent fires decimated the major metropolis and surrounding areas: buildings turned to ruins, hundreds of thousands of people left homeless, and a death toll exceeding 3,000. Today, as evidenced by the catastrophic 7.8-magnitude earthquake that struck Turkey in February 2023, these events still present a significant danger to life and economic security. To mitigate the potential devastation of future earthquakes and better prepare for these inevitable events, researchers are turning to high-performance computers to simulate the underlying geophysical processes and accurately quantify associated risks to critical infrastructure.

42 ENGINEERING↗
DOE OSTI2022

Residual Stress Modeling and Advanced Diffraction Measurements of 347H Steel Weldments

The 347H stainless steel is a primary high-temperature material for many energy and power generation industries. Stress relief cracking (SRC) has been a particular concern in welding of this material. The residual stress induced by welding and its evolution during post-welding heat treatment (PWHT) and subsequent operating and service conditions is one of the primary factors contributing to SRC. The lifetime of welded structure components is also controlled by the precipitation kinetics that accompanies PWHT, stress relaxation process, and long-term aging and complex synergistic factors. Various theories have been proposed in the past to explain SRC. However, a widely accepted approach to predicting the entire damage evolution and the resulting performance reduction is still lacking. This study is to demonstrate a reliable solution for two critical issues that affect the predictions. First, the residual stress distribution obtained from both simulation and neutron diffraction is compared, which increases the accuracy of mechanical analysis simulation model and therefore builds a solid basis for the lifetime prediction model. Second, the generation, evolution, and annihilation of precipitates are monitored by the synchrotron diffraction experiment. Preliminary results demonstrate the critical importance of precipitation kinetics on the residual stress distribution/redistribution during heat treatments.

Yang, Yi↗
NASA NTRS2005

Investigation of the Carbon Arc Source as an AM0 Solar Simulator for Use in Characterizing Multi-Junction Solar Cells

The operation of multi-junction solar cells used for production of space power is critically dependent on the spectral irradiance of the illuminating light source. Unlike single-junction cells where the spectral irradiance of the simulator and computational techniques may be used to optimized cell designs, optimization of multi-junction solar cell designs requires a solar simulator with a spectral irradiance that closely matches AM0.

Xu, Jianzeng↗
NASA NTRSDate not supplied

Flight Dynamics Modeling for a Rapid Conceptual Development Environment

This paper presents the integration of flight dynamics modeling (FDM) within an aircraft conceptual development framework. Within this framework, FDM and simulation are used for the preliminary analysis of the takeoff, landing, and critical loss of thrust performance of early-stage aircraft concepts. The aerodynamic and propulsive models are built using the XML-based DAVE-ML model exchange. Automated integration of these models is achieved using the Simulink® model generation available through the DAVEtools Java package in addition to Simulink®-based flight simulation blocks. Results demonstrating the established methods are presented for the analysis of a twin-engine light transport example aircraft. This work establishes methods that will be used to perform trade studies of aircraft designed for Regional Air Mobility operations. These aircraft will require short takeoff and landing distances in addition to steep climb and descent gradients, each of which can be estimated and applied as design space constraints using the current methods.

Flight Dynamics↗
DOE OSTI2021

Density Functional Tight-Binding Simulations Reveal the Presence of Surface Defects on the Quartz (101)–Water Interface

Understanding the structure and reactivity of quartz–water interfaces is critical for numerous applications in the geological, environmental, and biological sciences. However, disagreements on the atomic-level structure of the interfaces between experiments and simulations are hampering our ability to predict the surface reactivity. Here, we used density functional tight-binding (DFTB)-based molecular dynamics to simulate a series of quartz (101) surfaces having different types and densities of surface defects in water and compared them with the structures determined by X-ray reflectivity measurements. The DFTB simulations are able to reproduce previous classical and quantum mechanical predictions of the pristine quartz (101)–water interface that disagree with experimental observations. To remedy this situation, a set of defective quartz surfaces having various surface silicon (Si) vacancies were built as indicated by recent experimental studies. We found that the rotation of surface [SiO 4 ] tetrahedra near Si vacancies can lead to outward displacements of Si atoms similar to those observed in the experiments. The presence of additional surface Si vacancies caused inward relaxations of terminal oxygens through the formation of hydrogen bonds. Finally, the overall results indicate that the quartz (101)–water interface may include a mixture of geminal (≡Si–(OH) 2 )- and vicinal (≡Si–OH)-type silanol groups together with the presence of surface Si vacancies.

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

Critical density of topological defects upon a continuous phase transition

Using extensive Monte Carlo simulations, we test the hypothesis that the density of corresponding topological defects has a universal value at the temperature of a continuous phase transition. We consider several simple two-dimensional models where domain walls, vortices, so-called Z{sub 2} vortices or their combinations are presented. These topological defects are relevant correspondingly to an Ising second-order phase transition, a Berezinskii–Kosterlitz–Thouless transition and an explicit crossover. We compare results for square and triangular lattices as well as for the complicated situation when two types of defects are presented and two transitions occur separated in temperature. All considered cases demonstrate consentient results confirming the hypothesis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1998

Shock Wave Boundary Layer Interactions in High Mach Number Flows: A Critical Survey of Current Numerical Prediction Capabilities

The report assesses the capability for numerical simulation of 2-D and 3-D shock wave laminar and turbulent boundary layer interactions Three fundamental configurations are considered: single fin, double fin, and hollow cylinder flare. Thirteen separate cases were examined by a distinguished international group of researchers using the Reynolds-averaged Navier-Stokes (RANS) equations with a wide range of turbulence models from zero equation to full Reynolds Stress Equation formulations. The report presents an extensive comparison of computations and experimental data, summarizes the results, and makes recommendations for future research.

Doyle D. Knight↗
NASA NTRSDate not supplied

EMU Ventilation Loop Simulation and Assessment of Contamination of the EMU Sublimator Hydrophilic Coating

Controlling moisture is a critical function performed in the space suit during extra vehicular activity (EVA) missions. Currently, this function is carried out by a sublimator, which is coated with a hydrophilic material. Tests have shown that siloxane compounds can off gas from the helmet absorption pad extension (HAP-E) and helmet absorption band (HAB) when they are at elevated temperatures and initial testing indicated that siloxanes could damage the hydrophilic coating. Therefore, more representative testing was needed. Important test considerations were source and destination kinetics, system geometries, scrubbing efficacy by the activated charcoal from the Metox canister, and other benefits that derive from active condensation occurring on the hydrophilic coating during exposure. A ventilation loop that was originally constructed to test hardware for control of CO2 and trace contaminants under EVA conditions was modified to conduct the siloxane off gassing tests. A vacuum compatible acrylic chamber was installed to simulate the suit volume. A steel beaker located in the chamber contained the HAP-E and HAB materials and was heated with two band heaters. The flow exiting the chamber passed through a LiOH cannister, which contained carbon beds to simulate their function in the Metox. The flow then was directed into a smaller volume that contained a coupon, which was coated with the hydrophilic material and was maintained at 10°C (50°F). Moisture was injected upstream of the suit volume simulator to maintain the dew point between 12.7°C (55°F) and 18.3°C (65°F). The pressure was controlled to between 3.8 and 4.4 psia, and the flow was held at 6 acfm. The system was run for an equivalent of 25 ten-hour EVAs, where the HAP-E/HAB materials were changed out after each EVA. The effect of the potential off gassing on the coupon was measured after each EVA period.

PLSS↗