SEARCH · Engineering Papers

Results for “numerical continuation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,585 records · Page 88

All sources NASA NTRS DOE OSTI

NASA NTRS2003

Rapid Technology Assessment via Unified Deployment of Global Optical and Virtual Diagnostics

This paper discusses recent developments in rapid technology assessment resulting from an active collaboration between researchers at the Air Force Research Laboratory (AFRL) at Wright Patterson Air Force Base (WPAFB) and the NASA Langley Research Center (LaRC). This program targets the unified development and deployment of global measurement technologies coupled with a virtual diagnostic interface to enable the comparative evaluation of experimental and computational results. Continuing efforts focus on the development of seamless data translation methods to enable integration of data sets of disparate file format in a common platform. Results from a successful low-speed wind tunnel test at WPAFB in which global surface pressure distributions were acquired simultaneously with model deformation and geometry measurements are discussed and comparatively evaluated with numerical simulations. Intensity- and lifetime-based pressure-sensitive paint (PSP) and projection moire interferometry (PMI) results are presented within the context of rapid technology assessment to enable simulation-based R&D.

Jordan, Jeffrey D.↗

NASA NTRS2006

High-Resolution Analysis Products to Support Severe Weather and Cloud-to-Ground Lightning Threat Assessments over Florida

The Applied Meteorology Unit (AMU) located at the Kennedy Space Center (KSC)/Cape Canaveral Air Force Station (CCAFS) implemented an operational configuration of the Advanced Regional Prediction System (ARPS) Data Analysis System (ADAS), as well as the ARPS numerical weather prediction (NWP) model. Operational, high-resolution ADAS analyses have been produced from this configuration at the National Weather Service in Melbourne, FL (NWS MLB) and the Spaceflight Meteorology Group (SMG) over the past several years. Since that time, ADAS fields have become an integral part of forecast operations at both NWS MLB and SMG. To continue providing additional utility, the AMU has been tasked to implement visualization products to assess the potential for supercell thunderstorms and significant tornadoes, and to improve assessments of short-term cloud-to-ground (CG) lightning potential. This paper and presentation focuses on the visualization products developed by the AMU for the operational high-resolution ADAS and AR.PS at the NWS MLB and SMG. The two severe weather threat graphics implemented within ADAS/ARPS are the Supercell Composite Parameter (SCP) and Significant Tornado Parameter (SIP). The SCP was designed to identify areas with supercell thunderstorm potential through a combination of several instability and shear parameters. The SIP was designed to identify areas that favor supercells producing significant tornadoes (F2 or greater intensity) versus non-tornadic supercells. Both indices were developed by the NOAAINWS Storm Prediction Center (SPC) and were normalized by key threshold values based on previous studies. The indices apply only to discrete storms, not other convective modes. In a post-analysis mode, the AMU calculated SCP and SIP for graphical output using an ADAS configuration similar to the operational set-ups at NWS MLB and SMG. Graphical images from ADAS were generated every 15 minutes for 13 August 2004, the day that Hurricane Charley approached and made landfall on the Florida peninsula. Several tornadoes struck the interior of the Florida peninsula in advance of Hurricane Charley's landfall during the daylight hours of 13 August. Since SPC had previously examined this case using SCP and SIP graphics generated from output of the Rapid Update Cycle (RUC) model, this day served as a good benchmark to compare and validate the high-resolution ADAS graphics against the smoother RUC analyses, which serves as background fields to the ADAS analyses. The ADAS-generated SCP and STP graphics have been integrated into the suite of products examined operationally by NWS MLB forecasters and are used to provide additional guidance for assessment of the near-storm environment during convective situations.

Case, Jonathan↗

DOE OSTI2023

SCALE 6.2 Lattice Physics Performance Assessment

