Engineering Papers⌕ Search

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.

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

At least 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Dynamic optical spectroscopy and pyrometry of static targets under optical and x-ray laser heating at the European XFEL

Experiments accessing extreme conditions at x-ray free electron lasers (XFELs) involve rapidly evolving conditions of temperature. Here, we report time-resolved, direct measurements of temperature using spectral streaked optical pyrometry of x-ray and optical laser-heated states at the High Energy Density instrument of the European XFEL. This collection of typical experiments, coupled with numerical models, outlines the reliability, precision, and meaning of time dependent temperature measurements using optical emission at XFEL sources. Dynamic temperatures above 1500 K are measured continuously from spectrally- and temporally-resolved thermal emission at 450–850 nm, with time resolution down to 10–100 ns for 1–200 μs streak camera windows, using single shot and integrated modes. Targets include zero-pressure foils free-standing in air and in vacuo, and high-pressure samples compressed in diamond anvil cell multi-layer targets. Radiation sources used are 20-fs hard x-ray laser pulses at 17.8 keV, in single pulses or 2.26 MHz pulse trains of up to 30 pulses, and 250-ns infrared laser single pulses. A range of further possibilities for optical measurements of visible light in x-ray laser experiments using streak optical spectroscopy are also explored, including for the study of x-ray induced optical fluorescence, which often appears as background in thermal radiation measurements. We establish several scenarios where combined emissions from multiple sources are observed and discuss their interpretation. Challenges posed by using x-ray lasers as non-invasive probes of the sample state are addressed.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

TRIPOLI-4 neutronics calculations for IAEA-CRP benchmark of CEFR start-up tests using new libraries JEFF-3.3 and ENDF/B-VIII

Sodium-cooled Fast Reactors (SFR) are one of the promising Generation IV fast reactors. The China Experimental Fast Reactor (CEFR) is a 65 MWth pool-type SFR with a high neutron leakage core using high enriched uranium oxide fuel. The CEFR start-up tests in 2010 consist of series of neutronics experiments. Essential experimental data are available from the 2018 IAEA-CEFR-CRP benchmark document and useful for the validation of neutron transport codes and nuclear data libraries. The TRIPOLI-4 Monte Carlo transport code is a general-purpose neutronics code using continuous-energy nuclear data libraries. It has a rich validation database covering different computational and experimental benchmark data sets to assure the accuracy and credibility of numerical studies. Previous TRIPOLI-4 SFR core physics calculations used mainly MOX fuels. SFR control rod worth studies with TRIPOLI-4 were mainly on big-size cores. To test the modeling capability of TRIPOLI-4 on fuel loading patterns and on variable control rods positions for SFRs and to check recent nuclear data libraries, it is interesting to investigate the CEFR start-up tests in this work, including core states from subcritical to supercritical, control rods and rod group worth, point kinetics parameters, and radial foil activation measurements. TRIPOLI-4 calculation results using new data libraries JEFF-3.3, ENDF/B-VIII, and those of ENDF/B-VII.1 for different CEFR core states are reported here. By means of different options of TRIPOLI-4, control-rod SAs reactivity worth, point kinetics parameters, and radial {sup 237}Np(n,f) fission rate distribution were successfully evaluated. Both ENDF/B-VIII and JEFF-3.3 nuclear data libraries provided reliable results with a difference of 170 +/- 11 pcm in k{sub eff}.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Witness the Trinity test through Lab artifacts

July 16 marks the 77th anniversary of the Trinity test, conducted in a desert in New Mexico. The test subject, an atomic bomb called The Gadget, was successfully detonated from a 100-foot steel tower. This event marks the commencement of the Atomic Age, a new era where fission capabilities could be employed for national security purposes. Shortly after the Trinity test, two Los Alamos-created atomic weapons were released above Japan, helping to end the world’s bloodiest conflict just weeks later. “Trinity was one of the greatest scientific experiments ever,” said NSRC Senior Historian Alan Carr said. “Los Alamos scientists changed the world forever on that day. Not only was it the dawn of the Atomic Age, but also the beginning of the Lab’s eight decades of cutting-edge science and its national security charge.” To preserve this event, and to continue to learn more about this critical moment in history, the National Security Research Center (NSRC) curates a collection of photographs, films, notes, unclassified artifacts and numerous other materials related to the science of the test. Notably, the collection includes a novel material that formed at the site, trinitite, and artifacts from one of the intriguing scientists present at the test, Enrico Fermi.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2020

"Russians in the Grid" Enrichment Process

