Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “source code”

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,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Manufactured solutions for an electromagnetic slot model

The accurate modeling of electromagnetic penetration is an important topic in computational electromagnetics. Electromagnetic penetration occurs through intentional or inadvertent openings in an otherwise closed electromagnetic scatterer, which prevent the contents from being fully shielded from external fields. To efficiently model electromagnetic penetration, aperture or slot models can be used with surface integral equations to solve Maxwell's equations. A necessary step towards establishing the credibility of these models is to assess the correctness of the implementation of the underlying numerical methods through code verification. Surface integral equations and slot models yield multiple interacting sources of numerical error and other challenges, which render traditional code-verification approaches ineffective. In this paper, we provide approaches to separately measure the numerical errors arising from these different error sources for the method-of-moments implementation of the electric-field integral equation with a slot model. Finally, we demonstrate the effectiveness of these approaches for a variety of cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2001

Code Sharing and Collaboration: Experiences From the Scientist's Expert Assistant Project and Their Relevance to the Virtual Observatory

In the Virtual Observatory (VO), software tools will perform the functions that have traditionally been performed by physical observatories and their instruments. These tools will not be adjuncts to VO functionality but will make up the very core of the VO. Consequently, the tradition of observatory and system independent tools serving a small user base is not valid for the VO. For the VO to succeed, we must improve software collaboration and code sharing between projects and groups. A significant goal of the Scientist's Expert Assistant (SEA) project has been promoting effective collaboration and code sharing among groups. During the past three years, the SEA project has been developing prototypes for new observation planning software tools and strategies. Initially funded by the Next Generation Space Telescope, parts of the SEA code have since been adopted by the Space Telescope Science Institute. SEA has also supplied code for the SIRTF (Space Infrared Telescope Facility) planning tools, and the JSky Open Source Java library. The potential benefits of sharing code are clear. The recipient gains functionality for considerably less cost. The provider gains additional developers working with their code. If enough users groups adopt a set of common code and tools, de facto standards can emerge (as demonstrated by the success of the FITS standard). Code sharing also raises a number of challenges related to the management of the code. In this talk, we will review our experiences with SEA--both successes and failures, and offer some lessons learned that might promote further successes in collaboration and re-use.

Korathkar, Anuradha↗
NASA NTRS2007

The Continuous Intercomparison of Radiation Codes (CIRC): Phase I Cases

CIRC aspires to be the successor to ICRCCM (Intercomparison of Radiation Codes in Climate Models). It is envisioned as an evolving and regularly updated reference source for GCM-type radiative transfer (RT) code evaluation with the principle goal to contribute in the improvement of RT parameterizations. CIRC is jointly endorsed by DOE's Atmospheric Radiation Measurement (ARM) program and the GEWEX Radiation Panel (GRP). CIRC's goal is to provide test cases for which GCM RT algorithms should be performing at their best, i.e, well characterized clear-sky and homogeneous, overcast cloudy cases. What distinguishes CIRC from previous intercomparisons is that its pool of cases is based on observed datasets. The bulk of atmospheric and surface input as well as radiative fluxes come from ARM observations as documented in the Broadband Heating Rate Profile (BBHRP) product. BBHRP also provides reference calculations from AER's RRTM RT algorithms that can be used to select the most optimal set of cases and to provide a first-order estimate of our ability to achieve radiative flux closure given the limitations in our knowledge of the atmospheric state.

Oreopoulos, Lazaros↗
DOE OSTI2024

Code for Experiment in Publication “Sensitivity Analysis in the Presence of Intrinsic Stochasticity for Discrete Fracture Network Simulations”

