Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiscale modeling”

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 775 records · Page 43

The Method of Finite Averages

The Method of Finite Averages (MoFA) is a rigorous multiscale modeling methodology for efficiently modeling multi-physical phenomena in heterogeneous porous media. The code developed in this project aims to perform the numerical calculations required to formulate, implement, and verify MoFA models for Earth and Energy systems (i.e., model verification refers to performing fully-resolved simulations of the systems and comparing their results to those of the models). In general, MoFA transforms partial differential equations (PDEs) describing the fine-scale physics of a system into coupled ordinary differential equations (ODEs)---in time---that describe the coarse-scale---or "average"---physical behaviors of the system. This transformation significantly expedites system simulation, as the coarse-scale ODEs involve vastly fewer degrees of freedom than the fine-scale PDEs. The code developed under this project will allow users to 1.) generate system geometries and numerical meshes, 2.) solve the PDE and ODE systems required for MoFA model formulation and implementation, 3.) solve the PDE systems required to obtain fully-resolved simulation results for model verification, and 4.) compare and plot results (e.g., the model and fully-resolved simulation solutions, the error between the solutions, etc.).

Pietrzyk, KyleM [Lawrence Livermore National Labor↗

The Second Wind Forecast Improvement Project (WFIP 2): General Overview

In 2015 the U.S. Department of Energy initiated a four-year study, the second Wind Forecast Improvement Project (WFIP 2), to improve the representation of boundary-layer physics and related processes in mesoscale models for better wind and wind power forecasts. This goal challenges numerical weather prediction models in complex terrain due to inherent assumptions underlying their boundary-layer parameterizations. The WFIP 2 effort involved the wind industry, universities, NOAA, and DOE's national laboratories in an integrated observational and modeling study. Observations spanned 18-months to assure a full annual cycle of continuously recorded observations from remote-sensing and in situ measurement systems. The study area comprised the Columbia Basin of eastern Washington and Oregon, containing more than 4 GW of installed wind capacity; nests of observational systems captured important atmospheric scales from mesoscale to NWP sub-grid scale. Model improvements targeted NOAA's High-Resolution Rapid Refresh (HRRR) model to facilitate transfer of improvements to NWS operational forecast models. This paper describes the general WFIP 2 scope and objectives, the particular scientific challenges of improving wind forecasts in complex terrain, early successes of the project, and an integrated approach to archiving observations and model output. It provides an introduction for a set of more detailed BAMS papers addressing WFIP 2 observational science, modeling challenges and solutions, practical transfer of boundary layer advances to industry use, and advances in coupling improved mesoscale models to microscale models that can represent interactions between wind plants and the atmosphere.

WFIP2, boundary layer, complex terrain, forecastin↗

Increased Quasi Stationarity and Persistence of Winter Ural Blocking and Eurasian Extreme Cold Events in Response to Arctic Warming. Part II: A Theoretical Explanation

In Part I of this study, it was shown that the Eurasian cold anomalies related to Arctic warming depend strongly on the quasi stationarity and persistence of the Ural blocking (UB). The analysis here revealed that under weak mean westerly wind (MWW) and vertical shear (VS) (quasi barotropic) conditions with weak synoptic-scale eddies and a large planetary wave anomaly, the growth of UB is slow and its amplitude is small. For this case, a quasi-stationary and persistent UB is seen. However, under strong MWW and VS (quasi baroclinic) conditions, synoptic-scale eddies are stronger and the growth of UB is rapid; the resulting UB is less persistent and has large amplitude. In this case, a marked retrogression of the UB is observed. The dynamical mechanism behind the dependence of the movement and persistence of UB upon the background conditions is further examined using a nonlinear multiscale model. The results show that when the blocking has large amplitude under quasi-baroclinic conditions, the blocking-induced westward displacement greatly exceeds the strong mean zonal-wind-induced eastward movement and hence generates a marked retrogression of the blocking. By contrast, under quasi-barotropic conditions because the UB amplitude is relatively small the blocking-induced westward movement is less distinct, giving rise to a quasi-stationary and persistent blocking. It is further shown that the strong mid–high-latitude North Atlantic mean zonal wind is the quasi-barotropic condition that suppresses UB’s retrogression and thus is conducive to the quasi stationarity and persistence of the UB. The model results show that the blocking duration is longer when the mean zonal wind in the blocking region or eddy strength is weaker.

54 ENVIRONMENTAL SCIENCES↗

Finite-element analysis of top-casing electric source method for imaging hydraulically active fracture zones

Imaging hydraulically active fracture zones (HAFZ) is of paramount importance to subsurface resource extraction, geologic storage, and hazardous waste disposal. We have developed advanced 3D finite-element (FE) electrical imaging algorithms for HAFZ in the presence of a steel-cased well. The algorithms use tetrahedral FE meshes in the simulation domain and coarse rectangular finite-difference meshes in the imaging domain. This heterogeneous dual-mesh approach is well suited to modeling the multiscale earth model due to steel-cased wells. We find that the algorithms accurately and efficiently simulate surface electric field measurements over a 3D HAFZ at depth when one end point of a surface electric source is connected to a wellhead. For brevity, this configuration is called the top-casing electric source method. By replacing a hollow cased well with a solid prism, we improve our computational efficiency without affecting the solution accuracy. The sensitivity of the top-casing source method to HAFZ highly depends on the continuity of a steel-cased well because it makes currents preferentially flow to HAFZ. The sensitivity also depends on conductivity structures around the well because they control current leaking from the steel-cased well. We find that the method can image a localized HAFZ and detect changes in its width and height. Here, the imaging results are improved when a volume of the imaging domain is constrained from geomechanical perspectives. A primary advantage of the method is the fact that the sources and receivers are placed on the surface, thus not interrupting the well operation.

3D↗

Multiscale affinity maturation simulations to elicit broadly neutralizing antibodies against HIV

The design of vaccines against highly mutable pathogens, such as HIV and influenza, requires a detailed understanding of how the adaptive immune system responds to encountering multiple variant antigens (Ags). Here, we describe a multiscale model of B cell receptor (BCR) affinity maturation that employs actual BCR nucleotide sequences and treats BCR/Ag interactions in atomistic detail. We apply the model to simulate the maturation of a broadly neutralizing Ab (bnAb) against HIV. Starting from a germline precursor sequence of the VRC01 anti-HIV Ab, we simulate BCR evolution in response to different vaccination protocols and different Ags, which were previously designed by us. The simulation results provide qualitative guidelines for future vaccine design and reveal unique insights into bnAb evolution against the CD4 binding site of HIV. Our model makes possible direct comparisons of simulated BCR populations with results of deep sequencing data, which will be explored in future applications.

Conti, Simone (ORCID:000000017986567X)↗

Geospatial Information, Metadata, and Maps for Global River Corridor Science Focus Area Sites (v5)

This dataset provides geospatial information, metadata, and maps for the Pacific Northwest National Laboratory (PNNL) River Corridor Science Focus Area (RC-SFA; https://www.pnnl.gov/projects/river-corridor) sites. The RC-SFA works to transform understanding of spatial and temporal dynamics in river corridor hydrobiogeochemical functions from molecular reaction to watershed and basin scales. The knowledge we gain is used to formulate and test hypotheses and to improve mechanistic representation of river corridor processes and their response to disturbances in multiscale models of integrated hydrobiogeochemical function. The data provided includes Site ID, latitude, longitude, stream name, and common ID (COMID) for sites used across the RC-SFA. The COMID can be used to find and download data from NHDPlus (https://www.epa.gov/waterdata/nhdplus-national-hydrography-dataset-plus) and other platforms. The sites included are non-exhaustive. Sites (including past sites) will be added to this data package in the future. Data generated from the RC SFA can be accessed at https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA. This data package was originally published in April 2023. It was updated in June 2023 (v2; modified files), December 2023 (v3; modified files), January 2025 (v4; modified files), and December 2025 (v5; modified files). See the change history section in the readme for more details. This dataset is comprised of one main data folder. The data folder consists of (1) file-level metadata; (2) data dictionary; (3) readme; (4) methods codes; (5) geospatial information for all RC SFA sites including International Generic Sample Number (IGSN); (6) maps of all sites and sites in Washington State, USA; and (7) a subfolder with the shapefile of all sites. All files are .csv, .pdf, .shp, .cpg, .dbf, .prj, .qmd, or .shx. We thank the Confederated Tribes and Bands of the Yakama Nation for access to field locations where some data were collected in Washington state. We also thank the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate sample collection and optimization of data usage according to their values and worldview.

54 ENVIRONMENTAL SCIENCES↗

Characterization of microstructure and property evolution in advanced cladding and duct: Materials exposed to high dose and elevated temperature

Designing materials for performance in high-radiation fields can be accelerated through a carefully chosen combination of advanced multiscale modeling paired with appropriate experimental validation. Here, the studies reported in this work, the combined efforts of six universities working together as the Consortium on Cladding and Structural Materials, use that approach to focus on improving the scientific basis for the response of ferritic–martensitic steels to irradiation. A combination of modern modeling techniques with controlled experimentation has specifically focused on improving the understanding of radiation-induced segregation, precipitate formation and growth under radiation, the stability of oxide nanoclusters, and the development of dislocation networks under radiation. Experimental studies use both model and commercial alloys, irradiated with both ion beams and neutrons. Lastly, transmission electron microscopy and atom probe are combined with both first-principles and rate theory approaches to advance the understanding of ferritic–martensitic steels.

36 MATERIALS SCIENCE↗

Energy Frontier Research Center -- Center for Materials Science of Nuclear Fuels

The mission of the Center for Material Science of Nuclear Fuels (CMSNF) is to develop a first principles-based understanding of how complex defects affect thermal transport in irradiated nuclear fuel (specifically UO 2 ), and to develop an experimentally validated multiscale modeling capability for the predictive understanding of the microstructure dynamics and its impact on thermal transport in nuclear fuel.

36 MATERIALS SCIENCE↗

2169 Steel Waveform Experiments

In support of LLNL efforts to develop multiscale models of a variety of materials, we have performed a set of eight gas gun impact experiments on 2169 steel (21% Cr, 6% Ni, 9% Mn, balance predominantly Fe). These experiments provided carefully controlled shock, reshock and release velocimetry data, with initial shock stresses ranging from 10 to 50 GPa (particle velocities from 0.25 to 1.05 km/s). Both windowed and free-surface measurements were included in this experiment set to increase the utility of the data set, as were samples ranging in thickness from 1 to 5 mm. Target physical phenomena included the elastic/plastic transition (Hugoniot elastic limit), the Hugoniot, any phase transition phenomena, and the release path (windowed and free-surface). The Hugoniot was found to be nearly linear, with no indications of the Fe α – ε phase transition. Releases were non-hysteretic, and relatively consistent between 3- and 5-mm–thick samples (the 3 mm samples giving slightly lower wavespeeds on release). Reshock tests with explosively welded impactors produced clean results; those with glue bonds showed transient releases prior to the arrival of the reshock, reducing their usefulness for deriving strength information. The free-surface samples, which were steps on a single piece of steel, showed lower wavespeeds for thin (1 mm) samples than for thicker (2 or 4 mm) samples. A configuration used for the last three shots allows release information to be determined from these free surface samples. The sample strength appears to increase with stress from ~1 GPa to ~ 3 GPa over this range, consistent with other recent work but about 40% above the Steinberg model.

Furnish, Michael David↗

Final Technical Report for Award # ER64999

This report provides a summary of activities for Award # ER64999, a Genomes to Life Project funded by the Office of Science, Basic Energy Research. The project was entitled "Methanogenic archaea and the global carbon cycle: a systems biology approach to the study of Methanosarcina species". The long-term goal of this multi-investigator project was the creation of integrated, multiscale models that accurately and quantitatively predict the role of Methanosarcina species in the global carbon cycle under dynamic environmental conditions. To achieve these goals we pursed four specific aims: (1) genome sequencing of numerous members of the Order Methanosarcinales, (2) identification of genomic sources of phenotypic variation through in silico comparative genomics, (3) elucidation of the transcriptional networks of two Methanosarcina species, and (4) development of comprehensive metabolic network models for characterized strains to address the question of how metabolic models scale with genetic distance.

54 ENVIRONMENTAL SCIENCES↗

Energy Frontier Research Center (EFRC): Center for Materials Science of Nuclear Fuels (CMSNF) Major Accomplishments

The mission of the Center for Material Science of Nuclear Fuels (CMSNF) has been to develop a first-principles-based understanding of thermal transport in the most widely used nuclear fuel, UO 2 , in the presence of defect microstructure associated with radiation environments. The overarching goal within this mission was to develop an experimentally validated multiscale modeling capability directed toward a predictive understanding of the impact of radiation and fission-product induced defects and microstructure on thermal transport in nuclear fuel. Implementation of the mission was accomplished by integrating the physics of thermal transport in crystalline solids with microstructure science under irradiation through multi institutional experimental and computational materials theory teams from Idaho National Laboratory, Oak Ridge National Laboratory, Purdue University, the University of Florida, the University of Wisconsin, and the Colorado School of Mines. The Center’s research focused on five major areas: (i) The fundamental aspects of anharmonicity in UO 2 crystals and its impact on thermal transport; (ii) The effects of radiation microstructure on thermal transport in UO 2 ; (iii) The mechanisms of defect clustering in UO 2 under irradiation; (iv) The effect of temperature and oxygen environment on the stoichiometry of UO 2 ; and (v) The mechanisms of growth of dislocation loops and voids under irradiation. The Center has made important progress in each of these areas, as summarized below.

36 MATERIALS SCIENCE↗

High-Resolution Characterization of UMo Alloy Microstructure

This report highlights the capabilities and procedure for high-resolution characterization of UMo fuels in PNNL. Uranium-molybdenum (UMo) fuel processing steps, from casting to forming final fuel, directly affect the microstructure of the fuel, which in turn dictates the in-reactor performance of the fuel under irradiation. In order to understand the influence of processing on UMo microstructure, microstructure characterization techniques are necessary. Higher-resolution characterization techniques like transmission electron microscopy (TEM) and atom probe tomography (APT) are needed to interrogate the details of the microstructure. The findings from TEM and APT are also directly beneficial for developing predictive multiscale modeling tools that can predict the microstructure as a function of process parameters. This report provides background on focused-ion-beam–based TEM and APT sample preparation, TEM and APT analysis procedures, and the unique information achievable through such advanced characterization capabilities for UMo fuels, from a fuel fabrication capability viewpoint.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

TA [B] Predicting Microstructure-Creep Resistance Correlation in High Temperature Alloys over Multiple Time Scales

DoE-NETL partnered with Purdue University to predict the creep and associated microstructure evolution of tungsten-based refractory alloys. Researchers use grain boundary (GB) diagrams, a new concept, to establish time-dependent creep resistance and associated microstructure evolution of grain boundaries/intergranular films GB/IGF controlled creep as a function of load, environment, and temperature. The goal was to conduct a systematic study that includes the development of a theoretical framework, multiscale modeling, and experimental validation using W-based body-centered-cubic alloys, doped/alloyed with one or two of the following elements: nickel, palladium, cobalt, iron, and copper—typical refractory alloys. Prior work has already established and validated a basic theory for W-based binary and ternary alloys; the study conducted under this project extended this proven work. Based on interface diagrams phase field models were developed to predict long term microstructural evolution. In order to validate the models nanoindentation creep data was used to elucidate the role played by the interface properties in predicting long term creep strength and microstructure evolution.

20 FOSSIL-FUELED POWER PLANTS↗

Advanced Technology and Mitigation (ATDM) SPARC Re-Entry Code Fiscal Year 2017 Progress and Accomplishments for ECP.

The SPARC (Sandia Parallel Aerodynamics and Reentry Code) will provide nuclear weapon qualification evidence for the random vibration and thermal environments created by re-entry of a warhead into the earth’s atmosphere. SPARC incorporates the innovative approaches of ATDM projects on several fronts including: effective harnessing of heterogeneous compute nodes using Kokkos, exascale-ready parallel scalability through asynchronous multi-tasking, uncertainty quantification through Sacado integration, implementation of state-of-the-art reentry physics and multiscale models, use of advanced verification and validation methods, and enabling of improved workflows for users. SPARC is being developed primarily for the Department of Energy nuclear weapon program, with additional development and use of the code is being supported by the Department of Defense for conventional weapons programs.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

High Fidelity Fast Running Multi-scale Multi-physics Battery Pack Software

The project team comprising of General Motors LLC (GM), ANSYS, ESim, and the National Renewable Energy Laboratory (NREL) proposed to develop multiscale multiphysics battery pack engineering software that significantly surpassed the most advanced state-of-the-art software in computational speed and accuracy. The successful software was to run faster than current state-of-the-art software by a factor of 100 while resolving complex multiphysics interplays across the scales from materials to large automotive battery pack systems. The proposed project was to shift the paradigm to using software for battery pack system design and evaluation, potentially revolutionizing the standard development process for the broader electric vehicle industry. The project objectives were to: (a) enhance computation speed and stability of the pack level multiscale model while still resolving nonlinearities of the battery’s dynamic responses across the scales from particles to a battery pack and (b) enhance multiphysics integrity at the pack level model enabling model-based diagnostics and prognostics for life and safety.

25 ENERGY STORAGE↗

Theory and Simulation of Nanoscale Phenomena

The Center for Integrated Nanotechnologies (CINT) is a Department of Energy/Office of Science Nanoscale Science Research Center (NSRC), operating as a national user facility devoted to establishing the scientific principles that govern nanoscale integration. Nanoscale integration is defined as assembling diverse nanoscale materials across length scales to design and achieve new properties and functionality. The CINT Theory and Simulation of Nanoscale Phenomena thrust is the component of CINT dedicated to developing and applying theory to enable nanoscale integration. Our focus is on understanding and simulating the unique behavior of integrated materials and systems with nanoscale structure. This mission is achieved through collaborations with CINT Users, between thrust scientists, and with CINT scientists from other thrusts. Our research is focused on three science directions that together form the basis for integration at the nanoscale, namely (i) Hierarchical structure and dynamics in soft matter, (ii) Excitation and Transport in Nanostructured Systems, and (iii) Emergent phenomena at surfaces and interfaces. A broad spectrum of techniques is developed and applied including continuum fluid theory, atomistic and coarse-grained molecular dynamics simulations, static and dynamic electronic structure calculations, multiscale modeling, low-energy effective Hamiltonian methods, and perturbative and exact quantum many-body approaches. These tools are applied to physical systems of interest to CINT Users, the other CINT thrusts, and the general scientific community with the goals of understanding and controlling the interactions between nanoscale building blocks to assemble specific integrated structures, controlling energy transfer and other interactions over multiple length scales, and designing and exploiting the interactions within assembled structures to achieve new materials functionality.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Deformation Mechanisms at Grain Boundaries in Polycrystals

Grain boundaries (GBs) are important sources of deformation accommodation and strengthening in polycrystalline materials. Grain boundary sliding is a deformation mechanism that is known to activate at high temperatures, low strain rates, and for materials with small grain sizes. While crystalline slip has been extensively studied, investigations of grain boundary sliding has been mostly restricted to nanocrystalline materials, symmetric or synthetic GBs, and loading in high-temperature environments. This study is aimed at improving our understanding of the coupling between deformation mechanism relationships, namely between slip-GB interaction and grain boundary sliding, by means of a combined high-resolution experimental characterization in conjunction with multiscale modeling. The material of choice in the present study is high-purity aluminum with a through-the-specimen-thickness grain structure. Digital image correlation experiments are performed to measure 2D strains relative to the microstructural features. These experiments, performed at various load steps, can resolve individual slip bands and the onset of grain boundary sliding. Onset of grain boundary sliding is observed in the first loadstep (the first loadstep was collected just after macroscopic yield), which suggests that the mechanism is driven by high stress as opposed to only being active to accommodate large values of plastic strain. Further, grain boundary sliding is observed in cases where slip impingement is present, as well as at boundaries that do not interact with obvious slip events. A temperature-dependent crystal plasticity model is used to study the constitutive response, which is then compared to full-field strain response experimentally obtained at the microscale-level. From this analysis, it is determined that a combination of the resolved shear stress at the grain boundary plane and boundary straightness form a set of criteria for predicting sliding propensity. Further, the GBs characterized from the experiments are exactly reproduced in a molecular dynamics framework as bicrystal systems, and the underlying atomistic scale mechanisms for pure sliding are studied. It is seen that GBs in the present study undergo sliding similar to a fluid-like flow. To study slip and sliding interaction, notches are introduced to facilitate dislocation emission and eventual interaction with the sliding GB. It was seen that slipGB interaction causes local increases in stress concentration, and a subsequent increase in the local sliding rate. These atomistic scale insights are used to inform a modified 3 crystal plasticity model with explicit GB elements and additional degrees of freedom to include GB sliding. The modified model results in improved predictions of deformation near the GBs than a conventional crystal plasticity model. These results improve the sophistication and accuracy of deformation modeling by establishing a set of criteria to describe the grain boundary sliding mechanism and help better predict micromechanical quantities near the grain boundaries.

36 MATERIALS SCIENCE↗

Local-nonlocal coupling in Emu/PDMS

A technique called the splice method for coupling local to peridynamic subregions of a body is described. The method relies on ghost nodes, whose values of displacement are interpolated from nearby physical nodes, to make each subregion visible to the other. In each time step, the nodes in each subregion treat the nodes in the other subregion as boundary conditions. Adaptively changing the subregions is possible through the creation and deletion of ghost nodes. Example problems in 2D and 3D illustrate how the method is used to perform multiscale modeling of fracture and impact events within a larger structure.

42 ENGINEERING↗