Threat intelligence sources are numerous, and the quality of their intelligence is equally variable. The need for quality, actionable, and shareable intelligence is becoming increasingly important as cyber-threats continue to increase. In this paper, we describe the process we used to enrich U.S. CERT’s Technical Alert entitled “Russian Government Cyber-Activity Targeting Energy and Other Critical Infrastructure Sectors,” affectionately referred to as “Russians In the Grid,” from a series of separate but related reports to a final, actionable, unified Structured Threat Information Expression (STIX) v2 representation of the report. In the process, we detail the tools and procedures we used, some of which are automated and others highly manual. We identify some of the limitations of current STIX enrichment tools and techniques, and propose enhancements to allow for more effective enrichment and better overall threat intelligence to be used for more effective threat data-sharing, as well as machine learning-based predictive threat work being conducted at Idaho National Laboratory.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

Considerations for a Co-Located Aquaculture and Wave Energy Deployment: Quantifying Demand and Assessing Integration

The co-location of marine energy and aquaculture is a concept of increasing value and interest in the United States, as the desire for sustainably produced seafood and renewable energy continues to grow. With both industries being fairly nascent in their national development and facing challenges, the deployment of a wave-powered aquaculture farm requires numerous considerations. This report analyzes the energy demands of an integrated multi-trophic aquaculture system at various farm scales and discusses the preliminary technical, regulatory, and logistical considerations required for a co-located deployment. A pilot-scale deployment is outlined with potentially ideal characteristics, and an optimized planning and permitting approach is introduced.

16 TIDAL AND WAVE POWER↗
NASA NTRS1991

A computer code for multiphase all-speed transient flows in complex geometries. MAST version 1.0

The operation of the MAST code, which computes transient solutions to the multiphase flow equations applicable to all-speed flows, is described. Two-phase flows are formulated based on the Eulerian-Lagrange scheme in which the continuous phase is described by the Navier-Stokes equation (or Reynolds equations for turbulent flows). Dispersed phase is formulated by a Lagrangian tracking scheme. The numerical solution algorithms utilized for fluid flows is a newly developed pressure-implicit algorithm based on the operator-splitting technique in generalized nonorthogonal coordinates. This operator split allows separate operation on each of the variable fields to handle pressure-velocity coupling. The obtained pressure correction equation has the hyperbolic nature and is effective for Mach numbers ranging from the incompressible limit to supersonic flow regimes. The present code adopts a nonstaggered grid arrangement; thus, the velocity components and other dependent variables are collocated at the same grid. A sequence of benchmark-quality problems, including incompressible, subsonic, transonic, supersonic, gas-droplet two-phase flows, as well as spray-combustion problems, were performed to demonstrate the robustness and accuracy of the present code.

Chen, C. P.↗
NASA NTRS1992

Multigrid calculation of internal flows in complex geometries

The development, validation, and application of a general purpose multigrid solution algorithm and computer program for the computation of elliptic flows in complex geometries is presented. This computer program combines several desirable features including a curvilinear coordinate system, collocated arrangement of the variables, and Full Multi-Grid/Full Approximation Scheme (FMG/FAS). Provisions are made for the inclusion of embedded obstacles and baffles inside the flow domain. The momentum and continuity equations are solved in a decoupled manner and a pressure corrective equation is used to update the pressures such that the fluxes at the cell faces satisfy local mass continuity. Despite the computational overhead required in the restriction and prolongation phases of the multigrid cycling, the superior convergence results in reduced overall CPU time. The numerical scheme and selected results of several validation flows are presented. Finally, the procedure is applied to study the flowfield in a side-inlet dump combustor and twin jet impingement from a simulated aircraft fuselage.

Smith, K. M.↗
NASA NTRS2006

Diffusion Limited Supercritical Water Oxidation (SCWO) in Microgravity Environments

Tests designed to quantify the gravitational effects on thermal mixing and reactant injection in a Supercritical Water Oxidation (SCWO) reactor have recently been performed in the Zero Gravity Facility (ZGF) at NASA s Glenn Research Center. An artificial waste stream, comprising aqueous mixtures of methanol, was pressurized to approximately 250 atm and then heated to 450 C. After uniform temperatures in the reactor were verified, a controlled injection of air was initiated through a specially designed injector to simulate diffusion limited reactions typical in most continuous flow reactors. Results from a thermal mapping of the reaction zone in both 1-g and 0-g environments are compared. Additionally, results of a numerical model of the test configuration are presented to illustrate first order effects on reactant mixing and thermal transport in the absence of gravity.

Hicks, M. C.↗
NASA NTRS2015

Advanced Intensity-Modulation Continuous-Wave Lidar Techniques for ASCENDS O2 Column Measurements

