Search NASA⌕ Search

SEARCH · Search NASA

Results for “Physics Modeling”

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 757 records · Page 42

An Initial Exploration of Improved Numerics within the Guidelines of the Negative Spalart-Allmaras Turbulence Model

A simple modification to the negative Spalart-Allmaras turbulence model is suggested so that when the turbulence working variable, ~v, is negative, diagonal dominance is increased, as is the tendency for the time-advancement scheme to push e toward positive values. Owing to the fact that the modification is only active when ~v is less than zero, the physical model is left unchanged. Using the proposed modification with a strong implicit solver based on Newton's method, convergence rates can be somewhat improved, with typical reductions in iterations and computer time on the order of 15-50%. The benefits are realized primarily when second- or higher-order accuracy is used for discretizing the convective terms in the turbulence model because of large overshoots that can occur with these schemes at the edges of boundary layers and wakes. For flowfields with few regions of negative ~v, or on very fine meshes where ~v is always greater than zero, little or no benefits should be expected.

Anderson, William K.↗

Modeling of Optoelectronic Devices

Ultrafast modulation of semiconductor quantum well (QW) laser is of technological importance for information technology. Improvement by order(s) of magnitude in data transfer rate is possible as terahertz (THz) radiation is available for heating the laser at picosecond time scale. Optical gain modulation in the QW is achieved via temperature modulation of electron-hole plasma (EHP). Applications include free-space THz communication, optical switching, and pulse generation. The EHP in the semiconductor QW is described with a two-band model. Semiconductor Bloch equations with many-body effects are used to derive a hydrodynamical model for the active QW region. Because of ultrafast carrier-carrier scatterings in the order of 50 fs, EHP follows quasiequilibrium Fermi-Dirac distributions and THz field interacts incoherently with it. Carrier-longitudinal optical (LO) phonon scatterings and coherent laser-EHP interaction are treated microscopically in our physical model. A set of hydrodynamical equations for plasma density, temperature, and laser envelop amplitude are derived and Runge-Kutta method is adopted for numerical simulation. A typical 8 nm GaAs/Al(0.3)Ga(0.7) As single QW at 300 K is used. Additional information is contained in the original extended abstract.

Li, Jian-Zhong↗

ECLSS Four-Bed CO2 Scrubber Digital Twin

The ECLSS Digital Twin is a cloud-based simulation of the Four-Bed CO2 Scrubber, currently one of the primary means of removing carbon dioxide onboard the International Space Station, that operationalizes SME-developed multi-physics models and replicates conditions of the actual hardware in near real-time. In addition to providing insight into the system’s performance, it will enable prognostics, diagnostics, predictive maintenance, and simulated off-nominal scenarios. By leveraging digital representations, projects can save resources and gain a better understanding of their physical systems with the goal of building robust and reliable hardware for future deep space exploration. The use of data infrastructure in the cloud allows for streamlined analysis and visualization on a much larger scale than is possible with current tools.

Jared Wilson↗

Limitations in Representation of Physical Processes Prevents Successful Simulation of PM2.5 During KORUS-AQ

