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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Assessment of Static Delamination Propagation Capabilities in Commercial Finite Element Codes Using Benchmark Analysis

With capabilities for simulating delamination growth in composite materials becoming available, the need for benchmarking and assessing these capabilities is critical. In this study, benchmark analyses were performed to assess the delamination propagation simulation capabilities of the VCCT implementations in Marc TM and MD NastranTM. Benchmark delamination growth results for Double Cantilever Beam, Single Leg Bending and End Notched Flexure specimens were generated using a numerical approach. This numerical approach was developed previously, and involves comparing results from a series of analyses at different delamination lengths to a single analysis with automatic crack propagation. Specimens were analyzed with three-dimensional and two-dimensional models, and compared with previous analyses using Abaqus . The results demonstrated that the VCCT implementation in Marc TM and MD Nastran(TradeMark) was capable of accurately replicating the benchmark delamination growth results and that the use of the numerical benchmarks offers advantages over benchmarking using experimental and analytical results.

Orifici, Adrian C.↗
DOE OSTI2021

Designing high-performance superconductors with nanoparticle inclusions: Comparisons to strong pinning theory

One of the most promising routes for achieving high critical currents in superconductors is to incorporate dispersed, non-superconducting nanoparticles to control the dissipative motion of vortices. However, these inclusions reduce the overall superconducting volume and can strain the interlaying superconducting matrix, which can detrimentally reduce T c . Consequently, an optimal balance must be achieved between the nanoparticle density n p and size d. Determining this balance requires garnering a better understanding of vortex–nanoparticle interactions, described by strong pinning theory. Here, we map the dependence of the critical current on nanoparticle size and density in (Y 0.77 , Gd 0.23 )Ba 2 Cu 3 O 7-δ films in magnetic fields of up to 35 T and compare the trends to recent results from time-dependent Ginzburg–Landau simulations. We identify consistency between the field-dependent critical current J c (B) and expectations from strong pinning theory. Specifically, we find that J c ∝ B -α , where α decreases from 0.66 to 0.2 with increasing density of nanoparticles and increases roughly linearly with nanoparticle size d/ξ (normalized to the coherence length). At high fields, the critical current decays faster (~B -1 ), suggesting that each nanoparticle has captured a vortex. When nanoparticles capture more than one vortex, a small, high-field peak is expected in J c (B). Due to a spread in defect sizes, this novel peak effect remains unresolved here. Finally, we reveal that the dependence of the vortex creep rate S on nanoparticle size and density roughly mirrors that of α, and we compare our results to low-T nonlinearities in S(T) that are predicted by strong pinning theory.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Predicting hot electron generation in inertial confinement fusion with particle-in-cell simulations

A series of two-dimensional particle-in-cell simulations with speckled laser drivers was carried out to study hot electron generation in direct-drive inertial confinement fusion on OMEGA. Scaling laws were obtained for hot electron fraction and temperature as functions of laser/plasma conditions in the quarter-critical region. Using these scalings and conditions from hydro simulations, the temporal history of hot electron generation can be predicted. The scalings can be further improved to predict hard x-rays for a collection of OMEGA warm target implosions within experimental error bars. Furthermore, these scalings can be readily implemented into inertial confinement fusion design codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Effect of Illumination Area on the Ultrafast Temporal Response of MSM GaN Photodiodes

We investigate the influence of spatial illumination profiles on the temporal response of metal–semiconductor–metal GaN photodiodes. Using both simulation and experimental measurements, we compare the response curves under two scenarios: illumination confined to the active area between electrodes and extended illumination beyond the contacts. The results show that limiting the beam illumination to the active region significantly sharpens the response, reducing the long decay tails associated with slow carrier drift from peripheral regions. As a result, the experimental data closely match simulation predictions, confirming that illumination geometry plays a critical role in optimizing photodiode performance for ultrafast detection applications.

Carrier drift time↗
NASA NTRS2020

Modelling climate change impacts on maize yields under low nitrogen input conditions in sub‐Saharan Africa

