Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex 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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Computational modeling of microalgal biofilm growth in heterogeneous rotating algal biofilm reactors (RABRs) for wastewater treatment

Rotating algal biofilm reactors (RABRs) are innovative systems designed to cultivate microalgae biofilms efficiently. In this paper, we have developed a novel mathematical model to accurately capture the growth dynamics of algae biofilms within RABR. By considering the spatial heterogeneity of the RABR, we introduce a PDE-based model that addresses the spatial variations across the substratum, enabling a more accurate simulation of biofilm growth in RABRs. The photosynthesis process is modeled through reactive kinetics, driving the growth of the algae biofilm. To analyze the system's behavior, we employ finite difference numerical methods to solve the complex PDE model. We then conduct extensive numerical simulations to understand algae biofilm growth in the RABR environment under various operational factors and environmental conditions. One primary focus in these simulations is to investigate the impact of various harvesting strategies, harvesting frequencies, light intensity, and light exposure on the overall biomass productivity of the algae biofilm. The numerical results provide valuable insights into optimizing algae biofilm growth and designing harvesting techniques in RABR systems. Our proposed novel mathematical model provides an effective platform for the theoretical investigation and design of RABRs for wastewater treatment.

09 BIOMASS FUELS↗
DOE OSTI2024

Progressive Dynamics for Cloth and Shell Animation

We propose Progressive Dynamics, a coarse-to-fine, level-of-detail simulation method for the physics-based animation of complex frictionally contacting thin shell and cloth dynamics. Progressive Dynamics provides tight-matching consistency and progressive improvement across levels, with comparable quality and realism to high-fidelity, IPC-based shell simulations [Li et al. 2021] at finest resolutions. Together these features enable an efficient animation-design pipeline with predictive coarse-resolution previews providing rapid design iterations for a final, to-be-generated, high-resolution animation. In contrast, previously, to design such scenes with comparable dynamics would require prohibitively slow design iterations via repeated direct simulations on high-resolution meshes. We evaluate and demonstrate Progressive Dynamics's features over a wide range of challenging stress-tests, benchmarks, and animation design tasks. Here Progressive Dynamics efficiently computes consistent previews at costs comparable to coarsest-level direct simulations. Its matching progressive refinements across levels then generate rich, high-resolution animations with high-speed dynamics, impacts, and the complex detailing of the dynamic wrinkling, folding, and sliding of frictionally contacting thin shells and fabrics.

Computer Science↗
DOE OSTI2023

Simulations of Criticality Control Overpack Container Compaction at the Waste Isolation Pilot Plant

Criticality Control Overpack (CCO) containers are being considered for the disposal of defense-related nuclear waste at the Waste Isolation Pilot Plant (WIPP). At WIPP, these containers would be placed in underground disposal rooms, which will naturally close and compact the containers closer to one another over several centuries. This report details simulations to predict the final container configuration as an input to nuclear criticality assessments. Each container was discretely modeled, including the plywood and stainless steel pipe inside the 55-gallon drum, in order to capture its complex mechanical behavior. Although these high-fidelity simulations were computationally intensive, several different material models were considered in an attempt to reasonably bound the horizontal and vertical compaction percentages. When exceptionally strong materials were used for the containers, the horizontal and vertical closure respectively stabilized at 43:9 % and 93:7 %. At the other extreme, when the containers completely degraded and the clay seams between the salt layers were glued, the horizontal and vertical closure reached respective final values of 48:6 % and 100 %.

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

Physics-Guided, Physics-Informed, and Physics-Encoded Neural Networks and Operators in Scientific Computing: Fluid and Solid Mechanics

