Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion 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 217 records · Page 12

Multi-solid and multi-fluid diffuse interface model: Applications to dynamic fracture and fragmentation

We extend the model of diffuse solid–fluid interfaces developed earlier by authors of this paper to the case of arbitrary number of interacting hyperelastic solids. Plastic transformations of solids are taken into account through a Maxwell type model. The specific energy of each solid is given in separable form: it is the sum of a hydrodynamic part of the energy depending only on the density and the entropy, and an elastic part of the energy which is unaffected by the volume change. It allows us to naturally pass to the fluid description in the limit of vanishing shear modulus. In spite of a large number of governing equations, the model has a quite simple mathematical structure: it is a duplication of a single visco-elastic model. The model is well posed both mathematically and thermodynamically: it is hyperbolic and compatible with the second law of thermodynamics. The resulting model can be applied in the situations involving an arbitrary number of fluids and solids. In particular, we show the ability of the model to describe spallation and penetration phenomena occurring during high velocity impacts.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A comparison of probabilistic generative frameworks for molecular simulations

Generative artificial intelligence is now a widely used tool in molecular science. Despite the popularity of probabilistic generative models, numerical experiments benchmarking their performance on molecular data are lacking. Here, in this work, we introduce and explain several classes of generative models, broadly sorted into two categories: flow-based models and diffusion models. We select three representative models: neural spline flows, conditional flow matching, and denoising diffusion probabilistic models, and examine their accuracy, computational cost, and generation speed across datasets with tunable dimensionality, complexity, and modal asymmetry. Our findings are varied, with no one framework being the best for all purposes. In a nutshell, (i) neural spline flows do best at capturing mode asymmetry present in low-dimensional data, (ii) conditional flow matching outperforms other models for high-dimensional data with low complexity, and (iii) denoising diffusion probabilistic models appear the best for low-dimensional data with high complexity. Our datasets include a Gaussian mixture model and the dihedral torsion angle distribution of the Aib9 peptide, generated via a molecular dynamics simulation. We hope our taxonomy of probabilistic generative frameworks and numerical results may guide model selection for a wide range of molecular tasks.

Artificial intelligence↗

Forecasting residential solar photovoltaic deployment in California

Residential distributed photovoltaic (PV) deployment in the United States has experienced robust growth, and policy changes impacting the value of solar are likely to occur at the federal and state levels. To establish a credible baseline and evaluate impacts of potential new policies, this analysis employs multiple methods to forecast residential PV deployment in California, including a time-series forecasting model, a threshold heterogeneity diffusion model, a Bass diffusion model, and National Renewable Energy Laboratory's dSolar model. As a baseline, the residential PV market in California is modeled to peak in the early 2020s, with a peak annual installation of 1.5-2 GW across models. We then use the baseline results from the dSolar model and the threshold model to gauge the impact of the recent federal investment tax credit (ITC) extension, the newly approved California net energy metering (NEM) policy, and a hypothetical value-of-solar (VOS) compensation scheme. We find that the recent ITC extension may increase annual PV installations by 12%-18% (roughly 500 MW, MW) for the California residential sector in 2019-2020. The new NEM policy only has a negligible effect in California due to the relatively small new charges (< 100 MW in 2019-2020). Moreover, impacts of the VOS compensation scheme (0.12 cents per kilowatt-hour) are larger, reducing annual PV adoption by 32% (or 900-1300 MW) in 2019-2020.

14 SOLAR ENERGY↗

The Effect of Diffusional Transport and Surface Catalysis on the Aerothermodynamic Heating for Martian Atmospheric Entry

For the reentry heating of 70-degree blunt sphere-cones entering Mars at velocities in excess of 6 km/sec, a study was performed to determine the magnitude of both convective and catalytic heating on the cone forebody. Case studies of the peak heating conditions for the current NASA Mars Pathfinder entry probe, as well as the peak heating conditions of the proposed NASA Mars 2001 aerobraking orbiter mission were performed. Since the actual behavior of the chemical recombination of Mars atmosphere hypersonic shock layer species on heat shield materials is not known, and no experimental studies of any consequence have been done to determine surface reaction rates for the CO2 system, a parametric analysis of various reaction schemes and surface catalytic mechanistic models was performed. In many cases the actual Mars entry probe will have a heat shield composed of an ablative material which can either partially or completely preclude the existence of catalytic surface recombination phenomena. The extent of this blockage effect has not been examined at this time and is not considered in this effort. Instead only non-ablative computations are performed to separately evaluate the full extent of the surface catalysis and related diffusional phenomena. Parametric studies include peak heating point comparisons of non-catalytic, partially catalytic and fully catalytic total surface heat transfer for three difference surface recombination catalytic mechanisms (with and without simultaneous CO + O and O + O recombination) as well as with different species diffusion models. Diffusion model studies include constant Schmidt number (equal diffusivities) and the bifurcation, approximate multi-component diffusion model. Since the gas phase reaction kinetics for the Mars CO2 based atmosphere have also not been validated, the effect of two different (fast and slow) sets of gas kinetics on heat transfer is presented.