Smallholder farmers in sub‐Saharan Africa (SSA) currently grow rainfed maize with limited inputs including fertilizer. Climate change may exacerbate current production constraints. Crop models can help quantify the potential impact of climate change on maize yields, but a comprehensive multimodel assessment of simulation accuracy and uncertainty in these low‐input systems is currently lacking. We evaluated the impact of varying [CO2], temperature and rainfall conditions on maize yield, for different nitrogen (N) inputs (0, 80, 160 kg N/ha) for five environments in SSA, including cool subhumid Ethiopia, cool semi‐arid Rwanda, hot subhumid Ghana and hot semi‐arid Mali and Benin using an ensemble of 25 maize models. Models were calibrated with measured grain yield, plant biomass, plant N, leaf area index, harvest index and in‐season soil water content from 2‐year experiments in each country to assess their ability to simulate observed yield. Simulated responses to climate change factors were explored and compared between models. Calibrated models reproduced measured grain yield variations well with average relative root mean square error of 26%, although uncertainty in model prediction was substantial (CV = 28%). Model ensembles gave greater accuracy than any model taken at random. Nitrogen fertilization controlled the response to variations in [CO2], temperature and rainfall. Without N fertilizer input, maize (i) benefited less from an increase in atmospheric [CO2], (ii) was less affected by higher temperature or decreasing rainfall and (iii) was more affected by increased rainfall because N leaching was more critical. The model inter-comparison revealed that simulation of daily soil N supply and N leaching plays a crucial role in simulating climate change impacts for low-input systems. Climate change and N input interactions have strong implications for the design of robust adaptation practices across SSA, because the impact of climate change will be modified if farmers intensify maize production with more mineral fertilizer.

Crop simulation model↗
NASA NTRS2002

Land Surface Data Assimilation and the Gulf Coast Sea-Breeze

A technique has been developed for assimilating GOES-derived skin temperature tendencies and insolation into the surface energy budget equation of a mesoscale model so that the simulated rate of temperature change closely agrees with the satellite observations. A critical assumption of the technique is that the availability of moisture (either from the soil or vegetation) is the least known term in the model's surface energy budget. Therefore, the simulated latent heat flux, which is a function of surface moisture availability, is adjusted based upon differences between the modeled and satellite observed skin temperature tendencies. An advantage of this technique is that satellite temperature tendencies are assimilated in an energetically consistent manner that avoids energy imbalances and surface stability problems that arise from direct assimilation of surface shelter temperatures. The fact that the rate of change of the satellite skin temperature is used rather than the absolute temperature means that sensor calibration is not as critical. The focus of this paper is to examine how the satellite assimilation technique impacts simulations of near-surface meteorology on the 4- to 48-hour time scale. The technique is implemented within the PSU/NCAR MM5 V34 and applied on a 36-km CONUS domain and a 12-km nest centered over Oklahoma. Two Friday periods during 2000 were selected for study, one during the warm season (May-June) and the other in the cold season (Nov-Dec). Bulk verification statistics (BIAS and RMSE) of surface air temperature and dewpoint show that the assimilation technique can improve numerical simulations on both regional and continental scales. Comparison of the simulated surface energy and water fluxes with observations from Energy Balance Bowen Ratio measurements taken at CART/ARM sites in the Southern Great Plains showed the assimilation technique produces a realistic argument of the Bowen ratio in both seasons. Additional runs with the Oregon State University land surface scheme as implemented within MM5 will be made and comparisons with assimilation results will be presented at the meeting.

Lapenta, William M.↗
NASA NTRS2001

Assimilation of GOES Land Surface Data Within a Rapid Update Cycle Format: Impact on MM5 Warm Season QPF

A technique has been developed for assimilating GOES-derived skin temperature tendencies and insolation into the surface energy budget equation of a mesoscale model so that the simulated rate of temperature change closely agrees with the satellite observations. A critical assumption of the technique is that the availability of moisture (either from the soil or vegetation) is the least known term in the model's surface energy budget. Therefore, the simulated latent heat flux, which is a function of surface moisture availability, is adjusted based upon differences between the modeled and satellite-observed skin temperature tendencies. An advantage of this technique is that satellite temperature tendencies are assimilated in an energetically consistent manner that avoids energy imbalances and surface stability problems that arise from direct assimilation of surface shelter temperatures. The fact that the rate of change of the satellite skin temperature is used rather than the absolute temperature means that sensor calibration is not as critical. The focus of this paper is to examine how the satellite assimilation technique impacts simulations of near-surface meteorology on the 0-to 12-hour time scale when implemented within a local rapid update cycle (LRUC) format. The PSU/NCAR MM5 V34 is used and configured with a 36-km CONUS domain and a 12-km nest centered over the southeastern US. The LRUC format consists of a sequence of 12-hour forecasts initialized every hour between 12 and 18 UTC seven days a week. GOES skin temperature tendencies and solar insolation are assimilated in a 1-hour period prior to the start of each twelve-hour forecast. A unique aspect of the LRUC is the satellite assimilation and the continuous recycling of the adjusted moisture availability field from one forecast cycle to the next. Preliminary results for a seven-day trial period indicate that hourly LST tendencies assimilated in a 1 hour LRUC showed improved simulated air and dewpoint temperatures for all cycles on each day. The LRUC will be used during the 2001 summer months to identify the impact of the assimilation on warm season QPF Results will be presented at the meeting.