Abstract Advancements in computing power have recently made it possible to utilize machine learning and deep learning to push scientific computing forward in a range of disciplines, such as fluid mechanics, solid mechanics, materials science, etc. The incorporation of neural networks is particularly crucial in this hybridization process. Due to their intrinsic architecture, conventional neural networks cannot be successfully trained and scoped when data are sparse, which is the case in many scientific and engineering domains. Nonetheless, neural networks provide a solid foundation to respect physics-driven or knowledge-based constraints during training. Generally speaking, there are three distinct neural network frameworks to enforce the underlying physics: (i) physics-guided neural networks (PgNNs), (ii) physics-informed neural networks (PiNNs), and (iii) physics-encoded neural networks (PeNNs). These methods provide distinct advantages for accelerating the numerical modeling of complex multiscale multiphysics phenomena. In addition, the recent developments in neural operators (NOs) add another dimension to these new simulation paradigms, especially when the real-time prediction of complex multiphysics systems is required. All these models also come with their own unique drawbacks and limitations that call for further fundamental research. This study aims to present a review of the four neural network frameworks (i.e., PgNNs, PiNNs, PeNNs, and NOs) used in scientific computing research. The state-of-the-art architectures and their applications are reviewed, limitations are discussed, and future research opportunities are presented in terms of improving algorithms, considering causalities, expanding applications, and coupling scientific and deep learning solvers.

Computer Science↗
DOE OSTI2023

Probing DNA structural heterogeneity by identifying conformational subensembles of a bicovalently bound cyanine dye

DNA is a re-configurable, biological information-storage unit, and much remains to be learned about its heterogeneous structural dynamics. For example, while it is known that molecular dyes templated onto DNA exhibit increased photostability, the mechanism by which the structural dynamics of DNA affect the dye photophysics remains unknown. Here, we use femtosecond, two-dimensional electronic spectroscopy measurements of a cyanine dye, Cy5, to probe local conformations in samples of single-stranded DNA (ssDNA–Cy5), double-stranded DNA (dsDNA–Cy5), and Holliday junction DNA (HJ–DNA–Cy5). A line shape analysis of the 2D spectra reveals a strong excitation–emission correlation present in only the dsDNA–Cy5 complex, which is a signature of inhomogeneous broadening. Molecular dynamics simulations support the conclusion that this inhomogeneous broadening arises from a nearly degenerate conformer found only in the dsDNA–Cy5 complex. These insights will support future studies on DNA’s structural heterogeneity.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Probing exciton diffusion dynamics in photosynthetic supercomplexes via exciton–exciton annihilation

Photosynthesis converts solar energy into chemical energy through coordinated energy transfer between light-harvesting complexes and reaction centers (RCs). Understanding exciton motion, particularly the exciton diffusion length, is essential for optimizing energy efficiency in photosystems. In this work, we combine intensity-cycling transient absorption spectroscopy with kinetic Monte Carlo (kMC) simulation to investigate exciton motion in the C2S2 photosystem II supercomplex of spinach. Using exciton–exciton annihilation, revealed in the fifth-order response, we experimentally estimate an exciton diffusion length of 10.9 nm based on a 3D normal diffusion model, suggesting the ability of excitons to traverse the supercomplex. However, kMC simulations reveal that exciton motion is sub-diffusive because of spatial constraints and the strong RC traps. An anomalous diffusion model analysis of the experimental data yields a diffusion length of 9.7 nm, while the simulated diffusion length is 7.4 nm. The variable exciton residence time across subunits, partly influenced by their connectivity to the trap, indicates inhomogeneous annihilation probability and suggests how plants balance efficient light harvesting with photoprotection. We also explore the influence of specific assumptions in the annihilation simulation, which are challenging to access in more complex environments, such as the thylakoid membrane. Our study provides a framework for studying exciton dynamics using exciton–exciton annihilation, which can be extended to understand the light-harvesting efficiencies of larger, more complex photosynthetic assemblies.

Zhang, Kunyan (ORCID:000000026830409X)↗
DOE OSTI2025

CHARMM-GUI Bicelle Builder : An Extension of Membrane Builder for Modeling and Simulation of Bicelle Systems

