Search NASA⌕ Search

SEARCH · Search NASA

Results for “discrete phase model”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Sensitivity of Single-Scattering Properties to Precipitation Particle Meltwater Geometry for GPM Passive Microwave and Radar Remote Sensing Applications

The objective of the current work is to describe the microwave single-scattering properties of partially melted ice-phase precipitation particles using physically-based computational methods, taking into full account the varied geometries of the initially “dry” ice particles and the distributions of liquid water that develop during the melting process. Ultimately, the bulk properties of ensembles of these partially-melted particles will be included in “scattering tables” to support radar and combined radar-radiometer remote sensing of precipitation. A meshless Lagrangian melting procedure (snowMELT) based on smoothed-particle hydrodynamics is applied to both spherical and finely-structured ice particle models to describe the full evolution of the particles from dry ice particles to liquid drops. The melting of spherical ice particles using snowMELT is compared to an alternative continuum physics model to validate melt times and internal thermodynamics. The discrete dipole approximation is then utilized to calculate the single-scattering properties of different mixed-phase particles throughout the melting process. The sensitivities of particle single-scattering properties to thermodynamic and hydrodynamic assumptions in snowMELT are explored, and the implications for combined radar-radiometer precipitation remote sensing from GPM are discussed.

William S Olson↗
DOE OSTI2026

Nonequilibrium universality of the nonreciprocally coupled 𝑂⁡(𝑛 1 ) × 𝑂⁡(𝑛 2 ) model

Nonequilibrium dynamics play an important role in all contexts of physics, both classical and quantum as well as living and nonliving, so it is crucial to develop a foundational understanding of nonequilibrium phase transitions. In this work we investigate an important class of nonequilibrium dynamics in the form of nonreciprocal interactions. In particular we study how nonreciprocal coupling between two 𝑂⁡(𝑛𝑖) order parameters (with 𝑖 = 1,2) affects the universality at a multicritical point, extending the analysis of J. T. Young et al. [Phys. Rev. X 10, 011039 (2020)], which considered the case 𝑛 1 = 𝑛 2 = 1, i.e., a ℤ 2 × ℤ 2 model. We show that nonequilibrium fixed points (NEFPs) emerge for a broad range of 𝑛 1 ,𝑛 2 and exhibit intrinsically nonequilibrium critical phenomena, namely a violation of fluctuation-dissipation relations at all scales and underdamped oscillations near criticality in contrast to the overdamped relaxational dynamics of the corresponding equilibrium models. Furthermore, the NEFPs exhibit an emergent discrete scale invariance in certain physically relevant regimes of 𝑛 1 ,𝑛 2 , but not others, depending on whether the critical exponent 𝜈 is real or complex. The boundary between these two regions is described by an exceptional point in the renormalization group (RG) flow, leading to distinctive features in correlation functions and the phase diagram. Another contrast with the previous work is the number and stability of the NEFPs as well as the underlying topology of the RG flow. Lastly, we investigate an extreme form of nonreciprocity where one order parameter is independent of the other order parameter but not vice versa. Unlike the ℤ 2 × ℤ 2 model, which becomes nonperturbative in this case, we identify a distinct nonequilibrium universality class whose dependent field similarly violates fluctuation-dissipation relations but does not exhibit discrete scale invariance or underdamped oscillations near criticality.

Critical phenomena↗
DOE OSTI2022

Kinetic frustration by limited bond availability controls the LAT protein condensation phase transition on membranes

LAT is a membrane-linked scaffold protein that undergoes a phase transition to form a two-dimensional protein condensate on the membrane during T cell activation. Governed by tyrosine phosphorylation, LAT recruits various proteins that ultimately enable condensation through a percolation network of discrete and selective protein-protein interactions. Here, we describe detailed kinetic measurements of the phase transition, along with coarse-grained model simulations, that reveal that LAT condensation is kinetically frustrated by the availability of bonds to form the network. Unlike typical miscibility transitions in which compact domains may coexist at equilibrium, the LAT condensates are dynamically arrested in extended states, kinetically trapped out of equilibrium. Modeling identifies the structural basis for this kinetic arrest as the formation of spindle arrangements, favored by limited multivalent binding interactions along the flexible, intrinsically disordered LAT protein. These results reveal how local factors controlling the kinetics of LAT condensation enable formation of different, stable condensates, which may ultimately coexist within the cell.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