Henline, William D.↗

The small ice cap instability in diffusive climate models

Simple climate models which invoke diffusive heat transport and ice cap albedo feedback have equilibrium solutions with no stable ice cap smaller than a radius of about 20 deg on a great circle. Attention is presently given to a solution of this phenomenon which is physically appealing. The ice-free solution has a thermal minimum, and if the minimum temperature is just above the critical value for ice formation, then the artificial addition of a patch of ice leads to a widespread depression of the temperature below the critical freezing temperature. A second stable solution will then exist whose spatial extent is determined by the range of the influence function of a point sink of heat, due to the albedo shift in the patch.

North, G. R.↗

On modeling pressure diffusion in non-homogeneous shear flows

New models are proposed for the 'slow and 'rapid' parts of the pressure diffusive transport based on the examination of DNS databases for plane mixing layers and wakes. The model for the 'slow' part is non-local, but requires the distribution of the triple-velocity correlation as a local source. The latter can be computed accurately for the normal component from standard gradient diffusion models, but such models are inadequate for the cross component. More work is required to remedy this situation.

Demuren, A. O.↗

Experimental Validation of FE/BEM Dynamic Strain Model Under Diffuse Acoustic Field Loading

Structural finite element (FE) models naturally output displacement or acceleration response data. However, they can also be used to compute stress, internal forces, and strain response. When coupled with a boundary element model (BEM) of the fluid surrounding the structure, a fully coupled analysis can be performed. Modeling a diffuse acoustic field in the BEM fluid provides an excitation like that found when the structure is placed in a reverberation chamber. Fully coupling the structural FE model to the acoustic BEM model provides a means to predict not only the acceleration response of the panel to diffuse field loading, but also the ability to predict the dynamic stress and strain response. This type of model has been available with current predictive tools, but experimental validation of the prediction of dynamic stress or strain is difficult to find. An aluminum panel was instrumented with accelerometers and strain gages and hung in a reverberation room and subjected to a diffuse acoustic field. This paper presents the comparison of the experimental and predicted results.

Tsoi, W. Ben↗

Physically constrained 3D diffusion for inverse design of fiber-reinforced polymer composite materials

Designing fiber-reinforced polymer composites (FRPCs) with a tailored nonlinear stress-strain response is crucial for applications such as energy absorption in crash structures, flexible robotics, and impact-resistant protective gear. However, the inherent complexities of composite materials and the multitude of parameters involved, render traditional design and optimization methods inadequate for achieving effective inverse design of composites. In this paper, we present an AI-based inverse design framework that effectively and efficiently generates FRPCs with targeted nonlinear stress-strain responses. We introduce a physically constrained diffusion model (PC3D_Diffusion) capable of managing the complexities of composite materials and producing detailed, high-quality designs. We propose a loss-guided, learning-free approach to generate physically feasible microstructure designs by explicitly enforcing physical constraints during the generation process. For training purposes, 1.35 million FRPC samples were created, and their corresponding stress-strain curves were computed using established physics-based computational models. The results show that PC3D_Diffusion consistently generates high-quality designs with tailored mechanical behaviors, while guaranteeing compliance with the physical constraints. PC3D_Diffusion advances FRPC inverse design and may facilitate the inverse design of other 3D materials, offering potential applications in industries reliant on materials with custom mechanical properties.

Xu, Pei [Clemson Univ., SC (United States)]↗

A comparison of molecular dynamics and diffuse interface model predictions of Lennard-Jones fluid evaporation

The unsteady evaporation of a thin planar liquid film is studied by molecular dynamics simulations of Lennard-Jones fluid. The obtained results are compared with the predictions of a diffuse interface model in which capillary Korteweg contributions are added to hydrodynamic equations, in order to obtain a unified description of the liquid bulk, liquid-vapor interface and vapor region. Particular care has been taken in constructing a diffuse interface model matching the thermodynamic and transport properties of the Lennard-Jones fluid. The comparison of diffuse interface model and molecular dynamics results shows that, although good agreement is obtained in equilibrium conditions, remarkable deviations of diffuse interface model predictions from the reference molecular dynamics results are observed in the simulation of liquid film evaporation. It is also observed that molecular dynamics results are in good agreement with preliminary results obtained from a composite model which describes the liquid film by a standard hydrodynamic model and the vapor by the Boltzmann equation. The two mathematical model models are connected by kinetic boundary conditions assuming unit evaporation coefficient.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Importance Sampling Model-Based Diffusion for Trajectory Optimization

