Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solution regularity”

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 469 records · Page 26

Regular and in-plane skyrmions and antiskyrmions from boundary instabilities

We formulate a theory of skyrmion and antiskyrmion generation using magnetic field and charge current pulses. We show that the topological defect can be created at an edge of a system with Dzyaloshinskii-Moriya interaction (DMI) as well as at a boundary between regions with different DMI. We consider both perpendicular and in-plane (also known as magnetic bimerons) versions of skyrmions and antiskyrmions. We show that the magnetization twist in the vicinity of an edge or a boundary is described by a kink solution, the presence of which can instigate the generation of topological defects. We study the collective excitations of magnetization analytically and numerically, and demonstrate that under application of magnetic field and charge current pulses the magnon modes localized near boundaries can develop instabilities leading to the formation of skyrmions or antiskyrmions. As a result, due to the skyrmion and antiskyrmion Hall effects, a properly chosen current direction can push the topological defects away from the boundary, thus facilitating their generation.

36 MATERIALS SCIENCE↗

Verification of the REBUS Software

Ongoing design activities at Argonne National Laboratory are requiring a thorough verification of the Argonne Reactor Computation codes be performed. REBUS is central to this system. The driver for this effort requires the Triangular-Z and hexagonal-Z core geometry options of REBUS to be verified. Previous work identified the REBUS features required to be verified to support current design activities, features of which are generally applicable to hexagonal-Z fast reactor designs. The scope of this verification effort includes verifying REBUS’s ability to correctly intepret the user input model, verifying that the features identified yield the intended results, and verifying the correctness of the REBUS output tables. The REBUS software verification relies heavily upon the accuracy of the embedded DIF3D software, the verification of which was completed and documented elsewhere. Given that DIF3D produces an accurate solution, the primary focus of the verification in the REBUS software is to ensure that it properly uses the DIF3D solution and that the depletion system (Bateman equations) are correctly implemented. This manuscript reiterates the verification tasks and displays results with respect to the features needed for current design activities. Analytic solutions of the Batemen equations are displayed and the results calculated with REBUS are displayed demonstrating the accuracy. Since coupled Bateman and neutron diffusion/transport solutions are extremely difficult to obtain, much of the focus is placed on how REBUS uses a given DIF3D solution assuming the accuracy of the DIF3D solution. The verification effort identified no issues that are debilitating or otherwise impactful to the design usage of REBUS, and thus REBUS version 11.0, release 3012 is considered verified. It is important to note that several outputs of REBUS are identified to be inaccurate, such as burnup in MWD/MT. Most of the relevant ones for VTR are generally accurate with 10-20% errors which is not impactful as all regular REBUS users are aware of this issue and know how to hand calculate the results. The REBUS manual further makes it clear that these values are consistent with the methodology being used by REBUS and thus the “errors” are more of an inconsistent definition with respect to what a user would expect given a definition in literature. Other issues that were identified included unclear documentation and software bugs all of which were inconsequential to the final results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Explicit-implicit staggered procedure for multibody dynamics analysis

An explicit-implicit staggered time-integration procedure is presented for the solution of multibody dynamical equations involving large rotations and constraints. The algorithms adopts a two-stage modification of the central difference algorithm for integrating the translational coordinates and the angular velocity vector, and the midpoint implicit algorithm to solve the kinematical relation in terms of the Euler parameters for updating the angular orientations. The Lagrange multipliers to enforce the system constraints are obtained by implicitly integrating a parabolically regularized differential equation for the multipliers. The performance of the present procedure has been evaluated to applying the procedure to solve several sample problems. The results indicate that the procedure is robust in dealing with a variety of constraints and spatial kinematic motions, hence it is recommended for applications to general multibody dynamics analyses.

Park, K. C.↗

Impedance Eduction in a Duct Using the Linearized Euler Equations

An impedance eduction technique that is based upon a numerical solution to the linearized Euler equations coupled with the minimization of an objective function is presented. The linearized Euler equations are solved numerically using a cubic finite element method and the minimum of the objective function is obtained via a gradient based optimizer. The mean flow is allowed to have a gradient in two coordinate directions and the impedance of the liner may vary arbitrarily over the liner surface. A new regularized wall impedance boundary condition presented by Rienstra which contains the effects of the mean boundary layer thickness is implemented. Impedances are educed for a perforate over honeycomb liner using test data acquired in the Langley Grazing Flow Impedance Tube. The impedance of the liner is educed both with an upstream source and a downstream source. The method is validated by comparing the educed impedance to a benchmark method. The primary conclusion of the study is that the impedances educed for upstream and downstream sources are not well matched and that the primary effect of the boundary thickness is to reduce the resistance of the liner.

