Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lagrangian Particles”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A massively parallel and scalable multi-CPU material point method

Harnessing the power of modern multi-GPU architectures, we present a massively parallel simulation system based on the Material Point Method (MPM) for simulating physical behaviors of materials undergoing complex topological changes, self-collision, and large deformations. Our system makes three critical contributions. First, we introduce a new particle data structure that promotes coalesced memory access patterns on the GPU and eliminates the need for complex atomic operations on the memory hierarchy when writing particle data to the grid. Second, we propose a kernel fusion approach using a new Grid-to-Particles-to-Grid (G2P2G) scheme, which efficiently reduces GPU kernel launches, improves latency, and significantly reduces the amount of global memory needed to store particle data. Finally, we introduce optimized algorithmic designs that allow for efficient sparse grids in a shared memory context, enabling us to best utilize modern multi-GPU computational platforms for hybrid Lagrangian-Eulerian computational patterns. We demonstrate the effectiveness of our method with extensive benchmarks, evaluations, and dynamic simulations with elastoplasticity, granular media, and fluid dynamics. In comparisons against an open-source and heavily optimized CPU-based MPM codebase [Fang et al. 2019] on an elastic sphere colliding scene with particle counts ranging from 5 to 40 million, our GPU MPM achieves over 100x per-time-step speedup on a workstation with an Intel 8086K CPU and a single Quadro P6000 GPU, exposing exciting possibilities for future MPM simulations in computer graphics and computational science. Moreover, compared to the state-of-the-art GPU MPM method [Hu et al. 2019a], we not only achieve 2x acceleration on a single GPU but our kernel fusion strategy and Array-of-Structs-of-Array (AoSoA) data structure design also generalizes to multi-GPU systems. Our multi-GPU MPM exhibits near-perfect weak and strong scaling with 4 GPUs, enabling performant and large-scale simulations on a 10243 grid with close to 100 million particles with less than 4 minutes per frame on a single 4-GPU workstation and 134 million particles with less than 1 minute per frame on an 8-GPU workstation.

Wang, Xinlei↗
DOE OSTI2020

Solving gyrokinetic systems with higher-order time dependence

We discuss theoretical and numerical aspects of gyrokinetics as a Lagrangian field theory when the field perturbation is introduced into the symplectic part. A consequence is that the field equations and particle equations of motion in general depend on the time derivatives of the field. The most well-known example is when the parallel vector potential is introduced as a perturbation, where a time derivative of the field arises only in the equations of motion, so an explicit equation for the fields may still be written. We will consider the conceptually more problematic case where the time-dependent fields appear in both the field equations and equations of motion, but where the additional term in the field equations is formally small. The conceptual issues were described by Burby (J. Plasma Phys., vol. 82 (3), 2016, 905820304): these terms lead to apparent additional degrees of freedom to the problem, so that the electric field now requires an initial condition, which is not required in low-frequency (Darwin) Vlasov–Maxwell equations. Also, the small terms in the Euler–Lagrange equations are a singular perturbation, and these two issues are interlinked. For well-behaved problems the apparent additional degrees of freedom are spurious, and the physically relevant solution may be directly identified. Because we needed to assume that the system is well behaved for small perturbations when deriving gyrokinetic theory, we must continue to assume that when solving it, and the physical solutions are thus the regular ones. The spurious nature of the singular degrees of freedom may also be seen by changing coordinate systems so the varying field appears only in the Hamiltonian. We then describe how methods appropriate for singular perturbation theory may be used to solve these asymptotic equations numerically. We then describe a proof-of-principle implementation of these methods for an electrostatic strong-flow gyrokinetic system; two basic test cases are presented to illustrate code functionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

Gas dynamics for accretion disk simulations

The behavior of accretion disks can largely be understood in terms of the basic physical processes of mass, energy, and momentum conservation. Despite this, detailed modeling of these systems using modern computational techniques is challenging and controversial. Disturbing differences exist between methods used widely in astrophysics, namely Eulerian finite-difference techniques and particle codes such as SPH. Therefore neither technique is fully satisfactory for accretion disk simulations. This paper describes a new fully Lagrangian method designed to resolve these difficulties.

