Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CLUSTER 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 1,027 records · Page 57

SU-E-T-283: Research of the Irradiation Damage to the Skin Cell by the Contaminative Electron in External Radiotherapy

Purpose: To investigate the micro-damage mechanism of the contaminative electron to the skin cell in external radiotherapy, the cell damage yield was simulated. Methods: The physical interaction and the energy deposited events of contaminative electrons in the cell DNA were modeled based on Geant4-DNA low energy physical model. The densitybased cluster mining algorithm was used to analyze the micro-damage yield and obtain its detailed compositive information. By taking the irradiation sensitive parameter and the clinical feature dose threshold of the skin cell into consideration, the damage features of the low energy contaminative electron to the skin cell were studied. The DSB and SSB yield and ratio, the cluster size and the root mean square radius, the cell SF and the lethal coefficient ε of the complex cDSB were researched. Results: For some very low energy electrons such as 20keV and 100keV, the cluster size can be more than 5 SSBs. The irradiation protection will be more crucial for the later response and the high α tissue. The lethal coefficient ε of the complex cDSB will increase with the incident electron number and the accumulative dose increasing for the dose square term in LQ model. However, the ε increment will increase 3%∼15% when the incident electron increases up to 100000. Conclusion: The simulation of the direct physical damage of the cell can be adjusted by a probability parameter to offset the simulation of the indirect biochemical damage. Thus the micro-damage mechanism of the contaminative electron to the skin cell can be detected at a certain extent by Monte Carlo physical simulation. This damage model of the low energy electron to DNA and these simulated results could be used to evaluate the damage effect of the low energy contaminative electron to the skin cells in the external radiotherapy. Strategic Priority Research Program of the Chinese Academy of Sciences (XDA03040000), Fundamental Research Funds for the Central Universities (2012HGXJ0057,2013HGXJ0193), Seed Foundation of Hefei University of Technology (2012HGZY0007) and Hefei University of Technology undergraduate innovative experiment project(2012CXCY431)

60 APPLIED LIFE SCIENCES↗

The gamma-ray pulsar population of globular clusters: Implications for the GeV excess

In this study, it has been suggested that the GeV excess, observed from the region surrounding the Galactic Center, might originate from a population of millisecond pulsars that formed in globular clusters. With this in mind, we employ the publicly available Fermi data to study the gamma-ray emission from 157 globular clusters, identifying a statistically significant signal from 25 of these sources (ten of which are not found in existing gamma-ray catalogs). We combine these observations with the predicted pulsar formation rate based on the stellar encounter rate of each globular cluster to constrain the gamma-ray luminosity function of millisecond pulsars in the Milky Way's globular cluster system. We find that this pulsar population exhibits a luminosity function that is quite similar to those millisecond pulsars observed in the field of the Milky Way (i.e. the thick disk). After pulsars are expelled from a globular cluster, however, they continue to lose rotational kinetic energy and become less luminous, causing their luminosity function to depart from the steady-state distribution. Using this luminosity function and a model for the globular cluster disruption rate, we show that millisecond pulsars born in globular clusters can account for only a few percent or less of the observed GeV excess. Among other challenges, scenarios in which the entire GeV excess is generated from such pulsars are in conflict with the observed mass of the Milky Way's Central Stellar Cluster.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Energy Efficient Engine Low Pressure Subsystem Aerodynamic Analysis

The objective of this study was to demonstrate the capability to analyze the aerodynamic performance of the complete low pressure subsystem (LPS) of the Energy Efficient Engine (EEE). Detailed analyses were performed using three- dimensional Navier-Stokes numerical models employing advanced clustered processor computing platforms. The analysis evaluates the impact of steady aerodynamic interaction effects between the components of the LPS at design and off- design operating conditions. Mechanical coupling is provided by adjusting the rotational speed of common shaft-mounted components until a power balance is achieved. The Navier-Stokes modeling of the complete low pressure subsystem provides critical knowledge of component acro/mechanical interactions that previously were unknown to the designer until after hardware testing.

Hall, Edward J.↗

Energy Efficient Engine Low Pressure Subsystem Flow Analysis