Trajectory optimization for robotic systems remains a challenging problem. This is especially true for robotic systems featuring nonlinear dynamics and many degrees of freedom. Data-based or model-free diffusion has recently been popularized in the fields of artificial intelligence and trajectory optimization. Model-Based Diffusion provides a data-free method of trajectory optimization, trained at runtime on a system dynamics model, suitable for high-dimensional models. This paper examines how importance sampling can enhance the performance of Model-Based Diffusion for trajectory optimization. Here, we quantify the benefits of importance sampling across three long horizon planning tasks. These results show as much as a 13x improvement in sample efficiency depending on environment and optimization parameters.

Golembeski, Seth [Georgia Institute of Technology,↗

Contributions to Loss Across the Magnetopause During an Electron Dropout Event

Abstract Dropout events are dramatic decreases in radiation belt electron populations that can occur in as little as 30 minutes. Loss to magnetopause due to a combination of magnetopause shadowing and outward radial transport plays a significant role in these events. We examine the dropout of relativistic electron populations during the October 2012 geomagnetic storm using simulated electron phase space density, evaluating the contribution of different processes to losses across the magnetopause. We compare loss contribution from outward transport calculated using a standard empirical radial diffusion model that assumes a dipolar geomagnetic field to an event‐specific radial diffusion model evaluated with a non‐dipolar geomagnetic field. We additionally evaluate the contribution of Shabansky type 1 particles, which bounce along magnetic field lines with local equatorial maxima, to the loss calculated during this event. We find that the empirical radial diffusion model with a dipolar background field underestimates the contribution of radial diffusion to this dropout event by up to 10% when compared to the event‐specific, non‐dipolar radial diffusion model. We additionally find that including Shabansky type 1 particles in the initial electron phase space density, that is, allowing some magnetic field lines distorted from the typical single‐minima configuration in drift shell construction, increases the calculated loss by an average of 0.75%. This shows that the treatment of the geomagnetic field significantly impacts the calculation of electron losses to the magnetopause during dropout events, with the non‐dipolar treatment of radial diffusion being essential to accurately quantify the loss of outer radiation belt populations.

79 ASTRONOMY AND ASTROPHYSICS↗

Integration and validation of some modules for modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures

Three modules are integrated into the built-in OpenFOAM rhoCentralFoam solver towards accurate and efficient modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures within the OpenFOAM 10.0 framework. The first module is the mixture-averaged diffusion model. The second module is the built-in OpenFOAM Lagrangian solver coupled with optimised droplet drag coefficient and convective heat transfer coefficient sub-models. The last module is a sparse stiff chemistry solver based on dynamic adaptive hybrid integration (AHI-S). The optimised droplet sub-models are first verified in correct implementation for subsequent simulations in this work. Further, they show good accuracy against experimental and analytical data in the modelling of ammonia droplet acceleration and cooling in the flowing and/or low-temperature air. The accuracy and efficiency gains related to the mixture-averaged diffusion model and the AHI-S chemistry solver are examined by simulating 1-D detonation propagation in ammonia droplet-free/laden ammoniaoxygen mixtures. Numerical results of detonation propagation speed, gaseous temperature, density, and species distributions around the induction zone show good agreement with experimental data and analytical solutions. Compared to the built-in OpenFOAM diffusion model, the mixture-averaged diffusion model provides different numerical predictions of pulsating instabilities in detonation propagation. It shows better accuracy in depicting the detonation structure within the droplet-free section attributed to improved multi-component diffusion modelling. Compared to the built-in OpenFOAM solver EulerImplicit (backward Euler), the AHI-S chemistry solver reduces the computational cost by around 50%. It achieves satisfactory accuracy in calculating detonation propagation speed within the droplet-free section with the optimal efficiency when the safety factor, β, equals 0.5.

42 ENGINEERING↗

A diffuse radar scattering model from Martian surface rocks

Remote sensing of Mars has been done with a variety of instrumentation at various wavelengths. Many of these data sets can be reconciled with a surface model of bonded fines (or duricrust) which varies widely across the surface and a surface rock distribution which varies less so. A surface rock distribution map from -60 to +60 deg latitude has been generated by Christensen. Our objective is to model the diffuse component of radar reflection based on this surface distribution of rocks. The diffuse, rather than specular, scattering is modeled because the diffuse component arises due to scattering from rocks with sizes on the order of the wavelength of the radar beam. Scattering for radio waves of 12.5 cm is then indicative of the meter scale and smaller structure of the surface. The specular term is indicative of large scale surface undulations and should not be causally related to other surface physical properties. A simplified model of diffuse scattering is described along with two rock distribution models. The results of applying the models to a planet of uniform fractional rock coverage with values ranging from 5 to 20% are discussed.

Calvin, W. M.↗

Multi-scale diffuse interface modeling of multi-component two-phase flow with partial miscibility

In this paper, we introduce a diffuse interface model to simulate multi-component two-phase flow with partial miscibility based on a realistic equation of state (e.g. Peng–Robinson equation of state). Because of partial miscibility, thermodynamic relations are used to model not only interfacial properties but also bulk properties, including density, composition, pressure, and realistic viscosity. As far as we know, this effort is the first time to use diffuse interface modeling based on equation of state for modeling of multi-component two-phase flow with partial miscibility. In numerical simulation, the key issue is to resolve the high contrast of scales from the microscopic interface composition to macroscale bulk fluid motion since the interface has a nanoscale thickness only. To efficiently solve this challenging problem, we develop a multi-scale simulation method. At the microscopic scale, we deduce a reduced interfacial equation under reasonable assumptions, and then we propose a formulation of capillary pressure, which is consistent with macroscale flow equations. Moreover, we show that Young–Laplace equation is an approximation of this capillarity formulation, and this formulation is also consistent with the concept of Tolman length, which is a correction of Young–Laplace equation. At the macroscopical scale, the interfaces are treated as discontinuous surfaces separating two phases of fluids. Our approach differs from conventional sharp-interface two-phase flow model in that we use the capillary pressure directly instead of a combination of surface tension and Young–Laplace equation because capillarity can be calculated from our proposed capillarity formulation. A compatible condition is also derived for the pressure in flow equations. Furthermore, based on the proposed capillarity formulation, we design an efficient numerical method for directly computing the capillary pressure between two fluids composed of multiple components. Finally, numerical tests are carried out to verify the effectiveness of the proposed multi-scale method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Satellite ultraviolet measurements of nitric oxide fluorescence with a diffusive transport model.

Twilight measurements of fluorescence in the (1, 0) gamma band of nitric oxide were made from June 1967 to January 1969 by an ultraviolet scanning spectrometer on board the polar orbiting satellite Ogo 4. Nitric oxide vertical column emission rates were measured between solar zenith angles of 93 and 98 deg. Seasonal and latitudinal variations were found to be less than a factor of 1.3, the scatter and uncertainty in the data prohibiting more precise determinations from being made. Time independent chemical diffusion models for the vertical distribution of nitric oxide agree well with profiles measured from sounding rockets. The column emission rates calculated from the theoretical models are larger than the satellite measurements by a factor of 3.

Rusch, D. W.↗

The Implementation and Evaluation of the Emergency Response Dose Assessment System (ERDAS) at Cape Canaveral Air Station/Kennedy Space Center

The Emergency Response Dose Assessment System (ERDAS) is a system which combines the mesoscale meteorological prediction model RAMS with the diffusion models REEDM and HYPACT. Operators use a graphical user interface to run the models for emergency response and toxic hazard planning at CCAS/KCS. The Applied Meteorology Unit has been evaluating the ERDAS meteorological and diffusion models and obtained the following results: (1) RAMS adequately predicts the occurrence of the daily sea breeze during non-cloudy conditions for several cases. (2) RAMS shows a tendency to predict the sea breeze to occur slightly earlier and to move it further inland than observed. The sea breeze predictions could most likely be improved by better parameterizing the soil moisture and/or sea surface temperatures. (3) The HYPACT/REEDM/RAMS models accurately predict launch plume locations when RAMS winds are accurate and when the correct plume layer is modeled. (4) HYPACT does not adequately handle plume buoyancy for heated plumes since all plumes are presently treated as passive tracers. Enhancements should be incorporated into the ERDAS as it moves toward being a fully operational system and as computer workstations continue to increase in power and decrease in cost. These enhancements include the following: activate RAMS moisture physics; use finer RAMS grid resolution; add RAMS input parameters (e.g. soil moisture, radar, and/or satellite data); automate data quality control; implement four-dimensional data assimilation; modify HYPACT plume rise and deposition physics; and add cumulative dosage calculations in HYPACT.

Evans, Randolph J.↗