Membrane mimetics, such as detergent micelles, nanodiscs, and amphipol complexes, which can provide membrane-like environments while retaining small and soluble features, have been utilized to study membrane proteins. A bicelle, composed of varying lipids and detergents, is a useful membrane mimetic because the lipid-to-detergent ratio, the q-value, can be adjusted to alter the properties of the aggregate, including the thickness and size of the bicelle. However, building a bicelle model for modeling and simulation studies requires nontrivial efforts, even for experts. We introduce CHARMM-GUI Bicelle Builder, a web-based platform that can generate various all-atom bicelle systems via a graphical user interface with all available lipids and detergents in Membrane Builder. To illustrate and validate Bicelle Builder with practical systems, we have modeled and simulated pure bicelles consisting of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) lipids with 1,2-dihexanoyl-sn-glycero-3-phosphocholine (C6DHPC) detergents and protein–bicelle complexes, composed of DMPC with C6DHPC, foscholine-10 (FOS10), and lysophosphatidylcholine-12 (LPC12) detergents. Our simulation results indicate that Bicelle Builder can generate reliable and robust bicelle models with and without proteins that retain DMPC bilayer characteristics. Bicelle Builder is expected to help researchers better understand not only bicelles themselves but also atomistic-level structures of protein–bicelle complexes that are often difficult to access through experimental approaches.

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

FY2019 Energy Efficient Mobility Systems Annual Progress Report

EEMS Program activities during FY 2019 focused on analytical research to understand the impacts that new mobility technologies and services will have at the vehicle, traveler, and overall transportation system-level. This research included the development of vehicle and transportation system simulation models and tools to evaluate the complex interactions among the various actors within the mobility landscape, analysis of empirical data to characterize which solutions may provide the largest benefits, and development of new control systems and algorithms that use vehicle connectivity and automation to improve the performance and efficiency of individual vehicles as well as the overall traffic system. This document presents a brief overview of the EEMS Program and documents progress and results for projects within four of the five EEMS activity areas: (1) the SMART (Systems and Modeling for Accelerated Research in Transportation) Mobility Lab Consortium, (2) High Performance Computing and Big Data Solutions for Mobility Data, (3) Advanced R&D Projects conducted by industry and academia, and (4) Core Modeling, Simulation, and Evaluation, Similarly, the remaining EEMS activity area – (5) Living Labs (managed under VTO’s Technology Integration Program). Each of the individual progress reports provide a project overview and highlights of the technical results.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Editorial: Neuroscience, computing, performance, and benchmarks: Why it matters to neuroscience how fast we can compute

At the turn of the millennium the computational neuroscience community realized that neuroscience was in a software crisis: software development was no longer progressing as expected and reproducibility declined. The International Neuroinformatics Coordinating Facility (INCF) was inaugurated in 2007 as an initiative to improve this situation. The INCF has since pursued its mission to help the development of standards and best practices. In a community paper published this very same year, Brette et al. tried to assess the state of the field and to establish a scientific approach to simulation technology, addressing foundational topics, such as which simulation schemes are best suited for the types of models we see in neuroscience. In 2015, a Frontiers Research Topic “Python in neuroscience” by Muller et al. triggered and documented a revolution in the neuroscience community, namely in the usage of the scripting language Python as a common language for interfacing with simulation codes and connecting between applications. The review by Einevoll et al. documented that simulation tools have since further matured and become reliable research instruments used by many scientific groups for their respective questions. Open source and community standard simulators today allow research groups to focus on their scientific questions and leave the details of the computational work to the community of simulator developers. A parallel development has occurred, which has been barely visible in neuroscientific circles beyond the community of simulator developers: Supercomputers used for large and complex scientific calculations have increased their performance from ~10 TeraFLOPS (10 13 floating point operations per second) in the early 2000s to above 1 ExaFLOPS (10 18 floating point operations per second) in the year 2022. This represents a 100,000-fold increase in our computational capabilities, or almost 17 doublings of computational capability in 22 years. Moore's law (the observation that it is economically viable to double the number of transistors in an integrated circuit every other 18–24 months) explains a part of this; our ability and willingness to build and operate physically larger computers, explains another part. It should be clear, however, that such a technological advancement requires software adaptations and under the hood, simulators had to reinvent themselves and change substantially to embrace this technological opportunity. It actually is quite remarkable that—apart from the change in semantics for the parallelization—this has mostly happened without the users knowing. The current Research Topic was motivated by the wish to assemble an update on the state of neuroscientific software (mostly simulators) in 2022, to assess whether we can see more clearly which scientific questions can (or cannot) be asked due to our increased capability of simulation, and also to anticipate whether and for how long we can expect this increase of computational capabilities to continue.