zmodel-kokkos

The code implements the so-called ZModel on an unstructured mesh. This models two-phase (or feasibly multi-phase) incompressible flow by meshing and evolving only the interface between materials, and never directly discretizing the bulk flow.

May, Ian↗
NASA NTRS1990

Modelling and computationally efficient time domain linear equalisation of nonlinear bandlimited QPSK satellite channels

The problem of modeling and equalization of a nonlinear satellite channel is considered. The channel is assumed to be bandlimited and exhibits both amplitude and phase nonlinearities. In traditional models, computations are usually performed in the frequency domain and solutions are based on complex numerical techniques. A discrete time model is used to represent the satellite link with both uplink and downlink white Gaussian noise. Under conditions of practical interest, a simple and computationally efficient time-domain design technique for the minimum mean square error linear equalizer is presented. The efficiency of this technique is enhanced by the use of a fast and simple iterative algorithm for the computation of the autocorrelation coefficients of the output of the nonlinear channel. Numerical results on the evaluations of bit error probability and other relevant parameters needed in the design and analysis of a nonlinear bandlimited QPSK system demonstrate the simplicity and computational efficiency of the proposed approach.

Konstantinides, K.↗
NASA NTRSDate not supplied

Computational Study of Low-Speed Boundary-Layer Laminar-Turbulent Transition over a 30º Swept Plate in a Channel

Reliable boundary-layer laminar-turbulent transition predictions are critical for the design of aircraft that use natural or controlled laminar flow to achieve low drag and better fuel efficiency. During cruise flight through low turbulence, the three-dimensional boundary layer over swept surfaces can transition prematurely due to stationary crossflow instabilities induced by surface roughness. Parabolized stability equations analyses have been successfully used to explain the different phases of the transition process under well controlled perturbations such as steps or discrete roughness elements. However, methodologies are still lacking to model the boundary layer receptivity to distributed random roughness, and the following non-linear instability stages leading to transition. For engineering applications, several semi-empirical methods based on linear amplification factors or local helicity have found growing adoption, but their reliability is limited by their level of simplification and the sparsity of data they are based on. To address this gap, a new study of the effect of distributed random surface roughness on crossflow transition was started, using a 30◦ swept plate as a test vehicle. The plate, which is mounted in the 2 Foot × 3 Foot Low-Speed Boundary-Layer Channel at the NASA Langley Research Center, has been tested with micro-scale roughness of varied geometry near the leading edge. In order to analyze the flow in the test section and prepare Direct Numerical Simulations with resolved roughness, Reynolds-Averaged Navier-Stokes (RANS) simulations of the flow over the plate installed in the wind tunnel test section were performed using the Langtry-Menter transition model with stationary crossflow extension. This paper presents the results of these simulations, with comparison to the experiment.

ARMD↗
NASA NTRS1991

The zeroes of controlled structures - Sensor/actuator attributes and structural modelling

The effects of sensor/actuator design and modeling decisions on the zero frequencies of single input single output structural systems are explored. A zero definition based on the high gain closed loop character of the poles of a system is used to study parametrically zeroes as a function of sensor and actuator location. An infinite order model is adopted so that the effects of finite order modeling are removed. Examples focus on point or localized actuators and sensors and simple uniform structures. The sensitivity of the zero frequencies to changes in the sensor and actuator location for both collocated and non-collocated cases are presented. Non-minimum phase behavior is visible in only certain non-collocated situations. In addition to the sensor/actuator attributes, the zeroes are shown to be a function of the truncation and discretization used to model the structural response for collocated and dual systems. Methods for improving the zero predictions in the control bandwidth are suggested.

