Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation decomposition”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Investigating the role of edge neutrals in exciting tearing mode activity and achieving flat temperature profiles in LTX-β

Abstract We present observations, numerical simulations, and analysis from experiments in the Lithium Tokamak Experiment-Beta (LTX- β ) in which the electron temperature profile ( T e ( r )) shifts from flat to peaked and a tearing mode is also destabilized when the average density ( n e ave ) exceeds ∼10 19 m −3 . Flat T e ( r ) is obtained routinely in LTX- β , with a lithium coated, low-recycling first wall, once the external fueling is stopped and density decays [Boyle et al 2023 Nucl. Fusion 63 056020]. In the present experiment, flat T e profiles can be sustained while maintaining constant n e ave below a line averaged density threshold ( n e ave th ) of ∼10 19 m −3 . Above n e ave th , T e ( r ) shifts from flat to peaked and a tearing mode is destabilized. Due to low recycling, the achieved n e ave can be controlled precisely by external fueling and hence, a certain threshold of the edge neutral inventory from the external fueling is experimentally manifested through n e ave th . The goal of the present work is to investigate the role of edge neutrals in determining T e ( r ) and MHD stability in the unique low-recycling regime of LTX- β . Our hypothesis is that the peaking of T e ( r ) beyond n e ave th is due ultimately to the edge cooling by the cold neutrals beyond a critical fueling flux. At lower fueling flux, flat T e ( r ) results in broader pressure profile and lower resistivity, which in turn stabilizes the tearing mode. This hypothesis is supported by edge neutral density estimation by DEGAS 2 code. Mode analysis by singular value decomposition confirms the tearing mode structure to be m / n = 2/1 ( m and n being the poloidal and toroidal mode numbers). Linear tearing stability analysis with M3D-C1 predicts that plasmas with n e ave > 10 19 are highly susceptible to a n = 1 tearing mode. ORBIT simulations, however, confirmed that the tearing modes do not contribute to the loss of fast ions from neutral beam injection. This study shows for the first time that the neutral inventory at the edge could be one of the deciding factors for the achievability of the unique operation regime of flat T e ( r ) and the excitation of tearing activity that could be disruptive for the plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Root exudates induced coupled carbon and phosphorus cycling in a soil with low phosphorus availability

The amount and type of root exudates can influence P availability in the rhizosphere directly by desorption or dissolution of soil minerals, or indirectly by decomposition of soil organic matter (SOM). Here this study aimed to determine the mechanisms by which specific root exudates influence the distribution and availability of P in soils with low P availability. Water, glucose, alanine, and oxalate were delivered through a simulated root into soils for 15 days. Zymography and planar optodes were used to image potential phosphatase activity, and O2 and pH distribution, respectively. Soils were analyzed for resin extractable inorganic P (Pi), dissolved organic C (DOC), water soluble Fe, and Al, and microbial community structure. Characterization of SOM and P were conducted using ultra-high resolution mass spectrometry and 31P solution nuclear magnetic resonance (NMR), respectively. The addition of oxalate resulted in the greatest resin extractable Pi, DOC, and water-soluble Fe, and Al compared to the other exudates suggesting destabilization of mineral associated organic matter (MAOM) and release of organic P (Po). Both 31P solution NMR and ultra-high resolution mass spectrometry analysis provided evidence of mineralization of Po released from the destabilization of MAOM. The study demonstrates the important role microbial and plant-derived metal chelating ligands play in destabilizing MAOM, releasing SOM and importantly Po, that when mineralized may contribute to increasing Pi availability in soils with low P availability.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1989

Direct simulation of compressible turbulence

