Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model development”

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 991 records · Page 55

Controlled Environments Enable Adaptive Management in Aquatic Ecosystems Under Altered Environments

Ecosystems worldwide are impacted by altered environment conditions resulting from climate, drought, and land use changes. Gaps in the science knowledge base regarding plant community response to these novel and rapid changes limit both science understanding and management of ecosystems. We describe how CE Technologies have enabled the rapid supply of gap-filling science, development of ecosystem simulation models, and remote sensing assessment tools to provide science-informed, adaptive management methods in the impacted aquatic ecosystem of the California Sacramento-San Joaquin River Delta. The Delta is the hub for California's water, supplying Southern California agriculture and urban communities as well as the San Francisco Bay area. The changes in environmental conditions including temperature, light, and water quality and associated expansion of invasive aquatic plants negatively impact water distribution and ecology of the San Francisco Bay/Delta complex. CE technologies define changes in resource use efficiencies, photosynthetic productivity, evapotranspiration, phenology, reproductive strategies, and spectral reflectance modifications in native and invasive species in response to altered conditions. We will discuss how the CE technologies play an enabling role in filling knowledge gaps regarding plant response to altered environments, parameterization and validation of ecosystem models, development of satellite-based, remote sensing tools, and operational management strategies.

Remote Sensing↗

Towards a New Generation of Agricultural System Data, Models and Knowledge Products: Design and Improvement

This paper presents ideas for a new generation of agricultural system models that could meet the needs of a growing community of end-users exemplified by a set of Use Cases. We envision new data, models and knowledge products that could accelerate the innovation process that is needed to achieve the goal of achieving sustainable local, regional and global food security. We identify desirable features for models, and describe some of the potential advances that we envisage for model components and their integration. We propose an implementation strategy that would link a "pre-competitive" space for model development to a "competitive space" for knowledge product development and through private-public partnerships for new data infrastructure. Specific model improvements would be based on further testing and evaluation of existing models, the development and testing of modular model components and integration, and linkages of model integration platforms to new data management and visualization tools.

Agricultural systems↗

Chemistry and Transport In a Multi-Dimensional Model

The focus of the work funded under this proposal is the exchange between the stratosphere and the troposphere, and between the troposphere and the blaspheme. These two interfaces represent the frontiers of atmospheric chemistry. It is the combination of exchange processes at both interfaces that ultimately controls how the blaspheme (including human activities) affects the ozone layer. The modeling work was motivated by and attempts to integrate information obtained by aircraft, spacecraft, shuttle and oceanic measurements. The model development and research activities accomplished in the past three years provide a technical and intellectual basis for the research in this group. The innovative part of our research program is related to the IAV of ozone and the hydrological cycle. Other related but independently supported work include the study of isotopic fractionation of atmospheric species, e.g., N2O and CO2. Our theory suggests that we now have the ability to probe the middle atmosphere at a level of sensitivity where subtle details such as the isotopic composition of simple molecules can yield measurable systematic effects. This creates the possibility for probing the chemistry and dynamics of the middle atmosphere using all of the N2O and CO2 isotopologues. In the following we will briefly describe the model development and review the highlights of recent accomplishments.

Yung, Yuk L.↗

Distributed Macroscopic Traffic Simulation with Open Traffic Models

This paper presents OTM-MPI, an extension of the Open Traffic Models platform (OTM) for running macroscopic traffic simulations in high-performance computing environments. OTM-MPI represents the first open-source, distributed-memory, macroscopic simulation model developed for modern high performance parallel machines and large networks. Macroscopic simulations are appropriate for studying regional traffic scenarios when aggregate trends are of interest, rather than individual vehicle traces. They are also appropriate for studying the routing behavior of classes of vehicles, such as app-informed vehicles. The network partitioning was performed with METIS. Inter-process communication was done with MPI (message-passing interface). Results are provided for two networks: one realistic network which was obtained from Open Street Maps for Chattanooga, TN, and another larger synthetic grid network. The software recorded a speedups of 198x using 256 cores for Chattanooga, and 475x with 1,024 cores for the synthetic network.

macro-scopic traffic simulation↗

Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows

The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.

Source record↗

Muon Borehole Detector Design for Use in 4-D Density Overburden Monitoring