Lapenta, William M.↗
NASA NTRS2011

Spectral Bio-indicator Simulations for Tracking Photosynthetic Activities in a Corn Field

Accurate assessment of vegetation canopy optical properties plays a critical role in monitoring natural and managed ecosystems under environmental changes. In this context, radiative transfer (RT) models simulating vegetation canopy reflectance have been demonstrated to be a powerful tool for understanding and estimating spectral bio-indicators. In this study, two narrow band spectroradiometers were utilized to acquire observations over corn canopies for two summers. These in situ spectral data were then used to validate a two-layer Markov chain-based canopy reflectance model for simulating the Photochemical Reflectance Index (PRI), which has been widely used in recent vegetation photosynthetic light use efficiency (LUE) studies. The in situ PRI derived from narrow band hyperspectral reflectance exhibited clear responses to: 1) viewing geometry which affects the asset of light environment; and 2) seasonal variation corresponding to the growth stage. The RT model (ACRM) successfully simulated the responses to the variable viewing geometry. The best simulations were obtained when the model was set to run in the two layer mode using the sunlit leaves as the upper layer and shaded leaves as the lower layer. Simulated PRI values yielded much better correlations to in situ observations when the cornfield was dominated by green foliage during the early growth, vegetative and reproductive stages (r = 0.78 to 0.86) than in the later senescent stage (r = 0.65). Further sensitivity analyses were conducted to show the important influences of leaf area index (LAI) and the sunlit/shaded ratio on PRI observations.

Cheng, Yen-Ben↗
NASA NTRS1989

Model-free simulations of turbulent reactive flows

The current computational methods for solving transport equations of turbulent reacting single-phase flows are critically reviewed, with primary attention given to those methods that lead to model-free simulations. In particular, consideration is given to direct numerical simulations using spectral (Galerkin) and pseudospectral (collocation) methods, spectral element methods, and Lagrangian methods. The discussion also covers large eddy simulations and turbulence modeling.

Givi, Peyman↗
DOE OSTI2024

Aboveground Biomass Estimation Using NISAR Simulated ALOS-2 Time Series Data

Aboveground biomass (AGB) is a critical parameter to better understand the global carbon cycle and to develop sustainable forest management. However, a large uncertainty prevails. L-band SAR data have demonstrated strong potential to accurately retrieve AGB over low-biomass regions (<100 Mg ha-1). The upcoming NASA-ISRO Synthetic Aperture Radar mission will collect data at L- and S-band over earth’s landmass with a repeat period of 12 days, allowing us to have ample data for monitoring biomass and its dynamics. One of the key science requirements of the mission is to produce annual AGB maps at 1-ha resolution with RMS accuracy of 20 Mg/ha for 80 percentage of area over low-biomass regions in Calibration/Validation sites. The NISAR biomass algorithm will generate AGB maps based on the parameterization of semi-empirical model along with NISAR time-series dual pol data (HH and HV). To calibrate and validate the model for mission requirements, the mission will use reference estimates of AGB produced from ground inventory plots and airborne LiDAR data collected over selected sites distributed across different global ecoregions. This paper presents the initial results of the calibration/validation of the NISAR AGB retrieval algorithm over the Lenoir Landing (LENO), Alabama, USA site using NISAR simulated ALOS-2 time series data. Five multi-temporal dual-pol HH and HV NISAR Simulated ALOS 2 data collections were used as input to assess the performance of the model. The model AGB retrieval results shows that the NISAR model was able to achieve RMS accuracy within 20 Mg/ha.