Watson, Willie R.↗

Use of satellite telecommunications for telemedicine

The use of space telecommunications facilities in the interests of medicine should establish a reliable and timely connection between the specialists involved and access for each of them to information sources the need. Since a certain global-scale program should be implemented for this goal to be achieved, one of the top-priority tasks is to determine appropriate orbital facilities for telecommunications as well as appropriate areas and countries for them, outlining the territories where portable data receive/transport records can be placed and stations that are cheap and simple enough to operate can be stored. Eventually the satellite telecommunications system should encourage health workers to unite in a world community on the basis of a regular and efficient information exchange with the goal of protecting the health of man and the peace of Earth. This goal could only be achieved by gradually removing differences in the domain of medical education and information in different regions over the globe, by providing a possibility of rapid and coordinated access to the information on medical problems of interest on the global scale, as well as by coordinating the discussion and solution of all problems associated with the development and progress of the satellite communications system.

Tamkovich, Gennady↗

Results from the GPS Flight Experiment on the High Earth Orbit AMSAT OSCAR-40 Spacecraft

A GPS receiver flying on the High Earth Orbit (HEO) AMSAT-OSCAR 40 (AO-40) spacecraft has been returning GPS observations from high above the altitude of the GPS constellation. AO-40, an amateur radio satellite launched November 16, 2000, is currently in a low inclination, 1000 by 59000 lan altitude orbit. This low-cost experiment utilizes a mid 1990's era, 6-channel, CIA code receiver configured with high gain receiving antennas for tracking above the GPS constellation. The receiver has performed well, despite operating significantly outside of its original design environment. It has regularly returned GPS observations from points all around the orbit, with over ten weeks of GPS tracking data collected to date. Signal to noise levels as high as 48 B-Hz have been recorded near apogee, when the spacecraft was at an altitude of close to 60000 km. GPS side lobe signals have been tracked on several occasions, primarily from Block IIR GPS satellites. Although the receiver has not computed a solution in real-time, point solutions have been computed on the ground using simultaneous measurements from four satellites. This experiment has provided important experience dealing with the many challenges inherent to GPS tracking at high altitudes, and the measurements returned are providing valuable information about the characteristics of GPS signals available for future HE0 users.

Bauer, Frank↗

Equitable Urban Electric Vehicle Charging: Feasibility and Benefits of Streetlight Charging in Kansas City Right-of-Way

With an increasing global emphasis on sustainability, electric vehicles (EVs) play a crucial role in reducing urban pollution and carbon emissions. For EVs to be widely adopted, equitable and convenient access to charging infrastructure is essential. Equity in this paper refers to the proactive engagement with the community to ensure that the benefits of streetlight charging are distributed equitably across diverse neighborhoods in Kansas City, providing fair charging opportunities and resources to all community members. This research explores the potential of utilizing streetlights—ubiquitous elements of urban electrical infrastructure—as low-cost, equitable EV charging solutions. Compared with conventional chargers, streetlight charging offers several notable advantages, including proximity to roadways, potential boosts to the local economy, and easier usage due to city ownership. These chargers also leverage existing power setups to minimize costs and maximize efficiency by reutilizing existing structures. In this manuscript, we introduce a systematic framework to develop, analyze, and evaluate a scalable and cost-effective streetlight charging solution. Initially, we employ a two-tiered site selection framework that accounts for both charging demand and equity considerations to identify optimal locations for streetlight chargers. We then assess the feasibility of these chargers through observational data from 23 units installed in Kansas City, Missouri, comparing their performance with traditional chargers. Here, our evaluation extends to the environmental impact, comparing reduced gasoline consumption and greenhouse gas (GHG) emission reductions between streetlights and regular charging stations. Our findings reveal that streetlight charging stations offer faster charging speeds than conventional stations, validating their viability for EV charging. Additionally, the duration of stay at streetlight chargers tends to be shorter, likely due to parking constraints or associated costs. Furthermore, streetlight-powered chargers demonstrate considerable environmental advantages, achieving 11.94% greater avoided gasoline consumption and 11.24% higher GHG reductions compared with regular charging stations, highlighting their substantial environmental benefits.