Global atmospheric carbon dioxide (CO2) measurements for the NASA Active Sensing of CO2 Emissions over Nights, Days, and Seasons (ASCENDS) space mission are critical for improving our understanding of global CO2 sources and sinks. Advanced Intensity- Modulated Continuous-Wave (IM-CW) lidar techniques are investigated as a means of facilitating CO2 measurements from space to meet the ASCENDS measurement requirements. In recent numerical, laboratory and flight experiments we have successfully used the Binary Phase Shift Keying (BPSK) modulation technique to uniquely discriminate surface lidar returns from intermediate aerosol and cloud contamination. We demonstrate the utility of BPSK to eliminate sidelobes in the range profile as a means of making Integrated Path Differential Absorption (IPDA) column CO2 measurements in the presence of optically thin clouds, thereby eliminating the need to correct for sidelobe bias errors caused by the clouds. Furthermore, high accuracy and precision ranging to the surface as well as to the top of intermediate cloud layers, which is a requirement for the inversion of column CO2 number density measurements to column CO2 mixing ratios, has been demonstrated using new hyperfine interpolation techniques that takes advantage of the periodicity of the modulation waveforms. This approach works well for both BPSK and linear swept-frequency modulation techniques. The BPSK technique under investigation has excellent auto-correlation properties while possessing a finite bandwidth. A comparison of BPSK and linear swept-frequency is also discussed in this paper. These results are extended to include Richardson-Lucy deconvolution techniques to extend the resolution of the lidar beyond that implied by limit of the bandwidth of the modulation, where it is shown useful for making tree canopy measurements.

Campbell, Joel F.↗
NASA NTRS2019

initMIP-Antarctica: an Ice Sheet Model Initialization Experiment of ISMIP6

Ice sheet numerical modeling is an important tool to estimate the dynamic contribution of the Antarctic ice sheet to sea level rise over the coming centuries. The influence of initial conditions on ice sheet model simulations, however, is still unclear. To better understand this influence, an initial state intercomparison exercise (initMIP) has been developed to compare, evaluate, and improve initialization procedures and estimate their impact on century-scale simulations. initMIP is the first set of experiments of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), which is the primary Coupled Model Intercomparison Project Phase 6 (CMIP6) activity focusing on the Greenland and Antarctic ice sheets. Following initMIP-Greenland, initMIP-Antarctica has been designed to explore uncertainties associated with model initialization and spin-up and to evaluate the impact of changes in external forcings. Starting from the state of the Antarctic ice sheet at the end of the initialization procedure, three forward experiments are each run for 100 years: a control run, a run with a surface mass balance anomaly, and a run with a basal melting anomaly beneath floating ice. This study presents the results of initMIP-Antarctica from 25 simulations performed by 16 international modeling groups. The submitted results use different initial conditions and initialization methods, as well as ice flow model parameters and reference external forcings. We find a good agreement among model responses to the surface mass balance anomaly but large variations in responses to the basal melting anomaly. These variations can be attributed to differences in the extent of ice shelves and their upstream tributaries, the numerical treatment of grounding line, and the initial ocean conditions applied, suggesting that ongoing efforts to better represent ice shelves in continental-scale models should continue.

ice sheet model↗
NASA NTRS2022

Ten Years of VIIRS On-Orbit Geolocation Calibration and Performance

The first innovative Visible Infrared Imaging Radiometer Suite (VIIRS) sensor aboard the Suomi National Polar-orbiting Partnership (SNPP) satellite has been in operation for 10 years since its launch on 28 October 2011. The second VIIRS sensor aboard the first Join Polar Satellite System (JPSS-1) satellite has been in operation for 4 years since its launch on 18 November 2017, which became NOAA-20. Well-geolocated and radiometrically calibrated Level-1 sensor data records (SDRs) from VIIRS are crucial to numerical weather prediction (NWP) and Level-2+ environmental data record (EDR) algorithms and products. The high quality of Level-2+ EDRs is a requirement for the continuity of NASA Earth science data records (ESDRs) and climate data records (CDRs), one of the two objectives of the SNPP mission and one of the three elements in the JPSS mission objective. The other objective of the SNPP mission is risk reduction for the follow-on JPSS missions. This paper summarizes the on-orbit geolocation calibration and validation (Cal/Val) activities for both VIIRS sensors onboard SNPP and NOAA-20 in the past 10 years. These activities include nominal geolocation Cal/Val activities, risk reduction activities, and improvements for the on-orbit VIIRS sensor operations. After these activities, sub-pixel geolocation accuracy is achieved. Nadir equivalent geolocation uncertainty is generally within 75 m (1-σ), or 20% imagery band pixels, in either the along-scan or along-track direction for both SNPP and NOAA-20 VIIRS sensors. The worst 16-day measured geolocation errors (radial, 3-σ) are 280 m and 267 m, respectively, in the latest SNPP and NOAA-20 VIIRS data collections, which are better than the required accuracy of 375 m (radial, 3-σ). The risk reduction activities also improved VIIRS builds for JPSS-3 and JPSS-4 satellites, and provide lessons learned for other VIIRS-like sensor builds.