Whitehurst, R.↗
NASA NTRS2020

Thermal Impact of Cosmic Ray Interaction with an XRay Microcalorimeter Array

The X-ray Integral Field Unit (X-IFU) instrument on the Athena mission will be positioned at the Lagrangian point L2 and be subject to cosmic rays generated by astrophysics sources, primarily relativistic protons. Previous simulations have shown that particles of energy higher than 150 MeV will make it through the outer layers of the satellite. They will reach the detector wafer with a rate of 3 cts/sq. cm/s and a most probable energy deposited in the Si frame supporting the array at 150 keV. These events can affect the energy resolution of the detectors through the thermal fluctuations that they produce. This study assesses this potential problem and discusses two suggested design approaches to decrease the impact of cosmic ray in order to limit their effect to their allocation of 0.2 eV within the Athena/X-IFU energy-resolution budget. The first is the addition of a coating layer of high heat capacity material (e.g., Pd) and the second is the splitting of this coating into two thermal regions near the TES array to keep the heat away from the array. Implementing these two features is predicted to cause a decrease in the number of events above 1 μK by more than a factor 10 to ~ 1.5 cps when compared to an equivalent design without these features.

Antoine R. Miniussi↗
DOE OSTI2025

pnnl/Lagrangian-Droplets

The Lagrangian Droplet cloud microphysics scheme uses high-resolution data of air velocity, temperature, and water vapor mixing ratio to drive particle-based simulations of cloud microphysics. The model was first implemented using a high-resolution Large Eddy Simulation of the Pi Chamber at Michigan Tech, which is small-scale convection chamber the generates clouds in a laboratory setting