The increase of atmospheric CO2 concentrations has initiated research into carbon sequestration methods. One possibility is to store CO2 in subsurface porous reservoirs. The monitoring of the injected CO2 plume is vital because escaping CO2 poses health and environmental risks. Typically, seismic reflection methods are used to determine the change in density due to the replacement of brine by CO2 in the reservoir, but this is expensive and not continuous. A potential alternative is to use cosmic muon tomography to measure density changes in the reservoir as a function of time. This paper describes the development of a muon detector that will be capable of being deployed in boreholes. The detector has the required dimensions, an angular resolution of approximately 2 degrees, and is robust. The design is based on alternating layers of scintillating rods, which can provide 4-D tomographic data to detect small changes in density at depths up to approximately 2 km. Geant4, a Monte Carlo simulation code, is being used to develop models to guide the design of the physical configuration. Preliminary testing and measurements have been performed to validate the simulation predictions and optimize the performance of the system. The simulation and preliminary experimental results are presented here.

Borehole Detector, CO¬2 Sequestration, Geant4, Muo↗

The Solar Influencer Next Door: Predicting Low-Income Solar Referrals and Leads

Increasing the adoption of solar among low-to-moderate income (LMI) households remains an important policy goal because of its promise to simultaneously reduce energy burden and support the just distribution of benefits of renewable energy. However, scaling LMI solar remains challenging due to affordability and access issues. Most existing LMI adoption has occurred under public-funded programs, highlighting the importance of increasing the cost-effectiveness of these programs at scale. We develop a new household-level data set on LMI solar lead acquisition, referrals, and adoption to understand the processes through which LMI solar uptake has occurred in California. Then, we develop models to predict two sub-mechanisms in the solar adoption process: whether an otherwise qualified lead becomes "lost" i.e. non-responsive to outreach and, for existing clients, whether they refer solar to others. For the program analyzed, participants received their solar system at no cost, which deemphasizes economic drivers of solar adoption and could differ from other program experiences. Both models substantially improved the accuracy of prediction relative to a baseline. Overall, we find that peer effects and solar economics are important to predicting referrals, and household demographic factors in lead loss prediction. Finally, we find that referrals are both the highest quality and largest source of LMI solar leads, providing a promising mechanism to expand LMI programs further.

customer acquisition costs↗

Progress on HL-LHC Nb 3 Sn Magnets

The high-luminosity Large Hadron Collider (HL-LHC) project aims at allowing to increase the collisions in the LHC by a factor of ten in the decade 2025-2035. One essential element is the superconducting magnet around the interaction region points, where the large aperture magnets will be installed to allow to further reduce the beam size in the interaction point. The core of this upgrade is the Nb 3 Sn triplet, made up of 150-mm aperture quadrupoles in the range of 7-8 m. The project is being shared between the European Organization for Nuclear Research and the US Accelerator Upgrade Program, based on the same design, and on the two strand technologies. The project is ending the short model phase, and entering the prototype construction. We will report on the main results of the short model program, including the quench performance and field quality. A second important element is the 11 T dipole that replaces a standard dipole making space for additional collimators. The magnet is also ending the model development and entering the prototype phase. A critical point in the design of this magnet is the large current density, allowing increase of the field from 8 to 11 T with the same coil cross section as in the LHC dipoles. This is also the first two-in-one Nb 3 Sn magnet developed so far. We will report the main results on the test and the critical aspects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modeling personnel turnover in the parametric organization

A model is developed for simulating the dynamics of a newly formed organization, credible during all phases of organizational development. The model development process is broken down into the activities of determining the tasks required for parametric cost analysis (PCA), determining the skills required for each PCA task, determining the skills available in the applicant marketplace, determining the structure of the model, implementing the model, and testing it. The model, parameterized by the likelihood of job function transition, has demonstrated by the capability to represent the transition of personnel across functional boundaries within a parametric organization using a linear dynamical system, and the ability to predict required staffing profiles to meet functional needs at the desired time. The model can be extended by revisions of the state and transition structure to provide refinements in functional definition for the parametric and extended organization.

Dean, Edwin B.↗

EGS Reservoir Modeling for Developing Geothermal District Heating at Cornell University: Preprint

Cornell University is pursuing development of an enhanced geothermal system (EGS) for providing heating to its main campus in Upstate New York. A ~10,000 ft (~3 km) deep vertical observation well ("CUBO") was drilled in 2022 to characterize the subsurface using wellbore logging, borehole imaging, fluid sampling, mini-frac tests, coring and drill cutting analysis. Down-hole temperatures measured at 3 km depth are about 80 degrees C, sufficiently high for direct-use heating. The well drilled through generally low porosity and low permeability Paleozoic sedimentary formations and into metamorphic basement rock, encountered at about 9,400 ft depth. Leveraging subsurface data obtained through CUBO, we investigated technical feasibility and design requirements of a doublet well system with horizontal laterals connected to a fracture network created through hydraulic fracturing. The EGS reservoir is sized to provide a nominal heat output in the range 5 to 10 MWth of continuous heating over a 15 year-lifetime with limited thermal drawdown. We applied the Gringarten Multiple parallel fracture model, the Cornell Discrete Fracture Simulator FOXFEM and the commercial simulator ResFrac to estimate required heat transfer area and design a potential hydraulic stimulation treatment. Reservoir simulations indicate that, depending on fluid flow rate and injection temperature, 2 to 3 km2 of effective fracture heat transfer area is required to supply the target heat output of 5 to 10 MWth over 15 years.