The US Nuclear Regulatory Commission relies on the lattice physics analysis capabilities of the SCALE code system to perform confirmatory licensing analyses. Either SCALE lattice physics code—TRITON/NEWT or Polaris—can be used to generate cross section data used by the PARCS nodal core simulator for full-core neutronics calculations. This report presents an assessment of the accuracy of SCALE lattice physics codes for preparation of lattice physics data that are used to support simulator codes such as the NRC’s PARCS, for UO 2 -mixed oxide (MOX)/Zr fueled light water reactor (LWR) analyses. Due to the nature of lattice physics calculations, critical reactor experiment benchmarks cannot be modeled in explicit detail in a lattice physics code. However, this limitation does not mean that these measurement data are not usable for lattice physics studies. Therefore, either geometry approximations or axial buckling must be implemented to determine the critical water height. These modeling limitations have led to development of a three-phase assessment strategy. In the first phase, selected critical experiment benchmarks are modeled using the SCALE 3D continuous-energy (CE) Monte Carlo (MC) code KENO, which is the most rigorous neutron transport method available in SCALE, with no approximations in the spatial, angle, or energy treatments. Biases and statistical uncertainties in quantities of interest such as k eff and pin power distributions are determined by comparing CE KENO results to experimental data. This first phase in the accuracy assessment is to demonstrate and establish the use of CE KENO as a reference solution for the second phase. In the second phase, 14 numerical test suites are used to compare SCALE lattice physics calculations with CE KENO as a reference solution. In the third phase, SCALE lattice physics depletion calculations are performed, and the spent fuel isotopic results are compared with available radiochemical assay measurements. Isotopic measurement comparisons provide quantitative assessment of isotopic density distribution predictions with the depletion models in SCALE lattice physics codes. This report documents results for all test suites. The assessment was performed using standard production techniques unless otherwise noted. Both TRITON/NEWT and Polaris exhibited acceptable accuracy for most test cases. For the few test cases in which acceptable accuracy criteria were not met, further code and data development are planned. The computer codes used in this assessment are as follows: SCALE 6.2 rev19189 (pre-release of 6.2.1) was used for CE KENO, TRITON/NEWT and Polaris (PWR only) calculations; SCALE 6.2.2 was used for Polaris calculations to address several updates in support of support boiling water reactor (BWR) lattice geometries. For this work, there is no difference between SCALE 6.2.1 and SCALE 6.2.2 beyond the inclusion of the Polaris BWR analysis capability; The ENDF/B-VII.1 continuous energy and 252 group libraries deployed in SCALE 6.2 were generated with AMPX 6.2 that is distributed with SCALE.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Electric field control of magnetism: multiferroics and magnetoelectrics

This article is written on behalf of a large number of colleagues, collaborators, and researchers in the field of complex oxides as well as current and former students and postdocs who continue to enable and undertake cutting-edge research in the field of multiferroics, magnetoelectrics, and the pursuit of electric-field control of magnetism. What we present is something that is extremely exciting from both a fundamental science and applications perspective and has the potential to revolutionize our world. Needless to say, to realize this potential will require numerous new innovations, both in the fundamental science arena as well as translating these scientific discoveries into real applications. Thus, this article will attempt to bridge the gap between fundamental condensed-matter physics and the actual manifestations of the physical concepts into real-life applications. Here, we hope this article will help spur more translational research within the broad materials physics community.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

A phase-field diffraction model for thermo-hydro-mechanical propagating fractures

This paper introduces a novel diffraction based thermo-hydraulic–mechanical (THM) model for fracture propagation using a phase-field fracture (PFF) approach. The key innovation of the THM-PFF model lies in its integrated treatment of four solution variables—displacements, phase-field, pressure, and temperature—each governed by a combination of conservation of momentum (mechanics problem), a variational inequality (constrained minimization problem), mass conservation (pressure problem), and energy conservation (temperature problem). This leads to a new formulation of a coupled variational inequality system. A major advancement is the development of an extended fixed-stress algorithm, where displacements, phase-field, pressures, and temperatures are solved in a staggered sequence. An important aspect of this work is the global coupling of pressures and temperatures across the domain using diffraction systems, with diffraction coefficients defined by material parameters weighted by the diffusive phase-field variable. To ensure robust local mass conservation, we employ enriched Galerkin finite elements (EG) for both pressure and temperature diffraction equations. By enriching the continuous Galerkin basis functions with discontinuous piecewise constants, EG accurately represents solution and parameter discontinuities while preserving local mass and energy conservation—crucial aspects for THM problems and realistic behavior. Moreover, the use of a predictor–corrector local mesh adaptivity scheme is employed, allowing the model to handle small phase-field length-scale parameters while maintaining high numerical accuracy and reasonable computational cost. Furthermore, these new model and algorithmic developments represent significant advances in the field and have been substantiated through rigorous numerical tests.