42 ENGINEERING↗

Clean and High-Efficiency Natural Gas Furnace with Advanced Acidic Gas Trap Technology

Ccurrent natural gas condensing furnaces generate substantial acidic condensate and SOx/NOx/CO/HC/methane emissions, causing long-term environment issues related to air, water and soil. Acidic condensate occurs when flue gases within the furnaces cool down below the dew points of acidic gases, which upon condensation, combine with water vapor to produce acidic solutions. The acidic condensate causes serious corrosion and fouling problems in natural gas condensing furnaces. This paper presents a novel solution to use an advanced adsorption technology based on monolithic acidic gas trap (AGT) adsorbers for SOx trapping, NOx redox, and formic acid/CO/HC/methane oxidation, enabling a new natural gas furnace with utra-clean flue gas and neutral condensate. When the AGT adsorber is appropriately integrated with natural gas furnaces, the acidic gases of SOx, NOx and formic gas/CO/HC will be continuously adsorbed, redoxed or oxidized, respectively, from the flue gas during a regular heating season, except for occasional regeneration activities for SOx trapping. In the paper, we present results from a prototype furnace employing the novel AGT component, which enables a pH of the condensate =7, NOx emissions of 1-2 nanograms/joule and an annual fuel utilization efficiency (AFUE) of 96%, providing a cost-effective technology for condensing furnaces that is eco-friendly to the environment.

Gao, Zhiming↗

Grid Resiliency using Hydro-based Microgrids During Wildfires

The primary objective of this project is to assess the techno-economic viability of hydropower resources with supplementary investments, including microgrid controllers, storage solutions, load banks, additional generation like diesel generators, communication links, and transmission infrastructure. The aim is to enable these systems to operate autonomously (in islanded mode) and ensure grid resiliency during and after wildfire events. The project seeks to establish a comprehensive framework for setting up hydropower-based microgrids, from risk estimation to final investment decisions. These microgrids can provide critical support during emergencies and offer additional grid stability and revenue opportunities during regular operations and other system contingencies.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Investigating Biases in VLBI Clock Functions and Position Solutions

A systematic bias causing a statistically significant drift between group delay and phase delay measurements in geodetic Very Long Baseline Interferometry (VLBI)has long been present. The cause of this drift remains unknown, although it is readily reproduced in VLBI experiments including a short baseline such as those conducted with the twin 13-meter telescopes at the Wettzell observatory, given the station names WETTZ13N and WETTZELL. This statistical incompatibility has been an obstacle in adopting phase delays in routine VLBI processing. The aim of this study is to provide additional evidence to aid in identifying the source of this systematic bias. A series of four regular VLBI experiments are processed with both group delays and phase delays, and post fit residuals, estimated positions, and clock functions are presented with and without phase calibration applied. These results are also computed for individual intermediate frequency bands.

Geodesy↗

A Green’s function fast multipole method for computation of micromechanical fields in heterogeneous materials

Computation of micromechanical fields in heterogeneous materials is usually performed using either the finite element method or the Green’s function method based on FFTs. The finite element method allows for accurate discretization and for non-periodic boundary conditions but is computationally expensive. On the other hand, the FFT-based method is computationally efficient but requires discretization on a regular grid of hexahedral voxels. In this paper, a Green’s function method allowing for accurate discretization using tetrahedral elements and for non-periodic boundary conditions is proposed. The convolution is computed using the fast multipole method, which provides good accuracy even for low-order expansion due to the fast decay of interactions between elements. The proposed Green’s function fast multipole method is verified by comparison with analytical and FFT-based solutions. Furthermore, the computational time is analyzed and compared to the FFT-based method for non-periodic convolution. Finally, effective properties of an elastic polycrystalline microstructure containing thin intergranular cracks are computed and analyzed.

36 MATERIALS SCIENCE↗

Calibration Issues and Operating System Requirements for Electron-Probe Microanalysis