biophysically detailed models↗
DOE OSTI2023

Microstructure modeling of nuclear structural materials: Recent progress and future directions

Modeling and simulation of microstructures are essential to understand the complex responses and behaviors of nuclear materials in extreme environments. The needs to assess the extended life operation as well as the growing interest in accelerating nuclear materials development and qualification have stimulated the use of high-fidelity multiscale models aided by empirical and ab initio data. This paper reviews the role of various models across different length and time scales in investigating irradiation effects on microstructure evolution and degradation, in particular the embrittlement caused by radiation induced or enhanced formation of nanoscale chemical heterogeneities. The strength and limitations of these models, including classical rate theories, cluster dynamics, phase-field methods, and atomistic models informed by ab initio energies, are discussed with seminal examples. Challenges regarding the lack of thermo-kinetic data and theoretical treatments considering chemical complexities and magnetic excitations, as well as the stabilizing effect by excess point defects in nuclear structural materials are presented, along with potential solutions based on ab initio informed surrogate energy models and statistical sampling by Monte Carlo simulations. Further, the review then highlights the opportunities to leverage the advantages of different methods by establishing hybrid models by shared variables or coupled codes and applications. Finally, the review concludes with forward-looking remarks on how the use of physics-based models can aid the improvement of machine-learning models of property degradation and vice versa.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Thermodynamic Basis for the Stabilization of Helical Peptoids by Chiral Sidechains

Peptoids are a class of highly customizable biomimetic foldamers that retain properties from both proteins and polymers. It has been shown that peptoids can adopt peptide-like secondary structures through the careful selection of sidechain chemistries, but the underlying conformational landscapes that drive these assemblies at the molecular level remain poorly understood. Given the high flexibility of the peptoid backbone, it is essential that methods applied to study peptoid secondary structure formation possess the requisite sensitivity to discriminate between structurally similar yet energetically distinct microstates. In this work, a generalizable simulation scheme is used to robustly sample the complex folding landscape of various 12mer polypeptoids, resulting in a predictive model that links sidechain chemistry with preferential assembly into one of 12 accessible backbone motifs. Using a variant of the metadynamics sampling method, four peptoid dodecamers are simulated in water: sarcosine, N-(1-phenylmethyl)glycine (Npm), (S)-N-(1-phenylethyl)glycine (Nspe), and (R)-N-(1-phenylethyl)glycine (Nrpe)–to determine the underlying entropic and energetic impacts of hydrophobic and chiral peptoid sidechains on secondary structure formation. Our results indicate that the driving forces to assemble Nrpe and Nspe sequences into polyproline type-I helices in water are found to be enthalpically driven, with small benefits from an entropic gain for isomerization and steric strain due to the presence of the chiral center. The minor entropic gains from bulky chiral sidechains in Nrpe- and Nspe-containing peptoids can be explained through increased configurational entropy in the cis state. However, overall assembly into a helix is found to be overall entropically unfavorable. Furthermore, these results highlight the importance of considering the many various competing interactions in the rational design of peptoid secondary structure building blocks.

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

LES ARM Symbiotic Simulation and Observation Scenario for CACTI