Ramachandran, Naveen [Jet Propulsion Laboratory, C↗
NASA NTRS1993

Self-Organized Criticality and Phase Problem in Magnetoaocustic Emission Burst in Ferromagnets

Two types of ferromagnets, pure iron and steel with varying geometry and microstructural properties, were prepared for the present study of magnetoacoustic emission (MAE). The purpose was to separate the effects of structural property variations from changes caused by differences in the sample geometry. The position and shape of the leading MAE sub-burst and its variation among the samples are explained by magnetic anisotropy and the results of numerical simulations which utilized the concept of self-organized criticality (SOC). The amplitude and duration of the second sub-burst, which previously was thought to occur as a result of a complicated interaction between non-180 deg domain walls and lattice defects, can easily be explained by the results of our simulation.

Namkung, M.↗
NASA NTRS2008

Characterization and Glass Formation of JSC-1 Lunar and Martian Soil Simulants

The space exploration mission of NASA requires long duration presence of human being beyond the low earth orbit (LEO), especially on Moon and Mars. Developing a human habitat or colony on these planets would require a diverse range of materials, whose applications would range from structural foundations, (human) life support, (electric) power generation to components for scientific instrumentation. A reasonable and cost-effective approach for fabricating the materials needed for establishing a self-sufficient human outpost would be to primarily use local (in situ) resources on these planets. Since ancient times, glass and ceramics have been playing a vital role on human civilization. A long term project on studying the feasibility of developing glass and ceramic materials using Lunar and Martian soil simulants (JSC-1) as developed by Johnson Space Center has been undertaken. The first step in this on-going project requires developing a data base on results that fully characterize the simulants to be used for further investigations. The present paper reports characterization data of both JSC-1 Lunar and JSC Mars-1 simulants obtained up to this time via x-ray diffraction analysis, scanning electron microscopy, thermal analysis (DTA, TGA) and chemical analysis. The critical cooling rate for glass formation for the melts of the simulants was also measured in order to quantitatively assess the glass forming tendency of these melts. The importance of the glasses and ceramics developed using in-situ resources for constructing human habitats on Moon or Mars is discussed.

Sen, Subhayu↗
DOE OSTI2022

Capillary Pressure to Relative Permeability: Task 7.3.4 Relative Permeability Analysis Final Report

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability relationship can be determined. Relative permeability is a numerical description of the interplay between two or more fluids and the porous media they flow through. It is a critical parameter for various tools used to characterized subsurface multiphase flow systems, such as numerical simulation for oil and gas development, carbon sequestration, and groundwater contamination remediation. Therefore, it is critical to get a good statistic distribution of relative permeability in the porous media under study. Empirical relationships for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility required to match laboratory derive relative permeability relationships. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability relationships. Existing laboratory-measured relative permeability data along with mercury intrusion capillary (MICP) data coupled with interfacial tension and contact angle measurements were used to determine the efficacy of this approach to relative permeability curve creation. The relative permeability relationships determined with this method were fit to existing laboratory data to elucidate common fitting parameters that were then used to create relative permeability relationships from MICP data that does not have an associated laboratory-measured relative permeability relationship.

58 GEOSCIENCES↗
DOE OSTI2022

Accurate Fe–He machine learning potential for studying He effects in BCC-Fe

Nanostructured Ferritic Alloys (NFAs) containing oxide nanoparticles are prospective advanced structural materials in future fusion reactors. Lacking fusion neutron sources with an adequate flux, modeling is critical to understand radiation damage in the NFAs, including helium effects. Machine learning interatomic potentials (MLPs) have been demonstrated to be capable of describing atomic interactions in chemically complex systems with comparable accuracy to density functional theory (DFT) but with much less computational costs. Hence, MLPs are chosen to study He bubble accumulation in chemically complex NFAs. As a preliminary step, we develop a Fe–He potential, based on ~10,000 atomic configurations. The developed MLP accurately predicts the bulk properties of BCC-Fe compared to DFT, including the lattice constant (2.832 Å), elastic constants (c 11 = 271 GPa, c 12 = 141, c 44 = 93), and phonon frequencies (maximum error < 3.5%). The MLP predicts He binding energies in He n V bubbles and He n clusters with a mean absolute error (MAE) of only ~70 meV, which is ~3 to ~5 times smaller than the MAE of binding energies estimated using empirical potentials. The MLP also accurately predicts the migration barrier of interstitial He (64 meV). Furthermore, the MLP correctly predicts the maximum number of He atoms (critical size) at which a He n V bubble (He n cluster) can kick-out a self-interstitial atom (SIA), forming a bigger He n V 2 bubble (He n V bubble). Subsequently, we use the MLP to perform molecular dynamics simulations to investigate the effect of temperature on the critical size of He n V bubbles and He n clusters. We find that the critical size decreases with increasing temperature — indicating that smaller bubbles and clusters can kick-out an SIA at higher temperatures, thereby escalating the He bubble growth. We report the current Fe–He MLP can be further developed to include all chemical interactions in NFAs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Disentangling the Effects of Vapor Pressure Deficit and Soil Water Availability on Canopy Conductance in a Seasonal Tropical Forest During the 2015 El Niño Drought

In land models used in Earth system models, vegetation water stress is heuristically related to soil water matric potential, and root water uptake in soil layers is simply proportional to root density to meet the transpiration demand. The efficiency of roots utilizing deeper soil water under water stress cannot be well represented, but it is of critical importance for transpiration regulation and carbon flux. To improve modeling of tropical forest response to drought, a plant hydrodynamic scheme (HYDRO), adapted from a trait-based model in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), has been implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) to represent the hydraulic dynamics and root water uptake strategy. The adapted model conforms to the big-leaf photosynthesis model in ELM. Vegetation water stress in the model is a function of prognostic leaf water potential and is used to regulate stomatal conductance. Point simulations were performed at two tropical forest sites in the Amazon and Panama. Model sensitivity to the numerical solver, plant hydraulic traits, soil properties, and atmospheric forcing has been evaluated. Our results showed that, even though water stress factors calculated at our study sites using HYDRO are comparable to that simulated by the heuristic water stress model, soil moisture prediction is improved with soil water redistributions through roots. Simulated results using HYDRO is not as sensitive to plant hydraulic traits measured for different tree species as to the soil properties and atmospheric forcing. Simplified numerical solver may cause inaccuracy as to where roots uptake water. This study suggests that plant hydrodynamic representation is critical for improved simulations of soil moisture under water stresses.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

