Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulator study”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,495 records · Page 83

Mechanical behavior of ultralight nickel metamaterial

The mechanical response of ultralight kagomé structures consisting of hollow nickel (Ni) nanotubes and solid Ni nanorods to compression is studied using molecular dynamics simulations. In both kagomé architectures, 16[112] Shockley partial dislocations and twin formation are observed under compression. The structure made from solid nanorods shows deformation near both the nodes and beams of the kagomé lattice. The hollow kagomé architecture has a higher yield point than the solid kagomé structure. The deformation in the hollow nanotube structure is mostly localized in the nodal region for strains less than 11%. At higher strains, the deformation sets in all the struts and nodes of the hollow kagomé lattice. Owing to this two-stage deformation mechanism, the hollow Ni nanotube kagomé structure shows less bending and greater toughness than the solid Ni nanorod kagomé architecture.

Rajak, Pankaj (ORCID:0000000263446056)↗

Evaluation of the Large-Scale and Regional Climatic Response Across North Africa to Natural Variability in Oceanic Modes and Terrestrial Vegetation Among the CMIP5 Models

Hydrologic variability is a serious threat to the poverty-stricken regions of North Africa. Meanwhile, the scientific community struggles to attribute these extreme climatic episodes to specific oceanic and land drivers. Prior modeling studies have not assessed simulated feedbacks over North Africa against an established observational benchmark, so their results are considered as untested, model-specific findings. The Coupled Model Intercomparison Project Phase Five (CMIP5) archive represents the state-of-the-art in climate projections, with most CMIP5 models now containing interactive vegetation phenology; however, there have been no studies to date of their simulated vegetation feedbacks. The capability of CMIP5 models at accurately simulating the forcing of both sea-surface temperature (SST) and regional leaf area index (LAI) anomalies on North African climate needs to be a key consideration in assessing the models’ overall credibility and determining appropriate model weighting for developing climate projections. The Generalized Equilibrium Feedback Assessment (GEFA) is a promising statistical method, which can be applied to either model output or observations, for isolating the local and remote impacts of individual oceanic or terrestrial forcings on regional climate. We performed a combined observational and modeling assessment of land-ocean-atmosphere interactions across the distinct ecological and moisture gradients of North Africa. We evaluated and demonstrated the reliability of the GEFA statistical method over North Africa using the Community Earth System Model (CESM), applied GEFA to observational data to quantify the observed forcing of ocean basin SST anomalies and regional LAI anomalies on North African climate, evaluated the CMIP5 models’ performance in terms of representing these key observed feedbacks, and formulated CMIP5 feedback performance metrics for weighting North African climate projections. Our study represented the first attempt to separate the observed roles of oceanic and vegetation feedbacks across North Africa, the first systematic assessment and intercomparison of land-ocean-atmosphere feedbacks in CMIP5, and the first exploration of vegetation feedbacks among CMIP5 models. The following work was accomplished. (1) The team generated 14 scientific publications and 30 presentations based on the DOE-funded research. (2) A stepwise version of the GEFA statistical method was developed in which unimportant forcings were dropped to increase the reliability of the results. (3) SGEFA was successfully validated through experiments with CESM in terms of its ability to isolate the atmospheric responses to individual oceanic or land forcings. (4) Observational evidence was revealed for the Sahel’s positive vegetation-rainfall feedback on the seasonal to interannual time scale, and it was attributed to a moisture recycling mechanism rather than an albedo mechanism. (5) An approach with SGEFA was developed in which the individual contributions of soil moisture versus vegetation forcings could be separated, revealing that the former forcing outweighed the latter for sub-Saharan Africa. (6) Tropical ocean temperatures were found to be key regulators of pan-tropical vegetation variability, especially for arid and semi-arid regions, including sub-Saharan Africa. (7) The CMIP5 models largely underestimated the importance of land feedbacks across the Sahel. (8) The general consensus among CMIP5 models indicates, for the late 21st century, a diminished seasonal predictability of sub-Saharan African regional climate and an elevated role of the land surface compared to oceanic drivers in regulating regional climate variability. (9) Land-ocean-atmosphere interactions were demonstrated to be key contributors to the seasonal predictability of African wildfire activity. (10) The El Djouf was determined to be of greater importance than the Bodélé depression in terms of providing trans-Atlantic dust transport to the Americas.