Several direct simulations of 3-D homogeneous, compressible turbulence are presented with emphasis on the differences with incompressible turbulent simulations. A fully spectral collocation algorithm, periodic in all directions coupled with a 3rd order Runge-Kutta time discretization scheme is sufficient to produce well-resolved flows at Taylor Reynolds numbers below 40 on grids of 128x128x128. A Helmholtz decomposition of velocity is useful to differentiate between the purely compressible effects and those effects solely due to vorticity production. In the context of homogeneous flows, this decomposition in unique. Time-dependent energy and dissipation spectra of the compressible and solenoidal velocity components indicate the presence of localized small scale structures. These structures are strongly a function of the initial conditions. Researchers concentrate on a regime characterized by very small fluctuating Mach numbers Ma (on the order of 0.03) and density and temperature fluctuations much greater than sq Ma. This leads to a state in which more than 70 percent of the kinetic energy is contained in the so-called compressible component of the velocity. Furthermore, these conditions lead to the formation of curved weak shocks (or shocklets) which travel at approximately the sound speed across the physical domain. Various terms in the vorticity and divergence of velocity production equations are plotted versus time to gain some understanding of how small scales are actually formed. Possible links with Burger turbulence are examined. To visualize better the dynamics of the flow, new graphic visualization techniques have been developed. The 3-D structure of the shocks are visualized with the help of volume rendering algorithms developed in-house. A combination of stereographic projection and animation greatly increase the number of visual cues necessary to properly interpret the complex flow.

Zang, T. A.↗
NASA NTRSDate not supplied

TPSAS-NF1676L-12725-DND

Uncertainty analysis and robust design - increase confidence and consistency in aerospace vehicle safety predictions by developing improved methods for quantifying and managing uncertainty. Quantifying: uncertainty modeling (model uncertainty based on experimental data, simulations and/or expert opinion) and uncertainty propagation (given uncertainty models of a system’s inputs, how to propagate them through system models, to efficiently evaluate the corresponding system’s outputs?). Managing: robust design (generate designs that robustly accommodate uncertainty) and uncertainty decomposition (identify uncertainties that contribute the most to performance degradation and determine the parameters that should (not) be modeled as uncertain).

Sean P Kenny↗
DOE OSTI2022

The effect of elastic anisotropy on the symmetry selection of irradiation-induced void superlattices in cubic metals

We report self-organized microstructures and patterns have been widely observed in non-equilibrium physical systems. In particular, irradiation in metals creates far-from-equilibrium environments, in which the competing dynamics of defect production and annihilation can lead to unique self-organized superlattice structures, e.g., void and gas bubble superlattices. From a physical point of view, the superlattice structures are dictated by the intrinsic symmetry breaking in the metals, i.e., anisotropy caused by the breaking of continuous rotational symmetry. In the literature, two distinctive anisotropies, elastic anisotropy and diffusion anisotropy of interstitials, have been proposed to be the origins of superlattice formation. However, it is still unclear which anisotropy dominates the symmetry selection of superlattice structures. In this paper, we study elastic anisotropy and its effect on the symmetry of void superlattices. By using theoretical analyses and phase field simulations, we show that elastic anisotropy in cubic metals can lead to either face-centered cubic or simple cubic superlattices depending on the Zener anisotropy ratio. The superlattices formed under this elastic anisotropy mechanism must form under the influence of spinodal decomposition, as the mechanism requires perturbations in the vacancy concentration field to develop into spatially-static concentration waves. We compare to existing work on symmetry selection in superlattices via diffusion anisotropy and to experimental observations, and we suggest that concentration wave development under the influence of elastic anisotropy is not the mechanism for symmetry selection during the formation of irradiation-induced void superlattices, but that diffusion anisotropy could be the dominant mechanism.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Unraveling local structures of Salt-in-Water and Water-in-Salt electrolytes via ab initio molecular dynamics

Water-in-salt electrolytes (WiSE) are attractive for electrochemical energy storage applications owing to their wide electrochemical stability windows and inherent safety. The high concentration of salts is widely known to suppress the decomposition of water, which otherwise limits the stability of conventional aqueous salt-in-water electrolytes (SiWE). Nevertheless, the microstructural features of WiSE that lead to this enhanced stability have not yet been fully elucidated. In this work, ab initio molecular dynamics simulations were performed to study the energetic, structural, and spectroscopic properties of LiTFSI SiWE and WiSE solutions. A detailed mapping of water-water and water-anion hydrogen bonding and cation–anion electrostatic interactions is reported for both electrolytes. Structural features are presented in terms of both radial and spatial distribution functions. Analysis of IR and vibrational power spectra reveal key differences in the intermolecular interactions in SiWE and WiSE that arise from modified solvation shell structures. Finally, the results obtained herein reveal the most important structural and spectroscopic differences between the electrolytes in normal and superconcentrated concentrations.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Flow field Reconstruction for Inhomogeneous Turbulence using Data and Physics Driven Models