Fleming, Farla M.↗
NASA NTRS2016

Polarized Radiative Transfer of a Cirrus Cloud Consisting of Randomly Oriented Hexagonal Ice Crystals: The 3 x 3 Approximation for Non-Spherical Particles

The reflection and transmission of polarized light for a cirrus cloud consisting of randomly oriented hexagonal columns were calculated by two very different vector radiative transfer models. The forward peak of the phase function for the ensemble-averaged ice crystals has a value of order 6 x 10(exp 3) so a truncation procedure was used to help produce numerically efficient yet accurate results. One of these models, the Vectorized Line-by-Line Equivalent model (VLBLE), is based on the doubling- adding principle, while the other is based on a vector discrete ordinates method (VDISORT). A comparison shows that the two models provide very close although not entirely identical results, which can be explained by differences in treatment of single scattering and the representation of the scattering phase matrix. The relative differences in the reflected I and Q Stokes parameters are within 0.5 for I and within 1.5 for Q for all viewing angles. In 1971 Hansen showed that for scattering by spherical particles the 3 x 3 approximation is sufficient to produce accurate results for the reflected radiance I and the degree of polarization (DOP), and he conjectured that these results would hold also for non-spherical particles. Simulations were conducted to test Hansen's conjecture for the cirrus cloud particles considered in this study. It was found that the 3 x 3 approximation also gives accurate results for the transmitted light, and for Q and U in addition to I and DOP. For these non-spherical ice particles the 3 x 3 approximation leads to an absolute error 2 x 10(exp -6) for the reflected and transmitted I, Q and U Stokes parameters. Hence, it appears to be an excellent approximation, which significantly reduces the computational complexity and burden required for multiple scattering calculations.

Polarization↗
DOE OSTI2020

Discretization limits of lattice‐Boltzmann methods for studying immiscible two‐phase flow in porous media

Summary Digital images of porous media often include features approaching the image resolution length scale. The behavior of numerical methods at low resolution is therefore important even for well‐resolved systems. We study the behavior of the Shan‐Chen (SC) and Rothman‐Keller (RK) multicomponent lattice‐Boltzmann models in situations where the fluid‐fluid interfacial radius of curvature and/or the feature size of the medium approaches the discrete unit size of the computational grid. Various simple, small‐scale test geometries are considered, and a drainage test is also performed in a Bentheimer sandstone sample. We find that both RK and SC models show very high ultimate limits: in ideal conditions the models can simulate static fluid configuration with acceptable accuracy in tubes as small as three lattice units across for RK model (six lattice units for SC model) and with an interfacial radius of curvature of two lattice units for RK and SC models. However, the stability of the models is affected when operating in these extreme discrete limits: in certain circumstances the models exhibit behaviors ranging from loss of accuracy to numerical instability. We discuss the circumstances where these behaviors occur and the ramifications for larger‐scale fluid displacement simulations in porous media, along with strategies to mitigate the most severe effects. Overall we find that the RK model, with modern enhancements, exhibits fewer instabilities and is more suitable for systems of low fluid‐fluid miscibility. The shortcomings of the SC model seem to arise predominantly from the high, strongly pressure‐dependent miscibility of the two fluid components.

Li, Zhe↗
DOE OSTI2021

Validation of SAS4A/SASSYS-1 for predicting steady-state single-phase natural circulation

