Generic Semi-Analytical Simulation Tool (GenSAS)
A generic simulation tool for power grid and general dynamic systems based on semi-analytical solutions (SAS) technology.
SEARCH · Engineering Papers
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.
A generic simulation tool for power grid and general dynamic systems based on semi-analytical solutions (SAS) technology.
This study introduces an innovative approach to modeling meandering river morphology, integrating and investigating the effects of geometric characteristics and vegetation-induced channel coarsening. The developed comprehensive framework combines several advanced techniques: Genetic Programming for refining the scour factor of transverse bed slope, a Leaf Area Index (LAI)-enhanced analytical model for quantifying vegetative flow resistance, and an upstream-weighted moving average method for efficient approximation of the convolution integral in meander migration calculations. The model is validated against both an idealized Kinoshita meander and a natural bend of the Tumen River (China) in equilibrium, demonstrating its robustness across diverse scales and conditions. The model's ability to simulate the long-term evolution, including cutoff events, provides valuable insight for river management strategies. The current findings demonstrate that channel geometry, particularly width-to-depth ratio, plays a dominant role in meander evolution, with wider channels prone to more complex and rapid morphological changes. Vegetation effects are most pronounced in channels with moderate width-to-depth ratios, where they can significantly influence migration rates and bed topography. A combination of channel widening and deepening, coupled with strategic vegetation management, can effectively enhance navigability while maintaining channel stability in the studied Tumen River reach. Sensitivity analyses highlight the complex interplay between hydraulic conditions, sediment characteristics, and vegetation in shaping river morphology. This research advances understanding of the multifaceted nature of meandering river systems and offers practical tools for informed decision-making in river engineering and environmental management, particularly in the context of climate change and increasing anthropogenic pressures on fluvial ecosystems.
This paper presents a GPU-accelerated simulation package, TRED, for next-generation neutrino detectors with pixelated charge readout, leveraging community-driven software ecosystems to ensure adaptability and extensibility. We introduce two generic contributions: (i) an effective-charge representation based on Gaussian quadrature rules, in which the linear- interpolation factors for the field response inside each voxel are absorbed into the effective charge, and (ii) a sparse, block- binned tensor representation that enables efficient FFT-based computation of induced signals on readout electrodes for sparsely activated detector volumes. The former captures structure inside a voxel without dense sampling, while the latter achieves low memory usage and scalable runtime, as demonstrated in bench- mark studies. The underlying data representation is applicable to large-scale detectors and to other computational problems involving sparse activity.
Analytic studies and two-dimensional “clean” radiation-magnetohydrodynamic (rad-MHD) simulations employing dynamical similarity driver-target scaling prescriptions [Ruiz et al., Phys. Plasmas 30, 032708 (2023)] suggest that Magnetized Liner Inertial Fusion (MagLIF) target implosions can scale to > 10 MJ DT fusion yields when peak drive current is increased beyond 60 MA. We present results from three-dimensional (3D) rad-MHD simulations of similarity-scaled MagLIF target implosions at peak drive currents ranging from 15 to 40 MA. Simulations in this study suggest that magneto-Rayleigh–Taylor instability (MRTI) growth and feedthrough to the fuel region are more severe at higher drive current scales, which reduces the fusion yield compared to prior analytic and 2D clean simulation predictions. In contrast to standard MagLIF, simulations of current-scaled MagLIF target implosions driven by a dynamic screw pinch (DSP) demonstrate reduced MRTI feedthrough and greater fuel magnetization, resulting in improved thermonuclear performance and enhanced performance scaling with peak drive current. DSP drive enables additional scaling of the liner mass to increase liner radius but maintain implosion time, resulting in higher implosion velocities at the expense of increased susceptibility to MRTI. We present a current- and mass-scaled simulated DSP-MagLIF target implosion at the ∼ 40 MA peak current level that produces ignition scale performance, demonstrating a burn-averaged Lawson ignition parameter above unity and DT fusion yield above 1 MJ.
An analytical and empirical-based 1-D, non-isothermal, steady-state model for anion exchange membrane fuel cell capable of capturing two-phase phenomena is presented in this study. Coupled multi-physics including mass and charge transport, electrochemical reactions, heat transfer, and two-phase water transport are considered in the model and the simulated results are compared to experimental data. To better represent actual material properties and localized conditions, the model applies multilayer discretization in the gas diffusion electrode to enhance prediction accuracy. The model successfully predicts the baseline performance at 70 °C, 131 kPa abs., 92% RH with pure H 2 /O 2 gas as well as the limiting current at 10% H 2 . The robust simulation approach allows for simplistic and accurate estimation of cell performance without the complications of applying two-phase parameters and expensive computational need for numerical models. In addition, the results from the sensitivity studies of material properties and operating conditions provide valuable insights on water management strategies and optimal component design for advancing anion exchange membrane fuel cell technology.
We present an ongoing work in which a surrogate model is being developed to reproduce the response dynamics of the third-integer resonant extraction process in the Delivery Ring (DR) at Fermilab. This is in pursuit of smoothly extracting circulating beam to the Mu2e Experiment’s production target, whereby the goal is to extract a uniform slice of the circulating 1e12 protons in the DR over 25,000 turns (43 ms). The DR contains 3 harmonic sextupoles that excite a third-integer resonance and three fast, tune-ramping quadrupole magnets that drive the horizontal tune towards the 29/3 resonance. In our initial work, the surrogate model trains on a semi-analytical simulation provided in the same format as live data. Using Reinforcement Learning (and other potential ML methods), the trained surrogate acts as the “environment” in which a simple ML control agent could learn to dynamically adjust the quadrupole ramp at 430 break points within the 43 microsecond spill window. The controller will be hosted on a dedicated Arria 10 FPGA. In this work, we report the accuracy and fidelity of the surrogate model in comparison to the response dynamics of the physics simulator.
Not Available
Not provided.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
An improved electrical contact resistance (ECR) model for elastic rough electrode contact is proposed, incorporating the effects of asperity interactions and temperature rise by frictional and joule heating. The analytical simulation results show that the ECR decreases steeply at the beginning of the contact between Al and Cu. However, it becomes stabilized after reaching a specific contact force. It is also found that the longer elapsed sliding contact time, the higher ECR due to the increase in electrical resistivity of electrode materials by the frictional temperature rise at the interface. The effects of surface roughness parameters on ECR are studied through the 3 2 full-factorial design-of-experiment analysis. Based on the two representative roughness parameters, i.e., root-mean-square (rms) roughness and asperity radius, their individual and coupled effects on the saturated ECR are examined. The saturated ECR increases with the rms roughness for a rough machined surface condition, but it is hardly affected by the asperity radius. On the other hand, the saturated ECR increases with both the rms roughness and the asperity radius under a smooth thin film surface condition.
Previous analytical and simulation-based analyses suggest that deeper land surface models are needed to realistically simulate the terrestrial thermal state in climate models, with implications for land-atmosphere interactions. Analytical approaches mainly focused on the subsurface propagation of harmonics such as the annual temperature signal, and a direct comparison with climate-change model output has been elusive. Here, this study addresses the propagation of a harmonic pulse fitted to represent the timescale and amplitude of anthropogenic warming. Its comparison to land model simulations with stepwise increased bottom boundary depth leads to an agreement between the simulation-based and analytical frameworks for long-term climate trends. Any depth increase gradually decreases the relative error in the subsurface thermodynamics, and a minimum depth of 170 m is recommended to simulate the ground climate adequately. The approach provides an accurate estimate of the required land-model depth for climate-change simulations and assesses the relative bias in insufficiently deep land models.
The objective is to better predict the bulk-scale mechanical behavior of porous metals that have over pressurized cavities (e.g., irradiated metals with helium bubbles) by quantifying the complex coupling among cavity aspects (e.g., size distribution, inhomogeneous overpressure values, spatial arrangement) and metal properties (e.g., rate-dependency, crystallographic lattice). This requires up-scaling local mechanical fields from the single crystal scale and will be accomplished using a homogenization approach that combines full-field numerical simulations, analytical formalisms, and physics-informed machine learning to produce symbolically-defined constitutive equations (e.g., gauge functions). These equations will satisfy the objective because they enable computationally efficient predictions that approach the accuracy of computationally expensive full-field numerical simulations, abide by theoretical requirements (e.g., conservation of energy, work conjugacy), and retain the transparency of analytical models.
A slow extraction system is being developed for the Fermilab’s Delivery Ring to deliver protons to the Mu2e experiment. During the extraction, the beam on target experiences small intensity variations owing to many factors. Various adaptive learning algorithms will be employed for beam regulation to achieve the required spill quality. We discuss here preliminary results of the slow and fast regulation algorithms validation through the computer simulations before their implementation in the FPGA. Particle tracking with sextupole resonance was used to determine the fine shape of the spill profile. Fast semi-analytical simulation schemes and Machine Learning models were used to optimize the fast regulation loop.
Abstract The hot-tail mechanism of runaway electron (RE) production (Harvey et al 2000 Phys. Plasmas 7 4590) is the primary source of RE in the case of rapidly cooling tokamak plasma. Quantifying this mechanism is very important as it can provide most of the post-thermal-quench (TQ) current, or a seed current for the secondary source of RE through the avalanche mechanism. An analytic model which omits pitch-angle scattering is often used in literature for estimating the hot-tail RE density (Smith and Verwichte 2008 Phys. Plasmas 15 072502). In the present study, we use the CQL3D bounce-averaged Fokker–Planck code (Harvey and McCoy 1992 Proc. IAEA Technical Committee Meeting on Advances in Simulation and Modeling of Thermonuclear Plasmas p 527) to test the limits of validity of the model. In particular, we examine the cases of Z = 1 and Z = 18 ions, for sets of different initial temperature, density, electric field and the characteristic time of temperature decay. We show that for Z = 1 plasma, the ratio of RE density computed by CQL3D to that estimated from the model is within 0.6–6.0 in studied cases. For the Z = 18 case, this factor is systematically a much smaller number, typically 0.02–0.6. We suggest a simple correction to the model that narrows down the range of this ratio to 0.3–3.8 in all of the cases, including Z = 1 and Z = 18 plasmas.
Workforce planning deals with determining the number of employees and associated skills necessary to meet the future operational needs of an organization. A workforce system consists of six elements: recruitment, attrition, promotion, training, retention, and scheduling. Historically, several workforce modeling and analysis methodologies have been developed to capture these elements. This paper reviews the results of workforce and manpower models published within peer-reviewed literature between 1959 and 2021 to provide an in-depth analysis of current models. The focus of this review is on analytical, simulation, and empirical models found in literature that were collected based on a citation requirement and keyword search criteria. Results demonstrate the trends in workforce modeling research and discuss the common uses of each model type and the advantages/disadvantages related to each model. Based on the common attributes of workforce systems, the discussion focuses on the most frequently used model type for each element and the best use for each model. Lastly, recommendations are made for the development of workforce models that allow the most comprehensive view of the workforce systems of the future.
Operations of multi-reservoir systems are challenged in-part by the interplay of complex physical processes functioning within the watershed. The employment of intelligent systems can be of aid by linking environmental sensing, information technology, data analytics, simulation and decision support to achieve a data-to-decision flow of information. A further challenge is that watershed resources are managed for multiple purposes requiring some level of coordination among numerous resource managers, asset operators and users. System intelligence in this context relies on shared community platforms (data portals, community models), and coordinated communication between decision makers. Opportunities to enrich watershed intelligence has been the subject of a roadmapping exercise for the Department of Energy’s Water Power Technologies Office which has relied on broad stakeholder engagement. Initial phases of engagement involved personal interviews and a series of virtual group meetings, which focused on identifying opportunities to improve the intelligence of the physical infrastructure within our watersheds—examples of feedback include improved sensing of snowpack and runoff, data standards for facilitated data sharing, and better forecasting tools. The latter phase of engagement involved the conduct of a case study in the Upper Colorado River basin where key stakeholders were interviewed to map how their decisions are informed by intelligence from other basin stakeholders. Our presentation will highlight the interdisciplinary flow of information in complex watershed systems and identify physical and institutional opportunities toward the strategic operation of water infrastructure.
Grazing-incidence reflective optics are commonly used in synchrotron radiation and free-electron laser facilities to transport and focus the emitted X-ray beams. To preserve the imaging capability at the diffraction limit, the fabrication of these optics requires precise control of both the residual height and slope errors. However, all the surface figuring methods are height based, lacking the explicit control of surface slopes. Although our preliminary work demonstrated a one-dimensional (1D) slope-based figuring model, its 2D extension is not straightforward. In this study, a novel 2D slope-based figuring method is proposed, which employs an alternating objective optimization on the slopes in the x - and y -directions directly. An analytical simulation revealed that the slope-based method achieved smaller residual slope errors than the height-based method, while the height-based method achieved smaller residual height errors than the slope-based method. Therefore, a hybrid height and slope figuring method was proposed to further enable explicit control of both the height and slopes according to the final mirror specifications. An experiment to finish an elliptical-cylindrical mirror using the hybrid method with ion beam figuring was then performed. Both the residual height and slope errors converged below the specified threshold values, which verified the feasibility and effectiveness of the proposed ideas.