A methodology combining Large Eddy Simulation (LES) trained data and a physics driven wave packet model to obtain a reduced order reconstruction for broadband, three-dimensional, temporally stationary but spatially inhomogeneous, incompressible turbulence. Wake turbulence generated by an axisymmetric dragging disk with a turbulent co-flow serves as the benchmark test case. We begin by studying the proper-orthogonal decomposition of the turbulent fluctuations taken from a high-resolution LES to first identify whether the fields demonstrate a low-rank character. It is argued that the presence of the turbulent co-flow results in a largely broadband character lacking any tonal properties. This is especially true for Strouhal numbers greater than 1 and only a small fraction of energy is contained in the leading order Kelvin-Helmholtz modes. As such reconstructions and reduced order modeling purely relying on data from LES does not appear to be a lucrative solution - contrary to problems with strongly tonal character. To supplement the missing energy from a low order truncated mode expansion, we utilize a physics based super-resolution (enrichment) algorithm that relies on spatio-temporally localized Gabor wave packets whose time evolution is described using a set of ordinary differential equations. The reconstructed flow has single- and two-point correlations that are consistent with the reference high resolution simulation data.

SLS↗
DOE OSTI2022

Reduced order models for thermal radiative transfer problems based on moment equations and data-driven approximations of the Eddington tensor

Here a new group of structure and asymptotic preserving reduced-order models (ROMs) for multidimensional nonlinear thermal radiative transfer (TRT) problems is presented. They are formulated by means of the nonlinear projective approach and data compression techniques. The nonlinear projection is applied to the Boltzmann transport equation (BTE) to derive a hierarchy of low-order moment equations. Approximation of the Eddington tensor that provides exact closure for the system of moment equations is found with projection-based data-driven methodologies. These include the (i) proper orthogonal decomposition (POD), (ii) dynamic mode decomposition (DMD) and (iii) a variant of the DMD. A parameterization is derived for this ROM for the temperature of radiation incoming to the problem domain (the radiation drive temperature). This parameterization is informed from results of a dimensionless study of the TRT problem. Analysis of the ROMs is performed on the classical Fleck-Cummings TRT multigroup test problem in 2D geometry with a radiation-driven Marshak wave. Numerical results are presented to demonstrate the performance of these ROMs for the simulation of evolving radiation and heat waves. Results show these models to be sufficiently accurate for practical computations with rather low-rank representations of the Eddington tensor. As the rank of the approximation is increased, the errors of solutions generated by the ROMs gradually decreases.

42 ENGINEERING↗
DOE OSTI2020

Denoising atomic resolution 4D scanning transmission electron microscopy data with tensor singular value decomposition

Tensor singular value decomposition (SVD) is a method to find a low-dimensional representation of data with meaningful structure in three or more dimensions. Here, tensor SVD has been applied to denoise atomic-resolution 4D scanning transmission electron microscopy (4D STEM) data. On data simulated from a SrTiO 3 [100] perfect crystal and a Si [110] edge dislocation, tensor SVD achieved an average peak signal-to-noise ratio (PSNR) of ~40 dB, which matches or exceeds the performance of other denoising methods, with processing times at least 100 times shorter. On experimental data from SrTiO 3 [100] and LiZnSb [11 2 ¯ 0]/GaSb [110] samples, tensor SVD denoises multiple GB 4D STEM data sets in ten minutes on a typical personal computer. Denoising with tensor SVD improves both convergent beam electron diffraction patterns and virtual-aperture annular dark field images.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Nonsolvating Fluoroaromatic Cosolvent Enabled Long-Term Cycling of High-Voltage Lithium-Ion Batteries with Organosulfur Electrolytes