The validation of system analysis codes for nuclear reactor systems is required for the development and application of these computational tools. Designed as a comprehensive system analysis code for advanced nuclear reactors, SAS4A/SASSYS-1 requires validation of its physics model for capturing single-phase natural circulation behavior. To support the validation of SAS4A/SASSYS-1, high-precision experiments are performed capturing steady-state single-phase natural circulation on a scaled facility with comprehensive instrumentation. In this work, dedicated tests are performed quantifying the critical modeling parameters, and a single-phase natural circulation benchmark dataset is obtained with well-documented uncertainty and comprehensive facility description. The validation is then performed against the dataset examining the capability of SAS4A/SASSYS-1 in simulating steady-state single-phase natural circulation. The experimental facility is modeled in the candidate code. Solution verification is performed using Richardson-extrapolation-based estimators which quantify and restrict numerical errors from discretization. Input uncertainty provided by the benchmark dataset is forward propagated through the candidate code, quantifying the output uncertainty in a Monte Carlo approach. The composition of the output uncertainty is also quantified through a variance-based sensitivity analysis. With the uncertainty quantified for each individual condition, a detailed comparison between the simulation results and experimental data is performed covering the whole dataset. The results show consistent agreement for all primary parameters. The current validation activity provides a valuable benchmark dataset for the validation of system analysis codes in capturing single-phase natural circulation and demonstrates satisfactory prediction capability of SAS4A/SASSYS1 for steady-state single-phase natural circulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Validation of Numerical Models of the Offshore Wind Turbine From the Alpha Ventus Wind Farm Against Full-Scale Measurements Within OC5 Phase III

The main objective of the Offshore Code Comparison Collaboration Continuation, with Correlation (OC5) project is validation of aero-hydro-servo-elastic simulation tools for offshore wind turbines (OWTs) through comparison of simulated results to the response data of physical systems. Phase III of the OC5 project validates OWT models against the measurements recorded on a Senvion 5M wind turbine supported by the OWEC Quattropod from the alpha ventus offshore wind farm. The following operating conditions of the wind turbine were chosen for the validation: (1) idling below the cut-in wind speed, (2) rotor-nacelle assembly (RNA) rotation maneuver below the cut-in wind speed, (3) power production below and above the rated wind speed, and (4) shutdown. A number of validation load cases were defined based on these operating conditions. The following measurements were used for validation: (1) strains and accelerations recorded on the support structure and (2) pitch, yaw, and azimuth angles, generator speed, and electrical power recorded from the RNA. Strains were not directly available from the majority of the OWT simulation tools; therefore, strains were calculated based on out-of-plane bending moments, axial forces, and cross-sectional properties of the structural members. The simulation results and measurements were compared in terms of time series, discrete Fourier transforms, power spectral densities, and probability density functions of strains and accelerometers. A good match was achieved between the measurements and models setup by OC5 Phase III participants.

17 WIND ENERGY↗
DOE OSTI2025

Probabilistic Predictions for Fastener Failure in the Sandia Mechanics Challenge Using the Discrete-Direct Uncertainty Quantification Approach

This paper documents the blind and post-blind analysis predictions for the 2023 Sandia Mechanics Challenge (SMC), which involved predicting the behavior of a threaded fastener joint structure subjected to shock loading. Utilizing repeat sets of fastener calibration data from various experimental configurations including tension, double shear, and joint tension, we developed a library of calibrated models which were propagated through the application model using the Discrete-Direct (DD) uncertainty quantification (UQ) approach. Although the initial blind predictions did not incorporate spare-sample processing to quantify fastener failure probabilities, the analyses yielded reasonable conclusions aligned with experimental results. In the post-blind analysis phase, we focused on enhancing the fidelity of the aluminum constitutive model and innovating the DD approach to obtain probabilistic predictions for fastener failure, particularly when quantities of interest (QoIs) approach their bounds. The improved aluminum model captures the behavior of the cantilever under shock loading more accurately, predicting both partial and complete cracks, although it tends to underpredict failure propagation. The enhanced DD approach facilitates probabilistic predictions that reflect the interdependent failure mechanisms of the fasteners and the cantilever, revealing that while certain fasteners are more likely to fail, the failure does not necessarily follow a progressive pattern. Overall, the post-blind analyses significantly improved the predictive capabilities of the model, providing valuable insights into the SMC application and establishing a robust foundation for informed engineering decisions. The methodology demonstrates a cost-effective and extensible approach suitable for a wide range of applications, highlighting the importance of uncertainty quantification to provide context for engineering decision making.