High levels of fine particulate matter (PM2.5) pollution in East Asia often exceed local air quality standards. Observations from the Korea United States-Air Quality (KORUS-AQ) field campaign in May and June 2016 showed that development of extreme pollution (haze) occurred through a combination of long-range transport and favorable meteorological conditions that enhanced local production of PM2.5. Atmospheric models often have difficulty simulating PM2.5 chemical composition during haze, which is of concern for the development of successful control measures. We use observations from KORUS-AQ to examine the ability of the GEOS-Chem chemical transport model to simulate PM2.5 composition throughout the campaign and identify the mechanisms driving the pollution event. At the surface, the model underestimates sulfate by -64% but overestimates nitrate by +36%. The largest underestimate in sulfate occurs during the pollution event, where models typically struggle to generate elevated sulfate concentrations due to missing heterogeneous chemistry in aerosol liquid water in the polluted boundary layer. Hourly surface observations show that the model nitrate bias is driven by an overestimation of the nighttime peak. In the model, nitrate formation is limited by the supply of nitric acid, which is biased by +100% against aircraft observations. We hypothesize that this is due to a large missing sink, which we implement here as a factor of five increase in dry deposition. We show that the resulting increased deposition velocity is consistent with observations of total nitrate as a function of photochemical age. The model does not account for factors such as the urban heat island effect or the heterogeneity of the built-up urban landscape resulting in insufficient model turbulence and surface area over the study area that likely results in insufficient dry deposition. Other species such as NH3 could be similarly affected but were not measured during the campaign. Nighttime production of nitrate is driven by NO2 hydrolysis in the model, while observations show that unexpectedly elevated nighttime ozone (not present in the model) should result in N2O5 hydrolysis as the primary pathway. The model is unable to represent nighttime ozone due to an overly rapid collapse of the afternoon mixed layer and excessive titration by NO. We attribute this to missing nighttime heating driving deeper nocturnal mixing that would be expected to occur in a city like Seoul. This urban heating is not considered in air quality models run at large enough scales to treat both local chemistry and long-range transport. Key model failures in simulating nitrate, mainly overestimated daytime nitric acid, incorrect representation of nighttime chemistry, and an overly shallow and insufficiently turbulent nighttime mixed layer, exacerbate the model’s inability to simulate the buildup of PM2.5 during haze pollution. To address the underestimate in sulfate most evident during the haze event, heterogeneous aerosol uptake of SO2 is added to the model which previously only considered aqueous production of sulfate from SO2 in cloud water. Implementing a simple parameterization of this chemistry improves the model abundance of sulfate but degrades the SO2 simulation implying that emissions are underestimated. We find that improving model simulations of sulfate has direct relevance to determining local vs. transboundary contributions to PM2.5. During the haze pollution event, the inclusion of heterogeneous aerosol uptake of SO2 decreases the fraction of PM2.5 attributable to long-range transport from 66% to 54%. Locally-produced sulfate increased from 1% to 25% of locally-produced PM2.5, implying that local emissions controls would have a larger effect than previously thought. However, this additional uptake of SO2 is coupled to the model nitrate prediction which affects the aerosol liquid water abundance and chemistry driving sulfate-nitrate-ammonium partitioning. An additional simulation of the haze pollution with heterogeneous uptake of SO2 to aerosol and simple improvements to the model nitrate simulation results in 30% less sulfate due to 40% less nitrate and aerosol water, and results in an underestimate of sulfate during the haze event. Future studies need to better consider the impact of model physical processes such as dry deposition and nighttime boundary layer mixing on the simulation of nitrate and the effect of improved nitrate simulations on the overall simulation of secondary inorganic aerosol (sulfate+nitrate+ammonium) in East Asia. Foreign emissions are rapidly changing, increasing the need to understand the impact of local emissions on PM2.5 in South Korea to ensure continued air quality improvements.

Katherine R. Travis↗

Influence of Sub-grid-Scale Isentropic Transports on McRAS Evaluations using ARM-CART SCM Datasets

In GCM-physics evaluations with the currently available ARM-CART SCM datasets, McRAS produced very similar character of near surface errors of simulated temperature and humidity containing typically warm and moist biases near the surface and cold and dry biases aloft. We argued it must have a common cause presumably rooted in the model physics. Lack of vertical adjustment of horizontal transport was thought to be a plausible source. Clearly, debarring such a freedom would force the incoming air to diffuse into the grid-cell which would naturally bias the surface air to become warm and moist while the upper air becomes cold and dry, a characteristic feature of McRAS biases. Since, the errors were significantly larger in the two winter cases that contain potentially more intense episodes of cold and warm advective transports, it further reaffirmed our argument and provided additional motivation to introduce the corrections. When the horizontal advective transports were suitably modified to allow rising and/or sinking following isentropic pathways of subgrid scale motions, the outcome was to cool and dry (or warm and moisten) the lower (or upper) levels. Ever, crude approximations invoking such a correction reduced the temperature and humidity biases considerably. The tests were performed on all the available ARM-CART SCM cases with consistent outcome. With the isentropic corrections implemented through two different numerical approximations, virtually similar benefits were derived further confirming the robustness of our inferences. These results suggest the need for insentropic advective transport adjustment in a GCM due to subgrid scale motions.