Instrument purchase requirements and dialogue with manufacturers have established hardware parameters for alignment, stability, and reproducibility, which have helped improve the precision and accuracy of electron microprobe analysis (EPMA). The development of correction algorithms and the accurate solution to quantitative analysis problems requires the minimization of systematic errors and relies on internally consistent data sets. Improved hardware and computer systems have resulted in better automation of vacuum systems, stage and wavelength-dispersive spectrometer (WDS) mechanisms, and x-ray detector systems which have improved instrument stability and precision. Improved software now allows extended automated runs involving diverse setups and better integrates digital imaging and quantitative analysis. However, instrumental performance is not regularly maintained, as WDS are aligned and calibrated during installation but few laboratories appear to check and maintain this calibration. In particular, detector deadtime (DT) data is typically assumed rather than measured, due primarily to the difficulty and inconvenience of the measurement process. This is a source of fundamental systematic error in many microprobe laboratories and is unknown to the analyst, as the magnitude of DT correction is not listed in output by microprobe operating systems. The analyst must remain vigilant to deviations in instrumental alignment and calibration, and microprobe system software must conveniently verify the necessary parameters. Microanalysis of mission critical materials requires an ongoing demonstration of instrumental calibration. Possible approaches to improvements in instrument calibration, quality control, and accuracy will be discussed. Development of a set of core requirements based on discussions with users, researchers, and manufacturers can yield documents that improve and unify the methods by which instruments can be calibrated. These results can be used to continue improvements of EPMA.

Carpenter, P.↗

A browsing tool for the Internet Logical Library of the HPCC Software Exchange

As the quantity of information available on the Internet grows, locating a particular piece of information becomes more difficult. One possible solution is for a database of pointers to all available information to be maintained at a central site. Subject classifications for all the information could also be maintained in order to make searching possible. This paper describes one possible method of searching such an index. In particular a prototype browsing tool has been created using TCL/TK to demonstrate several possible features: rapidly scanning at any rank of the index, narrowing the index to any scope, regular-expression searching, and creation of a list of pointers answering to any set of index terms. The prototype browser is an easy-to-use independent X application designed for use in the Catalog of Repositories of the HPCC (High Performance Computing and Communications) Software Exchange.

Biro, Ross↗

Reconstruction of the 4D beam matrix

The widely used transverse parameters characterizing particle beams are the Twiss parameters. These parameters can be measured experimentally but they do not fully characterize the beam since they do not account for possible correlations in particle distribution between two transverse coordinates. These correlations may occur due to uncompensated magnetic field at the cathode or misalignment of focusing quadrupoles in the transport beamline. We test a novel diagnostic for diagnosing full 4D beam matrix which may be used to identify such imperfections. The diagnostic is based on transporting the beam through the beamline which includes a quadrupole and a skew quadrupole magnets and measuring the resulting 2D beam distribution at the screen downstream. Such a measurement can be viewed as measuring a 2D projection of the 4D distribution. Different settings of the quads provide measurements of different slices of the phase space. The reconstruction of the original beam matrix from a number of measurements is done using machine learning algorithm, which provides a fast and reliable way of reconstruction for an arbitrary configuration of the scanning beamline. In August 2024, we set up the diagnostic beamline to perform a quadrupole scan of the beam. The setup includes a skew quadrupole, a regular quadrupole, and a screen. The images on the screen were post-processed to remove experimental artifacts and enhance contrast by eliminating background noise outside the core of the distribution=. The rms parameters of the distribution were then calculated and used as inputs for the reconstruction algorithm. This algorithm attempts to determine the initial beam matrix that produces expected images on the screen closely matching the observed images across all quadrupole settings. The algorithm found a solution in which the expected rms parameters closely align with the observations. Validation of the results is planned for FY25.

43 PARTICLE ACCELERATORS↗

Elimination of secular terms from the differential equations for the elements of perturbed two-body motion