Here, the structure–activity relationships of nonsolvating cosolvents for organosulfur-based electrolyte systems were revealed. The performance of nonsolvating dilutant fluorobenzene (FB) was compared to various fluorinated ether dilutants in high-voltage electrolytes containing a concentration of 1.2 M LiPF 6 dissolved in fluoroethylene carbonate (FEC), ethyl methyl sulfone (EMS), and the dilutant. In a high-voltage and high-loading LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811) full cell configuration, the organosulfur-based electrolyte containing FB dilutant enabled superior electrochemical performance compared to the electrolytes using other nonsolvating fluorinated ether formulations. Moreover, the FB-containing electrolyte exhibited the highest ionic conductivity and lowest viscosity among all organosulfur-based electrolytes containing nonsolvating dilutant. These improvements are attributed to the enhanced physical properties of electrolyte and lithium-ion mobility. Furthermore, by employing first-principles simulations, the observed suppression of side reactions at high voltage is linked to FB’s lower reactivity toward singlet dioxygen, which is likely produced at the NMC interface. Overall, FB is considered an excellent diluent that does not impede cell operation by mass decomposition at the cathode.

25 ENERGY STORAGE↗
DOE OSTI2024

Dynamical Janus Interface Design for Reversible and Fast-Charging Zinc–Iodine Battery under Extreme Operating Conditions

Aqueous zinc (Zn) iodine (I 2 ) batteries have emerged as viable alternatives to conventional metal-ion batteries. However, undesirable Zn deposition and irreversible iodine conversion during cycling have impeded their progress. Here, to overcome these concerns, we report a dynamical interface design by cation chemistry that improves the reversibility of Zn deposition and four-electron iodine conversion. Due to this design, we demonstrate an excellent Zn-plating/-stripping behavior in Zn||Cu asymmetric cells over 1000 cycles with an average Coulombic efficiency (CE) of 99.95%. Moreover, the Zn||I 2 full cells achieve a high-rate capability (217.1 mA h g –1 at 40 A g –1 ; C rate of 189.5C) at room temperature and enable stable cycling with a CE of more than 99% at -50 °C at a current density of 0.05 A g –1 . In situ spectroscopic investigations and simulations reveal that introducing tetraethylammonium cations as ion sieves can dynamically modulate the electrode–electrolyte interface environment, forming the unique water-deficient and chloride ion (Cl – )-rich interface. Such Janus interface accounts for the suppression of side reactions, the prevention of ICl decomposition, and the enrichment of reactants, enhancing the reversibility of Zn-stripping/-plating and four-electron iodine chemistry. This fundamental understanding of the intrinsic interplay between the electrode–electrolyte interface and cations offers a rational standpoint for tuning the reversibility of iodine conversion.

25 ENERGY STORAGE↗
DOE OSTI2020

Permafrost thermal conditions are sensitive to shifts in snow timing

Abstract Changes in snow precipitation at high latitudes can significantly affect permafrost thermal conditions and thaw depth, potentially exposing more carbon-laden soil to microbial decomposition. A fully coupled process-based surface/subsurface thermal hydrology model with surface energy balance is used to analyze the impact of intra-annual variability in snow on permafrost thermal regime and the active layer thickness. In the four numerical scenarios considered, simulations were forced by the same meteorological data, except the snow precipitation, which was systematically altered to change timing of snowfall. The scenarios represent subtle shifts in snow timing, but the snow onset/melt days, the end of winter snowpack depth, and total annual snow precipitation are unchanged among scenarios. The simulations show that small shifts in the timing of snow accumulation can have significant effects on subsurface thermal conditions leading to active layer deepening and even talik formation when snowfall arrives earlier in the winter. The shifts in snow timing have a stronger impact on wetter regions, especially soil underneath small ponds, as compared to drained regions. This study highlights the importance of understanding potential changes in winter precipitation patterns for reliable projections of active-layer thickness in a changing Arctic climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Electronic Structure Theories of Singlet Fission and Multiple Exciton Generation (Final Report)