Sud, Y. C.↗

Effects of Forcing Differences and Initial Conditions on Inter-Model Agreement in the VolMIP Volc-Pinatubo-Full Experiment

This paper provides initial results from a multi-model ensemble analysis based on the volc-pinatubo-full experiment performed within the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP) as part of the sixth phase of the Coupled Model Intercomparison Project (CMIP6). The volc-pinatubo-full experiment is based on an ensemble of volcanic forcing-only climate simulations with the same volcanic aerosol dataset across the participating models (the 1991–1993 Pinatubo period from the CMIP6-GloSSAC dataset). The simulations are conducted within an idealized experimental design where initial states are sampled consistently across models from the CMIP6-piControl simulation providing unperturbed preindustrial background conditions. The multi-model ensemble includes output from an initial set of six participating Earth system models (CanESM5, GISS-E2.1-G, IPSL-CM6A-LR, MIROC-E2SL, MPI-ESM1.2-LR and UKESM1). The results show overall good agreement between the different models on the global and hemispheric scales concerning the surface climate responses, thus demonstrating the overall effectiveness of VolMIP's experimental design. However, small yet significant inter-model discrepancies are found in radiative fluxes, especially in the tropics, that preliminary analyses link with minor differences in forcing implementation; model physics, notably aerosol–radiation interactions; the simulation and sampling of El Niño–Southern Oscillation (ENSO); and, possibly, the simulation of climate feedbacks operating in the tropics. We discuss the volc-pinatubo-full protocol and highlight the advantages of volcanic forcing experiments defined within a carefully designed protocol with respect to emerging modelling approaches based on large ensemble transient simulations. We identify how the VolMIP strategy could be improved in future phases of the initiative to ensure a cleaner sampling protocol with greater focus on the evolving state of ENSO in the pre-eruption period.

metadata↗

Finite element modeling of electropolishing of radio frequency dipole Nb crab cavity in hydrofluoric-sulfuric acid electrolyte

The superior performance of niobium superconducting radio frequency (SRF) cavities is achieved through meticulous surface treatment, notably via chemical electropolishing, ensuring exceptionally smooth surfaces. While this technique has been extensively employed for cylindrically symmetric structures, addressing more intricate geometries poses significant challenges in achieving uniform polishing and controlled material removal, especially when moving away from retractable cathodes. This presents a multifaceted electrochemical, thermal, and fluid dynamics conundrum. A prime example is the 197 MHz radio frequency dipole (RFD) crabbing cavity proposed for the Electron Ion Collider (EIC) project, exemplifying such complex structures. Our groundbreaking work integrates the localized oxide thickness variation, considering its impact on current distribution and Joule heating, within the framework of multi-physics modeling using the COMSOL® simulation suite. This approach was applied to a comprehensive model of the RFD cavity, allowing us to investigate optimal external cooling water flow conditions necessary for achieving desirable outcomes. In conclusion, this illustrates the increasing utility of such multi-physics codes to greatly aid in designing solutions to challenging processing requirements for increasingly complex accelerator cavities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

New Third-Order Moments for the PBL

