Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled Simulation”

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,585 records · Page 88

Assessment of the Contribution of the Atmosphere to Uncertainties in Normalized Water-Leaving Radiance: A Combined Modeling and Data Analysis Approach

Our research for the Sensor Intercomparison and Merger for Biological and Interdisciplinary Oceanic Studies (SIMBIOS) program has been focused on modeling and simulation studies as well as the development of atmospheric correction algorithms. Based on our original proposal and discussions at the SIMBIOS team meetings and workshops, the objectives of our research can be summarized as follows: (1) Use our radiative transfer model for the coupled atmosphere-ocean system to simulate the radiation field at arbitrary levels and in any desired direction in the atmosphere and ocean so as to provide a firm connection between the signal received by the satellite sensor and by a sensor looking down into the water column just above the surface and just below it; (2) Use the simulations to quantify the influence of atmospheric aerosols on the water-leaving radiance, and to quantify the error in the water-leaving radiance as a function of uncertainties in the aerosol optical properties, mass loading and vertical extent; (3) Use the model simulations in conjunction with validation measurements taken by other SIMBIOS investigators (for satellite overpasses) to assess our understanding of the radiative transfer process in the coupled atmosphere-ocean column, and to examine the extent to which the model provides a realistic prediction of simultaneously measured in situ water-leaving radiance and the radiance received by the satellite sensor; (4) Modify and improve an existing atmospheric correction algorithm, based on the work above, as needed by constructing new look-up tables that include scattering by ocean particles; and (5) Carry out approaches to developing an atmospheric correction algorithm for ocean color imagery with strongly absorbing aerosols.

Stamnes, Knut↗

Synthetic, Realistic Transmission and Distribution Co-Simulation for Voltage Control Benchmarking

As distributed energy resources (DERs) play an increasing role in generation and grid services, control and market paradigms will need to be tested with models that incorporate transmission and distribution systems. The test systems will need to be realistic and accessible for researchers to form conclusions about realistic behaviors and maintain comparability and repeatability of results. These test systems also must have high resolution models of sufficient size for realistic distribution and transmission interaction. This paper demonstrates the tightly coupled HELICS-based co-simulation of realistic synthetic distribution and transmission models for the Austin metropolitan area and how this framework can be used to analyze resultant behaviors from distributed resources and their controls. These results are compared to the same distribution system run without co-simulation using a fixed feed-in voltage. Results show that even without advanced voltage controls, the feed-in voltage response when coupled to the transmission system results a lower level of voltage excursions, suggesting the importance of using this type of realistic, tightly coupled, transmission and distribution cosimulation when assessing DERs.

27 ARPA - Advanced Research Projects Agency-Energy↗

Accelerating the rate of discovery: toward high-repetition-rate HED science

As high-intensity short-pulse lasers that can operate at high-repetition-rate (HRR) (>10 H z ) come online around the world, the high energy density (HED) science they enable will experience a radical paradigm shift. The >10 3 increase in shot rate over today's shot-per-hour drivers translates into dramatically faster data acquisition, more experiments, and the ability to exploit machine learning, and thus the potential to significantly accelerate the advancement of HED science. A wide range of HED experiments, from opacity investigations to secondary source generation to plasma nuclear physics, will benefit from the increased statistics, precision, and exploration of phase space. Besides increasing the rate at which scientific experiments can be performed, HRR also allows for the rapid delivery of optimal experiments supported by simulations and modeling augmented by close coupling to empirical data. To fully realize such an HRR framework, numerous subsystems must be developed and brought together, including feedback laser control loops, high-throughput targetry and diagnostics, cognitive simulation, enhanced HED codes, and advanced data analytics. This paper describes the vision for an integrated HRR laser experimental HED system and outlines some of the major considerations and challenges for realizing it.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Unsteady transonic aerodynamics and aeroelastic calculations at low-supersonic freestreams