In 1961, Sperling linearized and regularized the differential equations of motion of the two-body problem by changing the independent variable from time to fictitious time by Sundman's transformation (r = dt/ds) and by embedding the two-body energy integral and the Laplace vector. In 1968, Burdet developed a perturbation theory which was uniformly valid for all types of orbits using a variation of parameters approach on the elements which appeared in Sperling's equations for the two-body solution. In 1973, Bond and Hanssen improved Burdet's set of differential equations by embedding the total energy (which is a constant when the potential function is explicitly dependent upon time.) The Jacobian constant was used as an element to replace the total energy in a reformulation of the differential equations of motion. In the process, another element which is proportional to a component of the angular momentum was introduced. Recently trajectories computed during numerical studies of atmospheric entry from circular orbits and low thrust beginning in near-circular orbits exhibited numerical instability when solved by the method of Bond and Gottlieb (1989) for long time intervals. It was found that this instability was due to secular terms which appear on the righthand sides of the differential equations of some of the elements. In this paper, this instability is removed by the introduction of another vector integral called the delta integral (which replaces the Laplace Vector) and another scalar integral which removes the secular terms. The introduction of these integrals requires a new derivation of the differential equations for most of the elements. For this rederivation, the Lagrange method of variation of parameters is used, making the development more concise. Numerical examples of this improvement are presented.

Bond, Victor R.↗

CFD analysis of turbopump volutes

An effort is underway to develop a procedure for the regular use of CFD analysis in the design of turbopump volutes. Airflow data to be taken at NASA Marshall will be used to validate the CFD code and overall procedure. Initial focus has been on preprocessing (geometry creation, translation, and grid generation). Volute geometries have been acquired electronically and imported into the CATIA CAD system and RAGGS (Rockwell Automated Grid Generation System) via the IGES standard. An initial grid topology has been identified and grids have been constructed for turbine inlet and discharge volutes. For CFD analysis of volutes to be used regularly, a procedure must be defined to meet engineering design needs in a timely manner. Thus, a compromise must be established between making geometric approximations, the selection of grid topologies, and possible CFD code enhancements. While the initial grid developed approximated the volute tongue with a zero thickness, final computations should more accurately account for the geometry in this region. Additionally, grid topologies will be explored to minimize skewness and high aspect ratio cells that can affect solution accuracy and slow code convergence. Finally, as appropriate, code modifications will be made to allow for new grid topologies in an effort to expedite the overall CFD analysis process.

Ascoli, Edward P.↗

A coupled discontinuous Galerkin-Finite Volume framework for solving gas dynamics over embedded geometries

Herein, we present a computational framework for solving the equations of inviscid gas dynamics using structured grids with embedded geometries. The novelty of the proposed approach is the use of high-order discontinuous Galerkin (dG) schemes and a shock-capturing Finite Volume (FV) scheme coupled via an hp adaptive mesh refinement (hp-AMR) strategy that offers high-order accurate resolution of the embedded geometries. The hp-AMR strategy is based on a multi-level block-structured domain partition in which each level is represented by block-structured Cartesian grids and the embedded geometry is represented implicitly by a level set function. The intersection of the embedded geometry with the grids produces the implicitly-defined mesh that consists of a collection of regular rectangular cells plus a relatively small number of irregular curved elements in the vicinity of the embedded boundaries. High-order quadrature rules for implicitly-defined domains enable high-order accuracy resolution of the curved elements with a cell-merging strategy to address the small-cell problem. The hp-AMR algorithm treats the system with a second-order finite volume scheme at the finest level to dynamically track the evolution of solution discontinuities while using dG schemes at coarser levels to provide high-order accuracy in smooth regions of the flow. On the dG levels, the methodology supports different orders of basis functions on different levels. The space-discretized governing equations are then advanced explicitly in time using high-order Runge-Kutta algorithms. Numerical tests are presented for two-dimensional and three-dimensional problems involving an ideal gas. The results are compared with both analytical solutions and experimental observations and demonstrate that the framework provides high-order accuracy for smooth flows and accurately captures solution discontinuities.

97 MATHEMATICS AND COMPUTING↗

Solution adaptivity using a triangular mesh

Solution adaptivity is discussed first from a general perspective and then from the specific viewpoint of triangular meshes. The use of a general connectivity triangular mesh is emphasized. The development of monitor surfaces and their geometric properties is discussed. Mesh point movement is addressed, as is the dynamic restructuring of the connectivity pattern among moving modes. Changes in the number of mesh nodes to obtain suitable refinement for a physical simulation is examined, and the use of locally regular structures to offset the data structure limitation of a general connectivity triangular mesh and thus to obtain an enhanced range of application is considered. To illustrate the basic features of the adaptive triangular mesh strategy, an application to the study of plasma equilibrium is briefly considered.

Eiseman, P. R.↗