Diffraction systems↗

DOE OSTI2023

Predictive scale-bridging simulations through active learning

Throughout computational science, there is a growing need to utilize the continual improvements in raw computational horsepower to achieve greater physical fidelity through scale-bridging over brute-force increases in the number of mesh elements. For instance, quantitative predictions of transport in nanoporous media, critical to hydrocarbon extraction from tight shale formations, are impossible without accounting for molecular-level interactions. Similarly, inertial confinement fusion simulations rely on numerical diffusion to simulate molecular effects such as non-local transport and mixing without truly accounting for molecular interactions. With these two disparate applications in mind, we develop a novel capability which uses an active learning approach to optimize the use of local fine-scale simulations for informing coarse-scale hydrodynamics. Our approach addresses three challenges: forecasting continuum coarse-scale trajectory to speculatively execute new fine-scale molecular dynamics calculations, dynamically updating coarse-scale from fine-scale calculations, and quantifying uncertainty in neural network models.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Pulse-by-pulse transient thermal deformation in crystal optics under high-repetition-rate FEL

Time-domain modeling of the thermal deformation of crystal optics can help define acceptable operational ranges across the pulse-energy repetition-rate phase space. In this paper, we have studied the transient thermal deformation of a water-cooled diamond crystal for a cavity-based X-ray free-electron laser (CBXFEL), either an X-ray free-electron laser oscillator (XFELO) or a regenerative amplifier X-ray free-electron laser (RAFEL), by numerical simulations including finite-element analysis and advanced data processing. Pulse-by-pulse transient thermal deformation of a 50 µm-thick diamond crystal has been performed with X-ray pulse repetition rates between 50 kHz and 1 MHz. Results for temperature and thermal deformation have been compared with the results of transient analysis using a continuous wave (CW) power loading. Temperature and thermal deformation results from pulse-by-pulse transient analysis vary with time about the results for the CW case for the same average power. The variation amplitude increases with pulse energy and decreases with repetition rate. When the repetition rate increases to infinity, both temperature and thermal deformation converge to the results for the CW case. Two critical time scales for the operation of crystal optics in a CBXFEL are (1) first-turn time, i.e. the time for the XFEL pulse to complete the first turn around the cavity so that the crystal sees the recirculated XFEL pulse, and (2) period-end time, i.e. the time that the next electron bunch arrives for the amplification, so that the crystal outcouples the amplified FEL power. For the same average power, simulation results show that the crystal thermal deformation seen by the XFEL beam decreases with repetition rate at the first-turn time of a 300 m-long cavity and increases with repetition rate at the period-end time. For the wavefront preservation requirement of the crystal optics, a pulse-energy versus repetition-rate phase space has been established. The upper bounds of the pulse energy at both first-turn and period-end times decreases with repetition rate, especially at the period-end time. The upper bound of the thermal deformation of the crystal at the period-end time for any repetition frequency can be estimated from the CW case. For a water-cooled diamond crystal of dimension 5 mm × 5 mm × 0.05 mm, the time to reach a quasi steady-state is about 50 ms for temperature and 50 µs for thermal deformation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Galvanostatic Intermittent Titration Technique Reinvented: Part I. A Critical Review