district heating↗

EGS Reservoir Modeling for Developing Geothermal District Heating at Cornell University

Cornell University is pursuing development of an enhanced geothermal system (EGS) for providing heating to its main campus in Upstate New York. A ~10,000 ft (~3 km) deep vertical observation well ('CUBO') was drilled in 2022 to characterize the subsurface using wellbore logging, borehole imaging, fluid sampling, mini-frac tests, coring and drill cutting analysis. Down-hole temperatures measured at 3 km depth are about 80 degrees C, sufficiently high for direct-use heating. The well drilled through generally low porosity and low permeability Paleozoic sedimentary formations and into metamorphic basement rock, encountered at about 9,400 ft depth. Leveraging subsurface data obtained through CUBO, we investigated technical feasibility and design requirements of a doublet well system with horizontal laterals connected to a fracture network created through hydraulic fracturing. The EGS reservoir is sized to provide a nominal heat output in the range 5 to 10 MWth of continuous heating over a 15 year-lifetime with limited thermal drawdown. We applied the Gringarten Multiple parallel fracture model, the Cornell Discrete Fracture Simulator FOXFEM and the commercial simulator ResFrac to estimate required heat transfer area and design a potential hydraulic stimulation treatment. Reservoir simulations indicate that, depending on fluid flow rate and injection temperature, 2 to 3 km2 of effective fracture heat transfer area is required to supply the target heat output of 5 to 10 MWth over 15 years.

district heating↗

Numerical simulation of spall in lead accounting for its occurrence in solid as well as liquid phases

The capabilities of a recently-developed continuum-level model developed for the purpose of predicting the evolution of spall, including the compaction of damaged material, and applicable to the formation of spall in the material that is fully solid as well as material that is locally liquid, are evaluated against empirical data for lead. Here, the model accounts for inertial forces, elastoplastic constitutive behavior, viscosity, and surface tension. Predictions of spall layer thickness are shown to match the empirical data reasonably well.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Magnetoresistance Hysteresis Evolution in the Granular Y–Ba–Cu–O High-Temperature Superconductor in a Wide Temperature Range

The temperature evolution of the magnetoresistance hysteresis in the granular YBa{sub 2}Cu{sub 3}O{sub 7-δ} high-temperature (T{sub C} ≈ 92 K) superconductor has been investigated. The measurements have been performed in the high-temperature region (78–90 K) and at the liquid helium temperature (4.2 K). The results obtained have been analyzed using the developed model of the behavior of transport properties of a granular high-temperature superconductor in an external magnetic field. Within the discussed model, the dissipation of the grain boundary subsystem is determined by the intergrain spacing-averaged effective field B{sub eff}, which is a superposition of external field H and the field induced by the magnetic moments of superconducting grains. Such a consideration yields the expression B{sub eff}(H) = H − 4πM(H) α for the effective field in the intergrain medium, where M(H) is the experimental hysteretic dependence of magnetization and α is the parameter of magnetic flux crowding in the intergrain medium. Here, the magnetoresistance is assumed to be proportional to the absolute value of the effective field: R(H) ~ |B{sub eff}(H)|. Analysis of the experimental R(H) and M(H) dependences obtained under the same conditions for the investigated high-temperature superconductor sample showed that in the high-temperature region this parameter is α ≈ 25. At the low temperature (4.2 K), we may state that the degree of flux crowding increases and the estimated α value is ~ 50. The estimates made are indicative of the strong effect of flux compression in the intergrain medium on the magnetotransport properties of the investigated granular high-temperature superconductor system. Possible reasons for a discrepancy between the developed model concepts and experimentally observed low-temperature R(H) hysteresis are analyzed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Implementation of A New Microwave Scattering Database and A Forward Model for Active Microwave Sensors in CRTM

