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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Concurrent Probabilistic Simulation of High Temperature Composite Structural Response

A computational structural/material analysis and design tool which would meet industry's future demand for expedience and reduced cost is presented. This unique software 'GENOA' is dedicated to parallel and high speed analysis to perform probabilistic evaluation of high temperature composite response of aerospace systems. The development is based on detailed integration and modification of diverse fields of specialized analysis techniques and mathematical models to combine their latest innovative capabilities into a commercially viable software package. The technique is specifically designed to exploit the availability of processors to perform computationally intense probabilistic analysis assessing uncertainties in structural reliability analysis and composite micromechanics. The primary objectives which were achieved in performing the development were: (1) Utilization of the power of parallel processing and static/dynamic load balancing optimization to make the complex simulation of structure, material and processing of high temperature composite affordable; (2) Computational integration and synchronization of probabilistic mathematics, structural/material mechanics and parallel computing; (3) Implementation of an innovative multi-level domain decomposition technique to identify the inherent parallelism, and increasing convergence rates through high- and low-level processor assignment; (4) Creating the framework for Portable Paralleled architecture for the machine independent Multi Instruction Multi Data, (MIMD), Single Instruction Multi Data (SIMD), hybrid and distributed workstation type of computers; and (5) Market evaluation. The results of Phase-2 effort provides a good basis for continuation and warrants Phase-3 government, and industry partnership.

Abdi, Frank↗
NASA NTRS1998

A 2-D Interface Element for Coupled Analysis of Independently Modeled 3-D Finite Element Subdomains

Over the past few years, the development of the interface technology has provided an analysis framework for embedding detailed finite element models within finite element models which are less refined. This development has enabled the use of cascading substructure domains without the constraint of coincident nodes along substructure boundaries. The approach used for the interface element is based on an alternate variational principle often used in deriving hybrid finite elements. The resulting system of equations exhibits a high degree of sparsity but gives rise to a non-positive definite system which causes difficulties with many of the equation solvers in general-purpose finite element codes. Hence the global system of equations is generally solved using, a decomposition procedure with pivoting. The research reported to-date for the interface element includes the one-dimensional line interface element and two-dimensional surface interface element. Several large-scale simulations, including geometrically nonlinear problems, have been reported using the one-dimensional interface element technology; however, only limited applications are available for the surface interface element. In the applications reported to-date, the geometry of the interfaced domains exactly match each other even though the spatial discretization within each domain may be different. As such, the spatial modeling of each domain, the interface elements and the assembled system is still laborious. The present research is focused on developing a rapid modeling procedure based on a parametric interface representation of independently defined subdomains which are also independently discretized.

Kandil, Osama A.↗
DOE OSTI2020

Recent Progress in Nuclear Lattice Simulations

We review several recent results on lattice simulations by the Nuclear Lattice Effective Field Theory Collaboration. In the first part we discuss the implementation of nuclear forces on the lattice using chiral effective field theory. The new development we highlight is the use of non-local lattice operators to achieve a simpler spin channel decomposition, in contrast with previous studies that considered only local interactions. In the second part, we present evidence that nuclear physics is close to a quantum phase transition. This development is also linked to the study of the differences between local and non-local interactions. In the final part we further explore the link between the nuclear forces and nuclear structure. We consider the simplest possible nuclear interaction which can accurately reproduce the ground state energies of neutron matter, light nuclei, and medium-mass nuclei. We discuss what these recent developments say about the emergence of nuclear structure from nuclear forces and the road ahead for nuclear lattice simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2016

Effects of Different Heave Motion Components on Pilot Pitch Control Behavior