The galvanostatic intermittent titration technique (GITT), introduced in 1977 by Weppner and Huggins, provided a readily accessible means to measuring the chemical diffusion coefficient of electrochemical electrode materials. The method continues to be widely used today, but the reported diffusivity values are highly inconsistent, ranging as much as four orders of magnitude for some Li layered oxide compositions. Even qualitative trends of diffusivity are inconsistent, suggesting significant flaws in the implementation of the method. Other variants of the GITT method also suffer from similar inconsistency problems. Here we identify numerous sources of significant error including composition-dependent reaction overpotentials, mathematical flaws in the relaxation analysis methods, finite-size and non-planar geometry effects, inter-particle inhomogeneity issues, early transient effects, and surface area uncertainties. We propose a simple relaxation analysis scheme using the time variable t relax + τ - t relax , where t relax is relaxation time and τ is the galvanostatic pulse duration. We also propose to use dense diffusion-limited samples to isolate the bulk-diffusion process in the time domain. Chemical diffusivity can be extracted much more reliably with this improved implementation of the GITT method.

25 ENERGY STORAGE↗

DOE OSTI2026

Roadmap to Advance Heliostat Technologies for High Temperature Solar-Thermal Systems

Since its establishment, the Heliostat Consortium (HelioCon) has made substantial progress toward closing many of the gaps in concentrating solar power (CSP) research. Numerous techno-economic studies have been performed, investigating topics ranging from the trade-off between size and temperature for industrial process heat applications to optimization of the heliostat design itself for various applications. Significant improvements have been made in optical metrology techniques, with first steps toward in situ measurement of heliostat fields. Several standards have been, and continue to be, developed with the coordination of an international group of CSP industry participants. Training programs have been developed, with universities including CSP in their engineering curricula, and many public webinars have been held to provide broad access to the latest CSP research. Improved CSP components such as mirror facets and wireless communication systems have been developed, and the solar tower at Sandia National Laboratories has been upgraded with a testbed for closed-loop controls research and development. Field deployment challenges involving heliostat foundations and sensitive wildlife habitats have been explored, with progress made toward methods for streamlining project development and permitting. Additional knowledge has been added to the body of work on wind behavior of heliostats and arrays of heliostats, with progress made toward a holistic understanding of wind design methods. Finally, techniques have been developed and demonstrated for assessing soiling conditions at a proposed project site, with predictive models for the soiling rate showing good results. Taking these results together, HelioCon has contributed greatly to the global CSP research and development effort over the past several years.

14 SOLAR ENERGY↗

NASA NTRS1982

Robustness of adaptive control algorithms in the presence of unmodeled dynamics

This paper reports the outcome of an exhaustive analytical and numerical investigation of stability and robustness properties of a wide class of adaptive control algorithms in the presence of unmodeled dynamics and output disturbances. The class of adaptive algorithms considered are those commonly referred to as model-reference adaptive control algorithms, self-tuning controllers, and dead-beat adaptive controllers; they have been developed for both continuous-time systems and discrete-time systems. The existing adaptive control algorithms have been proven to be globally asymptotically stable under certain assumptions, the key ones being (1) that the number of poles and zeroes of the unknown plant are known, and (2) that the primary performance criterion is related to good command following. These theoretical assumptions are too restrictive from an engineering point of view. Real plants always contain unmodeled high-frequency dynamics and small delays, and hence no upper bound on the number of the plant poles and zeroes exists. Also real plants are always subject to unmeasurable output additive disturbances, although these may be guide small. Hence, it is important to critically examine the stability robustness properties of the existing adaptive algorithms when some of the theoretical assumptions are removed; in particular, their stability and performance properties in the presence of unmodeled dynamics and output disturbances. Previously announced in STAR as N83-16061

Rohrs, C. E.↗

NASA NTRS1984

Investigations of flow field perturbations induced on slotted transonic-tunnel walls