54 ENVIRONMENTAL SCIENCES↗

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↗

Direct Numerical Simulation Database of High-Speed Flow over Parameterized Curved Walls

This study presents a direct numerical simulation (DNS) database of high-speed turbulent boundary layers (TBLs) subject to pressure gradients due to parametrically varied backward-facing and forward-facing wall curvatures, with an inflow Mach number of 4.9 and a friction Reynolds number of [Formula: see text] immediately before the onset of wall curvature. The Mach and Reynolds numbers are significantly higher than those reported in the literature for the DNS of pressure-gradient TBLs. The flow conditions and baseline wall geometries are representative of experiments in the high-speed blowdown wind tunnel at the National Aerothermochemistry Laboratory at Texas A&M University. The wall steepness of the baseline geometry for both the backward-facing and forward-facing walls was systematically varied to cause attached, incipiently separated, and fully separated flows. Precomputed flow statistics, including turbulent kinetic energy budgets, are available on the website of the Turbulence Modeling Resource of the NASA Langley Research Center, allowing other investigators to query any property of interest.

Engineering↗

Simulations of Grid-Responsive HVAC Cooling Measures via Ice/PCM storage

This study uses EnergyPlus building energy simulations to assess grid-responsive HVAC cooling measures during peak hours having high electricity price, including: lock the high capacity of a two-speed cooling coil during grid-response hours; and use ice/PCM energy storage to drive a chilled water coil and meet the cooling load, when shut off the main electric cooling coil. Here, the impacts on peak power reductions, energy and utility cost savings, and comfort levels in residential buildings of two U.S. cities were revealed.

Shen, Bo↗

Manipulation and statistical analysis of the fluid flow of polymer semiconductor solutions during meniscus-guided coating

Recent work in structure–processing relationships of polymer semiconductors have demonstrated the versatility and control of thin-film microstructure offered by meniscusguided coating (MGC) techniques. Here, we analyze the qualitative and quantitative aspects of solution shearing, a model MGC method, using coating blades augmented with arrays of pillars. The pillars induce local regions of high strain rates—both shear and extensional—not otherwise possible with unmodified blades, and we use fluid mechanical simulations to model and study a variety of pillar spacings and densities. We then perform a statistical analysis of 130 simulation variables to find correlations with three dependent variables of interest: thin-film degree of crystallinity and transistor field-effect mobilities for charge-transport parallel (Μ para ) and perpendicular (Μ perp ) to the coating direction. Our study suggests that simple fluid mechanical models can reproduce substantive correlations between the induced fluid flow and important performance metrics, providing a methodology for optimizing blade design.

36 MATERIALS SCIENCE↗

Simulations of neutron unbound physics for Geant4

The study of neutron unbound systems via the invariant mass technique is the primary focus of the MoNA Collaboration, which built and operates the Modular Neutron Array (MoNA) and the Large multi-Institutional Scintillator Array (LISA) at FRIB. Advancements in nuclear structure from theory and experiment along the neutron dripline have presented opportunities to understand the nature of unbound systems in higher mass nuclei. The GEometry ANd Tracking (Geant4) platform has been used in high-energy and nuclear physics to simulate particle interactions with as much detail as the user desires. Geant4 currently does not have a physics class to simulate neutron unbound systems. Given the advancement of accelerator facilities and active searches along the neutron dripline, detailed simulations to study the breakup of neutron unbound systems, are necessary. Furthermore, the implementation of the breakup of neutron unbound systems in Geant4 will be presented.

Neutron dripline↗

Isotropic seismic radiation from rock damage and dilatancy