The objective of this project is to provide the capability to analyze the aerodynamic performance of the complete low pressure subsystem (LPS) of the Energy Efficient Engine (EEE). The analyses were performed using three-dimensional Navier-Stokes numerical models employing advanced clustered processor computing platforms. The analysis evaluates the impact of steady aerodynamic interaction effects between the components of the LPS at design and off-design operating conditions. Mechanical coupling is provided by adjusting the rotational speed of common shaft-mounted components until a power balance is achieved. The Navier-Stokes modeling of the complete low pressure subsystem provides critical knowledge of component aero/mechanical interactions that previously were unknown to the designer until after hardware testing.

Hall, Edward J.↗

A Theoretical Investigation of the Plausibility of Reactions Between Ammonia and Carbonyl Species (Formaldehyde, Acetaldehyde, and Acetone) in Interstellar Ice Analogs at Ultracold Temperatures

We have reexamined the reaction between formaldehyde and ammonia, which was previously studied by us and other workers in modestly sized cluster calculations. Larger model systems with up to 12H2O were employed, and reactions of two more carbonyl species, acetaldehyde and acetone, were also carried out. Calculations were performed at the B3LYP/6-31+G** level with bulk solvent effects treated with a polarizable continuum model; limited MP2/6-31+G** calculations were also performed. We found that while the barrier for the concerted proton relay mechanism described in previous work remains modest, it is still prohibitively high for the reaction to occur under the ultracold conditions that prevail in dense interstellar clouds. However, a new pathway emerged in more realistic clusters that involves at least one barrierless step for two of the carbonyl species considered here: ammonia reacts with formaldehyde and acetaldehyde to form a partial charge transfer species in small clusters (4H2O) and a protonated hydroxyamino intermediate species in large clusters (9H2O, 12H2O); modest barriers that decrease sharply with cluster size are found for the analogous processes for the acetone-NH3 reaction. Furthermore, if a second ammonia replaces one of the water molecules in calculations in the 9H2O clusters, deprotonation can occur to yield the same neutral hydroxyamino species that is formed via the original concerted proton relay mechanism. In at least one position, deprotonation is barrierless when zero-point energy is included. In addition to describing the structures and energetics of the reactions between formaldehyde, acetaldehyde, and acetone with ammonia, we report spectroscopic predictions of the observable vibrational features that are expected to be present in ice mixtures of different composition.

Chen, Lina↗

Probing the Structural Evolution of the Hydrated Electron in Water Cluster Anions (H 2 O) n – , n ≤ 200, by Electronic Absorption Spectroscopy

Electronic absorption spectra of water cluster anions (H 2 O) n - , n ≤ 200, at T = 80 K are obtained by photodissociation spectroscopy and compared with simulations from literature and experimental data for bulk hydrated electrons. Two almost isoenergetic electron binding motifs are seen for cluster sizes 20 ≤ n ≤ 40, which are assigned to surface and partially embedded isomers. With increasing cluster size, the surface isomer becomes less populated, and for n ≥ 50, the partially embedded isomer prevails. The absorption shifts to the blue, reaching a plateau at n ≈ 100. In this size range, the absorption spectrum is similar to that of the bulk hydrated electron but is slightly red-shifted; spectral moment analysis indicates that these clusters are reasonable model systems for hydrated electrons near the liquid-vacuum interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Intermediate-age globular clusters in four galaxy merger remnants

We present the results of combining Hubble Space Telescope optical photometry with ground-based K{sub s} -band photometry from the Gemini imagers NIRI and FLAMINGOS-I to study the globular cluster (GC) populations in four early-type galaxies that are candidate remnants of recent mergers (NGC 1700, NGC 2865, NGC 4382, and NGC 7727). These galaxies were chosen based on their blue colors and fine structure, such as shells and ripples that are indicative of past interactions. We fit the combined VIK{sub s} GC data with simple toy models of mixed cluster populations that contain three subpopulations of different age and metallicity. The fits, done via chi-squared mapping of the parameter space, yield clear evidence for the presence of intermediate-age clusters in each galaxy. We find that the ages of ∼1-2 Gyr for these GC subpopulations are consistent with the previously estimated merger ages for the host galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