The free-stream interference caused by the flow through the slotted walls of the test sections of transonic wind tunnels has continuously a problem in transonic tunnel testing. The adaptive-wall transonic tunnel is designed to actively control the near-wall boundary conditions by sucking or blowing through the wall. In order to make the adaptive-wall concept work, parameters for computational boundary conditions must be known. These parameters must be measured with sufficient accuracy to allow numerical convergence of the flow field computations and must be measured in an inviscid region away from the model that is placed inside the wind tunnel. The near-wall flow field was mapped in detail using a five-port cone probe that was traversed in a plane transverse to the free-stream flow. The initial experiments were made using a single slot and recent measurements used multiple slots, all with the tunnel empty. The projection of the flow field velocity vectors on the transverse plane revealed the presence of a vortex-like flow with vorticity in the free stream. The current research involves the measurement of the flow field above a multislotted system with segmented plenums behind it, in which the flow is controlled through several plenums simultaneously. This system would be used to control a three-dimensional flow field.

Wu, J. M.↗

NASA NTRS1988

Cut-and-connect of two antiparallel vortex tubes

Motivated by an early conjecture that vortex cut-and-connect plays a key role in mixing and production of turbulence, helicity and aerodynamic noise, the cross-linking of two antiparallel viscous vortex tubes via direct numerical simulation is studied. The Navier-Stokes equations are solved by a dealiased pseudo-spectral method with 64 cubed grid points in a periodic domain for initial Reynolds numbers Re up to 1000. The vortex tubes are given an initial sinusoidal perturbation to induce a collision and keep the two tubes pressed against each other as annihilation continues. Cross-sectional and wire plots of various properties depict three stages of evolution: (1) Inviscid induction causing vortex cores to first approach and form a contact zone with a dipole cross-section, and then to flatten and stretch; (2) Vorticity annihilation in the contact zone accompanied by bridging between the two vortices at both ends of the contact zone due to a collection of cross-linked vortex lines, now orthogonal to the initial vortex tubes. The direction of dipole advection in the contact zone reverses; and (3) Threading of the remnants of the original vortices in between the bridges as they pull apart. The crucial stage 2 is shown to be a simple consequence of vorticity annihilation in the contact zone, link-up of the un-annihilated parts of vortex lines, and stretching and advection by the vortex tube swirl of the cross-linked lines, which accumulate at stagnation points in front of the annihilating vortex dipole. It is claimed that bridging is the essence of any vorticity cross-linking and that annihilation is sustained by stretching of the dipole by the bridges. Vortex reconnection details are found to be insensitive to asymmetry. Modeling of the reconnection process is briefly examined. The 3D spatial details of scalar transport (at unity Schmidt number), enstrophy production, dissipation and helicity are also examined.

Melander, Mogens V.↗

NASA NTRS1993

Numerical study of the effects of icing on viscous flow over wings

The progress made during the period 1 Jan. - 30 Jun. 1993 on the numerical study of the effects of icing on viscous flow over wings is summarized. The tasks performed were development of 3-D boundary layer methods for accurate estimates of the velocity field and surface heat transfer rates in the vicinity of the leading edge ice shape; studies of the effects of icing on 3-D highlift system performance; and continued improvement and validation of the 3-D Navier-Stokes solver. Results are discussed.

Sankar, L. N.↗

NASA NTRS1994

3-D simulation of high-latitude interaction regions: Comparison with Ulysses results

A three-dimensional (3-D) magnetohydrodynamic (MHD) numerical model is used to simulate the global evolution of a steady, tilted-dipole solar wind flow configuration similar to that prevalent in interplanetary space in 1993. Systematic latitudinal changes in the structure of a corotating interaction region (CIR) near 5 AU is shown to agree well with recent Ulysses observations. The abrupt disappearance of forward shocks and continued persistence of reverse shocks poleward of the latitude where Ulysses crossed the southern edge of the coronal streamer belt is explained as a natural consequence of the 3-D flow geometry.

Pizzo, V. J.↗

NASA NTRSDate not supplied

Observing and Projecting the Lasting Fate of the Hunga Eruption on Atmospheric Water Vapor and Hydroxyl Radical