VIIRS↗
NASA NTRS2022

Operational maize yield forecasts for Sub-Saharan Africa using Earth observation data and machine learning

Food insecurity continues to grow in Sub-Saharan Africa (SSA). In 2019, chronically malnourished people numbered nearly 240 million, or 20% of the population in SSA. Globally, numerous efforts have been made to anticipate potential droughts, crop conditions, and food shortages to foster improved food insecurity early warning and risk management. To support this goal, we develop an Earth Observation (EO) and machine-learning-based operational, subnational maize yield forecast system and evaluate its out-of-sample forecast skills during the growing seasons for Kenya, Somalia, Malawi, and Burkina Faso. In general, forecast skills improve substantially during the vegetative growth period (VP) and gradually during the reproductive development period (RP). Thus, mid-season assessment can provide effective early warning months before harvest. Skillful forecasts (Nash Sutcliffe Efficiency (NSE) > 0.6 and Mean Absolute Percentage Error (MAPE) < 20%) appear approximately two dekads after the VP; for example, they appear in May in Kenya and Somalia, January in Malawi, and July in Burkina Faso.

Donghoon Lee↗
NASA NTRS1987

A numerical method for the solution of three dimensional, incompressible, viscous flows over slender bodies

A marching iterative method for the solution of the three dimensional, incompressibhle, steady and parabolized Navier-Stokes equations is described. The equations are written in primitive variables and discretized in general axisymmetric orthogonal coordinate systems. The coupled set of finite-difference equations are solved without any splitting or factorization errors. Moreover, the continuity equation and the two crossflow momentum equations are exactly satisfied at every step of the iterative process. The solution scheme is equivalent to the solution of one Poisson equation by the Successive Plane Over Relaxation method and has good convergence properties. Other existing solution methods resemble a Jacobi-type iterative scheme and therefore are less efficient. Numerical experiments include the laminar, incompressible flow over prolate spheroids at incidence.

Rosenfeld, Moshe↗
DOE OSTI2023

A numerical solver for investigating the space charge effect on the electric field in liquid argon time projection chambers

This paper reports the development of a numerical solver aimed to simulate the interaction between the space charge (i.e. ions) distribution and the electric field in liquid argon time projection chamber (LArTPC) detectors. The ion transport equation is solved by a time-accurate, cell-centered finite volume method and the electric potential equation by a continuous finite element method. The electric potential equation updates the electric field which provides the drift velocity to the ion transport equation. The ion transport equation updates the space charge density distribution which appears as the source term in the electric potential equation. The interaction between the space charge distribution and the electric field is numerically simulated within each physical time step. The convective velocity in the ion transport equation can include the background flow velocity in addition to the electric drift velocity. The numerical solver has been parallelized using the Message Passing Interface (MPI) library. Numerical tests show and verify the capability and accuracy of the current numerical solver. It is planned that the developed numerical solver, together with a Computational Fluid Dynamics (CFD) package which provides the flow velocity field, can be used to investigate the space charge effect on the electric field in large-scale particle detectors.

42 ENGINEERING↗
NASA NTRS2011

Climate Change, Extreme Weather Events, and Fungal Disease Emergence and Spread

Empirical evidence from multiple sources show the Earth has been warming since the late 19th century. More recently, evidence for this warming trend is strongly supported by satellite data since the late 1970s from the cryosphere, atmosphere, oceans, and land that confirms increasing temperature trends and their consequences (e.g., reduced Arctic sea ice, rising sea level, ice sheet mass loss, etc.). At the same time, satellite observations of the Sun show remarkably stable solar cycles since the late 1970s, when direct observations of the Sun's total solar irradiance began. Numerical simulation models, driven in part by assimilated satellite data, suggest that future-warming trends will lead to not only a warmer planet, but also a wetter and drier climate depending upon location in a fashion consistent with large-scale atmospheric processes. Continued global warming poses new opportunities for the emergence and spread of fungal disease, as climate systems change at regional and global scales, and as animal and plant species move into new niches. Our contribution to this proceedings is organized thus: First, we review empirical evidence for a warming Earth. Second, we show the Sun is not responsible for the observed warming. Third, we review numerical simulation modeling results that project these trends into the future, describing the projected abiotic environment of our planet in the next 40 to 50 years. Fourth, we illustrate how Rift Valley fever outbreaks have been linked to climate, enabling a better understanding of the dynamics of these diseases, and how this has led to the development of an operational predictive outbreak model for this disease in Africa. Fifth, We project how this experience may be applicable to predicting outbreaks of fungal pathogens in a warming world. Lastly, we describe an example of changing species ranges due to climate change, resulting from recent warming in the Andes and associated glacier melt that has enabled amphibians to colonize higher elevation lakes, only to be followed shortly by the emergence of fungal disease in the new habitats.