The study described in this paper had two objectives. The first objective was to investigate if a different weighting of heave motion components decomposed at the center of gravity, allowing for a higher fidelity of individual components, would result in pilot manual pitch control behavior and performance closer to that observed with full aircraft motion. The second objective was to investigate if decomposing the heave components at the aircraft's instantaneous center of rotation rather than at the center of gravity could result in additional improvements in heave motion fidelity. Twenty-one general aviation pilots performed a pitch attitude control task in an experiment conducted on the Vertical Motion Simulator at NASA Ames under different hexapod motion conditions. The large motion capability of the Vertical Motion Simulator also allowed for a full aircraft motion condition, which served as a baseline. The controlled dynamics were of a transport category aircraft trimmed close to the stall point. When the ratio of center of gravity pitch heave to center of gravity heave increased in the hexapod motion conditions, pilot manual control behavior and performance became increasingly more similar to what is observed with full aircraft motion. Pilot visual and motion gains significantly increased, while the visual lead time constant decreased. The pilot visual and motion time delays remained approximately constant and decreased, respectively. The neuromuscular damping and frequency both decreased, with their values more similar to what is observed with real aircraft motion when there was an equal weighting of the heave of the center of gravity and heave due to rotations about the center of gravity. In terms of open- loop performance, the disturbance and target crossover frequency increased and decreased, respectively, and their corresponding phase margins remained constant and increased, respectively. The decomposition point of the heave components only had limited effects on pilot manual control behavior and performance.

manual control↗
NASA NTRS1990

Numerical simulation of spatially-evolving instability

A computational study of the spatial stability of plane Poiseuille flow is presented. The numerical scheme employs a time-splitting method to integrate the full Navier-Stokes equations using spectral collocation/finite-difference discretization on a non-staggered mesh. The eigenvalue decomposition procedure is applied for the solution of the Poisson equations using the capacitance matrix technique. The buffer domain method is incorporated for the outflow boundary conditions. The input perturbation velocities are obtained by solving the Orr-Sommerfeld equation for the nonlinear eigenvalue problem employing the companion matrix method. Computational results are compared with the linear theory for two-dimensional disturbances.

Danabasoglu, G.↗
NASA NTRSDate not supplied

Dynamic Analysis of Reynolds Number Effects on Trailing Edge Transonic Vortex Shedding and Its Impact on Turbine Blade Aerodynamic Performance

Time-resolved, high-speed self-aligned focusing schlieren images were acquired in the NASA Glenn Research Center Transonic Turbine Blade Cascade facility to help understand the aerodynamic behavior of high-pressure, thick trailing edge turbine blades. The trailing edge thickness of 9% of axial chord tested represents simulated ceramic matrix composite fabrication constraints, which was verified previously to possess a high-loss flow regime at high inlet turbulence conditions over a narrow range of Reynolds numbers and at a fixed design exit Mach number of 0.74. Our high-speed images, which were acquired at 10 distinct Reynolds numbers, show a significant increase in energy from flow oscillations due to transonic vortex shedding at Reynolds numbers corresponding to the high loss conditions. For those conditions, strong acoustic waves turn into shock waves. Spectral Proper Orthogonal Decomposition of the high-speed images shows acoustic waves from trailing edge vortex shedding at all conditions, with increased spectral energy at the high-loss conditions and slightly increasing frequency (about 6%) as a function of Reynolds number. Analysis of potential feedback timing is performed using velocity fields from a previous LES simulation, considering different feedback mechanisms. Most noteworthy is the acoustic/shock-boundary layer interaction mechanism on the suction surface at the blade geometric throat, which likely plays an important role in realistic curved blade passages.

Trailing Edge↗
DOE OSTI2025

Unveiling the Mechanism of Mn Dissolution Through a Dynamic Cathode‐Electrolyte Interphase on LiMn2O4

Abstract Understanding the formation and evolution of the cathode‐electrolyte interphase (CEI), which forms at the interface between the cathode and electrolyte, is crucial for revealing degradation mechanisms in cathode materials, especially for developing strategies to stabilize the interphase in the strongly oxidizing conditions that evolve at high operating voltages in next‐generation Li‐ion batteries. However, The present understanding of the CEI is challenged by its complex and dynamic nature. In this work, near‐edge X‐ray absorption fine structure spectroscopy, electrochemical characterization, and reactive molecular dynamics simulations are combined to reveal a mechanism for CEI formation and evolution above model LiMn 2 O 4 (LMO) thin‐film electrodes in contact with conventional carbonate‐based electrolytes. It is found that Mn dissolution from LMO can be understood in terms of repetitive Mn 3 O 4 formation and dissolution behavior during cycling, which is closely connected to electrolyte decomposition and a key aspect of the CEI formation and growth. The behavior of the CEI in this model system offers detailed insight into the dynamic chemistry of the interphase, underscoring the important role of electrolyte composition and cathode surface structure in interphase degradation.