The Large-Eddy Simulation (LES) Atmospheric Radiation Measurement (ARM) Symbiotic Simulation and Observation (LASSO) deep-convection scenario for the Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign. The LASSO-CACTI scenario consists of a library of mesoscale and LES simulations spanning grid spacings from 7.5 km to 100 m for convection near the Sierras de Córdoba mountain range, roughly centered on the ARM Mobile Facility. Supporting skill scores and observation comparisons complement the simulations in addition to subsets of the model output for various categories of model output variables. More information can be found at https://www.arm.gov/capabilities/modeling/lasso.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Predicting the Dissolution Rate of Borosilicate Glasses using QSPR analysis based on molecular dynamics simulations

Quantitative Structure Property Relationship (QSPR) analysis based on molecular dynamics (MD) simulations is a promising approach for establishing the composition-property relationships of glass and other materials with complex structures. A series of 20 borosilicate and boroaluminosilicate glasses have been modeled by using MD simulations with recently developed effective potentials. Short- and medium-range structures of these glasses were analyzed and, based on these structural information, QSPR analysis of the initial dissolution rate (r0) was made and compared with measured r0 at 90°C and pH 9 using various structural descriptors such as percentage of bridging oxygen species, network connectivity and average ring size. The structural descriptors, Fnet, containing energetic information such as single bond strength and other structural information were also used. It was found that overall network connectivity, average ring size and Fnet give reasonable predictions of the r0 of studied glasses, given the conditions that the glasses are homogeneous and dissolve congruently. Modifying glass compositions to account preferential release of modifiers gives a better prediction for incongruently dissolving glasses. The results were compared with our recent work of predicting glass dissolution behavior from compositions using the topological-constraints-based models.

Du, Jincheng↗
DOE OSTI2021

The Simons Observatory: HoloSim-ML: machine learning applied to the efficient analysis of radio holography measurements of complex optical systems

Near-field radio holography is a common method for measuring and aligning mirror surfaces for millimeter and sub-millimeter telescopes. In instruments with more than a single mirror, degeneracies arise in the holography measurement, requiring multiple measurements and new fitting methods. We present HoloSim-ML, a Python code for beam simulation and analysis of radio holography data from complex optical systems. This code uses machine learning to efficiently determine the position of hundreds of mirror adjusters on multiple mirrors with few micron accuracy. We apply this approach to the example of the Simons Observatory 6m telescope.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Updated Thermofluid Performance of the Simplified Flat Variant of the HEMJ

Our group has recently developed and studied “finger”-type divertors that are a simplified version of the helium-cooled modular divertor with multiple jets (HEMJ) using coupled computational fluid dynamics and thermal stress simulations. Such a simplified geometry could reduce complexity and cost given the large number of fingers required to cover the total divertor target area. Previous experimental studies for this simplified flat design reported lower heat transfer coefficients and higher pressure drops than the HEMJ, contrary to numerical predictions. Subsequent measurements determined that the original test section had significant dimensional variations in the jet exit holes. A new test section was therefore manufactured and tested in the Georgia Tech (GT) helium loop. The experimental results presented here for this test section at maximum heat flux of 7.1 MW/m 2 are in good agreement with numerical predictions. Correlations developed from these experimental data are extrapolated to predict the maximum heat flux that can be accommodated by the flat design and the coolant pumping power requirements under prototypical conditions. Lastly, numerical simulations are used to estimate the sensitivity of the flat design to geometric variations typical of manufacturing tolerances and variations in the gap width.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Influence of NOx chemistry on the prediction of natural gas end-gas autoignition in CFD engine simulations