The effect of galaxy triaxiality on globular clusters

Evidence of triaxiality and central mass concentration in the Galaxy and in M31 and M87 is examined. It is proposed that globular clusters on box orbits in the inner parts of these systems pass close enough to the center to be destroyed by tidal shocks. Remaining clusters will preferentially populate tube orbits with relatively high angular momentum. The process is used to explain the cluster distribution in M87 reported by Lauer and Kormendy (1986). Models are presented for cluster destruction by massive black holes in M87's dark halo. Consideration is given to techniques for testing the suggestion that clusters form only on box orbits.

Ostriker, J. P.↗

TRACING EMBEDDED STELLAR POPULATIONS IN CLUSTERS AND GALAXIES USING MOLECULAR EMISSION: METHANOL AS A SIGNATURE OF THE LOW-MASS END OF THE IMF

Most low-mass protostars form in clusters, in particular high-mass clusters; however, how low-mass stars form in high-mass clusters and what the mass distribution is are still open questions both in our own Galaxy and elsewhere. To access the population of forming embedded low-mass protostars observationally, we propose using molecular outflows as tracers. Because the outflow emission scales with mass, the effective contrast between low-mass protostars and their high-mass cousins is greatly lowered. In particular, maps of methanol emission at 338.4 GHz (J = 7{sub 0}–6{sub 0} A{sup +}) in low-mass clusters illustrate that this transition is an excellent probe of the low-mass population. We present here a model of a forming cluster where methanol emission is assigned to every embedded low-mass protostar. The resulting model image of methanol emission is compared to recent ALMA observations toward a high-mass cluster and the similarity is striking: the toy model reproduces observations to better than a factor of two and suggests that approximately 50% of the total flux originates in low-mass outflows. Future fine-tuning of the model will eventually make it a tool for interpreting the embedded low-mass population of distant regions within our own Galaxy and ultimately higher-redshift starburst galaxies, not just for methanol emission but also water and high-J CO.

79 ASTRONOMY AND ASTROPHYSICS↗

Tidal alignment of galaxies

We develop an analytic model for galaxy intrinsic alignments (IA) based on the theory of tidal alignment. We calculate all relevant nonlinear corrections at one-loop order, including effects from nonlinear density evolution, galaxy biasing, and source density weighting. Contributions from density weighting are found to be particularly important and lead to bias dependence of the IA amplitude, even on large scales. This effect may be responsible for much of the luminosity dependence in IA observations. The increase in IA amplitude for more highly biased galaxies reflects their locations in regions with large tidal fields. We also consider the impact of smoothing the tidal field on halo scales. We compare the performance of this consistent nonlinear model in describing the observed alignment of luminous red galaxies with the linear model as well as the frequently used 'nonlinear alignment model,' finding a significant improvement on small and intermediate scales. We also show that the cross-correlation between density and IA (the 'GI' term) can be effectively separated into source alignment and source clustering, and we accurately model the observed alignment down to the one-halo regime using the tidal field from the fully nonlinear halo-matter cross correlation. Inside the one-halo regime, the average alignment of galaxies with density tracers no longer follows the tidal alignment prediction, likely reflecting nonlinear processes that must be considered when modeling IA on these scales. Finally, we discuss tidal alignment in the context of cosmic shear measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

Synchronous parallel spatially resolved stochastic cluster dynamics

In this work, a spatially resolved stochastic cluster dynamics (SRSCD) model for radiation damage accumulation in metals is implemented using a synchronous parallel kinetic Monte Carlo algorithm. The parallel algorithm is shown to significantly increase the size of representative volumes achievable in SRSCD simulations of radiation damage accumulation. Additionally, weak scaling performance of the method is tested in two cases: (1) an idealized case of Frenkel pair diffusion and annihilation, and (2) a characteristic example problem including defect cluster formation and growth in α-Fe. For the latter case, weak scaling is tested using both Frenkel pair and displacement cascade damage. To improve scaling of simulations with cascade damage, an explicit cascade implantation scheme is developed for cases in which fast-moving defects are created in displacement cascades. For the first time, simulation of radiation damage accumulation in nanopolycrystals can be achieved with a three dimensional rendition of the microstructure, allowing demonstration of the effect of grain size on defect accumulation in Frenkel pair-irradiated α-Fe.