Following the Open Research requirements for AGU journals, we must release the code used to perform the experiment described in our recent publication, posted at (https://arxiv.org/abs/2312.04722). This code fits a joint emulator to data from a Discrete Fracture Network (DFN) simulation, performed using the open-source software DFNworks (https://dfnworks.lanl.gov/). All code to be released implements existing methods; there are no novel algorithms nor any major innovations to existing software.

Murph, Alexander↗
DOE OSTI2024

Code for Experiment in Publication “Determining the dominant factors for carbon-1 mineralization in three-dimensional fracture networks”

We plan to release the code used to perform the experiment described in our recent publication, entitled “Determining the dominant factors for carbon-1 mineralization in three-dimensional fracture networks.” This code fits a joint emulator to data from several Discrete Fracture Network (DFN) simulations, performed using the open-source software DFNworks (https://dfnworks.lanl.gov/). All code to be released implements existing methods; there are no novel algorithms nor any major innovations to existing software.

Murph, Alexander↗
DOE OSTI2021

Seismic Source Modeling Software Enhancements (FY21)

Seismic source modeling allows researchers both to simulate how a source that induces seismic waves interacts with the Earth to produce observed seismograms and, inversely, to infer what the time histories, sizes, and force distributions were for a seismic source given observed seismograms. In this report, we discuss improvements made in FY21 to our software as applies to both the forward and inverse seismic source modeling problems. For the forward portion of the problem, we have added the ability to use full 3-D nonlinear simulations by implementing 3-D time varying boundary conditions within Sandia’s linear seismic code Parelasti. Secondly, on the inverse source modeling side, we have developed software that allows us to invert seismic gradiometer-derived observations in conjunction with standard translational motion seismic data to infer properties of the source that may improve characterization in certain circumstances. First, we describe the basic theory behind each software enhancement and then demonstrate the software in action with some simple examples.

58 GEOSCIENCES↗
NASA NTRS1992

Entropy-Based Bounds On Redundancies Of Huffman Codes

Report presents extension of theory of redundancy of binary prefix code of Huffman type which includes derivation of variety of bounds expressed in terms of entropy of source and size of alphabet. Recent developments yielded bounds on redundancy of Huffman code in terms of probabilities of various components in source alphabet. In practice, redundancies of optimal prefix codes often closer to 0 than to 1.

Smyth, Padhraic J.↗
DOE OSTI2013

Explicit estimation of higher order modes in fission source distribution of Monte-Carlo calculation

Magnitude of higher order modes in fission source distribution of a multi-group Monte-Carlo calculation is estimated using the orthogonal property of forward and adjoint fission source distributions. Calculation capability of the forward and adjoint fission source distributions for fundamental and higher order modes are implemented in the AEGIS code, which is a two-dimensional transport code based on the method of characteristics. With the calculation results of the AEGIS code, magnitudes of the first to fifth higher order modes in fission source distribution obtained by the multi-group Monte-Carlo code GMVP are estimated. There are two contributions in the present study - (1) establishment of a surrogate model, which represents convergence of fission source distribution taking into account the inherent statistical 'noise' of higher order modes of Monte-Carlo calculations and (2) independent confirmation of the estimated dominance ratio in a Monte-Carlo calculation. The surrogate model would contribute to studies of the inter-cycle correlation and estimation of sufficient number of inactive/active cycles. (authors)

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2005

Computer Simulation of the VASIMR Engine

The goal of this project is to develop a magneto-hydrodynamic (MHD) computer code for simulation of the VASIMR engine. This code is designed be easy to modify and use. We achieve this using the Cactus framework, a system originally developed for research in numerical relativity. Since its release, Cactus has become an extremely powerful and flexible open source framework. The development of the code will be done in stages, starting with a basic fluid dynamic simulation and working towards a more complex MHD code. Once developed, this code can be used by students and researchers in order to further test and improve the VASIMR engine.

Garrison, David↗
DOE OSTI2024

3D modeling of deep borehole electromagnetic measurements with energized casing source for fracture mapping at the Utah Frontier Observatory for Research in Geothermal Energy

Here, we present a 3D numerical modelling analysis evaluating the deployment of a borehole electromagnetic measurement tool to detect and image a stimulated zone at the Utah Frontier Observatory for Research in Geothermal Energy geothermal site. As the depth to the geothermal reservoir is several kilometres and the size of the stimulated zone is limited to several 100 m, surface-based controlled-source electromagnetic measurements lack the sensitivity for detecting changes in electrical resistivity caused by the stimulation. To overcome the limitation, the study evaluates the feasibility of using a three-component borehole magnetic receiver system at the Frontier Observatory for Research in Geothermal Energy site. To provide sufficient currents inside and around the enhanced geothermal reservoir, we use an injection well as an energized casing source. To efficiently simulate energizing the injection well in a realistic 3D resistivity model, we introduce a novel modelling workflow that leverages the strengths of both 3D cylindrical-mesh-based electromagnetic modelling code and 3D tetrahedral-mesh-based electromagnetic modelling code. The former is particularly well-suited for modelling hollow cylindrical objects like casings, whereas the latter excels at representing more complex 3D geological structures. In this workflow, our initial step involves computing current densities along a vertical steel-cased well using a 3D cylindrical electromagnetic modelling code. Subsequently, we distribute a series of equivalent current sources along the well's trajectory within a complex 3D resistivity model. We then discretize this model using a tetrahedral mesh and simulate the borehole electromagnetic responses excited by the casing source using a 3D finite-element electromagnetic code. This multi-step approach enables us to simulate 3D casing source electromagnetic responses within a complex 3D resistivity model, without the need for explicit discretization of the well using an excessive number of fine cells. We discuss the applicability and limitations of this proposed workflow within an electromagnetic modelling scenario where an energized well is deviated, such as at the Frontier Observatory for Research in Geothermal Energy site. Using the workflow, we demonstrate that the combined use of the energized casing source and the borehole electromagnetic receiver system offer measurable magnetic field amplitudes and sensitivity to the deep localized stimulated zone. The measurements can also distinguish between parallel-fracture anisotropic reservoirs and isotropic cases, providing valuable insights into the fracture system of the stimulated zone. Besides the magnetic field measurements, vertical electric field measurements in the open well sections are also highly sensitive to the stimulated zone and can be used as additional data for detecting and imaging the target. We can also acquire additional multiple-source data by grounding the surface electrode at various locations and repeating borehole electromagnetic measurements. This approach can increase the number of monitoring data by several factors, providing a more comprehensive dataset for analysing the deep-localized stimulated zone. The numerical analysis indicates that it is feasible to use the combination of the energized casing and downhole electromagnetic measurements in monitoring localized stimulated zone at large depths.

58 GEOSCIENCES↗
DOE OSTI2016

A time-dependent neutron transport method of characteristics formulation with time derivative propagation

A new time-dependent Method of Characteristics (MOC) formulation for nuclear reactor kinetics was developed utilizing angular flux time-derivative propagation. This method avoids the requirement of storing the angular flux at previous points in time to represent a discretized time derivative; instead, an equation for the angular flux time derivative along 1D spatial characteristics is derived and solved concurrently with the 1D transport characteristic equation. This approach allows the angular flux time derivative to be recast principally in terms of the neutron source time derivatives, which are approximated to high-order accuracy using the backward differentiation formula (BDF). This approach, called Source Derivative Propagation (SDP), drastically reduces the memory requirements of time-dependent MOC relative to methods that require storing the angular flux. An SDP method was developed for 2D and 3D applications and implemented in the computer code DeCART in 2D. DeCART was used to model two reactor transient benchmarks: a modified TWIGL problem and a C5G7 transient. The SDP method accurately and efficiently replicated the solution of the conventional time-dependent MOC method using two orders of magnitude less memory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

EMPulse, a new 3-D simulation code for EMP formation and propagation

EMPulse is a comprehensive and modern 3-D simulation code for electro-magnetic pulse (EMP) formation and propagation studies, being developed at LLNL under a Laboratory Directed R&D (LDRD) initiative that was launched in the fall of 2014. Furthermore, EMPulse builds upon the open-source Warp particle-in-cell code framework which was developed by members of this team and collaborators at other institutions. Our goal of this endeavor is a new tool, using methods different from those in common use, and enabling the detailed and self-consistent study of multidimensional effects in geometries that have typically been treated only approximately. We describe the underlying Warp framework and its capabilities; the models that have been or are being newly developed and incorporated into EMPulse, tests of these models, and number of approaches to increased computational efficiency that we are studying.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Modifications of impurity transport and divertor sources by lithium wall conditioning in the National Spherical Torus Experiment

In the National Spherical Torus Experiment (NSTX), lithium coatings are evaporated on graphite plasma facing components (PFCs) for wall conditioning. In lithium-conditioned H-mode discharges, carbon accumulation is observed with core concentrations <10%, leading to a lack of density control, while lithium ions have concentrations <0.1%. In this thesis, modifications of carbon and lithium divertor sources as well as scrape-off layer (SOL) and core transport due to lithium conditioning are studied. Spectroscopic impurity influxes (measured by filtered cameras) and 2D multi-fluid edge transport simulations via the UEDGE code are employed to study divertor impurity sources and SOL transport, respectively. Core transport of carbon and lithium is analyzed using the impurity transport code MIST and the neoclassical transport codes NEO and NCLASS. A reduction of the carbon sputtering yield in the lower divertor is observed with lithium evaporation. However, weaker divertor impurity retention resulting from reduced recycling (inferred from UEDGE simulations) and the possible importance of wall sources can counteract this reduction in divertor carbon influxes. The suppression of edge-localized-modes (ELMs) is the primary cause of the increased carbon inventories in lithium-conditioned discharges, leading to lack of density control. Deviations from neoclassical predictions for carbon transport are observed at the pedestal top in lithium-conditioned discharges, indicating the presence of anomalous outward convection. While the lithium sputtering yield from lithium-coated graphite in the divertor is consistent with physical and temperature-enhanced sputtering, a strong reduction in ionized lithium influxes is observed, possibly due to prompt re-deposition. The different poloidal source distribution and the stronger divertor retention for lithium (inferred from UEDGE simulations) contribute to a lower edge lithium source with respect to carbon. The latter is due to the shorter lithium ionization mean-free-path and weaker classical parallel forces. The high neoclassical lithium diffusivity, due to collisions with carbon ions, is partially responsible for the low lithium core densities. Furthermore, the higher neoclassical levels make lithium transport closer to neoclassical expectations. The NSTX-Upgrade baseline scenario includes lithium coatings on graphite PFCs for deuterium particle control. While the remarkably good lithium divertor retention makes it an attractive choice, the suppression of ELMs leads to carbon accumulation and therefore to the need to develop impurity mitigation techniques.

impurity trasport↗
DOE OSTI2022

Combined Thermal-Water Microgrid Model

SAND2022-3874 O The Combined Thermal-Water Microgrid Model is a set of three models combined with controls algorithms for a Water Thermal Storage (WTS), Ground Source Heat Pump (GSHP), and Thermal Bridge Loop (TBL) that serves a residential neighborhood. The model allows a parametric investigation of volume of water storage versus number of ground heat exchanger bore-holes versus seasonal energy efficiency measures. The code is being released as open source so that the entire study can be replicated per an open access journal publication. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Villa, Daniel↗
DOE OSTI2014

A 2D Particle in Cell model for ion extraction and focusing in electrostatic accelerators

Negative ions are fundamental to produce intense and high energy neutral beams used to heat the plasma in fusion devices. The processes regulating the ion extraction involve the formation of a sheath on a scale comparable to the Debye length of the plasma. On the other hand, the ion acceleration as a beam is obtained on distances greater than λ{sub D}. The paper presents a model for both the phases of ion extraction and acceleration of the ions and its implementation in a numerical code. The space charge of particles is deposited following usual Particle in Cell codes technique, while the field is solved with finite element methods. Some hypotheses on the beam plasma transition are described, allowing to model both regions at the same time. The code was tested with the geometry of the NIO1 negative ions source, and the results are compared with existing ray tracing codes and discussed.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

Modelling low Mach number stellar hydrodynamics with MAESTROeX

Modelling long-time convective flows in the interiors of stars is extremely challenging using conventional compressible hydrodynamics codes due to the acoustic timestep limitation. Many of these flows are in the low Mach number regime, which allows us to exploit the relationship between acoustic and advective time scales to develop a more computationally efficient approach. MAESTROeX is an open source low Mach number stellar hydrodynamics code that allows much larger timesteps to be taken, therefore enabling systems to be modelled for much longer periods of time. This is particularly important for the problem of convection in the cores of rotating massive stars prior to core collapse. To fully capture the dynamics, it is necessary to model these systems in three dimensions at high resolution over many rotational periods. We present an overview of MAESTROeX's current capabilities, describe ongoing work to incorporate the effects of rotation and discuss how we are optimising the code to run on GPUs.

Harpole, A.↗
DOE OSTI2019

Gauge structures: From stabilizer codes to continuum models

Stabilizer codes are a powerful method for implementing fault-tolerant quantum memory and in the case of topological codes, they form useful models for topological phases of matter. In this paper, we discuss the theory of stabilizer codes as a discrete version of a linear gauge structure, a concept we introduce here. A linear gauge structure captures all the familiar ingredients of U(1) gauge theory including a generalization of charge conservation, Maxwell equations and topological terms. Using this connection, we prove some important results for stabilizer codes which can be used to characterize their error-correction properties. However, this perspective does not depend on any particular Hamiltonian or Lagrangian, which as a consequence is agnostic to the source of the gauge redundancy. Based upon the connection to stabilizer codes, we view the source of the gauge redundancy as an ambiguity in the tensor product structure of the Hilbert space. That is, the gauge provides a set of equivalent but arbitrary ways to factorize the Hilbert space. We apply this formalism to the d=2 and 3 toric code as well as the paradigm fracton models, X-Cube and Haah’s cubic code. From this perspective, we are able to map all these models to some continuum version which captures all or some of the fractonic and topological features of their parent discrete models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1994

Bit-wise arithmetic coding for data compression

This article examines the problem of compressing a uniformly quantized independent and identically distributed (IID) source. We present a new compression technique, bit-wise arithmetic coding, that assigns fixed-length codewords to the quantizer output and uses arithmetic coding to compress the codewords, treating the codeword bits as independent. We examine the performance of this method and evaluate the overhead required when used block-adaptively. Simulation results are presented for Gaussian and Laplacian sources. This new technique could be used as the entropy coder in a transform or subband coding system.

Kiely, A. B.↗