A computational procedure is presented to simulate transonic unsteady flows and corresponding aeroelasticity of wings at low-supersonic freestreams. The flow is modeled by using the transonic small-perturbation theory. The structural equations of motions are modeled using modal equations of motion directly coupled with aerodynamics. Supersonic freestreams are simulated by properly accounting for the boundary conditions based on pressure waves along the flow characteristics in streamwise planes. The flow equations are solved using the time-accurate, alternating-direction implicit finite-difference scheme. The coupled aeroelastic equations of motion are solved by an integration procedure based on the time-accurate, linear-acceleration method. The flow modeling is verified by comparing calculations with experiments for both steady and unsteady flows at supersonic freestreams. The unsteady computations are made for oscillating wings. Comparisons of computed results with experiments show good agreement. Aeroelastic responses are computed for a rectangular wing at Mach numbers ranging from subtransonic to upper-transonic (supersonic) freestreams. The extension of the transonic dip into the upper transonic regime is illustrated.

Guruswamy, Guru P.↗

A TOpographic Mapping (ATOM) Method to Design Magnetic Cores

Wireless power transfer offers safe, convenient, and efficient way of charging electric vehicles. Ongoing research is targeting wireless charging pad design optimization; designing the magnetic component is the most important part of the coupler design because the magnetic part determines the coupling factor and efficiency. Optimizing the coil layout and geometry as well as ferrite design requires finite elements analysis based modeling and simulation for minimized core losses, maximized magnetic coupling, and minimized material use for cost-effectiveness. Although parametric finite element analysis or emerging artificial intelligence methods can generate very accurate results, simulation times are extremely long. To address this issue, this study proposes a simple, effective core design called A TOpographic Mapping (ATOM). The proposed design is based on the design of magnetic core by using the magnetic flux distribution. The thickness of the core increases with increasing magnetic flux density, forming a variable thickness core design with less material and minimized core losses compared to conventional designs. A superimposing method is used to create an optimal design for a rotational magnetic field-based system. According to simulation results, the ATOM design reduces the required material volume by 13.19% and yields the lowest core loss and highest mutual inductance compared to other designs. In addition, misalignment, electromagnetic interference, and thermal performance were evaluated for the proposed design.