BISON microstructure-based pulverization criterion in high burnup structure

To improve the economics of commercial nuclear power production, utilities are seeking to increase the allowable burnup limit of UO$_2$ fuel. One of the main factors that contributes to the current burnup limit of 62 GWd/MTU in commercial light water reactors (LWRs) is the risk of fine fragmentation or pulverization during a loss of coolant accident (LOCA). Pulverization primarily occurs at high burnups, especially when the high burnup structure (HBS) has formed. To allow the industry to pursue increased burnup and develop mitigation strategies, it is essential to have improved capability to predict the onset of pulverization. However, the mechanism of pulverization is not well understood, and the existing predictive capabilities implemented in the BISON fuel performance code are empirical in nature. In this report, mesoscale simulations are used to improve understanding of the formation mechanism of the HBS and how it responds during a LOCA transient, and inform development of a BISON pulverization criterion. A phase-field model was used to simulate the evolution of bubble pressure as a result of HBS formation. The simulations showed that gas atoms diffuse from grain interiors to the new grain boundaries created during HBS formation, and diffuse rapidly along these grain boundaries to reach existing bubbles. This causes an increase in bubble pressure in existing bubbles, leading to bubble growth during steady-state operation. To simulate the response of HBS bubbles to a LOCA transient, a newly developed phase-field model was used; in agreement with preliminary results from FY20, bubble size did not change significantly during the duration of the transient. A phase-field fracture model was used to study fragmentation patterns in the HBS, including using input from the phase-field model as initial conditions. Phase-field fracture simulations were used to determine a pulverization criterion for BISON. A function for the critical pressure for pulverization to occur was fit to data from the phase-field fracture simulations, and this function was implemented as a material property in BISON. For comparison, an analytical criterion for pulverization was developed and implemented within the same material property.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Mechanism of Peptide Agonist Binding in CXCR4 Chemokine Receptor

Chemokine receptors are key G-protein-coupled receptors (GPCRs) that control cell migration in immune system responses, development of cardiovascular and central nervous systems, and numerous diseases. In particular, the CXCR4 chemokine receptor promotes metastasis, tumor growth and angiogenesis in cancers. CXCR4 is also used as one of the two co-receptors for T-tropic HIV-1 entry into host cells. Therefore, CXCR4 serves as an important therapeutic target for treating cancers and HIV infection. Apart from the CXCL12 endogenous peptide agonist, previous studies suggested that the first 17 amino acids of CXCL12 are sufficient to activate CXCR4. Two 17-residue peptides with positions 1–4 mutated to RSVM and ASLW functioned as super and partial agonists of CXCR4, respectively. However, the mechanism of peptide agonist binding in CXCR4 remains unclear. Here, we have investigated this mechanism through all-atom simulations using a novel Peptide Gaussian accelerated molecular dynamics (Pep-GaMD) method. The Pep-GaMD simulations have allowed us to explore representative binding conformations of each peptide and identify critical low-energy states of CXCR4 activated by the super versus partial peptide agonists. Our simulations have provided important mechanistic insights into peptide agonist binding in CXCR4, which are expected to facilitate rational design of new peptide modulators of CXCR4 and other chemokine receptors.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1997

Visualization in the Design of Modern Aircraft Aerodynamics

Modem aircraft design involves study of airflow through both windtunnel testing and computer simulation. These computer simulations result in often very large and complex sets of numbers, which contain information critical to the aircrafts performance. This talk will describe how visualization is used to understand these simulations, using a variety of techniques including low-level analysis such as simulated particles, high-level feature detection, and virtual-reality-based techniques for exploration. We will focus on the challenges of extremely large data sets, interactive performance, and information extraction. The talk will close with a vision of the future including the integration of simulation and visualization.

Bryson, Steve↗