Radiative transfer models are extensively used for the assimilation of satellite observations into NWP models as well as retrieving geophysical products from satellite measurements. CRTM is a community model developed by NOAA JCSDA and widely used for different purposes requiring RT calculations. CRTM requires bulk scattering lookup tables in order to perform all-sky RT calculations. However, the current CRTM lookup tables for microwave frequencies were generated based on the Mie theory by assuming spherical frozen particles. The scattering lookup tables generated using the DDA technique has shown to largely improve the RT scattering calculations in the MW region. This presentation targets (i) the implementation and validation of a DDA database that was originally developed for the ARTS RT model into CRTM, and (ii) developing the CRTM active sensor module that takes advantage of the backscattering coefficients computed using the DDA method. The DDA database only provides single scattering properties of different habits, while CRTM requires bulk scattering properties. The CRTM cloud coefficients were previously generated based on the effective radius for representing the size of the particles. However, effective radius is neither measurable nor provided by the NWP models, thus need to be estimated from other geophysical variables such as water content. Therefore, in addition to calculating the CRTM bulk scattering properties from the DDA single scattering database, the CRTM was also largely modified to use cloud water content (kg.m-3), instead of effective radius, for performing the interpolation over size/mass of the particles. CRTM already requires water content as input, thus no extra variables are required for performing scattering calculations using the new ARTS DDA database. The CRTM scattering modules search for effective radius in cloud coefficient files and will use the cloud water content if the effective radius dimension is not found in the cloud coefficient files. Figure 1 shows the CRTM simulated brightness temperatures computed using different cloud coefficients versus ATMS observed values over Hurricane Irma on September 7, 2017 at 18:00 UTC. We used all the cloud water content values included in ERA5 with default CRTM/DDA habits for water, rain, snow, ice, hail, and graupel. ERA5 does not provide separate water content values for ice, hail, and graupel, thus the ice water content values were divided between ice, hail, and graupel clouds similar to what was explained in the previous section. In channels with a frequency lower than 90 GHz, emission from water and rain clouds can compensate for cloud scattering so that cloud contaminated Tbs are larger than corresponding clear sky Tbs. The DDA simulations for channels 1-7 largely perform better than the Mie simulations. The DDA simulations show a mix of small negative and positive simulated minus observed values, while the Mie results show large negative biases. The weighting functions for some of the ATMS temperature sounding channels (channels 9-15) peak mostly above the clouds, therefore the measured Tbs become less sensitive to clouds so that the results of both Mie and DDA become very similar. The Mie lookup tables generate excessive scattering for channel 16, but not enough scattering for the water vapor channels. In the specific case of Hurricane Maria, the DDA lookup tables do not generate enough scattering for channel 16, but the DDA results are much more consistent with observations for water vapor channels than for channel 16. It should be noted that the results may vary if we use other habits to represent snow, hail, and graupel in the DDA simulations. Although these results clearly show the advantage of the DDA database over the Mie dataset, different error sources such as error in the observations, displacement of clouds in the ERA5 reanalysis, and also lack of convective clouds or in general errors in the input atmospheric and cloud profiles contribute to the differences between the simulated and observed values. Aside from the improvements in the simulations, a major advantage of the new dataset is a large number of habits that can be used to tune the data assimilation systems to perform well in different weather conditions.

Isaac Moradi↗

Sodium-Cooled Fast Reactor Reference Plant Model

This report details the progress of Idaho National Laboratory (INL) in creating a reference plant multiphysics model for the Advanced Burner Test Reactor (ABTR). This model was developed under Task 13 of the U.S. Nuclear Regulatory Commission project “Development and Modeling Support for Advanced Non-Light Water Reactors,” and is an extension of the reference plant model developed in Task 4b, which was improved upon in the following ways. (1) The discrete ordinates method was used in lieu of the super-homogenization (SPH)-corrected diffusion approximation in order to better capture the anisotropic scattering contribution and the neutron leakage change due to thermal expansion. (2) The novel neutronic spatial discretization approach, termed the ring-heterogeneous (RH) approximation, was conceptualized and introduced to capture the differential expansion of the materials in the core. This new technique proved capable of preserving fission rates and maintaining the eigenvalue within 2.5% and 266 pcm with 9 neutron energy groups, respectively. Separating the different materials in the core enables the differential expansion of materials to be explicitly accounted for, eliminating the need for problem-specific cross-section functionalization techniques. (3) The SAM model for the core and system thermal-hydraulics analysis was updated to include 61 channels instead of just four representative ones. This enables users to obtain improved spatial resolution for sodium temperature and density scalar fields. (4) All the mesh files were created via the Multiphysics Object-Oriented Simulation Environment (MOOSE) Reactor module, eliminating all reliance on external tools for mesh creation. (5) Finally, the fuel axial expansion now leverages the HT9 and UPuZr material properties that have been validated against experimental data. The reference plant model was used to perform a full-core unprotected loss of flow (ULOF) transient calculation, including neutronics, thermal and mechanical feedback mechanisms. Future work will be devoted to further enhancements of the model. Potential improvements to the model include the addition of the control rod driveline expansion feedback and the upgrading of the support plate model so as to explicitly include 3D effects. Additionally, a Nuclear Energy Advanced Modeling and Simulation funded parallel effort has completely automated the creation of the ring-heterogeneous (RH) mesh from the fully heterogeneous (FH) geometry, thus maximizing user friendliness for the sodium fast reactor sodium-cooled fast reactor (SFR) workflow and will be incorporated in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Ignition and flame characteristics of cryogenic hydrogen releases