SUMMARY We present semi-analytical and numerical results on isotropic radiation expected to be generated by seismic sources with damaged dilatant materials that sustain during failure an overall shear deformation. The results are derived for inclusions with different stiffness components, including those responsible for shear-volumetric coupling in dilatant material, and different geometries, in a surrounding elastic solid. The obtained relations indicate that the isotropic seismic potency and seismic moment components are proportional to the shear potency components multiplied by the aspect ratio between the thickness and length of the source zone. For a typical rupture zone with a thin tabular geometry, the expected small aspect ratio explains the general seismological observations of small isotropic radiation. However, isotropic radiation induced by dynamic dilation for short and thick sources is more pronounced and may reach for spherical sources sustaining shear deformation 18 per cent of the shear components. Sources that include also tensile components and explosions can produce larger isotropic radiation. The results provide rigorous relations that can be used to estimate different source properties in numerical simulations, experiments and field studies.

Geochemistry & Geophysics↗

Short-term impacts of 2017 western North American wildfires on meteorology, the atmosphere’s energy budget, and premature mortality

Abstract Western North American fires have been increasing in magnitude and severity over the last few decades. The complex coupling of fires with the atmospheric energy budget and meteorology creates short-term feedbacks on regional weather altering the amount of pollution to which Americans are exposed. Using a combination of model simulations and observations, this study shows that the severe fires in the summer of 2017 increased atmospheric aerosol concentrations leading to a cooling of the air at the surface, reductions in sensible heat fluxes, and a lowering of the planetary boundary layer height over land. This combination of lower-boundary layer height and increased aerosol pollution from the fires reduces air quality. We estimate that from start of August to end of October 2017, ∼400 premature deaths occurred within the western US as a result of short-term exposure to elevated PM 2.5 from fire smoke. As North America confronts a warming climate with more fires the short-term climate and pollution impacts of increased fire activity should be assessed within policy aimed to minimize impacts of climate change on society.

Bernstein, Diana N. (ORCID:0000000308279336)↗

Injection simulations of space charge dominated proton beams in IOTA

Injection simulations of space charge dominated proton beams in IOTA ABSTRACT = A 2.5 MeV proton injector has been constructed for the IOTA ring at Fermilab to study the interaction of nonlinear integrable optics (NIO) with high space charge beams. Space charge in the transport line from the RFQ to the injection location has a significant current dependent effect on the phase space. Simulation studies to support efficient injection of intense bunches into IOTA are presented, included schemes to inject directly into NIO lattices.

Wieland, J. [Fermilab]↗

How Well Does the DOE Global Storm Resolving Model Simulate Clouds and Precipitation Over the Amazon?

This study assesses a 40-day 3.25-km global simulation of the Simple Cloud-Resolving E3SM Model (SCREAMv0) using high-resolution ground-based observations from the Atmospheric Radiation Measurement (ARM) Green Ocean Amazon (GoAmazon) field campaign. SCREAMv0 reasonably captures the diurnal timing of boundary layer clouds yet underestimates the boundary layer cloud fraction and mid-level congestus. SCREAMv0 well replicates the precipitation diurnal cycle, however it exhibits biases in the precipitation cluster size distribution compared to scanning radar observations. Specifically, SCREAMv0 overproduces clusters smaller than 128 km, and does not form enough large clusters. Such biases suggest an inhibition of convective upscale growth, preventing isolated deep convective clusters from evolving into larger mesoscale systems. This model bias is partially attributed to the misrepresentation of land-atmosphere coupling. This study highlights the potential use of high-resolution ground-based observations to diagnose convective processes in global storm resolving model simulations, identify key model deficiencies, and guide future process-oriented model sensitivity tests and detailed analyses.

54 ENVIRONMENTAL SCIENCES↗

Orientation and Conformation of Proteins at the Air–Water Interface Determined from Integrative Molecular Dynamics Simulations and Sum Frequency Generation Spectroscopy

Understanding the assembly of proteins at the air-water interface (AWI) informs the formation of protein films, emulsion properties, and protein aggregation. Determination of protein conformation and orientation at an interface is difficult to resolve with a single experimental or simulation technique alone. To date, the interfacial structure of even one of the most widely studied proteins, lysozyme, at the AWI remains unresolved. Here in this study, molecular dynamics (MD) simulations are used to determine if the protein adopts a side-on, head-on, or axial orientation at the AWI with two different forcefields, GROMOS-53a6 + SPC/E and a99SB-disp + TIP4P-D. Vibrational sum frequency generation (SFG) spectroscopy experiments and spectral SFG calculations validate consistency between the structure determined from MD and experiments. Overall, we show with strong agreement that lysozyme adopts an axial conformation at pH 7. Further, we provide molecular-level insight as to how pH influences the binding domains of lysozyme resulting in side-on adsorption near the isoelectric point of the lysozyme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Resolving root causes of experiment discrepancies guided by machine learning