Turbulent convection is inherently non-local and a primary condition for a successful treatment of the PBL is a reliable model of non-locality. In the dynamic equations governing the convective flux, turbulent kinetic energy, etc, non-locality enters through the third-order moments, TOMs. Since the simplest form, the so-called down gradient approximation (DGA , severely underestimates the TOMs (by up to an order of magnitude), a more physical model is needed. In 1994, an analytical model was presented which was derived directly from the dynamical equations for the TOMs. It considerably improved the DGA but was a bit cumbersome to use. Here, we present a new analytic expression for the TOMs which is considerably simpler than the 1994 expression and which at the same time yields a much better fit to the LES data.

Canuto, V. M.↗

A High-Fidelity Model of the Peach Bottom 2 Turbine-Trip Benchmark Using VERA

This work presents a high-fidelity simulation of the Peach Bottom turbine trip (PBTT) benchmark using the Virtual Environment for Reactor Applications (VERA), a multiphysics reactor modeling tool developed by the U.S. Department of Energy’s Consortium for Advanced Simulation of Light Water Reactors energy innovation hub. The PBTT benchmark, based on a 1977 transient event at the end of cycle 2 in a General Electric Type-4 boiling water reactor (BWR), is a critical test case for validating core physics models with thermal feedback during rapid reactivity events. VERA was employed to perform end-to-end, pin-resolved simulations from conditions at the beginning of cycle 1 through the turbine-trip transient, incorporating detailed neutron transport, fuel depletion, and subchannel thermal hydraulics. The simulation reproduced key benchmark observables with high accuracy: the peak power excursion occurred at 0.75 s, matching the scram time and closely aligning with the benchmark average of 0.742 s; the simulated maximum power spike was approximately 7600 MW, which is within 3% of the benchmark average of 7400 MW; and void-collapse dynamics were consistent with benchmark expectations. Reactivity predictions during cycles 1 and 2 remained within 1500 pcm and 400 pcm of criticality, respectively. These results confirm VERA’s ability to model complex coupled neutronic and thermal hydraulic behavior in a BWR turbine-trip transient, which will support its use in future studies of modeling dryout, fuel performance, and uncertainty quantification for transients of this type.

BWR↗

New Third-Order Moments for the CBL

Turbulent convection is an inherently non-local phenomenon and a primary condition for a successful treatment of the CBL (convective boundary layer) is a reliable model of non-locality. In the dynamic equations governing the convective flux, the turbulent kinetic energy, etc., non-locality is represented by the third-order moments, TOMs. Since the simplest form, the so-called down gradient approximation (DGA), severely underestimates the TOMs (up to an order of magnitude), a more physical model is needed. In 1994, an analytical model was presented which was derived directly from the dynamical equations for the TOMs. It considerably improved the DGA but was a bit cumbersome to use and, more importantly, it was based on the quasi-normal (QN) approximation for the fourth-order moments. Here, we present a new analytic expression for the TOMs which is structurally simpler than the 1994 expression and which avoids the QN approximation. The resulting fit to the LES data is superior to that of the 1994 model.

Canuto, V. M.↗

Clumped-Isotope Constraint on Upper-Tropospheric Cooling During the Last Glacial Maximum

Ice cores and other paleotemperature proxies, together with general circulation models, have provided information on past surface temperatures and the atmosphere's composition in different climates. Little is known, however, about past temperatures at high altitudes, which play a crucial role in Earth's radiative energy budget. Paleoclimate records at high-altitude sites are sparse, and the few that are available show poor agreement with climate model predictions. These disagreements could be due to insufficient spatial coverage, spatiotemporal biases, or model physics; new records that can mitigate or avoid these uncertainties are needed. Here, we constrain the change in upper-tropospheric temperature at the global scale during the Last Glacial Maximum (LGM) using the clumped-isotope composition of molecular oxygen trapped in polar ice cores. Aided by global three-dimensional chemical transport modeling, we exploit the intrinsic temperature sensitivity of the clumped-isotope composition of atmospheric oxygen to infer that the upper troposphere (effective mean altitude 10–11 km) was 6–9°C cooler during the LGM than during the late preindustrial Holocene. A complementary energy balance approach supports a minor or negligible steepening of atmospheric lapse rates during the LGM, which is consistent with a range of climate model simulations. Proxy-model disagreements with other high-altitude records may stem from inaccuracies in regional hydroclimate simulation, possibly related to land-atmosphere feedbacks.

clumped isotope↗

Simulated lightning test shuttle .03 scale model

Lightning Attach Point tests were conducted for the space shuttle launch configuration (Orbiter, External Tank and Solid Rocket Boosters). A series of 250 long spark tests (15 to 20 foot sparks) determined that the orbiter may be struck on the nose, windshield brow, tail and wingtips during launch but not on the main engine nozzles which have been shown to be vulnerable to lightning damage. The orbiter main engine and SRB exhaust plumes were simulated electrically with physical models coated with graded resistance paints. The tests showed that the exhaust plumes from the SRB provide additional protection for the main engine nozzles. However, the tests showed that the Orbiter Thermal Protection System (TPS), which has also been shown to be vulnerable to lightning damage, may be struck during launch. Therefore further work is indicated in the areas of swept stroke studies on the model and on TPS panels. Further attach point testing is also indicated on the free-flying orbiter. Photographs of the test setup are shown.

Clifford, D. W.↗

Numerical Modeling Of Two-Phase, Reactive Flows

Capabilities of several computational schemes combined and enhanced. Advanced Rocket Injector-Combustor Code is computer program for numerical simulation of two-phase flows of chemically reacting or nonreacting components. Applied to sprays of fuel in injectors and combustors, flows in heat exchangers, evaporators, and mass-diffusion systems; and two-dimensional or axisymmetric flows containing liquid or solid suspensions. Used in research for development and verification of various physical models.

Liang, Pak-Yan↗

The Electric Propulsion Interactions Code (EPIC): A Member of the NASA Space Environment and Effects Program (SEE) Toolset

Science Applications International Corporation is currently developing the Electric Propulsion Interactions Code, EPIC, as part of a project sponsored by the Space Environments and Effects Program at NASA Marshall Space Flight Center. Now in its second year of development, EPIC is an interactive computer toolset that allows the construction of a 3-D spacecraft model, and the assessment of a variety of interactions between its subsystems and the plume from an electric thruster. This paper reports on the progress of EPZC including the recently added ability to exchange results the NASA Charging Analyzer Program, Nascap-2k. The capability greatly enhances EPIC's range of applicability. Expansion of the toolset's various physics models proceeds in parallel with the overall development of the software. Also presented are recent upgrades of the elastic scattering algorithm in the electric propulsion Plume Tool. These upgrades are motivated by the need to assess the effects of elastically scattered ions on the SIC for ion beam energies that exceed loo0 eV. Such energy levels are expected in future high-power (>10 kW) ion propulsion systems empowered by nuclear sources.

Mikellides, Ioannis G.↗

A Thermal Runaway Failure Model for Low-Voltage BME Ceramic Capacitors with Defects

Reliability of base metal electrode (BME) multilayer ceramic capacitors (MLCCs) that until recently were used mostly in commercial applications, have been improved substantially by using new materials and processes. Currently, the inception of intrinsic wear-out failures in high quality capacitors became much greater than the mission duration in most high-reliability applications. However, in capacitors with defects degradation processes might accelerate substantially and cause infant mortality failures. In this work, a physical model that relates the presence of defects to reduction of breakdown voltages and decreasing times to failure has been suggested. The effect of the defect size has been analyzed using a thermal runaway model of failures. Adequacy of highly accelerated life testing (HALT) to predict reliability at normal operating conditions and limitations of voltage acceleration are considered. The applicability of the model to BME capacitors with cracks is discussed and validated experimentally.

reliability↗

Electron number density profiles for the Aeroassist Flight Experiment

The basic features of a Microwave Reflectometer Ionization Sensor (MRIS) as designed for utilization on the Aeroassist Flight Experiment are described. A parametric study of the effects of trajectory and unknowns in the thermochemical nonequilibrium models for translational and vibrational-electronic energy exchange rates, reaction rates, and the average electronic excitation level of atoms is performed to illustrate how the MRIS data may be employed for code validation. This study, implemented with program LAURA, shows a particular sensitivity of the onset and severity of an electron avalanche phenomena associated with changes in these physical models.

Greendyke, Robert B.↗

Investigation for improving Global Positioning System (GPS) orbits using a discrete sequential estimator and stochastic models of selected physical processes

GEODYNII is a conventional batch least-squares differential corrector computer program with deterministic models of the physical environment. Conventional algorithms were used to process differenced phase and pseudorange data to determine eight-day Global Positioning system (GPS) orbits with several meter accuracy. However, random physical processes drive the errors whose magnitudes prevent improving the GPS orbit accuracy. To improve the orbit accuracy, these random processes should be modeled stochastically. The conventional batch least-squares algorithm cannot accommodate stochastic models, only a stochastic estimation algorithm is suitable, such as a sequential filter/smoother. Also, GEODYNII cannot currently model the correlation among data values. Differenced pseudorange, and especially differenced phase, are precise data types that can be used to improve the GPS orbit precision. To overcome these limitations and improve the accuracy of GPS orbits computed using GEODYNII, we proposed to develop a sequential stochastic filter/smoother processor by using GEODYNII as a type of trajectory preprocessor. Our proposed processor is now completed. It contains a correlated double difference range processing capability, first order Gauss Markov models for the solar radiation pressure scale coefficient and y-bias acceleration, and a random walk model for the tropospheric refraction correction. The development approach was to interface the standard GEODYNII output files (measurement partials and variationals) with software modules containing the stochastic estimator, the stochastic models, and a double differenced phase range processing routine. Thus, no modifications to the original GEODYNII software were required. A schematic of the development is shown. The observational data are edited in the preprocessor and the data are passed to GEODYNII as one of its standard data types. A reference orbit is determined using GEODYNII as a batch least-squares processor and the GEODYNII measurement partial (FTN90) and variational (FTN80, V-matrix) files are generated. These two files along with a control statement file and a satellite identification and mass file are passed to the filter/smoother to estimate time-varying parameter states at each epoch, improved satellite initial elements, and improved estimates of constant parameters.

Goad, Clyde C.↗

Efficient First-Order Algorithms for Large-Scale, Non-Smooth Maximum Entropy Models with Application to Wildfire Science

Maximum entropy (MaxEnt) models are a class of statistical models that use the maximum entropy principle to estimate probability distributions from data. Due to the size of modern data sets, MaxEnt models need efficient optimization algorithms to scale well for big data applications. State-of-the-art algorithms for MaxEnt models, however, were not originally designed to handle big data sets; these algorithms either rely on technical devices that may yield unreliable numerical results, scale poorly, or require smoothness assumptions that many practical MaxEnt models lack. In this paper, we present novel optimization algorithms that overcome the shortcomings of state-of-the-art algorithms for training large-scale, non-smooth MaxEnt models. Our proposed first-order algorithms leverage the Kullback–Leibler divergence to train large-scale and non-smooth MaxEnt models efficiently. For MaxEnt models with discrete probability distribution of n elements built from samples, each containing m features, the stepsize parameter estimation and iterations in our algorithms scale on the order of O(mn) operations and can be trivially parallelized. Moreover, the strong ℓ1 convexity of the Kullback–Leibler divergence allows for larger stepsize parameters, thereby speeding up the convergence rate of our algorithms. To illustrate the efficiency of our novel algorithms, we consider the problem of estimating probabilities of fire occurrences as a function of ecological features in the Western US MTBS-Interagency wildfire data set. Our numerical results show that our algorithms outperform the state of the art by one order of magnitude and yield results that agree with physical models of wildfire occurrence and previous statistical analyses of wildfire drivers.

Physics↗