Ou, Wenhan↗
DOE OSTI2023

A Hybrid Domain Overlapping Method for Coupling System Thermal Hydraulics and CFD Codes

A hybrid multiscale coupling methodology based on a domain overlapping approach has been developed for coupling System Thermal Hydraulics (STH) and Computational Fluid Dynamics (CFD) codes. The method has been implemented between the modern STH code SAM and the CFD code NekRS, using the coupling tool Cardinal. The coupling aims to extend the STH code's applicability to scenarios where local momentum and energy transfers are important yet difficult for STH codes to capture, such as three-dimensional mixing. Two coupling strategies are implemented and compared: a hybrid domain overlapping method and the conventional domain decomposition method. Here, the strategies are applied to two closed-loop applications, and the present method shows superior stability behavior when compared to the domain decomposition method. Then, the present coupling method is validated against experimental data from a double T-junction experiment. The present STH/CFD coupling shows improved agreement with experimental data when compared to STH standalone simulations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Phase transformation mechanism in irradiation-induced superlattice formation

Atomic kinetic Monte Carlo simulations were used to model void superlattice formation under irradiation in molybdenum, driven by anisotropic diffusion of self-interstitial atoms. A change in the phase transformation mechanism from nucleation and growth to spinodal decomposition occurred with increasing dose rate, with both mechanisms leading to superlattice formation. Analysis of a rate-theory based analytical model showed that an observed change in the kinetics of vacancy accumulation, the appearance of a region of positive second derivative in the plot of average vacancy concentration versus time, was caused by the onset of spinodal instability. Further, the analytical model showed that for molybdenum and several other metals where void superlattice formation is commonly observed, the phase transformation likely occurs by nucleation and growth. However, nickel may offer the possibility of experimental observation of the transition between phase transformation mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Bundling measures for food systems transformation: a global, multimodel assessment

Background Current food systems leave one in ten individuals at risk of hunger while driving unsustainable environmental impacts. Inaction risks further exacerbating negative impacts on both human and planetary health. These challenges emerge from complex system interactions, requiring approaches that engage with this complexity and consider how transformation measures interact across food systems. We aimed to quantify the magnitude and uncertainty of the impacts of key food systems transformation measures both individually and in a bundle using an ensemble of global economic models. Methods In this global multimodel assessment, we applied an ensemble of ten state-of-the-art global economic models to evaluate the potential of four key measures in transforming food systems: increasing agricultural productivity, halving food loss and waste, shifting towards healthier diets, and economy-wide climate mitigation policies aligned with limiting warming to 1·5°C. The scenarios used a middle-of-the-road shared socioeconomic pathway for population and gross domestic product growth, climate impact data from Jägermeyr and colleagues, Thornton and colleagues, and Nelson and colleagues, and dietary targets based on the EAT–Lancet healthy reference diet, with model simulations conducted from 2020 to 2050. We then assessed the effect of these measures in isolation and in combination in a bundled scenario. To further understand the interactions between these measures, we conducted a decomposition analysis that distinguishes between the individual effects of a measure (effect when implemented alone), total effects (its contribution within the bundle), and interaction effects (the difference between total and individual effects). This approach aimed to show complementarities and trade-offs that emerge when multiple measures are implemented simultaneously. Findings Our analysis showed that individual measures in isolation are insufficient to achieve high-level environmental objectives and might generate unintended consequences. In contrast, bundling measures produces co-benefits: avoiding 50% of projected agricultural greenhouse gas emissions by 2050 and almost 20% of anticipated land conversion, while moderating food price increases associated with ambitious climate change mitigation policies. Our decomposition analysis further shows that measures can have varying effects across different dimensions. Although dietary shifts and climate mitigation policies are the largest drivers of environmental benefits (each contributing to a median decline of >10 percentage points in non-CO 2 emissions and 5 percentage points in agricultural land use globally), productivity improvements and reducing food loss and waste play essential roles in moderating price increases (each contributing to a median decline of >5 percentage points in average prices). Interpretation This study highlights the importance of implementing coordinated approaches to food system transformation and climate change mitigation rather than relying on isolated interventions. Comprehensive transformation requires understanding how supply-side and demand-side changes can interact with climate mitigation policies, enabling policy makers to design intervention packages that maximise benefits while minimising trade-offs across environmental, economic, and social dimensions.