36 MATERIALS SCIENCE↗

Visualizing Clustering and Uncertainty Analysis with Multivariate Longitudinal Data

Multivariate time-series datasets are intrinsic to the study of dynamic, naturalistic behavior, such as in the applications of finance and motion video analysis. Statistical models provide the ability to identify event patterns in these data under conditions of uncertainty, but researchers must be able to evaluate how well a model uses available information in a dataset for clustering decisions and for uncertainty information. The Hidden Markov Model (HMM) is an established method for clustering time-series data, where the hidden states of the HMM are the clusters. We develop novel methods for quantifying the uncertainty of the performance of and for visualizing the clustering performance and uncertainty of fitting a HMM to multivariate time-series data. We explain the usefulness of uncertainty quantification and visualization with evaluating the performance of clustering models, as well as how information exploitation of time-series datasets can be enhanced. We implement our methods to cluster patterns of scanpaths from raw eye tracking data.

97 MATHEMATICS AND COMPUTING↗

The DESI One-Percent survey: exploring the Halo Occupation Distribution of Emission Line Galaxies with AbacusSummit simulations

The One-Percent survey of the Dark Energy Spectroscopic Instrument collected ~ 270k emission line galaxies (ELGs) at 0.8 < z < 1.6. The high completeness of the sample allowed the clustering to be measured down to scales never probed before, 0.04 Mpc/h in rp for the projected 2-point correlation function (2PCF) and 0.17 Mpc/h in galaxy pair separation s for the 2PCF monopole and quadrupole. The most striking feature of the measurements is a strong signal at the smallest scales, below 0.2 Mpc/h in rp and 1 Mpc/h in s. We analyse these data in the halo occupation distribution framework. We consider different distributions for central galaxies, a standard power law for satellites with no condition on the presence of a central galaxy and explore several extensions of these models. For all considered models, the mean halo mass of the sample is found to be log10 $\langle$Mh$\rangle$ ~ 11.9. We obtain a satellite mean occupation function which agrees with physically motivated ELG models only if we introduce central-satellite conformity, meaning that the satellite occupation is conditioned by the presence of central galaxies of the same type. To achieve in addition a good modelling of the clustering between 0.1 and 1 Mpc/h in rp, we allow for ELG positioning outside of the halo virial radius and find 0.5% of ELGs residing in the outskirts of halos. Furthermore, the satellite velocity dispersion inside halos is found to be ~ 30% larger than that of the halo dark matter particles. These are the main findings of our work. Here we investigate assembly bias as a function of halo concentration, local density or local density anisotropies and observe no significant change in our results. We split the data sample in two redshift bins and report no significant evolution with redshift. Lastly, changing the cosmology in the modelling impacts only slightly our results.

79 ASTRONOMY AND ASTROPHYSICS↗

Velocity correlations of galaxy clusters

We determine the velocity correlation function, pairwise peculiar velocity difference, and rms pairwise peculiar velocity dispersion of rich clusters of galaxies, as a function of pair separation, for three cosmological models: Omega = 1 and Omega = 0.3 cold dark matter (CDM), and Omega = 0.3 primeval baryonic isocurvature (PBI) models (all flat and Cosmic Background Explorer (COBE)-normalized). We find that close cluster pairs, with separation r is less than or equal to 10/h Mpc, exhibit strong attractive peculiar velocities in all models; the cluster pairwise velocities depend sensitively on the model. The mean pairwise attractive velocity of clusters on 5/h Mpc scale ranges from approximately 1700 km/s for Omega = 1 CDM to approximately 1000 km/s for PBI to approximately 700 km/s for Omega = 0.3 CDM. The small-scale pairwise velocities depend also on cluster mass: richer, more massive clusters exhibit stronger attractive velocities than less massive clusters. On large scales, from approximately 20 to 200/h Mpc, the cluster peculiar velocities are increasingly dominated by bulk and random motions; they are independent of cluster mass. The cluster velocity correlation function is negative on small scales for Omega = 1 and Omega = 0.3 CDM, indicating strong pairwise motion relative to bulk motion on small scales; PBI exhibits relatively larger bulk motions. The cluster velocity correlation function is positive on very large scales, from r approximately 10/h Mpc to r approximately 200/h Mpc, for all models. These positive correlations, which decrease monotonically with scale, indicate significant bulk motions of clusters up to approximately 200/h Mpc. The strong dependence of the cluster velocity functions on models, especially at small separations, makes them useful tools in constraining cosmological models when compared with observations.