42 ENGINEERING↗
DOE OSTI2026

Residual symmetries and scalar multiplet vacuum alignment in non-Abelian flavour models

We demonstrate that, upon minimizing a renormalizable, single-scalar potential invariant under a non-Abelian symmetry, special orientations in the associated vacuum alignment of the scalar multiplet correspond to the preservation of a discrete residual flavour symmetry in the broken phase of the theory. Conversely, we show that these special scalar alignments are perturbed when additional Lagrangian operators (e.g. renormalizable, multi-flavon operators and/or effective, higher-dimensional operators) are present that break said residual symmetry, leading to a vacuum reorientation and phenomenological consequences. We therefore construct a one-to-one correspondence principle between broken residual symmetries and vacuum alignment corrections, providing a mechanism to identify (and correct) a subtle but persistent form of phenomenologically relevant fine-tuning embedded in — but often ignored by — most successful non-Abelian flavour models. We first establish this correspondence in a set of toy models based on the S4 permutation symmetry, and then apply the lessons learned to the more realistic A4 Altarelli-Feruglio and ∆(27) Universal Texture Zero models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Quantification of Energy Savings and Demand Reduction for a Heat Pump Integrated with Thermal Energy Storage (Final Report)

To mitigate the variation in demand on the electric grid, thermal energy storage (TES) is an alternative to electric batteries or installing new peaking power plants. Stakeholders and policy makers across the United States have expressed interests in promoting TES, as demonstrated by the US Department of Energy’s Grid-Interactive Efficient Buildings program and the efforts of various state legislatures. However, the cost value provided by TES are unclear. If reliable cost benefits were determined, stakeholders would have a clearer picture of the financial returns that can be gained from their investment in TES. The study in this report is conducted by ORNL with collaboration with Emerson the Helix Innovation Center. In the first part of this report, EnergyPlus was used to perform whole-building simulations for two residential buildings in Indianapolis and Atlanta. The HVAC system in both buildings were equipped with phase change material TES. The TES tank was charged in off-peak hours and discharged in peak hours to perform load shifting. First, the economic value implied by existing time-of-use (TOU) rates offered by utility companies was analyzed via whole-building simulation. Second, existing demand reduction (DR) incentives sourced from 3 different electrical grid administrators (i.e., California, Texas, and New England region) were surveyed to determine their implied value. The study suggests that the traditional value analysis that focuses on ROI for the building owner significantly undervalues TES technology making economic viability difficult. A more comprehensive value analysis that includes peak demand management and deferred capital for peaking power plants shows that TES should be economically viable but here the value is greater for the utility and requires large market penetration and aggregation to fully realize the benefits. Therefore, to facilitate commercialization, new business models are needed that include a broader range of stakeholders and distribute the value of TES proportionally. In the 2 nd part of this project, the benefits of a novel phase change material (PCM) integrated heat pump configuration were evaluated via detailed component based simulation. A one-dimensional PCM heat exchanger model which discretizes the PCM tank and refrigerant tubes into small control volumes is developed. Each control volume can have different PCM temperatures, PCM properties, and heat transfer coefficients. The PCM tank is charged by a wrapped tank condenser and discharged by an internal refrigerant coil. The PCM heat exchanger model is integrated into DOE/ORNL Heat Pump Design Model for heat pump system simulation. To demonstrate the performance of the PCM integrated heat pump, a case study in Chicago was performed. A Time-of-Use utility structure-based control strategy is implemented to schedule the PCM tank charging and discharging mode switching. Compared with a conventional electric heat pump, the PCM integrated heat pump shows superior performance on load shifting and utility cost reduction. As a result, the proposed system demonstrates 24.6% utility saving for cooling application and 25.8% utility saving for heating application.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS2011