Natural gas (NG) represents a promising low-cost/low-emission alternative to diesel fuel when used in high-efficiency internal combustion engines. Advanced combustion strategies utilizing high EGR rates and controlled end-gas autoignition can be implemented with NG to achieve diesel-like efficiencies; however, to support the design of these next-generation NG ICEs, computational tools, including single- and multi-dimensional simulation packages will need to account for the complex chemistry that can occur between the reactive species found in EGR (including NOx) and the fuel. Research has shown that NOx plays an important role in the promotion/inhibition of large hydrocarbon autoignition and when accounted for in CFD engine simulations, can significantly improve the prediction of end-gas autoignition for these fuels. However, reduced NOx-enabled NG mechanisms for use in CFD engine simulations are lacking, and as a result, the influence of NOx chemistry on NG engine operation remains unknown. Here, we analyze the effects of NOx chemistry on the prediction of NG/oxidizer/EGR autoignition and generate a reduced mechanism of a suitable size to be used in engine simulations. Results indicate that NG ignition is sensitive to NOx chemistry, where it was observed that the addition of EGR, which included NOx, promoted NG autoignition. The modified mechanism captured well all trends and closely matched experimentally measured ignition delay times for a wide range of EGR rates and NG compositions. Here, the importance of C2-C3 chemistry is noted, especially for wet NG compositions containing high fractions of ethane and propane. Finally, when utilized in CFD simulations of a Cooperative Fuels Research (CFR) engine, the new reduced mechanism was able to predict the knock onset crank angle (KOCA) to within one crank angle degree of experimental data, a significant improvement compared to previous simulations without NOx chemistry.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Using Co-Simulation to Model Interconnect-Scale Power Systems from Loads to Generators

Co-simulation is a modeling technique that allows analysts to combine simulation tools and their corresponding models to exchange data during run-time, allowing the creation of larger and more complex models across heterogeneous domains. HELICS is a co-simulation platform developed over the past six years that has been shown to be effective for these multi-domain analysis. Recently, a HELICS-based analysis was completed where the ERCOT electrical interconnect in the United States was modeled in high detail from bulk power system generation to individual customer loads. This model was used to evaluate a flat-rate and transactive energy tariff with integrated wholesale and retail real-time and day-ahead energy markets. This modeling allows detailed analysis showing how the operations of the power system under these tariffs impact all actors in the power system, from individual customers to bulk power system operators.

co-simulation, HELICS, transactive energy system, ↗
DOE OSTI2024

Integration of Decentralized Graph-Based Multi-Agent Reinforcement Learning with Digital Twin for Traffic Signal Optimization

Machine learning (ML) methods, particularly Reinforcement Learning (RL), have gained widespread attention for optimizing traffic signal control in intelligent transportation systems. However, existing ML approaches often exhibit limitations in scalability and adaptability, particularly within large traffic networks. This paper introduces an innovative solution by integrating decentralized graph-based multi-agent reinforcement learning (DGMARL) with a Digital Twin to enhance traffic signal optimization, targeting the reduction of traffic congestion and network-wide fuel consumption associated with vehicle stops and stop delays. In this approach, DGMARL agents are employed to learn traffic state patterns and make informed decisions regarding traffic signal control. The integration with a Digital Twin module further facilitates this process by simulating and replicating the real-time asymmetric traffic behaviors of a complex traffic network. The evaluation of this proposed methodology utilized PTV-Vissim, a traffic simulation software, which also serves as the simulation engine for the Digital Twin. The study focused on the Martin Luther King (MLK) Smart Corridor in Chattanooga, Tennessee, USA, by considering symmetric and asymmetric road layouts and traffic conditions. Comparative analysis against an actuated signal control baseline approach revealed significant improvements. Experiment results demonstrate a remarkable 55.38% reduction in Eco_PI, a developed performance measure capturing the cumulative impact of stops and penalized stop delays on fuel consumption, over a 24 h scenario. In a PM-peak-hour scenario, the average reduction in Eco_PI reached 38.94%, indicating the substantial improvement achieved in optimizing traffic flow and reducing fuel consumption during high-demand periods. These findings underscore the effectiveness of the integrated DGMARL and Digital Twin approach in optimizing traffic signals, contributing to a more sustainable and efficient traffic management system.

42 ENGINEERING↗