On 15th Jan. 2022 the submarine Hunga volcano erupted injecting approximately 0.5 Tg of sulfur dioxide into the stratosphere, but more significantly added 150-170 Tg of water vapor to the stratospheric background (over a 10% perturbation) in a matter of several hours. The sulfur dioxide rapidly converted to sulfate aerosol and along with water vapor, was transported around the Southern Hemisphere sub-tropics into midlatitudes and beyond with some transport into the Northern Hemisphere. With a much longer lifetime than sulfate aerosol, measurable water vapor mass anomalies have persisted with only small losses over more than 2 years since the eruption and are likely to continue above background for the remainder of the decade. Satellite measurements from limb and nadir viewing observing instruments provide the information needed to reasonably initialize the Hunga eruption in the Goddard Earth Observing System (GEOS) model using the “replay” framework coupled to the Global Modeling Initiative (GMI) stratosphere-troposphere chemical mechanism for the recent past and have continued the simulations into the future with the free running chemistry climate model (CCM). Using a number of model ensemble members together with the satellite observations, we quantify how the Hunga eruption is perturbing atmospheric composition and climate and projecting the influences to come as the enhanced water vapor continues in the stratosphere with only very slow removal mechanisms. Modeling suggests that the water rich eruption enhanced stratospheric hydroxyl radiacal (OH) and is supported by the speed of the observed conversion of SO2 into sulfate aerosol. These simulations also suggest a significant reduction in tropospheric OH over the mid-high latitude Southern Hemisphere, which has implications for numerous species that are impacted by this radical. The Hunga eruption provides a useful test of chemistry climate models and an opportunity for observation-based process understanding and in making new connections.

Hunga volcano↗

DOE OSTI2024

A simplified all-frequency stable formulation with an implicit Coulomb gauge

A potential-based finite element formulation has been developed in the past to circumvent the low-frequency breakdown issues commonly encountered in electromagnetic (EM) simulations of low-frequency and multiscale problems. In this formulation, the magnetic vector and electric scalar potentials have been employed to express the electric field and magnetic flux, leading to a set of two equations that represents all four Maxwell’s equations and the current continuity equation. To enforce the Coulomb gauge, an auxiliary potential has been introduced, which results in a total of three equations to be solved simultaneously. To reduce the number of equations and unknowns needed in a simulation, a simplified formulation is proposed in this paper to enforce the Coulomb gauge implicitly without the need for the auxiliary potential. A numerical example is given to demonstrate the accuracy of the proposed formulation.

Mekonnen, Minyechil↗

NASA NTRS1987

Development and evaluation of improved numerical schemes for recirculating flows

The paper examines the performance of the flux-spline scheme for convection-diffusion. Computations are presented for a number of test cases, both linear and nonlinear. It is shown that in all cases the flux-spline scheme yields results which are superior to those obtained with the lower-order formulations such as hybrid differencing. In order to improve the computational efficiency, the flux-spline scheme has been combined with a direct solution algorithm for the continuity and momentum equations. Such an approach eliminates the need for an equation for pressure or pressure correction and is found to be rapidly convergent.

Patankar, S. V.↗

NASA NTRS1989

LEO high voltage solar array arcing response model, continuation 5

The modeling of the Debye Approximation electron sheaths in the edge and strip geometries was completed. Electrostatic potentials in these sheaths were compared to NASCAP/LEO solutions for similar geometries. Velocity fields, charge densities and particle fluxes to the biased surfaces were calculated for all cases. The major conclusion to be drawn from the comparisons of our Debye Approximation calculations with NASCAP-LEO output is that, where comparable biased structures can be defined and sufficient resolution obtained, these results are in general agreement. Numerical models for the Child-Langmuir, high-voltage electron sheaths in the edge and strip geometries were constructed. Electrostatic potentials were calculated for several cases in each of both geometries. Velocity fields and particle fluxes were calculated. The self-consistent solution process was carried through one cycle and output electrostatic potentials compared to NASCAP-type input potentials.

Metz, Roger N.↗