Abstract Scientists rely on accurate experimental data to explain nature and then harness this knowledge for applications addressing human needs. However, discrepancies between experiments of the same observable can impede scientific progress if one does not understand the underlying causes. Here, we developed a process that unravels data discrepancies by first using Bayesian machine learning to relate discrepancies to few of many, potentially biasing metadata features that encode experiment procedures. This machine learning output guides human experts to study discrepancy causes by simulating suspicious aspects of historical experiments or designing modern ones to address open questions. The study findings then lead to rejecting or correcting historical data on firm scientific bases. This process is demonstrated for the energy spectrum of neutrons emitted promptly (<1 ns) after fission of 252 Cf, a trusted nuclear physics Standard. It reduces the spread in experimental 252 Cf spectra by up to a factor of 6.

Neudecker, D. (ORCID:0000000339200627)↗

Conservation of angular momentum in ultrafast spin dynamics

The total angular momentum of a closed system is a conserved quantity, which should remain constant in time for any excitation experiment once the pumping signal has been extinguished. Such conservation, however, is never satisfied in practice in any real-time first-principles description of the demagnetization process. Furthermore, there is growing experimental evidence that the same takes place in experiments. The missing angular momentum is usually associated to lattice vibrations, which are not measured experimentally and are never considered in real-time simulations. Here, in this study, we critically analyze the issue and conclude that current state-of-the-art simulations violate angular momentum conservation already at the electronic level of description. This shortcoming originates from an oversimplified description of the spin-orbit coupling, which includes atomic contributions but neglects completely that of itinerant electrons. We corroborate our findings with time-dependent simulations using model tight-binding Hamiltonians, and show that indeed such conservation can be reintroduced by an appropriate choice of spin-orbit coupling. The consequences of our findings on recent experiments are also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Closing the Gap Between Modeling and Experiments in the Self-assembly of Biomolecules at Interfaces and in Solution

Molecular self-assembly is a powerful tool in materials design, wherein non-covalent interactions like electrostatic, hydrophobic, hydrogen bonding, and van der Waals can be exploited to produce supramolecular nanostructures that are functional and highly tunable. Biomolecules are attractive building blocks, as they are biocompatible, biodegradable and adopt a wide array of higher order structures. Moreover, naturally occurring protein systems display a manifold of structures and interactions that can be replicated in synthetic biomolecules. In this perspective, we highlight advances in multiscale simulation techniques across broad spatiotemporal scales that can aid in characterizing self-assembly of hybrid and hierarchical bionanomaterial systems, with an emphasis on physics-based simulation approaches currently employed to study biomolecules at mineral interfaces. The power of these approaches is highlighted across a few recent areas where molecular simulations have advanced our understanding of self-assembly spanning peptides to protein self-assembly. Looking forward, we discuss how in the near future emerging methods in statistical and machine learning will advance this research field in all areas from expanding the capabilities of physics-based simulation methods to enabling new analyses of high throughput experiments. These advances will pave the way for understanding the molecular recognition patterns in systems that are dictated by self-assembly - biomineralizing peptides, hierarchical peptoids, and large protein assemblies, and will aid in the development of a new synthesis science for achieving precise molecular control in materials design

Sampath, Janani↗

Comparison of tungsten versus molybdenum for double shell capsules using machine learning design optimization