Sundiang, Marina [Cornell Univ., Ithaca, NY (Unite↗
DOE OSTI2022

Decentralized modular hybrid supervisory control for the formation of unmanned helicopters

Abstract Formation control of Unmanned Aerial Vehicles (UAVs) requires them to tightly cooperate to reach and keep the formation, while avoiding collision. This paper proposes a novel decentralized hybrid supervisory control approach for the formation control of multiple UAVs. This is achieved by developing a symbolic motion planning technique to polarly partition the motion space resulting in a finite state discrete event model for the motion dynamics of each UAV. Then, a modular discrete supervisor is designed for different components of the formation mission including reaching the formation, keeping the formation, and collision avoidance. Further, for the collision avoidance mechanism, a novel top‐down decomposition‐based approach is developed to design local supervisors decentralizedly. It is formally proved that with the proposed top‐down decomposition‐based approach, the (locally) supervised UAVs, as a whole, can cooperatively satisfy the desired (global) collision avoidance specification. The proposed decentralized supervisory control algorithm is also verified through a hardware‐in‐the‐loop simulator for the formation control of unmanned helicopters.

Karimoddini, Ali↗
DOE OSTI2024

Beam correction for multi-pass arcs in FFA@CEBAF: status update

This work examines the multi-pass steering of six electron beams in an FFA arc ranging from approximately 10.5 GeV to 22 GeV. Shown here is an algorithm based on singular value decomposition (SVD) to successfully steer all six beams through the arc given precise knowledge of all beam positions at each of one hundred and one diagnostic locations with one hundred individual corrector magnets: that is successive application of SVD to different 100 × 101 response matrices—one for each beam energy. Further, a machine learning scheme is developed which only requires knowledge of the energy-averaged beam position at each location to provide equivalent steering. Extension of this scheme to other beam optics quantities as well as transverse and longitudinal coupling is explored.

Accelerator Physics↗
DOE OSTI2023

Development of Building Design Optimization Methodology: Residential Building Applications

Building design optimization is a highly complex problem, requiring long computational running processes because of the many options that exist when a building is being designed. This paper introduces an integrated approach through which to perform this optimization within an acceptable time frame. The approach includes the methods of variable selection, model simplification, and a sequential optimization process. Using singular value decomposition, a large number of design variables is reduced to a smaller subset that can be solved more quickly through the optimization algorithm. To expedite the variable selection process, a modeling approach that quickly simulates annual energy consumption was developed to replace full annual energy simulations. The developed methodology was applied to two residential buildings in the US, and the results are discussed herein. To assess the accuracy of the integrated optimization methodology, the optimized life cycle costs are compaa variables demonstrating the strongest contributions in the optimization study were identified. The proposed methodology significantly shortened the time requirements for the optimization processes of the two case studies by 74% and 84%; the optimized life cycle costs were within 0.05% and 0.06%, respectively, of the optimum point.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1994

Optical Implementation of Matching Pursuit for Image Representation

We have developed a technique for image analysis, representation, and decomposition. This technique was motivated by Stephane Mallat's matching-pursuit algorithm. We've altered and simplified the mechanics of his algorithm to enable an extremely fast implementation via optical processing. Initial computer simulations show that our algorithm is capable of decomposing and representing a 2-D image as a linear combination of basis images with both high speed and high fidelity.

matching pursuit adaptive wavelet transform image ↗
DOE OSTI2023

Preparations for quantum simulations of quantum chromodynamics in 1+1 dimensions. I. Axial gauge

Here, tools necessary for quantum simulations of 1+1 dimensional quantum chromodynamics are developed. When formulated in axial gauge and with two flavors of quarks, this system requires 12 qubits per spatial site with the gauge fields included via nonlocal interactions. Classical computations and D-wave’s quantum annealer advantage are used to determine the hadronic spectrum, enabling a decomposition of the masses and a study of quark entanglement. Color “edge” states confined within a screening length of the end of the lattice are found. IBM’s seven-qubit quantum computers, ibmq_jakarta and ibm_perth, are used to compute dynamics from the trivial vacuum in one-flavor QCD with one spatial site. More generally, the Hamiltonian and quantum circuits for time evolution of 1+1 dimensional SU(N c ) gauge theory with Nf flavors of quarks are developed, and the resource requirements for large-scale quantum simulations are estimated.

1+1 dimensions↗
DOE OSTI2025

A Parametric, Data-Driven, Non-Intrusive Reduced-Order Model Framework for Crystal Plasticity Simulations of Voids

The influence of the internal structure at micrometer length scales on the deformation of polycrystalline materials can be effectively captured using crystal plasticity finite element methods (CPFEM). However, the complexity and nonlinearity of the deformation equations CPFEM solves demand significant computational power and resources to achieve accurate predictions, limiting its broader application. To address this challenge, we have identified a reduced-order representation of the complex data in order to establish a computationally efficient reduced-order models (ROM) and drastically reduce the computational expense of CPFEM. Specifically, in this work, we developed a parametric, data-driven, and non-intrusive ROM framework for CPFEM using proper orthogonal decomposition (POD) and sparse variational Gaussian process (SVGP) regression for single-crystal microstructures under tensile loading conditions. The developed protocol enables one to compress field into a latent/low-dimensional space described by principal component analysis (PCA) via the singular value decomposition (SVD) algorithm. As a result, the high-dimensional data are reduced to a significantly smaller amount of dimensions with POD bases and POD coefficients. Furthermore, we deployed an ensemble of SVGPs—extended from the classical Gaussian process (GP) regression for scalability and handling big data—in a massively parallel manner to train and predict latent POD coefficients using known POD bases from a set of previously obtained simulations results. Lastly, using the predicted POD coefficients, we reconstructed the full-field results and showed reasonable agreement compared with the true values obtained from running CPFEM. The developed framework is validated with a set of CPFEM simulations of a single embedded void in single-crystal aluminum alloy. While the framework is broadly applicable, this work specifically focuses on single-crystal microstructures, a single load case (e.g., tensile), and a specific void geometry (spherical).

Anisotropy↗
DOE OSTI2022

Investigating coastal backwater effects and flooding in the coastal zone using a global river transport model on an unstructured mesh

Abstract. Coastal backwater effects are caused by the downstream water level increase as a result of elevated sea level, high river discharge and their compounding influence. Such effects have crucial impacts on floods in densely populated regions but have not been well represented in large-scale river models used in Earth system models (ESMs), partly due to model mesh deficiency and oversimplifications of river hydrodynamics. Using two mid-Atlantic river basins as a testbed, we perform the first attempt to simulate the backwater effects comprehensively over a coastal region using the MOSART river transport model under an ESM framework, i.e., Energy Exascale Earth System Model (E3SM) configured on a regionally refined unstructured mesh, with a focus on understanding the backwater drivers and their long-term variations. By including sea level variations at the river downstream boundary, the model performance in capturing backwaters is greatly improved. We also propose a new flood event selection scheme to facilitate the decomposition of backwater drivers into different components. Our results show that while storm surge is a key driver, the influence of extreme discharge cannot be neglected, particularly when the river drains to a narrow river-like estuary. Compound flooding, while not necessarily increasing the flood peaks, exacerbates the flood risk by extending the duration of multiple coastal and fluvial processes. Furthermore, our simulations and analysis highlight the increasing strength of backwater effects due to sea level rise and more frequent storm surge during 1990–2019. Thus, backwaters need to be properly represented in ESMs to improve the predictive understanding of coastal flooding.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Modeling diffusion and types I-V sorption of water vapor in heterogeneous systems

Here, a sorption and diffusion model is developed using a mobile-immobile decomposition of the vapor concentration. Henry’s mode represents the mobile species which diffuses through the solid, while Langmuir and pooling modes represent possible immobilization processes that cause nonlinearities in sorption capacity. The model can simulate the five classical sorption types and can incorporate additional dynamics compared to models based on thermodynamic equilibrium. Additionally, the framework allows for discontinuities in material properties encountered in multi-material systems. The resulting non-linear coupled equations are solved by employing a finite element method to discretize in space, and a backward Euler method to discretize in time. The discretized system of equations is solved via a Newton-type iteration scheme at each time-step. Four different materials are parameterized for the model, which are then used in examples to demonstrate the model’s ability to capture complex sorption processes in multi-material systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