Tucker, Compton J.↗
DOE OSTI2022

Electronic energy loss and ion velocity correlation effects in track production in swift-ion-irradiated LiNbO3: A quantitative assessment between structural damage morphology and energy deposition

The primary motivation for studying how irradiation modifies the structures and properties of solid materials involves the understanding of undesirable phenomena, including irradiation-induced degradation of components in nuclear reactors and space exploration, and beneficial applications, including material performance tailoring through ion beam modification and defect engineering. In this work, the formation mechanism of latent tracks with different damage morphologies in LiNbO 3 crystals under 0.09–6.17 MeV/u ion irradiation with an electronic energy loss from 2.6–13.2 keV/nm is analyzed by experimental characterizations and numerical calculations. Irradiation-induced damage is preliminarily evaluated via the prism coupling technique to analyze the correlation between the dark-mode spectra and energy loss profiles of irradiated regions. Under the irradiation conditions of different ion velocities and electronic energy losses, different damage morphologies, from individual spherical defects to discontinuous and continuous tracks, are experimentally characterized. During ion penetration process, the ion velocity determines the spatiotemporal distribution of deposited irradiation energy induced by electronic energy loss, meaning that the two essential factors including electronic energy loss and ion velocity co-affect the track damage. The inelastic thermal spike model is used to numerically calculate the spatiotemporal evolutions of energy deposition and the corresponding atomic temperature under different irradiation conditions, and a quantitative relationship is proposed by comparison with corresponding experimentally observed track damage morphologies. Additionally, the obtained quantitative relationship between irradiation conditions and track damage provides deep insight and guidance for understanding the damage behavior of crystal materials in extreme radiation environments and selecting irradiation parameters, including ion species and energies, for ion beam technique application in atomic-level defect manipulation, material modification, and micro/nanofabrication.

36 MATERIALS SCIENCE↗
NASA NTRS2000

Electromagnetic Field Penetration Studies

A numerical method is presented to determine electromagnetic shielding effectiveness of rectangular enclosure with apertures on its wall used for input and output connections, control panels, visual-access windows, ventilation panels, etc. Expressing EM fields in terms of cavity Green's function inside the enclosure and the free space Green's function outside the enclosure, integral equations with aperture tangential electric fields as unknown variables are obtained by enforcing the continuity of tangential electric and magnetic fields across the apertures. Using the Method of Moments, the integral equations are solved for unknown aperture fields. From these aperture fields, the EM field inside a rectangular enclosure due to external electromagnetic sources are determined. Numerical results on electric field shielding of a rectangular cavity with a thin rectangular slot obtained using the present method are compared with the results obtained using simple transmission line technique for code validation. The present technique is applied to determine field penetration inside a Boeing-757 by approximating its passenger cabin as a rectangular cavity filled with a homogeneous medium and its passenger windows by rectangular apertures. Preliminary results for, two windows, one on each side of fuselage were considered. Numerical results for Boeing-757 at frequencies 26 MHz, 171-175 MHz, and 428-432 MHz are presented.

Deshpande, M.D.↗
NASA NTRS2009

Numerical Evaluation of the "Dual-Kernel Counter-flow" Matric Convolution Integral that Arises in Discrete/Continuous (D/C) Control Theory

Discrete/Continuous (D/C) control theory is a new generalized theory of discrete-time control that expands the concept of conventional (exact) discrete-time control to create a framework for design and implementation of discretetime control systems that include a continuous-time command function generator so that actuator commands need not be constant between control decisions, but can be more generally defined and implemented as functions that vary with time across sample period. Because the plant/control system construct contains two linear subsystems arranged in tandem, a novel dual-kernel counter-flow convolution integral appears in the formulation. As part of the D/C system design and implementation process, numerical evaluation of that integral over the sample period is required. Three fundamentally different evaluation methods and associated algorithms are derived for the constant-coefficient case. Numerical results are matched against three available examples that have closed-form solutions.

Nixon, Douglas D.↗