Cen, Renyue↗

Coupled Modeling of Hydrodynamics and Sound in Coastal Ocean for Renewable Ocean Energy Development

An underwater sound model was developed to simulate sound propagation from marine and hydrokinetic energy (MHK) devices or offshore wind (OSW) energy platforms. Finite difference methods were developed to solve the 3D Helmholtz equation for sound propagation in the coastal environment. A 3D sparse matrix solver with complex coefficients was formed for solving the resulting acoustic pressure field. The Complex Shifted Laplacian Preconditioner (CSLP) method was applied to solve the matrix system iteratively with MPI parallelization using a high performance cluster. The sound model was then coupled with the Finite Volume Community Ocean Model (FVCOM) for simulating sound propagation generated by human activities, such as construction of OSW turbines or tidal stream turbine operations, in a range-dependent setting. As a proof of concept, initial validation of the solver is presented for two coastal wedge problems. This sound model can be useful for evaluating impacts on marine mammals due to deployment of MHK devices and OSW energy platforms.

Long, Wen↗

Lyman-alpha clouds as a relic of primordial density fluctuations

Primordial density fluctuations are studied using a CDM model and primordial clouds some of which are expanding, driven by pressure gradients created when the medium is photionized, and some of which are massive enough to continue collapsing in spite of the pressure. Normalization of CDM models to the clustering properties on large scales are used to predict the parameters of collapsing clouds of subgalactic mass at early epochs. It is shown that the abundance and dimensions of these clouds are comparable to those of the Lyman-alpha systems. The evolutionary history of the clouds is computed, utilizing a spherically symmetric hydrodynamics code with the dark matter treated as a collisionless fluid, and the H I column density distribution is evaluated as a function of N(H I) and redshift. The observed cloud parameters come out naturally in the CDM model and suggest that Lyman-alpha clouds are the missing link between primordial density fluctuations and the formation of galaxies.

Bond, J. R.↗

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗

Numerical Study of Turbulent Separation Bubbles with Varying Pressure Gradient and Reynolds Number

A family of cases each containing a small separation bubble is treated by direct numerical simulation (DNS), varying two parameters: the severity of the pressure gradients, generated by suction and blowing across the opposite boundary, and the Reynolds number. Each flow contains a well-developed entry region with essentially zero pressure gradient, and all are adjusted to have the same value for the momentum thickness, extrapolated from the entry region to the centre of the separation bubble. Combined with fully defined boundary conditions this will make comparisons with other simulations and turbulence models rigorous; we present results for a set of eight Reynolds-averaged Navier–Stokes turbulence models. Even though the largest Reynolds number is approximately 5.5 times higher than in a similar DNS study we presented in 1997, the models have difficulties matching the DNS skin friction very closely even in the zero pressure gradient, which complicates their assessment. In the rest of the domain, the separation location per se is not particularly difficult to predict, and the most definite disagreement between DNS and models is near reattachment. Curiously, the better models tend to cluster together in their predictions of pressure and skin friction even when they deviate from the DNS, although their eddy-viscosity levels are widely different in the outer region near the bubble (or they do not rely on an eddy viscosity). Stratford’s square-root law is satisfied by the velocity profiles, both at separation and reattachment. The Reynolds-number range covers a factor of two, with the Reynolds number based on the extrapolated momentum thickness equal to approximately 1500 and 3000. This allows tentative estimates of the improvements that even higher values will bring to the model comparisons. The solutions are used to assess models through pressure, skin friction and other measures; the flow fields are also used to produce effective eddy-viscosity targets for the models, thus guiding turbulence-modelling work in each region of the flow.

Coleman, G. N.↗