This DOE grant supported the development of electronic structure theory, algorithms, and computer programs for simulating singlet fission and multiple excision generation processes from Sept. 2013 to Aug. 2018 (which was terminated due to the PI’s departure from Northwestern University to industry). One of the main achievements is the development of the so-called active space decomposition (ASD) method hat describes exciton dynamics in organic semiconductors and spin dynamics in strongly correlated materials. The resulting software was interfaced with the Brilliantly Advanced General Electronic-structure Library (BAGEL) package and widely disseminated to the research community.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1990

Statistical analysis of near-wall structures in turbulent channel flow

The two-point correlation tensor is used to analyze near-wall structures in turbulent channel flow. Simulation results have been utilized to obtain the correlation tensor with sufficient spatial resolution to permit accurate differentiation for computing the vorticity correlation, and the resolution of sharp features such as the shear layer and its corresponding region of spanwise vorticity. Proper orthogonal decomposition and stochastic estimation were observed to yield similar results when appropriate conditions were employed for the estimation. Resulting structures were such that they would be detected by standard conditional sampling methods. The results of decomposing the vorticity field and the velocity field were significantly different, indicating that if a coherent structure is found that dominates the velocity fields, the curl of that velocity structure will not dominate the vorticity field.

Moser, Robert D.↗
NASA NTRS1996

Development of Multiobjective Optimization Techniques for Sonic Boom Minimization

A discrete, semi-analytical sensitivity analysis procedure has been developed for calculating aerodynamic design sensitivities. The sensitivities of the flow variables and the grid coordinates are numerically calculated using direct differentiation of the respective discretized governing equations. The sensitivity analysis techniques are adapted within a parabolized Navier Stokes equations solver. Aerodynamic design sensitivities for high speed wing-body configurations are calculated using the semi-analytical sensitivity analysis procedures. Representative results obtained compare well with those obtained using the finite difference approach and establish the computational efficiency and accuracy of the semi-analytical procedures. Multidisciplinary design optimization procedures have been developed for aerospace applications namely, gas turbine blades and high speed wing-body configurations. In complex applications, the coupled optimization problems are decomposed into sublevels using multilevel decomposition techniques. In cases with multiple objective functions, formal multiobjective formulation such as the Kreisselmeier-Steinhauser function approach and the modified global criteria approach have been used. Nonlinear programming techniques for continuous design variables and a hybrid optimization technique, based on a simulated annealing algorithm, for discrete design variables have been used for solving the optimization problems. The optimization procedure for gas turbine blades improves the aerodynamic and heat transfer characteristics of the blades. The two-dimensional, blade-to-blade aerodynamic analysis is performed using a panel code. The blade heat transfer analysis is performed using an in-house developed finite element procedure. The optimization procedure yields blade shapes with significantly improved velocity and temperature distributions. The multidisciplinary design optimization procedures for high speed wing-body configurations simultaneously improve the aerodynamic, the sonic boom and the structural characteristics of the aircraft. The flow solution is obtained using a comprehensive parabolized Navier Stokes solver. Sonic boom analysis is performed using an extrapolation procedure. The aircraft wing load carrying member is modeled as either an isotropic or a composite box beam. The isotropic box beam is analyzed using thin wall theory. The composite box beam is analyzed using a finite element procedure. The developed optimization procedures yield significant improvements in all the performance criteria and provide interesting design trade-offs. The semi-analytical sensitivity analysis techniques offer significant computational savings and allow the use of comprehensive analysis procedures within design optimization studies.

Chattopadhyay, Aditi↗
NASA NTRS1998

Design and Performance Analysis of a Massively Parallel Atmospheric General Circulation Model

In the 1990's computer manufacturers are increasingly turning to the development of parallel processor machines to meet the high performance needs of their customers. Simultaneously, atmospheric scientists study weather and climate phenomena ranging from hurricanes to El Nino to global warming that require increasingly fine resolution models. Here, implementation of a parallel atmospheric general circulation model (GCM) which exploits the power of massively parallel machines is described. Using the horizontal data domain decomposition methodology, this FORTRAN 90 model is able to integrate a 0.6 deg. longitude by 0.5 deg. latitude problem at a rate of 19 Gigaflops on 512 processors of a Cray T3E 600; corresponding to 280 seconds of wall-clock time per simulated model day. At this resolution, the model has 64 times as many degrees of freedom and performs 400 times as many floating point operations per simulated day as the model it replaces.