Aydin, Emrullah [Oak Ridge National Laboratory (OR↗

EVALUATION OF HOT CHANNEL FACTOR FOR SODIUM-COOLED FAST REACTORS WITH MULTI-PHYSICS TOOLKIT

The evaluation of hot channel factor (HCF) is of great significance to the quantification of safety margins for reactor designs. In this paper, HCFs for a sodium-cooled fast reactor (SFR) are evaluated with the Simulation-based High-efficiency Advanced Reactor Prototyping (SHARP) toolkit, which is developed under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) Campaign of DOE for multi-physics reactor performance and safety simulations. The high-fidelity neutronics and thermal hydraulics solvers PROTEUS and Nek5000 in the SHARP toolkit are coupled to perform the multi-physics simulations for HCF evaluation. The HCFs induced by cladding manufacturing tolerance, fissile content mal-distribution, wire orientation and uncertainties on the cladding, coolant, and fuel properties are evaluated for a reference core SFR design (AFR-100). The HCFs calculated with the SHARP toolkit are compared to legacy HCFs for similar reactor types. The comparison demonstrates the reduction or elimination of modeling uncertainties in the calculation of HCFs using high fidelity advanced modeling and simulation tools without the need of expensive experiments. Moreover, the reduction of the uncertainties on HCFs evaluation allows an increase in nominal parameters and safety margin, which in turn improves the economic competitiveness of the SFR.

Hot channel factor (HCF)↗

Evolution of microstructures in radiation fields using a coupled binary-collision Monte Carlo phase field approach

The simulation of radiation effects in materials broadly falls into two categories. At short time and length scales lies the modeling of primary radiation damage, such as point defect creation, energy deposition, and ballistic mixing. This is followed by the modeling at longer time scales of thermally activated microstructure evolution and defect reactions, such as recombination, clustering, and coarsening. The binary collision Monte Carlo method is an established, numerically efficient method for the computation of primary radiation damage. Conversely, the phase field method is a state of the art method for microstructure evolution on longer time and length scales. Here we present a concurrent coupling of these two methods, overcoming the difference between the discrete object Monte Carlo paradigm for primary radiation damage and the continuum field variable approach for microstructure evolution. The coupling is bidirectional, in which the microstructure evolution in the MOOSE finite element frame- work provides the spatial scattering data set for the charged particle trans- port and receives point defect, mass transport, and heat source terms from the simulated collision cascades that contribute to the field variable evolution. The concurrent coupling scheme is implemented in the code Magpie and demonstrated by investigating patterning for a model irradiated immiscible binary alloy. The results from the coupled binary collision Monte Carlo/phase field simulations reproduce the results of analytical models for phase separation, phase mixing, and patterning, supporting the approach and indicating its utility for modeling real materials systems.

36 MATERIALS SCIENCE↗

Powersheds

Powersheds is an open scientific software project for simulating river–reservoir cascades. It combines the performance of Rust with a friendly Python interface to model storage, pool elevation, head, releases, spills, routing lags, and power generation at hourly resolution. Designed for coupling with power-system models, simulations are driven by plant-level target power schedules and report realized generation after accounting for hydrologic and operational constraints.

Turner, Sean [Oak Ridge National Laboratory (ORNL)↗

Improving Climate Projections Using "Intelligent" Ensembles

Recent changes in the climate system have led to growing concern, especially in communities which are highly vulnerable to resource shortages and weather extremes. There is an urgent need for better climate information to develop solutions and strategies for adapting to a changing climate. Climate models provide excellent tools for studying the current state of climate and making future projections. However, these models are subject to biases created by structural uncertainties. Performance metrics-or the systematic determination of model biases-succinctly quantify aspects of climate model behavior. Efforts to standardize climate model experiments and collect simulation data-such as the Coupled Model Intercomparison Project (CMIP)-provide the means to directly compare and assess model performance. Performance metrics have been used to show that some models reproduce present-day climate better than others. Simulation data from multiple models are often used to add value to projections by creating a consensus projection from the model ensemble, in which each model is given an equal weight. It has been shown that the ensemble mean generally outperforms any single model. It is possible to use unequal weights to produce ensemble means, in which models are weighted based on performance (called "intelligent" ensembles). Can performance metrics be used to improve climate projections? Previous work introduced a framework for comparing the utility of model performance metrics, showing that the best metrics are related to the variance of top-of-atmosphere outgoing longwave radiation. These metrics improve present-day climate simulations of Earth's energy budget using the "intelligent" ensemble method. The current project identifies several approaches for testing whether performance metrics can be applied to future simulations to create "intelligent" ensemble-mean climate projections. It is shown that certain performance metrics test key climate processes in the models, and that these metrics can be used to evaluate model quality in both current and future climate states. This information will be used to produce new consensus projections and provide communities with improved climate projections for urgent decision-making.

Baker, Noel C.↗

Improving Climate Projections Using "Intelligent" Ensembles

Recent changes in the climate system have led to growing concern, especially in communities which are highly vulnerable to resource shortages and weather extremes. There is an urgent need for better climate information to develop solutions and strategies for adapting to a changing climate. Climate models provide excellent tools for studying the current state of climate and making future projections. However, these models are subject to biases created by structural uncertainties. Performance metrics-or the systematic determination of model biases-succinctly quantify aspects of climate model behavior. Efforts to standardize climate model experiments and collect simulation data-such as the Coupled Model Intercomparison Project (CMIP)-provide the means to directly compare and assess model performance. Performance metrics have been used to show that some models reproduce present-day climate better than others. Simulation data from multiple models are often used to add value to projections by creating a consensus projection from the model ensemble, in which each model is given an equal weight. It has been shown that the ensemble mean generally outperforms any single model. It is possible to use unequal weights to produce ensemble means, in which models are weighted based on performance (called "intelligent" ensembles). Can performance metrics be used to improve climate projections? Previous work introduced a framework for comparing the utility of model performance metrics, showing that the best metrics are related to the variance of top-of-atmosphere outgoing longwave radiation. These metrics improve present-day climate simulations of Earth's energy budget using the "intelligent" ensemble method. The current project identifies several approaches for testing whether performance metrics can be applied to future simulations to create "intelligent" ensemble-mean climate projections. It is shown that certain performance metrics test key climate processes in the models, and that these metrics can be used to evaluate model quality in both current and future climate states. This information will be used to produce new consensus projections and provide communities with improved climate projections for urgent decision-making.

Baker, Noel C.↗

Self-consistent modelling of the polar thermosphere and ionosphere to magnetospheric convection and precipitation (invited review)

It has recently been demonstrated that the dramatic effects of plasma precipitation and convection on the composition and dynamics of the polar thermosphere and ionosphere include a number of strong interactive, or feedback, processes. To aid the evaluation of these feedback processes, a joint three dimensional time dependent global model of the Earth's thermosphere and ionosphere was developed in a collaboration between University College London and Sheffield University. This model includes self consistent coupling between the thermosphere and the ionosphere in the polar regions. Some of the major features in the polar ionosphere, which the initial simulations indicate are due to the strong coupling of ions and neutrals in the presence of strong electric fields and energetic electron precipitation are reviewed. The model is also able to simulate seasonal and Universal time variations in the polar thermosphere and ionospheric regions which are due to the variations of solar photoionization in specific geomagnetic regions such as the cusp and polar cap.

Rees, D.↗

Tree tensor network hierarchical equations of motion based on time-dependent variational principle for efficient open quantum dynamics in structured thermal environments

In this work, we introduce an efficient method, TTN-HEOM, for exactly calculating the open quantum dynamics for driven quantum systems interacting with highly structured bosonic baths by combining the tree tensor network (TTN) decomposition scheme with the bexcitonic generalization of the numerically exact hierarchical equations of motion (HEOM). The method yields a series of quantum master equations for all core tensors in the TTN that efficiently and accurately capture the open quantum dynamics for non-Markovian environments to all orders in the system–bath interaction. These master equations are constructed based on the time-dependent Dirac–Frenkel variational principle, which isolates the optimal dynamics for the core tensors given the TTN ansatz. The dynamics converges to the HEOM when increasing the rank of the core tensors, a limit in which the TTN ansatz becomes exact. We introduce TENSO, tensor equations for non-Markovian structured open systems, as a general-purpose Python code to propagate the TTN-HEOM dynamics. We implement three general propagators for the coupled master equations: two fixed-rank methods that require a constant memory footprint during the dynamics and one adaptive-rank method with a variable memory footprint controlled by the target level of computational error. We exemplify the utility of these methods by simulating a two-level system coupled to a structured bath containing one Drude–Lorentz component and eight Brownian oscillators, which is beyond what can presently be computed using the standard HEOM. Our results show that the TTN-HEOM is capable of simulating both dephasing and relaxation dynamics of driven quantum systems interacting with structured baths, even those of chemical complexity, with an affordable computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Should Sea-Ice Modeling Tools Designed for Climate Research Be Used for Short-Term Forecasting?

In theory, the same sea-ice models could be used for both research and operations, but in practice, differences in scientific and software requirements and computational and human resources complicate the matter. Although sea-ice modeling tools developed for climate studies and other research applications produce output of interest to operational forecast users, such as ice motion, convergence, and internal ice pressure, the relevant spatial and temporal scales may not be sufficiently resolved. For instance, sea-ice research codes are typically run with horizontal resolution of more than 3 km, while mariners need information on scales less than 300 m. Certain sea-ice processes and coupled feedbacks that are critical to simulating the Earth system may not be relevant on these scales; and therefore, the most important model upgrades for improving sea-ice predictions might be made in the atmosphere and ocean components of coupled models or in their coupling mechanisms, rather than in the sea-ice model itself. This paper discusses some of the challenges in applying sea-ice modeling tools developed for research purposes for operational forecasting on short time scales, and highlights promising new directions in sea-ice modeling.

54 ENVIRONMENTAL SCIENCES↗

Preparation for quantum simulation of the ( 1 + 1 ) -dimensional O(3) nonlinear σ model using cold atoms

The 1+1D O(3) non-linear σ-model is a model system for future quantum lattice simulations of other asymptotically-free theories, such as non-Abelian gauge theories. We find that utilizing dimensional reduction can make efficient use of two-dimensional layouts presently available on cold atom quantum simulators. A new definition of the renormalized coupling is introduced, which is applicable to systems with open boundary conditions and can be measured using analog quantum simulators. Monte Carlo and tensor network calculations are performed to determine the quantum resources required to reproduce perturbative short-distance observables. In particular, we show that a rectangular array of 48 Rydberg atoms with existing quantum hardware capabilities should be able to adiabatically prepare low-energy states of the perturbatively-matched theory. Furthermore, these states can then be used to simulate non-perturbative observables in the continuum limit that lie beyond the reach of classical computers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fully gauge-fixed SU(2) Hamiltonian for quantum simulations

Here, we demonstrate how to construct a fully gauge-fixed lattice Hamiltonian for a pure SU(2) gauge theory. Our work extends upon previous work, where a formulation of an SU(2) lattice gauge theory was developed that is efficient to simulate at all values of the gauge coupling. That formulation utilized maximal-tree gauge, where all local gauge symmetries are fixed and a residual global gauge symmetry remains. By using the geometric picture of an SU(2) lattice gauge theory as a system of rotating rods, we demonstrate how to fix the remaining global gauge symmetry. In particular, the quantum numbers associated with total charge can be isolated by rotating between the lab and body frames using the three Euler angles. The Hilbert space in this new “sequestered” basis partitions cleanly into sectors with differing total angular momentum, which makes gauge-fixing to a particular total charge sector trivial, particularly for the charge-zero sector. In addition to this sequestered basis inheriting the property of being efficient at all values of the coupling, we show that, despite the global nature of the final gauge-fixing procedure, this Hamiltonian can be simulated using quantum resources scaling only polynomially with the lattice volume.

Lattice gauge theory↗

Direct numerical simulations of a reacting mixing layer with chemical heat release

In order to study the coupling between chemical heat release and fluid dynamics, direct numerical simulations of a chemically reacting mixing layer with heat release are performed. The fully compressible equations as well as an approximate set of equations that is asymptotically valid for low-Mach-number flows are treated. These latter equations have the computational advantage that high-frequency acoustic waves have been filtered out, allowing much larger time steps to be taken in the numerical solution procedure. A detailed derivation of these equations along with an outline of the numerical solution technique is given. Simulation results indicate that the rate of chemical product formed, the thickness of the mixing layer, and the amount of mass entrained into the layer all decrease with increasing rates of heat release.

Mcmurtry, P. A.↗

Evaluating LSM-Based Water Budgets Over a West African Basin Assisted with a River Routing Scheme

Within the framework of the African Monsoon Multidisciplinary Analysis (AMMA) Land Surface Model Intercomparison Project phase 2 (ALMIP-2), this study evaluates the water balance simulated by the Interactions between Soil, Biosphere, and Atmosphere (ISBA) over the upper Oum River basin, in Benin, using a mesoscale river routing scheme (RRS). The RRS is based on the nonlinear Muskingum Cunge method coupled with two linear reservoirs that simulate the time delay of both surface runoff and base flow that are produced by land surface models. On the basis of the evidence of a deep water-table recharge in that region,a reservoir representing the deep-water infiltration (DWI) is introduced. The hydrological processes of the basin are simulated for the 2005-08 AMMA field campaign period during which rainfall and stream flow data were intensively collected over the study area. Optimal RRS parameter sets were determined for three optimization experiments that were performed using daily stream flow at five gauges within the basin. Results demonstrate that the RRS simulates stream flow at all gauges with relative errors varying from -22% to 3% and Nash-Sutcliffe coefficients varying from 0.62 to 0.90. DWI varies from 24% to 67% of the base flow as a function of the sub-basin. The relatively simple reservoir DWI approach is quite robust, and further improvements would likely necessitate more complex solutions (e.g., considering seasonality and soil type in ISBA); thus, such modifications are recommended for future studies. Although the evaluation shows that the simulated stream flows are generally satisfactory, further field investigations are necessary to confirm some of the model assumptions.

Multidisciplinary↗

Augmenting a Simulation Campaign for Hybrid Computer Model and Field Data Experiments

The Kennedy and O’Hagan (KOH) calibration framework uses coupled Gaussian processes (GPs) to meta-model an expensive simulator (first GP), tune its “knobs” (calibration inputs) to best match observations from a real physical/field experiment and correct for any modeling bias (second GP) when predicting under new field conditions (design inputs). There are well-established methods for placement of design inputs for data-efficient planning of a simulation campaign in isolation, that is, without field data: space-filling, or via criterion like minimum integrated mean-squared prediction error (IMSPE). Analogues within the coupled GP KOH framework are mostly absent from the literature. Here, in this study, we derive a closed form IMSPE criterion for sequentially acquiring new simulator data for KOH. We illustrate how acquisitions space-fill in design space, but concentrate in calibration space. Closed form IMSPE precipitates a closed-form gradient for efficient numerical optimization. We demonstrate that our KOH-IMSPE strategy leads to a more efficient simulation campaign on benchmark problems, and conclude with a showcase on an application to equilibrium concentrations of rare earth elements for a liquid–liquid extraction reaction.

97 MATHEMATICS AND COMPUTING↗