Fierce, Laura [Pacific Northwest National Laborato↗
NASA NTRS1993

Large-scale isentropic mixing properties of the Antarctic polar vortex from analyzed winds

Winds derived from analyzed geopotential height fields are used to study quasi-horizontal mixing by the large-scale flow in the lower stratosphere during austral spring. This is the period when the Antarctic ozone hole appears and disappears. Trajectories are computed for large ensembles of particles initially inside and outside the main polar vortex. Mixing and transport are diagnosed through estimates of finite time Lyapunov exponents and Lagrangian dispersion statistics of the tracer trajectories. At 450 K and above prior to the vortex breakdown: Lyapunov exponents are a factor of 2 smaller inside the vortex than outside; diffusion coefficients are an order of magnitude smaller inside than outside the vortex; and the trajectories reveal little exchange of air across the vortex boundary. At lower levels (425 and 400 K) mixing is greater, and there is substantial exchange of air across the vortex boundary. In some years there are large wave events that expel small amounts of vortex air into the mid-latitudes. At the end of the spring season during the vortex breakdown there is rapid mixing of air across the vortex boundary, which is evident in the mixing diagnostics and the tracer trajectories.

Bowman, Kenneth P.↗
NASA NTRS2017

Modeling Solar Energetic Particle Events Using ENLIL Heliosphere Simulations

Solar energetic particle (SEP) event modeling (SEPMOD) has gained renewed attention in part because of the availability of a decade of multipoint measurements from STEREO (Solar TErrestrial RElations Observatory) and L1 (Lagrangian point 1) spacecraft at 1 AU (Astronomical Unit). These observations are coupled with improving simulations of the geometry and strength of heliospheric shocks obtained by using coronagraph images to send erupted material into realistic solar wind backgrounds. The STEREO and ACE (Aerosol, Cloud systems, ocean Ecosystems) measurements in particular have highlighted the sometimes surprisingly widespread nature of SEP events. It is thus an opportune time for testing SEP models, which typically focus on protons approximately 1-100 megaelectronvolts, toward both physical insight to these observations and potentially useful space radiation environment forecasting tools. Some approaches emphasize the concept of particle acceleration and propagation from close to the Sun, while others emphasize the local field line connection to a traveling, evolving shock source. Among the latter is the previously introduced SEPMOD treatment, based on the widely accessible and well-exercised WSA-ENLIL (Wang-Sheeley-Arge-ENLIL)-cone model. SEPMOD produces SEP proton time profiles at any location within the ENLIL domain. Here we demonstrate a SEPMOD version that accommodates multiple, concurrent shock sources occurring over periods of several weeks. The results illustrate the importance of considering longer-duration time periods and multiple CME (Coronal Mass Ejection) contributions in analyzing, modeling, and forecasting SEP events.

Luhmann, J. G.↗
DOE OSTI2024

Noncommutative gauge symmetry in the fractional quantum Hall effect

Abstract We show that a system of particles on the lowest Landau level can be coupled to a probe U(1) gauge field$$ \mathcal{A} $$ A μ in such a way that the theory is invariant under a noncommutative U(1) gauge symmetry. While the temporal component$$ \mathcal{A} $$ A 0 of the probe field is coupled to the projected density operator, the spatial components$$ \mathcal{A} $$ A i are best interpreted as quantum displacements, which distort the interaction potential between the particles. We develop a Seiberg-Witten-type map from the noncommutative U(1) gauge symmetry to a simpler version, which we call “baby noncommutative” gauge symmetry, where the Moyal brackets are replaced by the Poisson brackets. The latter symmetry group is isomorphic to the group of volume preserving diffeomorphisms. By using this map, we resolve the apparent contradiction between the noncommutative gauge symmetry, on the one hand, and the particle-hole symmetry of the half-filled Landau level and the presence of the mixed Chern-Simons terms in the effective Lagrangian of the fractional quantum Hall states, on the other hand. We outline the general procedure which can be used to write down effective field theories which respect the noncommutative U(1) symmetry.

Physics↗
NASA NTRS1990

Lidar and aircraft studies of deep Cirrus systems from the 1986 FIRE IFO

Several NCAR King Air flight missions were conducted during the Wisconsin FIRE IFO experiment in support of the University of Utah polarization lidar observations of deep cirrus cloud systems at the Wausau ground site. Data collected from four cirrus systems are included in this analysis, including those of 22 and 28 October, and 1 and 2 November. Lidar data were generally obtained at 2 min intervals in the zenith direction over observation periods that ranged from approximately 4 to 10 h, bracketing the aircraft missions. The data were processed to yield height-time (HTI) displays of lidar linear depolarization ratio sigma and relative range-normalized return power P. King Air operations consisted of a combination of rapid profiling and Lagrangian spiral descents and stacked racetrack patterns in the vicinity of the field site. From the spiral descents are constructed vertical profiles of ice particle concentration N(sub i) and ice mass content IWC derived from PMS 2-D probe imagery and, when detected, FSSP cloud droplet concentration N(sub W) and liquid water content, LWC. Aircraft flight leg data are presented for the vertical velocity W and the same ice and water cloud content parameters. In addition, aerosol particle concentrations obtained with the ASAS probe are examined, and photographs of ice particles collected in-situ on oil-coated slides are presented to illustrate ice particle habit.

Sassen, Kenneth↗
NASA NTRS2018

Towards an Enhanced Droplet Activation Scheme for Multi-Moment Bulk Microphysics Schemes

Initial droplet spectra produced upon activation impact the ensuing chain of microphysical processes andtherefore play a crucial role in cloud evolution. This work re-examines dependencies of newly formed clouddroplet size distribution (CDSD) characteristics on environmental and aerosol properties via parcel model simulationsthat serve as the basis for a multi-moment bulk microphysics droplet activation scheme suitable for acloud-resolving model (CRM). It is found that applying a fixed size threshold to define activated droplets versusemploying physical considerations can lead to erroneous activation and overly broad CDSDs for high aerosolconcentration and weak updraft conditions. Aerosol distributions characterized by larger median sizes and/orincreased solubility can result in greater activated droplet numbers, whereas impacts of these parameters onCDSD spectral width depend on both aerosol number concentration and updraft velocity. An expansion of theactivation scheme to include CDSD spectral width is proposed to aid efforts to extend high-order momentprediction to cloud droplet categories in CRMs as well as better represent variability in the activation process onthe cloud scale.simulations to investigate the regime dependence of the relative dispersion(d)1 of newly activated CDSDs, where d is the ratio of dropletradius standard deviation (σ) to the mean radius (r ). C16 demonstratedthat increasing Na resulted in increasing (decreasing) d values via reducedcondensational narrowing (spectral broadening) rates within theAL (UL) regime, with d values peaking in the TR regime. Their findingssuggest a similar regime dependence for d as R09 noted for Nc and helpexplain reportedly conflicting relationships between Na and CDSDspectral characteristics (cf. Hudson and Noble, 2014; Liu et al., 2014),although the applicability of these results within bulk microphysicalschemes was not addressed.Simulating aerosol-cloud interactions with CRMs employing bulkmicrophysics requires that the scheme minimally predict two CDSDparameters, namely mass and number concentrations, and represent thedroplet activation process. Various activation schemes aim to determineNc from aerosol and environmental properties and include analyticalexpressions (e.g., Abdul-Razzak et al., 1998; Morrison et al., 2005) aswell as lookup tables (LUTs) based on detailed parcel model calculations(e.g., Saleeby and Cotton, 2004, hereafter SC04; Segal and Khain,2006; Thompson and Eidhammer, 2014). Expressions to diagnose CDSDspectral width from Nc (Grabowski, 1998; Liu et al., 2006; Morrison andGrabowski, 2007) or cloud water content (Geoffroy et al., 2010) havealso been developed, although more robust methods to obtain CDSDspectral width upon activation are presently lacking. This latter point isrelevant for triple-moment (3 M) bulk microphysics that aim to predictdistribution spectral width alongside number and mass concentrations(e.g., Loftus et al., 2014; Milbrandt and Yau, 2005).The current work extends the findings of C16 to the current LUTbasedaerosol activation scheme used in the Regional AtmosphericModeling System (RAMS) (Cotton et al., 2003; SC04; Saleeby and vanden Heever, 2013, hereafter SvdH13) and additionally examinesaerosol size and solubility impacts on newly activated CDSD properties.Because early cloud development processes such as condensationalgrowth, evaporation, and droplet self-collection depend on and impactCDSD spectral width (Hudson and Yum, 1997; Seifert and Beheng 2001;Lu and Seinfeld, 2006; Igel and van den Heever, 2017), an expansion ofthe activation LUTs to include CDSD spectral width is proposed as apreliminary step for extending 3M prediction to CDSDs in CRMs forimproved simulations of aerosol-cloud interactions.2. MethodologyThe current RAMS two-moment microphysics module determinesthe fractional number of aerosol particles that activate to cloud dropletsfrom five-dimensional LUTs based on model predicted air temperature(T), w, Na, and the geometric median radius (rg) and soluble fraction (ε)of the aerosol size distribution (SvdH13). These LUTs are created offlineusing a one-dimensional Lagrangian adiabatic parcel model (Feingoldand Heymsfield, 1992; Heymsfield and Sabin, 1989; SC04) to simulateexplicit droplet activation and initial CDSD growth for a range of ambientatmospheric conditions [T, w] and binned lognormal aerosol sizedistributions given by= ⎡⎣ ⎢− ⎤⎦ ⎥N r Nr π σr rσ( )2 lnexp[ln( / )]2(ln )aggg22 (1)where r is the dry aerosol particle bin radius and σg is the geometricstandard deviation of the distribution. As the parcel model simulationsfocus on the activation process, other processes such as coalescence,sedimentation, and mixing are not considered. Details of the parcelmodel can be found in SC04 and SvdH13, and only a brief description isprovided here. At the onset of parcel model calculations, the initiallydry aerosol particles in all bins first deliquesce and reach theirequilibrium diameters in a sub-saturated environment based on theKöhler equation for solution droplets. The parcel is then lifted at a fixedupward velocity w, and particle growth by vapor diffusion, along withconcurrent changes in the ambient environment, are iteratively computedusing the Variable-coefficient Ordinary Differential Equation(VODE) solver (Brown et al., 1989). The time resolution of these calculationsis determined within the VODE solver, and the frequency atwhich the solver is called is controlled by a longer model time stepbased on fixed upward parcel displacement increments (Δz) at thespecified w (Δt=Δz/w). Model calculations proceed until the parcelreaches a height 50m beyond the level of maximum saturation ratio(Smax) or total parcel displacement exceeds 2 km. Upon model termination,Smax and the fractional number of aerosols (factv) resulting innewly formed cloud droplets, defined as particles having diameters of atleast 2 μm, are cataloged in the LUTs according to the specified T, w, Na,rg, and ε parameter values.A critical point regarding the creation of these LUTs is the use of afixed minimum diameter (Dmin) to define cloud droplets in the parcelmodel, which can produce erroneous CDSD characteristics, particularlywithin the UL regime. For aerosol distributions with large rg valuesunder low SS conditions, for example, deliquesced aerosols within thelarge tail of the distribution can exceed 2 μm in diameter yet remainunactivated as ‘haze’ particles (Levin and Cotton, 2009; McFigganset al., 2006). For this study, aerosol particles activate to cloud dropletsbased on the critical diameter Dcrit as a function of parcel supersaturationratio (Sr) as in R09:D = σ MS RTρ83 ln( ) critsol wr w (2)where σsol is the surface tension of a solution droplet, Mw and ρw are themolar mass and density of liquid water, respectively, and R is theuniversal gas constant. Additionally, at relatively large w values withinthe AL regime, Nc stabilizes shortly after reaching supersaturation.However, parcel ascent and condensational growth continue beyondthe level of Smax, potentially causing additional narrowing of the CDSD.In the current work, model calculations terminate upon reaching Smaxas changes in Nc are negligible with continued ascent (Peng et al., 2007;R09).Parcel model simulations are performed to examine the sensitivitiesof CDSD characteristics to w, Na, rg, and ε, with the ranges for theseparameters listed in Table 1. Aerosols are assumed to be a mix of solubleand insoluble material of equal density, specified by ε, where fullysoluble aerosols correspond to ammonium sulfate with hygroscopicityparameter κ=0.61 (Petters and Kreidenweis, 2007). FollowingSvdH13, aerosol geometric standard deviation is fixed at σg=1.8, andaerosol distributions (Eq. 1) are partitioned into 100 logarithmicallyspacedbins spanning a size range specific to each rg value. For all simulations,Δz=1 m, and initial values of relative humidity, air temperatureand pressure are set to RH=0.99, T=10 °C and p=900 hPa,respectively.

Loftus, Adrian M.↗
DOE OSTI2021

Neutron stars in $$f(\mathtt {R,L_m})$$ gravity with realistic equations of state: joint-constrains with GW170817, massive pulsars, and the PSR J0030+0451 mass-radius from NICER data

Abstract In this work, we investigate neutron stars (NS) in $$f(\mathtt {R,L_m})$$ f ( R , L m ) theory of gravity for the case $$f(\mathtt {R,L_m})= \mathtt {R}+ \mathtt {L_m}+ \sigma \mathtt {R}\mathtt {L_m}$$ f ( R , L m ) = R + L m + σ R L m , where $$\mathtt {R}$$ R is the Ricci scalar and $$\mathtt {L_m}$$ L m the Lagrangian matter density. In the term $$\sigma \mathtt {R}\mathtt {L_m}$$ σ R L m , $$\sigma $$ σ represents the coupling between the gravitational and particles fields. For the first time the hydrostatic equilibrium equations in the theory are solved considering realistic equations of state and NS masses and radii obtained are subject to joint constrains from massive pulsars, the gravitational wave event GW170817 and from the PSR J0030+0451 mass-radius from NASA’s Neutron Star Interior Composition Explorer ( NICER ) data. We show that in this theory of gravity, the mass-radius results can accommodate massive pulsars, while the general theory of relativity can hardly do it. The theory also can explain the observed NS within the radius region constrained by the GW170817 and PSR J0030+0451 observations for masses around $$1.4~M_{\odot }$$ 1.4 M ⊙ .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1994

Design of orbital debris shields for oblique hypervelocity impact

A new impact debris propagation code was written to link CTH simulations of space debris shield perforation to the Lagrangian finite element code DYNA3D, for space structure wall impact simulations. This software (DC3D) simulates debris cloud evolution using a nonlinear elastic-plastic deformable particle dynamics model, and renders computationally tractable the supercomputer simulation of oblique impacts on Whipple shield protected structures. Comparison of three dimensional, oblique impact simulations with experimental data shows good agreement over a range of velocities of interest in the design of orbital debris shielding. Source code developed during this research is provided on the enclosed floppy disk. An abstract based on the work described was submitted to the 1994 Hypervelocity Impact Symposium.

Fahrenthold, Eric P.↗
NASA NTRS1996

Simulations of Active Vortex Generators

We are interested in the study, via numerical simulations, of active vortex generators. Vortex generators may be used to modify the inner part of the boundary layer or to control separation thus enhancing the performance and maneuverability of aerodynamic configurations. We consider generators that consist of a surface cavity elongated in the streamwise direction and partially covered with a moving lid that at rest lies flush with the boundary. Streamwise voracity is generated and ejected due to the oscillatory motion of the lid. The present simulations c Implement relevant experimental investigations of active vortex generators that have been conducted at NASA Ames Research Center and Stanford University. Jacobson and Reynolds used a piezoelectric device in water, allowing for small amplitude high frequency oscillations. They placed the lid asymmetrically on the cavity and observed a strong outward velocity at the small gap of the cavity. Saddoughi used a larger mechanically driven device in air to investigate this flow and observed a jet emerging from the wide gap of the configuration, contrary to the findings of Jacobson and Reynolds. More recently, Lachowiez and Wlezien are investigating the flow generated by an electro-mechanically driven lid to be used for assertion control in aerodynamic applications. We are simulating the flows generated by these devices and we are conducting a parametric study that would help us elucidate the physical mechanisms present in the flow. Conventional computational schemes encounter difficulties when simulating flows around complex configurations undergoing arbitrary motions. Here we present a formulation that achieves this task on a purely Lagrangian frame by extending the formulation presented by Koumoutsakos, Leonard and Pepin. The viscous effects are taken into account by modifying the strength of the particles, whereas fast multipole schemes employing hundreds of thousands ol'particle's allow for high resolution simulations. We shall present simulation results of an oscillating plate at various Reynolds numbers and Strouhal frequencies. Estimates of the forces needed to drive the devices will also be presented.

Mansour, N. N.↗
NASA NTRSDate not supplied

Modeling Self-Pressurization and Spray Bar Pressure Control of a Cryogenic Storage Tank in Normal Gravity

This paper presents Computational Fluid Dynamics (CFD) models for simulating self-pressurization and spray bar cooling processes in a large-scale liquid hydrogen storage tank under normal gravity. The self-pressurization model utilizes the kinetics-based Schrage equation with the Volume-Of-Fluid (VOF) model to account for the evaporative and condensing interfacial mass fluxes. It uses a laminar approach to model natural convection and compares the pressures and temperatures predicted by the CFD model with experimental data from the Multipurpose Hydrogen Test Bed (MHTB) self-pressurization experiment. Additionally, the predicted interfacial mass flux is presented. A CFD model for simulating pressure control using a spray bar is also presented. An Eulerian-Lagrangian approach models the interactions between the discrete droplets and the continuous ullage phase. The spray model is coupled with the VOF model by tracking particles in the ullage, removing particles from the ullage when they reach the interface, and then adding their contributions to the liquid. The T-sat model is used to calculate the droplet-ullage heat and mass transfer. In this model, droplets warm up to the saturation temperature corresponding to the ullage vapor pressure thereafter evaporating while remaining at the saturation temperature. The tank pressure, vapor, and liquid temperature evolutions are validated against data from the MHTB spray bar mixing experiment.

Self-Pressurization↗
DOE OSTI2023

QES-Plume v1.0: a Lagrangian dispersion model

Low-cost simulations providing accurate predictions of transport of airborne material in urban areas, vegetative canopies, and complex terrain are demanding because of the small-scale heterogeneity of the features influencing the mean flow and turbulence fields. Common models used to predict turbulent transport of passive scalars are based on the Lagrangian stochastic dispersion model. The Quick Environmental Simulation (QES) tool is a low-computational-cost framework developed to provide high-resolution wind and concentration fields in a variety of complex atmospheric-boundary-layer environments. Part of the framework, QES-Plume, is a Lagrangian dispersion code that uses a time-implicit integration scheme to solve the generalized Langevin equations which require mean flow and turbulence fields. Here, QES-Plume is driven by QES-Winds, a 3D fast-response model that computes mass-consistent wind fields around buildings, vegetation, and hills using empirical parameterizations, and QES-Turb, a local-mixing-length turbulence model. In this paper, the particle dispersion model is presented and validated against analytical solutions to examine QES-Plume’s performance under idealized conditions. In particular, QES-Plume is evaluated against a classical Gaussian plume model for an elevated continuous point-source release in uniform flow, the Lagrangian scaling of dispersion in isotropic turbulence, and a non-Gaussian plume model for an elevated continuous point-source release in a power-law boundary-layer flow. In these cases, QES-Plume yields a maximum relative error below 6 % when compared with analytical solutions. In addition, the model is tested against wind-tunnel data for a uniform array of cubical buildings. QES-Plume exhibits good agreement with the experiment with 99 % of matched zeros and 59 % of the predicted concentrations falling within a factor of 2 of the experimental concentrations. Furthermore, results also emphasize the importance of using high-quality turbulence models for particle dispersion in complex environments. Finally, QES-Plume demonstrates excellent computational performance.

58 GEOSCIENCES↗
DOE OSTI2022

A method for generating moving, orthogonal, area preserving polygonal meshes

A new method for generating locally orthogonal polygonal meshes from a set of generator points is presented in which polygon areas are a constraint. The area constraint property is particularly useful for particle methods where moving polygons track a discrete portion of material. Because Voronoi polygon meshes have some very attractive mathematical and numerical properties for numerical computation, a generalization of Voronoi polygon meshes was formulated that enforces a polygon area constraint. Area constrained moving polygonal meshes allow one to develop hybrid particle-mesh numerical methods that display some of the most attractive features of each approach. It is shown that this mesh construction method can continuously reconnect a moving, unstructured polygonal mesh in a pseudo-Lagrangian fashion without change in cell area/volume, and the method's ability to simulate various physical scenarios is shown. Overall, the advantages are identified for incompressible fluid flow calculations, with demonstration cases that include material discontinuities of all three phases of matter and large density jumps.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A particle finite element-based model for droplet spreading analysis

Here, a particle finite element method-based model is proposed to analyze droplet dynamics problems, particularly droplet spreading on solid substrates (wetting). The model uses an updated Lagrangian framework to formulate the governing equations of the liquid. The curvature of the liquid surface is tracked accurately using a deforming boundary mesh. In order to predict the spreading rate of the droplet on the solid substrate and track the corresponding contact angle evolution, dissipative forces at the contact line are included in the formulation in addition to the Navier-slip boundary conditions at the solid–liquid interface. The inclusion of these boundary conditions makes it possible to account for the induced Young’s stress at the contact line and for the viscous dissipation along the solid–liquid interfacial region. These are found to be essential to obtain a mesh-independent physical solution. The temporal evolution of the contact angle and the contact line velocity of the proposed model are compared with spreading droplets and micro-sessile droplet injection experiments and are shown to be in good agreement.

42 ENGINEERING↗