Schaffer, Daniel S.↗
DOE OSTI2023

The Intrinsic 150‐Day Periodicity of the Southern Hemisphere Extratropical Large‐Scale Atmospheric Circulation

Abstract The variability of the Southern Hemisphere (SH) extratropical large‐scale circulation is dominated by the Southern Annular Mode (SAM), whose timescale is extensively used as a key metric in evaluating state‐of‐the‐art climate models. Past observational and theoretical studies suggest that the SAM lacks any internally generated (intrinsic) periodicity. Here, we show, using observations and a climate model hierarchy, that the SAM has an intrinsic 150‐day periodicity. This periodicity is robustly detectable in the power spectra and principal oscillation patterns (aka dynamical mode decomposition) of the zonal‐mean circulation, and in hemispheric‐scale precipitation and ocean surface wind stress. The 150‐day period is consistent with the predictions of a new reduced‐order model for the SAM, which suggests that this periodicity is associated with a complex interaction of turbulent eddies and zonal wind anomalies, as the latter propagate from low to high latitudes. These findings present a rare example of periodic oscillations arising from the internal dynamics of the extratropical turbulent circulations. Based on these findings, we further propose a new metric for evaluating climate models, and show that some of the previously reported shortcomings and improvements in simulating SAM's variability connect to the models' ability in reproducing this periodicity. We argue that this periodicity should be considered in evaluating climate models and understanding the past, current, and projected Southern Hemisphere climate variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Evaluating the Persistence of Soil Carbon After a 30-Year Nitrogen Fertilization Experiment at the Fernow Experimental Forest, WV

We sampled soils in a previously nitrogen (N)-fertilized watershed (Watershed 3), which received 35 kg N per hectare annually for 30 years, and a reference watershed (Watershed 7), 4 years after the end of a 30-year N fertilization experiment at the Fernow Experimental Forest in West Virginia. As N deposition has reduced soil pH and plant carbon (C) investments into the rhizosphere, we compared the extent to which removing these potential limitations to microbial decomposition by increasing soil pH, adding artificial root exudates, or elevating soil temperature would increase microbial decomposition (and soil C losses) in soils that have and have not received excess N inputs. As such, we incubated soils with and without the addition of artificial root exudates and dolomitic lime (2 amendment treatments and 1 untreated control), at three different temperatures (15, 20, and 25 C) to simulate warming. We measured soil respiration, enzyme activity (acid phosphatase, N -acetyl-β-glucosaminidase, beta-glucosidase, phenol oxidase, and peroxidase), and microbial biomass carbon in a 15-week microcosm experiment. All of the datasheets include the metadata for all of the treatment combinations from this incubation experiment (ID, incubation temperature, treatment, watershed).The "Soil_Chemistry_Data.csv" data includes (1) the ratio of dry soil weight (g) by wet (or fresh) soil weight (g) after 72 hours at 65 C, (2) the wet (or fresh) weight of the soil samples used to determine soil pH, and (3) soil pH, which was measured in CaCl2. Detailed methods for the measurements of soil pH are described in the datasheet itself.The "Respiration_Data.csv" data includes (1) respiration day for the measurement, and (2) cumulative CO2 respiration (mmol) measured in the microcosms on specific days.The "Microbial_Biomass_Data.csv" data includes (1) mmol microbial biomass carbon per g (dry weight) of soil in each microcosm.The "Enzyme_Data.csv" data includes (1) all extracellular soil enzyme activity (N-acetyl-glucosaminidase (NAG), acid phosphatase (AP), β-glucosidase (BG), phenol oxidase and peroxidase) that we measured in our microcosms.

54 ENVIRONMENTAL SCIENCES↗