In this work, under-expanded cryogenic hydrogen jets were investigated experimentally for their ignition and flame characteristics. The test facility described herein, was designed and constructed to release hydrogen at a constant temperature and pressure, to study the dispersion and thermo-physical properties of cryogenic hydrogen releases and flames. In this study, a non-intrusive laser spark focused on the jet axis was used to measure the maximum ignition distance. The radiative power emitted by the corresponding jet flames was also measured for a range of release scenarios from 37 K to 295 K, 2–6 bar abs through nozzles with diameters from 0.75 to 1.25 mm. The maximum ignition distance scales linearly with the effective jet diameter (which scales as the square root of the stagnant fluid density). A 1-dimensional (stream-wise) cryogenic hydrogen release model developed previously at Sandia National Laboratories (although this model is not yet validated for cryogenic hydrogen) was exercised to predict that the mean mole fraction at the maximum ignition distance is approximately 0.14, and is not dependent on the release conditions. The flame length and width were extracted from visible and infra-red flame images for several test cases. The flame length and width both scale as the square root of jet exit Reynolds number, as reported in the literature for flames from atmospheric temperature hydrogen. As shown in previous studies for ignited atmospheric temperature hydrogen, the radiative power from the jet flames of cold hydrogen scales as a logarithmic function of the global flame residence time. The radiative heat flux from jet flames of cold hydrogen is higher than the jet flames of atmospheric temperature hydrogen, for a given mass flow rate, due to the lower choked flow velocity of low-temperature hydrogen. Lastly, this study provides critical information with regard to the development of models to inform the safety codes and standards of hydrogen infrastructure.

08 HYDROGEN↗

Graph Theory Roots of Spatial Operators for Kinematics and Dynamics

Spatial operators have been used to analyze the dynamics of robotic multibody systems and to develop novel computational dynamics algorithms. Mass matrix factorization, inversion, diagonalization, and linearization are among several new insights obtained using such operators. While initially developed for serial rigid body manipulators, the spatial operators and the related mathematical analysis have been shown to extend very broadly including to tree and closed topology systems, to systems with flexible joints, links, etc. This work uses concepts from graph theory to explore the mathematical foundations of spatial operators. The goal is to study and characterize the properties of the spatial operators at an abstract level so that they can be applied to a broader range of dynamics problems. The rich mathematical properties of the kinematics and dynamics of robotic multibody systems has been an area of strong research interest for several decades. These properties are important to understand the inherent physical behavior of systems, for stability and control analysis, for the development of computational algorithms, and for model development of faithful models. Recurring patterns in spatial operators leads one to ask the more abstract question about the properties and characteristics of spatial operators that make them so broadly applicable. The idea is to step back from the specific application systems, and understand more deeply the generic requirements and properties of spatial operators, so that the insights and techniques are readily available across different kinematics and dynamics problems. In this work, techniques from graph theory were used to explore the abstract basis for the spatial operators. The close relationship between the mathematical properties of adjacency matrices for graphs and those of spatial operators and their kernels were established. The connections hold across very basic requirements on the system topology, the nature of the component bodies, the indexing schemes, etc. The relationship of the underlying structure is intimately connected with efficient, recursive computational algorithms. The results provide the foundational groundwork for a much broader look at the key problems in kinematics and dynamics. The properties of general graphs and trees of nodes and edge were examined, as well as the properties of adjacency matrices that are used to describe graph connectivity. The nilpotency property of such matrices for directed trees was reviewed, and the adjacency matrices were generalized to the notion of block weighted adjacency matrices that support block matrix elements. This leads us to the development of the notion of Spatial Kernel Operator SKO kernels. These kernels provide the basis for the development of SKO resolvent operators.

Jain, Abhinandan↗