Kennedy Space Center Orion Processing Team Planning for Ground Operations

Topics in this presentation are: Constellation Ares I/Orion/Ground Ops Elements Orion Ground Operations Flow Orion Operations Planning Process and Toolset Overview, including: 1 Orion Concept of Operations by Phase 2 Ops Analysis Capabilities Overview 3 Operations Planning Evolution 4 Functional Flow Block Diagrams 5 Operations Timeline Development 6 Discrete Event Simulation (DES) Modeling 7 Ground Operations Planning Document Database (GOPDb) Using Operations Planning Tools for Operability Improvements includes: 1 Kaizen/Lean Events 2 Mockups 3 Human Factors Analysis

Letchworth, Gary↗
DOE OSTI2024

Frank–Kasper Phases of Diblock Copolymer Melts: Self-Consistent Field Results of Two Commonly Used Models

We constructed phase diagrams of conformationally asymmetric diblock copolymer A-B melts using the polymer self-consistent field (SCF) calculations of both the dissipative particle dynamics chain (DPDC) model (i.e., compressible melts of discrete Gaussian chains with the DPD non-bonded potential) and the “standard” model (i.e., incompressible melts of continuous Gaussian chains with the Dirac δ-function non-bonded potential) in the χN-ε plane, where χN and ε characterize, respectively, the repulsion and conformational asymmetry between the A and B blocks, at the A-block volume fraction f = 0.2 and 0.3. Consistent with previous SCF calculations of the “standard” model, σ and A15 are the only stable Frank–Kasper (FK) phases among the five FK (i.e., σ, A15, C14, C15 and Z) phases considered. The stability of σ and A15 is due to their delicate balance between the energetic and entropic contributions to the Helmholtz free energy per chain of the system, which, within our parameter range, increases in the order of σ/A15, Z, and C14/C15. While in general the SCF phase diagrams of these two models are qualitatively consistent, A15 is not stable for the DPDC model at the copolymer chain length N = 10 and f = 0.3; any differences in the SCF phase diagrams are solely due to the differences between these two models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Validation of SAS4A/SASSYS-1 for Steady-State Single-Phase Natural Circulation

The development of numerical software for engineering applications requires the validation of the code against experimental benchmark datasets. To support the validation of SAS4A/SASSYS-1 in simulating the single-phase natural circulation within nuclear reactor systems, high-precision steady-state experiments are performed to capture single-phase natural circulation phenomena on an existing scaled facility with comprehensive instrumentation. Forced convection tests are designed precisely capturing the facility’s critical thermal-hydraulic parameters required for one-dimensional modeling, and a comprehensive single-phase natural circulation dataset is obtained with well-documented experimental uncertainty and facility description. Analysis and discussion based on the finalized dataset confirm the dataset’s ability in capturing dominant phenomena and address important physical interpretation of parameters. The dataset reported in this project provides a valuable benchmark resource for the validation of system analysis codes under single-phase natural circulation inside nuclear reactors. Validation of SAS4A/SASSYS-1 is then performed against the obtained benchmark dataset to confirm the capability of its physics model in simulating the steady-state single-phase natural circulation. The experimental facility is modeled in the candidate code. Solution verification is performed using Richardson-extrapolation-based estimators to quantify and restrict numerical errors from discretization. Numerical uncertainty originating from finite maximum pseudo-transient time is also quantified and restricted. Input uncertainty provided by the benchmark dataset is forward propagated through the candidate code, directly quantifying the simulation output in a Monte Carlo approach. The composition of the output uncertainty is further determined by the estimators of Sobol’ indices through a variance-based sensitivity analysis. With uncertainty quantified for each individual condition, detailed comparison between the simulation results and experimental data is performed covering the whole dataset, which shows consistent agreement for all important quantities. The validation activity reported in this project demonstrates that SAS4A/SASSYS-1 can predict the primary parameters with satisfactorily accuracy under steady-state single-phase natural circulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