Double shell targets are an alternative ignition platform for inertial confinement fusion. One design consideration for double shell targets is the choice of inner shell material to help trap radiation emitted by the hot fuel to aid ignition. Materials such as molybdenum and tungsten are of interest for the inner shell layer of the targets. While molybdenum has a lower density that could inhibit instability growth and allow for radiography and code benchmarking, tungsten has a higher density that could provide more compression and confinement. These tradeoffs have been explored using optimized designs for each material. Our previous work [Vazirani et al., “Coupling 1D xRAGE simulations with machine learning for graded inner shell design optimization in double shell capsules,” Phys. Plasmas 28, 122709 (2021); Vazirani et al., “Coupling multi-fidelity xRAGE with machine learning for graded inner shell design optimization in double shell capsules,” Phys. Plasmas 30, 062704 (2023); and Vazirani et al., “Bayesian batch optimization for molybdenum versus tungsten inertial confinement fusion double shell target design,” Stat. Anal. Data Min. 17, e11698 (2024)] resulted in a multi-fidelity Bayesian optimization framework to find yield-optimized double shell target geometries. By leveraging simulations of varying fidelities (one-dimensional and two-dimensional) to inform one another, the multi-fidelity optimization was able to optimize a design in the highest fidelity with significantly fewer simulations than would be used in a systematic parameter scan. In this work, we apply the multi-fidelity Bayesian optimization to explore the optimized designs of double shell targets with molybdenum and tungsten inner shells as well as the physics producing the high performing implosions. A physics exploration of all the simulations used in this study shows trends in designs that contribute to high yields, ion temperatures, and fuel areal densities. Comparison of molybdenum and tungsten simulations shows that they can produce similar implosion conditions with different geometries, which would be important to study in experiments. Graded density layers produce varying performances with the two materials but continue to be of interest for future studies along with studies of doped inner shell materials and applied surface roughness.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Deep State Space Model for Rainfall‐Runoff Simulations

The classical way of studying the rainfall‐runoff processes in the water cycle relies on conceptual or physically‐based hydrologic models. Deep learning (DL) has recently emerged as an alternative and blossomed in the hydrology community for rainfall‐runoff simulations. However, the decades‐old Long Short‐Term Memory (LSTM) network remains the benchmark for this task, outperforming newer architectures like Transformers. In this work, we propose a State Space Model (SSM), specifically the Frequency Tuned Diagonal State Space Sequence (S4D‐FT) model, for rainfall‐runoff simulations. The proposed S4D‐FT is benchmarked against the established LSTM and a physically‐based Sacramento Soil Moisture Accounting model under in‐sample and out‐of‐sample simulation setups across 531 watersheds in the contiguous United States (CONUS). Results show that S4D‐FT is able to outperform the LSTM model across diverse regions under both simulation setups, especially for regions that feature snowmelt‐driven or intermittent flow regimes. In contrast, S4D‐FT tends to underperform in flashier, high‐magnitude flow regimes, likely due to its global state‐space convolution computation that emphasizes slow, storage‐driven dynamics, which makes it less effective at picking up short bursts and noisy spikes in the data. In summary, our pioneering introduction of the S4D‐FT for rainfall‐runoff simulations challenges the dominance of LSTM in the hydrology community and expands the arsenal of DL tools available for hydrological modeling.

Wang, Yihan [Univ. of Oklahoma, Norman, OK (United↗

Sampling Rare Events in Aqueous Systems Using Molecular Simulations

Birth of a new distinct phase is a phenomenon encountered in a myriad of processes, and has wide ranging consequences in material processing, biological self-assembly, separations and several other processes. Several phase transitions are nucleation driven. The nucleation events occur over nanosecond timescales and involve hundreds to thousands of molecules. These length and timescales are difficult to access in experiments, thereby making experimental studies of nucleation challenging. On the other hand, molecular simulations sample the nanosecond and nanometer scales making them ideal to study nucleation. However, nucleation is a rare event, meaning that the waiting time to observe one nucleation event is significant. This makes simulation studies of rare events challenging. The project focused on a multi-pronged approach to address such challenges to develop the next generation rare event sampling methods for molecular simulations. The key outcomes of our work include developing more effective methods for sampling rare events, utilizing machine learning to better elucidate nucleation mechanisms, development of software for easy implementation of the methodologies, and applications of the methods to realistic systems to push the method applicability beyond model systems. Overall, this work has enabled pushing the frontiers of molecular simulations to study rare events with a focus on nucleation in aqueous solutions.

36 MATERIALS SCIENCE↗